According to Hanne, this is needed to turn 1d harmonic map heat flow into true
curve-shortening flow.

Hardwiring the weighting here is very hacky.  Something more flexible and
general will need to be implemented eventually.

This are needed by the gradientflow code to implement the implicit Euler method
for an L2 gradient flow.