Hello, I am new to Fenics and I was wondering if there is a way to evaluate the source and boundary condition functions exactly at the quadrature points rather than interpolating from the nodes as the ‘Expression’ functionality does. I thought that quadrature elements might be the solution, but I read that is just an optimization strategy whereby the values are interpolated and stored at the quadrature points.

Thank you in advance.