| 27f02ce8 | 23-May-2019 |
Matthew G. Knepley <knepley@gmail.com> |
Plex+FE: Updates for hybrid integration - Added support for hybrid Jacobians - Added PetscDSHasBdJacobian() - Added boundary Jacobian preconditioner - Needed a path to integrate the hybrid Jacobian p
Plex+FE: Updates for hybrid integration - Added support for hybrid Jacobians - Added PetscDSHasBdJacobian() - Added boundary Jacobian preconditioner - Needed a path to integrate the hybrid Jacobian preconditioner - Fixed indexing in hybrid element mat update - Be careful about the difference between intrinsic dimension and embedding dimension for hybrid integration - Need to use rectangular transformation of gradient in case they are different - Fix use of embedding dimension - Make sure FEGeom struct has the dimensions set - Pointwise function derivatives are in the embedding space - Use PetscArrayzero() - Change TransitiveClosure() for hybrid cells - We get the closures of the endcaps, and then the hybrid faces, this allows degenerate cells that we need in PyLith
show more ...
|