Handling the Finite Element Functions

Legacy DOLFIN
1

Hi all,

So if I have a mesh on my domain, and I’ve computed my stiffness and mass matrices, is there a way to get an expression for the finite element functions themselves at each mesh point so that I can manipulate them (add, subtract, integrate, etc.)? I’d like to replace, by manually computing, some entries of the stiffness matrix - in order to do this, I need to get a hold of the finite element functions themselves.

Thanks!

2

You can evaluate the basis functions at any point in your domain with evalute_basis

See for instance Computing shape function derivatives for each cell - #8 by dokken

3

Thank you!

Would you mind providing me with the arguments that this function takes? I am having difficulty finding documentation on it.