Expressions obtained via subclassing Expression in Python are interpolated not evaluated on

I have written code in Python using fenics, dolfin etc. to solve the three dimensional Schroedinger Equation. For my choice of potential and weight function I observer slow or non-existing convergence to theoretically known values. After counting the number of calls to the python class implementing the potential I realized, that this does not depend on the interpolation order chosen via parameters["form_compiler"]["quadraturedegree"] but is consistent with being called at all nodes of the interpolation. Is there a way to change a parameter such that the Expression is evaluated at all integration grid points, of course at the cost of much larger CPU time? I will not be able to use fenics for publishable results if this cannot be done/fixed

Comments (6)