[petsc-dev] attempt at understanding PetscFEM

Jed Brown jedbrown at mcs.anl.gov
Mon Nov 18 18:08:50 CST 2013


Geoffrey Irving <irving at naml.us> writes:

> Here's an attempt at the tensor definitions of all the functions in
> PetscFEM.  Could you check if it looks correct, especially the
> Jacobian section?

Yes, this looks right.  I write it without the index notation in my
paper (section 3.2 and 3.3).

http://59A2.org/na/Brown-EfficientNonlinearSolversNodalHighOrder3D-2010.pdf


FWIW, I do not agree with Matt's factoring into arrays of function
pointers, especially since material models should be reused.  I think
he's willing to add support for having one function pointer that fills
up some tensors.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20131118/ded898d5/attachment.sig>


More information about the petsc-dev mailing list