1509     PetscReal    elemDiff = 0.0;  in DMPlexComputeL2DiffLocal()  local
1566           elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q];  in DMPlexComputeL2DiffLocal()
1573 …PetscCall(PetscPrintf(PETSC_COMM_SELF, "  elem %" PetscInt_FMT " diff %g\n", c, (double)elemDiff));  in DMPlexComputeL2DiffLocal()
1574     localDiff += elemDiff;  in DMPlexComputeL2DiffLocal()
1625     PetscReal    elemDiff = 0.0;  in DMComputeL2GradientDiff_Plex()  local
1680           elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q];  in DMComputeL2GradientDiff_Plex()
1687 …PetscCall(PetscPrintf(PETSC_COMM_SELF, "  elem %" PetscInt_FMT " diff %g\n", c, (double)elemDiff));  in DMComputeL2GradientDiff_Plex()
1688     localDiff += elemDiff;  in DMComputeL2GradientDiff_Plex()
1775         PetscReal    elemDiff = 0.0;  in DMComputeL2FieldDiff_Plex()  local
1830 …          elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q];  in DMComputeL2FieldDiff_Plex()
1835         localDiff[fields[f]] += elemDiff;  in DMComputeL2FieldDiff_Plex()
1920     PetscScalar  elemDiff = 0.0;  in DMPlexComputeL2DiffVec()  local
1963 …          elemDiff += PetscSqr(PetscRealPart(interpolant[fc] - funcVal[fc])) * wt * fegeom.detJ[q];  in DMPlexComputeL2DiffVec()
1971     PetscCall(VecSetValue(D, c - cStart, elemDiff, INSERT_VALUES));  in DMPlexComputeL2DiffVec()