Handling the Finite Element Functions

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