| 863027ab | 12-Jan-2023 |
Jed Brown <jed@jedbrown.org> |
DMBlockingType: add support for field node blocking
We used point blocking in the past, meaning that all dofs on a topological point (cell, face, edge, vertex) were densely coupled in the point-bloc
DMBlockingType: add support for field node blocking
We used point blocking in the past, meaning that all dofs on a topological point (cell, face, edge, vertex) were densely coupled in the point-block diagonal. This is expensive and often unnecessary. Field node blocking reduces PCVPBJACOBI setup and apply cost for cubic and higher elements.
show more ...
|
| 90612307 | 12-Jan-2023 |
Jed Brown <jed@jedbrown.org> |
DMProjectCoordinates: fix projection for isoperiodicity
The former code relied on DMCreateInterpolation, which generates a matrix, but a matrix is not sufficient for mapping from coordinate global s
DMProjectCoordinates: fix projection for isoperiodicity
The former code relied on DMCreateInterpolation, which generates a matrix, but a matrix is not sufficient for mapping from coordinate global space to global space since it doesn't contain any information about the affine effect on elements in the support of isoperiodic points. This code directly evaluates on the new dual functions using DMProjectFieldLocal.
DMSetCoordinateDM: allow passing NULL to clear the coordinate DEM
show more ...
|