petsc - use PetscFECreateLagrange

Solids - Fix default number of faces in the box mesh

solids - refactor to separate problem formulations

Mooney-Rivlin model (#794)

* refactored MR to work with current version of main. Updated doc with dS for MR

* fixed documentation formatting

* update dS for MR in documentation

* cleaning u

Mooney-Rivlin model (#794)

* refactored MR to work with current version of main. Updated doc with dS for MR

* fixed documentation formatting

* update dS for MR in documentation

* cleaning up MR dS equations and wording

* minor changes in doc

* MR dS draft 1"

* MR coded, jacobian diverging

* added FSInitial-MR1.h

* fixed call option

* NH no longer needs mu_1 option

* processphysics and upated naming scheme

* WIP: Notes on using finite differencing for dS

* added print strain at every increment cl option

* MR cleaned. started CSV file but not quite there

* examples/solids: write csv file

can be read using pandas:

In [2]: pd.read_csv('test.csv')
Out[2]:
increment energy
0 0.000000 0.000000
1 0.066667 0.002402
[...]

* examples/solids: add targets for CSV logging and Altair chart

The MR data is currently hacked because it isn't converging as I'd expect.

* fixed style

* solids/doc: added derivation of dS for MR

* added analytical dS in Jacobian evaluation

* fixed style

* updated doc and energy plot for MR

* MR almost done

* started removing integer division

* got rid of all integer divisions

* missed some 1/3

* remove extra print statement

* added coulped MR FSInitial-MRc qfunction

* MR coupled version doc update

* refactored MR to work with current version of main. Updated doc with dS for MR

* fixed documentation formatting

* update dS for MR in documentation

* cleaning up MR dS equations and wording

* minor changes in doc

* MR dS draft 1"

* MR coded, jacobian diverging

* added FSInitial-MR1.h

* fixed call option

* NH no longer needs mu_1 option

* processphysics and upated naming scheme

* added print strain at every increment cl option

* MR cleaned. started CSV file but not quite there

* examples/solids: write csv file

can be read using pandas:

In [2]: pd.read_csv('test.csv')
Out[2]:
increment energy
0 0.000000 0.000000
1 0.066667 0.002402
[...]

* examples/solids: add targets for CSV logging and Altair chart

The MR data is currently hacked because it isn't converging as I'd expect.

* fixed style

* solids/doc: added derivation of dS for MR

* added analytical dS in Jacobian evaluation

* fixed style

* updated doc and energy plot for MR

* MR almost done

* started removing integer division

* got rid of all integer divisions

* missed some 1/3

* fixed some typo

* rerunning plot files

* Mooney-Rivlin coupled model

* some minor changes

* Mooney-Rivlin coupled model

* some clean-up and update

* clean-up, rebased

* examples/solids: revise docs/derivation and fix Newton

The (fixed) Newton linearization yields an indefinite operator causing
Newton to fail (so output/MR-strain1.csv is broken). This can be handled
by switching to minres, but a better solution would be nice.

* format: added newlines and removed unused/commented code blocks

* doc: removed GP model dropdown.

* examples/solids: Mooney-Rivlin docs/notation

* examples/solids: Mooney-Rivlin notation and experimental configuration

* examples/solids: process -dm_plex_box_faces with PETSc main

* make style

* examples/solids: syntax fix and minor cleanup

* make style

* doc: fixing Holzapfel ref in Mooney-Rivlin

* added TESTARGS for Mooney-Rivlin

* fixed TESTARGS for MR

* tests/tap.sh: make test for solids examples with names

* WIP: adding test option. added -expect_final_strain_energy and RegressionTests_solids function but not sure how to actually pull the reference strain energy from the reference file

* examples/solids: add regression tests using expected strain energy

* examples/solids: error on MMS with model that does not provide manufactured solutions

* examples/solids: fix HIP via CEED_QFUNCTION_HELPER

* update in-line comments

* changed as Jeremy suggested

* examples/solids: make style

Co-authored-by: Rezgar Shakeri <Rezgar.Shakeri@colorado.edu>
Co-authored-by: Jed Brown <jed@jedbrown.org>

show more ...

examples/solids: rewrite Cauchy-Green (b-I) in terms of grad(u)/Green-Euler strain for better stability (#769)

* FSCurrent-NH1.h: rewrite (b-I3) in terms of Grad(u)

* fixed the error

* rewrite

examples/solids: rewrite Cauchy-Green (b-I) in terms of grad(u)/Green-Euler strain for better stability (#769)

* FSCurrent-NH1.h: rewrite (b-I3) in terms of Grad(u)

* fixed the error

* rewrite (b-I3) in terms of Grad(u) in FSCurrent-NH2.h

* update documentation

* doc: defined Green_Euler strain tensor

* Update examples/solids/index.rst

Co-authored-by: Jed Brown <jed@jedbrown.org>

* added eq numbers in tip

* Update examples/solids/index.rst

Co-authored-by: Jed Brown <jed@jedbrown.org>

Co-authored-by: Jed Brown <jed@jedbrown.org>

show more ...

petsc - new compatibility macro for DMAddBoundary

style - updating C variable style

Install install backend headers under include/ceed/

This makes it possible to distribute source plugins that provide
additional backends. It's also used in MFEM, perhaps temporarily.

Deprecate ceed

Install install backend headers under include/ceed/

This makes it possible to distribute source plugins that provide
additional backends. It's also used in MFEM, perhaps temporarily.

Deprecate ceed-backend.h, which was not previously installed, but some
users accessed it from an in-place build.

Also install CUDA and HIP headers that allow users to provide CUfunction
and hipFunction_t.

Co-authored-by: Jeremy L. Thompson <jeremy.thompson@colorado.edu>
Requested-by: Andrew T. Barker <barker29@llnl.gov>

show more ...

examples/solids: Add residual and Jacobian in current configuration

There are some benefits to evaluating residuals and/or Jacobians in
current configuration. This includes simpler representation of

examples/solids: Add residual and Jacobian in current configuration

There are some benefits to evaluating residuals and/or Jacobians in
current configuration. This includes simpler representation of body
forces and contact, as well as lower storage and lower arithmetic in the
Jacobian action. We add these implementations as well documentation to
derive their representations via linearize-then-pushforward and
pushforward-then-linearize.

Co-authored-by: Arash Mehraban <arashm81@gmail.com>
Co-authored-by: Jed Brown <jed@jedbrown.org>
Co-authored-by: Jeremy L. Thompson <jeremy.thompson@colorado.edu>

show more ...

Command line option for viz in Fluids and Solids examples

examples/fluids: change the cl option -of to -output_dir for the output folder

examples/solids: add cl option -output_dir for the viz outpu

Command line option for viz in Fluids and Solids examples

examples/fluids: change the cl option -of to -output_dir for the output folder

examples/solids: add cl option -output_dir for the viz output

change: output folder -> output directory

Co-authored-by: Jed Brown <jed@jedbrown.org>

examples: change variable name outputfolder -> outputdir

examples: update the README files

show more ...

WIP: Solids-GivensBC (#702)

* Givens boundary condition added.

Co-authored with: @jedbrown and @jeremylt

* Comments first letter capitalized

examples/petsc: use VecGetArrayAndMemType() to support CUDA/HIP/Kokkos

We resolve default VECCUDA or VECHIP based on /gpu/cuda or /gpu/hip. For
others, one should pass -dm_vec_type kokkos or -dm_vec

examples/petsc: use VecGetArrayAndMemType() to support CUDA/HIP/Kokkos

We resolve default VECCUDA or VECHIP based on /gpu/cuda or /gpu/hip. For
others, one should pass -dm_vec_type kokkos or -dm_vec_type standard.

There is no longer a "libCEED User Requested MemType" field in the
output.

Enable CI testing using PETSc built locally with HIP. This should be
converted to building a specified version once some bugs are fixed and
the branch merges.

Exclude automatic selection of VECHIP for /gpu/hip/occa due to SEGV:
https://github.com/CEED/libCEED/issues/678

show more ...

solids - traction BCs

solids - fix missing ()s in clamp BCs

QFunction Context Data Object (#596)

* Ctx - create context object for QFunction context data

* Context - rename UserContext -> QFunctionContext

* Ctx - add lcov markers

* Ctx - fix leak in

QFunction Context Data Object (#596)

* Ctx - create context object for QFunction context data

* Context - rename UserContext -> QFunctionContext

* Ctx - add lcov markers

* Ctx - fix leak in identity QFunctions

* Hip/Cuda - rename sync functions for vector/context

* Tests - lcov marker update

* QFunction - drop unused function

* Python - fix copy-paste errors

* Ctx - update notes for Fortran usage

* Fortran - drop unneeded cast

Co-authored-by: Jed Brown <jed@jedbrown.org>

* Interface - use void* for SetData interfaces

* Make - use call quiet for NVCC

* Interface - use void* for GetData interfaces

* Make - add quiet call option for examples

* Makefile - create common makefile to reduce duplication/complexity in example makefiles

Co-authored-by: Jed Brown <jed@jedbrown.org>

show more ...

PETSc examples: Update DMPlexCreateSphereMesh() and DMAddBoundary() calls after changes in PETSc's API (#600)

* PETSc examples: Update DMPlexCreateSphereMesh call after change in API

* PETSc bps:

PETSc examples: Update DMPlexCreateSphereMesh() and DMAddBoundary() calls after changes in PETSc's API (#600)

* PETSc examples: Update DMPlexCreateSphereMesh call after change in API

* PETSc bps: Update DMAddBoundary() call after change in API

* Add compatibility macros

* solids/ and fluids/ : Update DMAddBoundary() call and add compatibility macros

* Travis: Update PETSc ref with latest hash from master

show more ...

Op - add interface for multigrid level creation (#579)

* Op - add interface for multigrid level creation

* Op - add implementation for OperatorMultigridLevelCreate

* make style

* make tidy

Op - add interface for multigrid level creation (#579)

* Op - add interface for multigrid level creation

* Op - add implementation for OperatorMultigridLevelCreate

* make style

* make tidy

* Op - add test t550, fix errors

* Tests - add Fortran version of t550

* Tests - add t511 for testing tensor basis multigrid level setup

* make style and tidy

* Tests - fix t55* memory leaks

* Tests - add t552 for non-tensor basis multigrid levels

* CUDA - use CeedIntMax in shared CUDA backend

* Tests - add OCCA test exception for t55*

* Op - add lvector global prolongation multiplicity, simplifies user interface

* Solids - convert example to new interface

* make style

* Tests - convert t550 to multicomponent

* Solids - drop unused ceed_fine

* Python - add new multigrid level interface

* Python - fix operator wrap, use ceed python obj rather than ceed pointer

* Gallery - update comment slightly

* Tests - remove accidental duplicate test

* Multigrid - add size=2 case as well

* Operator - drop unneeded inline

* QFunction - simplify context ownership to match vector

* make style

* Python - update multigrid function signature

* Operator - refactor prolong/restrict qfunctions as scaling qfunctions

* Vector - add testing for reciprocal and add to Fortran/Python interfaces

* CUDA - add VectorReciprocal on device

* Gallery - drop specalized versions for 'Scale', wil fix performance hit later

* Hip - add vector reciprocal

* Operator - add more flexible prologation basis creation interface

* Vec - make sure data is set for VectorReciprocal

* Tests - drop ncomp for t550/1 so kernel is not too large for Magma backend

* Tests - add missing lcov markers

* make style

* Travis - allow ARM job to fail

* Travis - fix intel install

* Travis - try different install dir name for inteloneapi

* Travis - add ifort, ipp packages

* Tests - add missing lcov marker

show more ...

Vec - add CeedVectorTakeArray

Solids - add comment explaining use of Set/Sync array for debugging

Solids - add missing MatShellSetVecType calls

Solids - add missing address-of operator

make style

Solids - enable use of GPU vectors

Op - change LinearAssemble* to accept CeedVector instead of pointer to CeedVector, allows for less memory movement and interfaces with parent code better

style - rename AssembleLinear* to LinearAssemble*