norm_square        68 src/grid_deformation.cpp   const double norm_square
norm_square        72 src/grid_deformation.cpp   const double factor = sqrt( norm_square - 2 * Esfem::Grid::grid_dim() * t);