Edit each elemental matrix of variational problem

Hello everyone.
I am triying to implement the calculus of curvature with ESFEM as i show below in the variatonal form

This problem is neccesary to be implemented with elements higher order than linear elements, but it is not possible because the reading of the mesh in Fenics for FEM in surfaces. The issue is that i want to acces to the elemental matrix for each element and edit it or define this for the user, is this possible?

This issue is a extension of this post Curvature by finite elements - variational formulation - FEniCS Project where is explained better

Thanks

You can use assemble_local, Element stiffness matrix - #2 by bhaveshshrimali

1 Like