Geoffrey Irving <[email protected]> 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.
pgpGy4W2JTM2I.pgp
Description: PGP signature
