308   PetscReal   fmoments[3]; // \int \hat f, \int x \hat f, \int r^2 \hat f  in TestParticlesToField()  local
317   PetscCall(computeFieldMoments(dm, fhat, fmoments, user));  in TestParticlesToField()
319 …n mass: %20.10e, x-momentum: %20.10e, energy: %20.10e.\n", fmoments[0], fmoments[1], fmoments[2]));  in TestParticlesToField()
322 …if (PetscAbsReal((fmoments[m] - pmoments[m]) / fmoments[m]) > user->mtol) PetscCall(PetscPrintf(Pe…  in TestParticlesToField()
324 …fmoments[m] - pmoments[m]) / fmoments[m]) <= user->mtol, PetscObjectComm((PetscObject)dm), PETSC_E…  in TestParticlesToField()
336   PetscReal   fmoments[3]; // \int \hat f, \int x \hat f, \int r^2 \hat f  in TestFieldToParticles()  local
344   PetscCall(computeFieldMoments(dm, fhat, fmoments, user));  in TestFieldToParticles()
345 …n mass: %20.10e, x-momentum: %20.10e, energy: %20.10e.\n", fmoments[0], fmoments[1], fmoments[2]));  in TestFieldToParticles()
348 …if (PetscAbsReal((fmoments[m] - pmoments[m]) / fmoments[m]) > user->mtol) PetscCall(PetscPrintf(Pe…  in TestFieldToParticles()
350 …fmoments[m] - pmoments[m]) / fmoments[m]) <= user->mtol, PetscObjectComm((PetscObject)dm), PETSC_E…  in TestFieldToParticles()
362   PetscReal   fmoments[3]; // \int \hat f, \int x \hat f, \int r^2 \hat f  in TestParticlesToGradientField()  local
371   PetscCall(computeFieldMoments(dm, fhat, fmoments, user));  in TestParticlesToGradientField()