Computing the material phase at a particular point may involve
a call to a Python function.  That is usually quite costly.
Therefore, make sure that this call is really made only once
per quadrature point.