How to get the first order derivative of the solution and store it in a numpy array

[dokken]

@raboynics

ChatGPT suggests to compute the derivative at ‘‘quadrature points’’ within each cell, rather than at the vertices or edges. It says those are the points where integration is performed and that the derivatives are well defined at those points. Does it make any sense?

Yes, that makes sense.

@raboynics

If so, is there a way to implement this with Fenicsx?

See for instance: Numerical values from ufl.SpatialCoordinate - #2 by dokken