| ac9d17c7 | 17-Feb-2025 |
Matthew G. Knepley <knepley@gmail.com> |
Plex+PetscFE: Fixes to allow embedded manifolds
- Add PetscFEGeomMode to describe geometry
- Reorganized PetscFEGeom
- PetscFEGeomCreate(), DMFieldCreateFEGeom() and DMSNESGetFEGeom() now take the Pe
Plex+PetscFE: Fixes to allow embedded manifolds
- Add PetscFEGeomMode to describe geometry
- Reorganized PetscFEGeom
- PetscFEGeomCreate(), DMFieldCreateFEGeom() and DMSNESGetFEGeom() now take the PetscFEGeomMode
- J and invJ have dE dimensions
- f1, g1, g2, g3 have dE dimensions
- FIX the dim argument of all point function is really dE because it is used to index the gradients
- Clarified behavior of PetscFEGeomGetCellPoint()
- Now -dm_plex_coordinate_dim can change the coordinate dimension of a mesh created by Plex
- Implicit cell geometry must be used when dim != cdim
show more ...
|
| d5def619 | 17-Feb-2025 |
Jonas Heinzmann <jheinzmann@ethz.ch> |
SNESLineSearch: fix computation of directional derivative in the case of the VI solvers
The two line search algorithms bisection and critical point work with the directional derivative, assuming tha
SNESLineSearch: fix computation of directional derivative in the case of the VI solvers
The two line search algorithms bisection and critical point work with the directional derivative, assuming that F(x) = grad G(x).
However, when working with the variational inequality solvers, the constraints were not correctly considered when computing the directional derivative.
This commit introduces a new linesearch->ops->vidirderiv function which correctly considers the constraints (similar to linesearch->ops->viproject and linesearch->ops->vinorm).
show more ...
|