ENH: runTime calculation of (dp/dt)^2 - issue #224

- implemented using magSqr() instead of sqr().
  For scalar fields they are the same, but can be useful
  if this function object is extended for more field types.