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.

Because such statements bleed into all the files that include
the headers, and this may lead to compilation errors later.