450 …PetscInt         dim, coordDim, numFields, numComponents = 0, qNc, Nq, cStart, cEnd, vStart, vEnd,…  in InterpolateGradient()  local
456   PetscCall(DMGetCoordinateDim(dm, &coordDim));  in InterpolateGradient()
457   fegeom.dimEmbed = coordDim;  in InterpolateGradient()
482 …coordDim * numComponents * 2, &gradsum, coordDim * numComponents, &interpolant, coordDim * Nq, &co…  in InterpolateGradient()
489     PetscCall(PetscArrayzero(gradsum, coordDim * numComponents));  in InterpolateGradient()
493       PetscScalar   *grad = &gradsum[coordDim * numComponents];  in InterpolateGradient()
504         PetscCall(PetscArrayzero(grad, coordDim * numComponents));  in InterpolateGradient()
521 …for (d = 0; d < coordDim; ++d) grad[fc * coordDim + d] += interpolant[fc * dim + d] * wt * fegeom.…  in InterpolateGradient()
528         for (d = 0; d < coordDim; ++d) gradsum[fc * coordDim + d] += grad[fc * coordDim + d];  in InterpolateGradient()
533           for (d = 0; d < coordDim; ++d) {  in InterpolateGradient()
535 …     PetscCall(PetscPrintf(PETSC_COMM_SELF, "%g", (double)PetscRealPart(grad[fc * coordDim + d])));  in InterpolateGradient()