Alternate minimization.
Using PETSc constrained non-linear solver for the phase field.
Few updates on variable names and cleanups.

- Adds issue report on CI log
- Adds full html report as artifact
- Adds code metric lizard_issues (currently 87)

A generalized property type for (almost) arbitrary functions. Based on
the exprtk library also used in the Function-type parameters.

In this initial implementation all *scalar* variables from the variable
array passed for the value and derivative evaluation can be used. Return
values are either scalar, 2- or 3-vectors, and 2x2- or 3x3 matrices.

The evaluation is somewhat slower than the dedicated implementation; in
the case of component transport process a bilinear property was replaced
by exactly the same function-type property and the total run-time for
the ctest has increased by around 6%.

Allowing not only a zero value to be applied to the inactive
nodes but also heterogeneous, time-dependent values.

 to WaterVapourLatentHeatWithCriticalTemperature

and changed saturation notation