1519       PetscInt     Nb, Nc, q, fc;  in DMPlexComputeL2DiffLocal()  local
1561         for (fc = 0; fc < Nc; ++fc) {  in DMPlexComputeL2DiffLocal()
1562           const PetscReal wt = quadWeights[q * qNc + (qNc == 1 ? 0 : qc + fc)];  in DMPlexComputeL2DiffLocal()
1564 …cInt_FMT ",%" PetscInt_FMT " point %g %g %g diff %g (%g, %g)\n", c, field, fc, (double)(coordDim >…  in DMPlexComputeL2DiffLocal()
1565 …lPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q]), (double)PetscRealPart(interpolant[fc…  in DMPlexComputeL2DiffLocal()
1566           elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q];  in DMPlexComputeL2DiffLocal()
1634       PetscInt    Nb, Nc, q, fc;  in DMComputeL2GradientDiff_Plex()  local
1670           for (fc = 0; fc < Nc; ++fc) {  in DMComputeL2GradientDiff_Plex()
1673             for (d = 0; d < dim; ++d) sum += interpolant[fc * dim + d] * n[d];  in DMComputeL2GradientDiff_Plex()
1674             interpolant[fc] = sum;  in DMComputeL2GradientDiff_Plex()
1677         for (fc = 0; fc < Nc; ++fc) {  in DMComputeL2GradientDiff_Plex()
1678           const PetscReal wt = quadWeights[q * qNc + (qNc == 1 ? 0 : qc + fc)];  in DMComputeL2GradientDiff_Plex()
1679 …T ",%" PetscInt_FMT " diff %g\n", c, field, fc, (double)(PetscSqr(PetscRealPart(interpolant[fc] - …  in DMComputeL2GradientDiff_Plex()
1680           elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q];  in DMComputeL2GradientDiff_Plex()
1774         PetscInt     Nb, Nc, q, fc;  in DMComputeL2FieldDiff_Plex()  local
1813           for (fc = 0; fc < Nc; ++fc) funcVal[fc] = 0.;  in DMComputeL2FieldDiff_Plex()
1825           for (fc = 0; fc < Nc; ++fc) {  in DMComputeL2FieldDiff_Plex()
1826             const PetscReal wt = quadWeights[q * qNc + (qNc == 1 ? 0 : qc + fc)];  in DMComputeL2FieldDiff_Plex()
1828 …tscInt_FMT ",%" PetscInt_FMT " point %g %g %g diff %g\n", cell, fields[f], fc, (double)(dE > 0 ? c…  in DMComputeL2FieldDiff_Plex()
1829 …          (double)(PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q])));  in DMComputeL2FieldDiff_Plex()
1830 …          elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q];  in DMComputeL2FieldDiff_Plex()
1930       PetscInt     Nb, Nc, q, fc;  in DMPlexComputeL2DiffVec()  local
1961           for (fc = 0; fc < Nc; ++fc) {  in DMPlexComputeL2DiffVec()
1962             const PetscReal wt = quadWeights[q * qNc + (qNc == 1 ? 0 : qc + fc)];  in DMPlexComputeL2DiffVec()
1963 …          elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q];  in DMPlexComputeL2DiffVec()
2062       for (PetscInt fc = 0; fc < mNc; ++fc) {  in DMPlexComputeL2FluxDiffVecLocal()  local
2063         const PetscReal wt = quadWeights[q * qNc + (qNc == 1 ? 0 : fc)];  in DMPlexComputeL2FluxDiffVecLocal()
2064 … fluxElemDiff += PetscSqr(PetscRealPart(interpolant[fc] - minterpolant[fc])) * wt * fegeom.detJ[q];  in DMPlexComputeL2FluxDiffVecLocal()
2194     PetscInt    starSize, st, fc;  in DMPlexComputeClementInterpolant()  local
2234           for (fc = 0; fc < fNc; ++fc) val[foff + fc] += interpolant[fc] * wt;  in DMPlexComputeClementInterpolant()
2240       for (fc = 0; fc < Nc; ++fc) valsum[fc] += val[fc];  in DMPlexComputeClementInterpolant()
2244         for (fc = 0; fc < Nc; ++fc) {  in DMPlexComputeClementInterpolant()
2245           if (fc) PetscCall(PetscPrintf(PETSC_COMM_SELF, ", "));  in DMPlexComputeClementInterpolant()
2246           PetscCall(PetscPrintf(PETSC_COMM_SELF, "%g", (double)PetscRealPart(val[fc])));  in DMPlexComputeClementInterpolant()
2251     for (fc = 0; fc < Nc; ++fc) valsum[fc] /= volsum;  in DMPlexComputeClementInterpolant()
2327     PetscInt    starSize, st, d, fc;  in DMPlexComputeGradientClementInterpolant()  local
2365           for (fc = 0; fc < Nc; ++fc) {  in DMPlexComputeGradientClementInterpolant()
2368 …for (d = 0; d < coordDim; ++d) grad[fc * coordDim + d] += interpolant[fc * dim + d] * wt * fegeom.…  in DMPlexComputeGradientClementInterpolant()
2376       for (fc = 0; fc < numComponents; ++fc) {  in DMPlexComputeGradientClementInterpolant()
2377         for (d = 0; d < coordDim; ++d) gradsum[fc * coordDim + d] += grad[fc * coordDim + d];  in DMPlexComputeGradientClementInterpolant()
2382         for (fc = 0; fc < numComponents; ++fc) {  in DMPlexComputeGradientClementInterpolant()
2384             if (fc || d > 0) PetscCall(PetscPrintf(PETSC_COMM_SELF, ", "));  in DMPlexComputeGradientClementInterpolant()
2385 …     PetscCall(PetscPrintf(PETSC_COMM_SELF, "%g", (double)PetscRealPart(grad[fc * coordDim + d])));  in DMPlexComputeGradientClementInterpolant()
2391     for (fc = 0; fc < numComponents; ++fc) {  in DMPlexComputeGradientClementInterpolant()
2392       for (d = 0; d < coordDim; ++d) gradsum[fc * coordDim + d] /= volsum;  in DMPlexComputeGradientClementInterpolant()
4758         PetscInt eOffset = 0, Nc, fc, foff;  in DMPlexComputeJacobian_Patch_Internal()  local
4764           for (fc = 0; fc < Nc; ++fc) {  in DMPlexComputeJacobian_Patch_Internal()
4765             const PetscInt i = foff + fc;  in DMPlexComputeJacobian_Patch_Internal()
6310     PetscInt     offsetI, NcI, NbI = 1, fc, f;  in DMPlexComputeJacobianByKey()  local
6325         for (fc = 0; fc < NcI; ++fc) {  in DMPlexComputeJacobianByKey()
6327             const PetscInt i = offsetI + f * NcI + fc;  in DMPlexComputeJacobianByKey()
6548         for (PetscInt fc = 0; fc < NcI; ++fc) {  in DMPlexComputeJacobianByKeyGeneral()  local
6550             const PetscInt i = offsetI + f * NcI + fc;  in DMPlexComputeJacobianByKeyGeneral()