Changed gradient type from double[dims] to Point<dims, double> because preferred by Paraview and because Point has already a function norm, s.t. getting the gradient magnitude doesn't require an extra function.