1469 …PetscInt         dim, coordDim, numFields, numComponents = 0, qNc, Nq, cellHeight, cStart, cEnd, c…  in DMPlexComputeL2DiffLocal()  local
1474   PetscCall(DMGetCoordinateDim(dm, &coordDim));  in DMPlexComputeL2DiffLocal()
1475   fegeom.dimEmbed = coordDim;  in DMPlexComputeL2DiffLocal()
1504 …mponents, &interpolant, coordDim * (Nq + 1), &coords, Nq, &fegeom.detJ, coordDim * coordDim * Nq, …  in DMPlexComputeL2DiffLocal()
1540         qgeom.J        = &fegeom.J[q * coordDim * coordDim];  in DMPlexComputeL2DiffLocal()
1541         qgeom.invJ     = &fegeom.invJ[q * coordDim * coordDim];  in DMPlexComputeL2DiffLocal()
1545           gcoords = &coords[coordDim * Nq];  in DMPlexComputeL2DiffLocal()
1546 …MPlexBasisTransformApplyReal_Internal(dm, &coords[coordDim * q], PETSC_TRUE, coordDim, &coords[coo…  in DMPlexComputeL2DiffLocal()
1548           gcoords = &coords[coordDim * q];  in DMPlexComputeL2DiffLocal()
1551         ierr = (*funcs[field])(coordDim, time, gcoords, Nc, funcVal, ctx);  in DMPlexComputeL2DiffLocal()
1557 …if (transform) PetscCall(DMPlexBasisTransformApply_Internal(dm, &coords[coordDim * q], PETSC_FALSE…  in DMPlexComputeL2DiffLocal()
1564 …double)(coordDim > 0 ? coords[coordDim * q] : 0), (double)(coordDim > 1 ? coords[coordDim * q + 1]…  in DMPlexComputeL2DiffLocal()
1594 …PetscInt         dim, coordDim, qNc = 0, Nq = 0, numFields, numComponents = 0, cStart, cEnd, c, fi…  in DMComputeL2GradientDiff_Plex()  local
1599   PetscCall(DMGetCoordinateDim(dm, &coordDim));  in DMComputeL2GradientDiff_Plex()
1600   fegeom.dimEmbed = coordDim;  in DMComputeL2GradientDiff_Plex()
1621 … &funcVal, coordDim * (Nq + 1), &coords, coordDim * coordDim * Nq, &fegeom.J, coordDim * coordDim …  in DMComputeL2GradientDiff_Plex()
1649         qgeom.J        = &fegeom.J[q * coordDim * coordDim];  in DMComputeL2GradientDiff_Plex()
1650         qgeom.invJ     = &fegeom.invJ[q * coordDim * coordDim];  in DMComputeL2GradientDiff_Plex()
1654           gcoords = &coords[coordDim * Nq];  in DMComputeL2GradientDiff_Plex()
1655 …MPlexBasisTransformApplyReal_Internal(dm, &coords[coordDim * q], PETSC_TRUE, coordDim, &coords[coo…  in DMComputeL2GradientDiff_Plex()
1657           gcoords = &coords[coordDim * q];  in DMComputeL2GradientDiff_Plex()
1660         ierr = (*funcs[field])(coordDim, time, gcoords, n, Nc, funcVal, ctx);  in DMComputeL2GradientDiff_Plex()
1666 …if (transform) PetscCall(DMPlexBasisTransformApply_Internal(dm, &coords[coordDim * q], PETSC_FALSE…  in DMComputeL2GradientDiff_Plex()
1882 …PetscInt         dim, coordDim, numFields, numComponents = 0, qNc, Nq, cStart, cEnd, c, field, fie…  in DMPlexComputeL2DiffVec()  local
1887   PetscCall(DMGetCoordinateDim(dm, &coordDim));  in DMPlexComputeL2DiffVec()
1916 …mComponents, &interpolant, coordDim * Nq, &coords, Nq, &fegeom.detJ, coordDim * coordDim * Nq, &fe…  in DMPlexComputeL2DiffVec()
1946           qgeom.J        = &fegeom.J[q * coordDim * coordDim];  in DMPlexComputeL2DiffVec()
1947           qgeom.invJ     = &fegeom.invJ[q * coordDim * coordDim];  in DMPlexComputeL2DiffVec()
1950           PetscCall((*funcs[field])(coordDim, time, &coords[q * coordDim], Nc, funcVal, ctx));  in DMPlexComputeL2DiffVec()
2288 …PetscInt         dim, coordDim, numFields, numComponents = 0, qNc, Nq, cStart, cEnd, vStart, vEnd,…  in DMPlexComputeGradientClementInterpolant()  local
2295   PetscCall(DMGetCoordinateDim(dm, &coordDim));  in DMPlexComputeGradientClementInterpolant()
2296   fegeom.dimEmbed = coordDim;  in DMPlexComputeGradientClementInterpolant()
2321 …coordDim * numComponents * 2, &gradsum, coordDim * numComponents, &interpolant, coordDim * Nq, &co…  in DMPlexComputeGradientClementInterpolant()
2329     PetscCall(PetscArrayzero(gradsum, coordDim * numComponents));  in DMPlexComputeGradientClementInterpolant()
2333       PetscScalar   *grad = &gradsum[coordDim * numComponents];  in DMPlexComputeGradientClementInterpolant()
2345         PetscCall(PetscArrayzero(grad, coordDim * numComponents));  in DMPlexComputeGradientClementInterpolant()
2359           qgeom.J        = &fegeom.J[q * coordDim * coordDim];  in DMPlexComputeGradientClementInterpolant()
2360           qgeom.invJ     = &fegeom.invJ[q * coordDim * coordDim];  in DMPlexComputeGradientClementInterpolant()
2368 …for (d = 0; d < coordDim; ++d) grad[fc * coordDim + d] += interpolant[fc * dim + d] * wt * fegeom.…  in DMPlexComputeGradientClementInterpolant()
2377         for (d = 0; d < coordDim; ++d) gradsum[fc * coordDim + d] += grad[fc * coordDim + d];  in DMPlexComputeGradientClementInterpolant()
2383           for (d = 0; d < coordDim; ++d) {  in DMPlexComputeGradientClementInterpolant()
2385 …     PetscCall(PetscPrintf(PETSC_COMM_SELF, "%g", (double)PetscRealPart(grad[fc * coordDim + d])));  in DMPlexComputeGradientClementInterpolant()
2392       for (d = 0; d < coordDim; ++d) gradsum[fc * coordDim + d] /= volsum;  in DMPlexComputeGradientClementInterpolant()