| c9ece94a | 13-Apr-2018 |
Matthew G. Knepley <knepley@gmail.com> |
Merge remote-tracking branch 'origin/knepley/fix-fe-bd-integral'
* origin/knepley/fix-fe-bd-integral: Plex+SNES: Fix BdResidual - Use DMPlexGetSubpoint() for aux closure, in case it is a boundary
Merge remote-tracking branch 'origin/knepley/fix-fe-bd-integral'
* origin/knepley/fix-fe-bd-integral: Plex+SNES: Fix BdResidual - Use DMPlexGetSubpoint() for aux closure, in case it is a boundary mesh - Be careful about Plex conversion FE: Fixed boundary integration - Allow auxiliary dm to be a different dimension, need dimAux and different tabulation FE: Small confusion about topological vs coordinate dimension PetscDS: The u_x array must have size dimEmbed because the derivatives are a function of the coordinate dim, not topological dim Plex: Added DMPlexGetSubpoint() - This returns the point in a submesh corresponding to some point in the original mesh - This replaces custom code in projection
show more ...
|
| 32c2c4bc | 14-Mar-2018 |
Toby Isaac <tisaac@cc.gatech.edu> |
PetscFEGeomComplete(): correct normal vector in 2D.
We get the normal vector from the unit basis vector added to complete the Jacobian in a way that makes the determinant positive. In 3D, this is f
PetscFEGeomComplete(): correct normal vector in 2D.
We get the normal vector from the unit basis vector added to complete the Jacobian in a way that makes the determinant positive. In 3D, this is fine, because by convention the facets are oriented with outward normal, but in 2D they are not, so we have to invert the sign.
show more ...
|