296       f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7);  in FormFunctionGradient()
297       f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4);  in FormFunctionGradient()
298       f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8);  in FormFunctionGradient()
299       f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2);  in FormFunctionGradient()
300       f5 = PetscSqrtScalar(1.0 + d2 * d2 + d5 * d5);  in FormFunctionGradient()
301       f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6);  in FormFunctionGradient()
321       ft = ft + PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2);  in FormFunctionGradient()
328       ft = ft + PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2);  in FormFunctionGradient()
336       ft = ft + PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4);  in FormFunctionGradient()
343       ft = ft + PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4);  in FormFunctionGradient()
350     ft += PetscSqrtScalar(1.0 + d1 * d1 + d2 * d2);  in FormFunctionGradient()
355     ft += PetscSqrtScalar(1.0 + d1 * d1 + d2 * d2);  in FormFunctionGradient()
503       f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7);  in FormHessian()
504       f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4);  in FormHessian()
505       f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8);  in FormHessian()
506       f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2);  in FormHessian()
507       f5 = PetscSqrtScalar(1.0 + d2 * d2 + d5 * d5);  in FormHessian()
508       f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6);  in FormHessian()
669         fnorm = PetscSqrtScalar(nf1 * nf1 + nf2 * nf2);  in MSA_BoundaryConditions()