108   Mat       M_p;  in createMp()  local
129   PetscCall(DMCreateMassMatrix(sw, dm, &M_p));  in createMp()
132   *Mp_out = M_p;  in createMp()
136 …t DM dm, DM sw, const PetscInt a_tid, const PetscInt dim, const PetscReal a_wp[], Vec rho, Mat M_p)  in particlesToGrid()  argument
151   PetscCall(MatMultTranspose(M_p, ff, rho));  in particlesToGrid()
159 PetscErrorCode gridToParticles(const DM dm, DM sw, const Vec rhs, Vec work_ferhs, Mat M_p, Mat Mass)  in gridToParticles()  argument
181     PetscCall(MatGetLocalSize(M_p, &M, &N));  in gridToParticles()
189       PetscCall(MatCreateVecs(M_p, &matshellctx->uu, &matshellctx->ff));  in gridToParticles()
191         PetscCall(MatTransposeMatMult(M_p, M_p, MAT_INITIAL_MATRIX, 4, &MtM));  in gridToParticles()
193         PetscCall(PetscInfo(M_p, "createMtM KSP with explicit Mp'Mp\n"));  in gridToParticles()
197         PetscCall(MatTranspose(M_p, MAT_INITIAL_MATRIX, &matshellctx->MpTrans));  in gridToParticles()
198         matshellctx->Mp = M_p;  in gridToParticles()
229 …PetscCall(PetscInfo(M_p, "createMtMKSP Have %" PetscInt_FMT " eqs, nzl = %" PetscInt_FMT "\n", N, …  in gridToParticles()
232           PetscCall(MatViewFromOptions(M_p, NULL, "-ftop2_Mp_mat_view"));  in gridToParticles()
246       PetscCall(KSPSetOperators(ksp, M_p, PM_p));  in gridToParticles()
248       PetscCall(KSPSetOperators(ksp, M_p, M_p));  in gridToParticles()
250     PetscCall(MatViewFromOptions(M_p, NULL, "-ftop2_Mp_mat_view"));  in gridToParticles()
256     PetscCall(MatMult(M_p, matshellctx->uu, ff));  in gridToParticles()