Sounds good. Note that if you use direct solvers, that extra node is basically
free because the vertex separators are unchanged. It's a marginal cost in the
storage of assembled matrices and the length of state vectors. And in 3D, even
less significant.
Susanne Claus writes:
> Dear Matthew
Dear Matthew and Jed,
Brilliant. Thank you so much!
Your changes work like a charm Matthew (I tested your branch on the gmsh
file I sent) and thank you so much for your advice Jed. The loss of one
order of convergence for an inf-sup stable pressure discretization seems
indeed a very high
Jed is right about the numerics. However, this does not look hard. Here is
my try at it:
https://gitlab.com/petsc/petsc/-/merge_requests/4838
Please tell me if this works and I will make a test and merge.
Thanks,
Matt
On Thu, Feb 10, 2022 at 6:47 PM Jed Brown wrote:
> Susanne, do
Susanne, do you want PetscFE to make the serendipity (8-node) finite element
space or do you just want to read these meshes? I.e., would it be okay with you
if the coordinates were placed in a Q_2 (9-node, biquadratic) finite element
space?
This won't matter if you're traversing the dofs per
On Thu, Feb 10, 2022 at 10:12 AM Susanne Claus
wrote:
> Hello,
>
> I am using DMPlex for the mesh structure of a solid mechanics finite
> element code. I mainly use gmsh as input file format. When I try to read in
> 8-noded Quadrilaterals (Element type 16 in gmsh) DMPlex tells me that this
>