178   DM           dmP;          // The DM for momentum density fields  member
272   options->dmP                    = NULL;  in ProcessOptions()
496     PetscCall(DMGetNamedGlobalVector(ctx->dmP, "p", &p));  in SetupContext()
498     PetscCall(DMRestoreNamedGlobalVector(ctx->dmP, "p", &p));  in SetupContext()
534     PetscCall(DMGetNamedGlobalVector(ctx->dmP, "pres", &pres));  in SetupContext()
536     PetscCall(DMRestoreNamedGlobalVector(ctx->dmP, "pres", &pres));  in SetupContext()
949   PetscCall(DMGetGlobalVector(ctx->dmP, &prhs));  in ComputeMomentFields()
950   PetscCall(DMGetGlobalVector(ctx->dmP, &pflux));  in ComputeMomentFields()
952   PetscCall(DMGetNamedGlobalVector(ctx->dmP, "p", &p));  in ComputeMomentFields()
1002   PetscCall(DMRestoreGlobalVector(ctx->dmP, &prhs));  in ComputeMomentFields()
1026   PetscCall(DMGetNamedGlobalVector(ctx->dmP, "pres", &pres));  in ComputeMomentFields()
1038   PetscCall(DMRestoreNamedGlobalVector(ctx->dmP, "pres", &pres));  in ComputeMomentFields()
1062   PetscCall(DMRestoreGlobalVector(ctx->dmP, &pflux));  in ComputeMomentFields()
1065   PetscCall(DMRestoreNamedGlobalVector(ctx->dmP, "p", &p));  in ComputeMomentFields()
1085   PetscCall(DMGetNamedGlobalVector(ctx->dmP, "p", &p));  in MonitorMomentFields()
1087   PetscCall(DMRestoreNamedGlobalVector(ctx->dmP, "p", &p));  in MonitorMomentFields()
1097   PetscCall(DMGetNamedGlobalVector(ctx->dmP, "pres", &pres));  in MonitorMomentFields()
1099   PetscCall(DMRestoreNamedGlobalVector(ctx->dmP, "pres", &pres));  in MonitorMomentFields()
1692   PetscCall(DMCreateSubDM(ctx->dmMom, 1, &field, &ctx->isP, &ctx->dmP));  in CreateMomentFields()
1693   PetscCall(DMCreateMassMatrix(ctx->dmP, ctx->dmP, &ctx->MP));  in CreateMomentFields()
3061   PetscCall(DMDestroy(&ctx.dmP));  in main()