| 4f9ab2b4 | 25-Feb-2022 |
Jed Brown <jed@jedbrown.org> |
DMPlex: skip coordinate projection when FE space matches modulo quadrature
This requires checking that dual spaces are equal so that an legit projection is not required. For example, high order Gmsh
DMPlex: skip coordinate projection when FE space matches modulo quadrature
This requires checking that dual spaces are equal so that an legit projection is not required. For example, high order Gmsh meshes have equally spaced vertices while we prefer Gauss-Lobatto spacing, so that's a case where a projection is required. (It could be made faster, but that's for another day.)
show more ...
|
| 5b8243e1 | 08-Mar-2022 |
Jed Brown <jed@jedbrown.org> |
DM: set local and global vector sizes based on gcd
This enables having block size 3 for elasticity with Q3 elements, for example. The alternative reduces to bs=1 in such cases because the block size
DM: set local and global vector sizes based on gcd
This enables having block size 3 for elasticity with Q3 elements, for example. The alternative reduces to bs=1 in such cases because the block sizes don't match. Also use gcd for DMGetLocalToGlobalMapping.
show more ...
|