Hello, NGSolve community,
Does mesh.SplitElements_Alfeld() work for 3D?
When I uses the following code:
print ("volels: ", mesh.GetNE(VOL))
print ("bndels: ", mesh.GetNE(BND))
before using mesh.SplitElements_Alfeld() and after using mesh.SplitElements_Alfeld() , the number of volume elements seems to be the same.
Many thanks
Best regards,
Lina
I think it’s not.
A more general interface to a wider class of splittings is under construction.
Joachim
Thank you so much, Joachim, for your classification.
Best regards,
Lina
While I don’t think Alfeld split is supported by native NGSolve in 3D, it is supported via ngsPETSc (I’ve been wanting to try implementing that for a while!).
Have a look here: PETSc DMPlex — ngsPETSc 0.0.1 documentation