Dear FEniCS family,
I was trying to implement C1 continuity for quadrilateral surface elements. I tried penalty method based on biharmonic example, but, that doesn’t give robust technique to work. The results are highly fluctuating with penalty term for quadrilateral mesh with CG (2nd order) elements.
I wanted to implemented custom Hermite shape functions (second order) in dolfinx where I use the same global assembly and symbolic computation. Can you please provide any resource where such is implemented and any robust method to apply C1 continuity.
I am using gmsh .msh file to input mesh data, so isogeometric doesn’t work for me.
Thanks