220       std::vector<moab::EntityHandle> connp;  in DMCreateInterpolation_Moab()  local
221       merr = dmbp->hierarchy->get_connectivity(parent, dmbp->hlevel, connp);  in DMCreateInterpolation_Moab()
224       for (unsigned ic = 0; ic < connp.size(); ++ic) {  in DMCreateInterpolation_Moab()
227 …if (found.find(connp[ic]) != found.end()) continue;                    /* make sure we don't doubl…  in DMCreateInterpolation_Moab()
228 …if (dmbp->vghost->find(connp[ic]) != dmbp->vghost->end()) onz[ldof]++; /* update out-of-proc onz */  in DMCreateInterpolation_Moab()
230         found.insert(connp[ic]);  in DMCreateInterpolation_Moab()
267   std::vector<moab::EntityHandle> connp, connc;  in DMCreateInterpolation_Moab()  local
278     merr = dmbp->hierarchy->get_connectivity(ehandle, dmbp->hlevel, connp);  in DMCreateInterpolation_Moab()
284     pcoords.resize(connp.size() * 3);  in DMCreateInterpolation_Moab()
286     values_phi.resize(connp.size() * connc.size());  in DMCreateInterpolation_Moab()
288     merr = dmbp->hierarchy->get_coordinates(&connp[0], connp.size(), dmbp->hlevel, &pcoords[0]);  in DMCreateInterpolation_Moab()
298 …PetscCall(DMMoabPToRMapping(dim, connp.size(), &pcoords[0], &ccoords[offset], &natparam[offset], &…  in DMCreateInterpolation_Moab()
319     merr = dmbp->hierarchy->get_connectivity(ehandle, dmbp->hlevel, connp);  in DMCreateInterpolation_Moab()
324     pcoords.resize(connp.size() * 3);  in DMCreateInterpolation_Moab()
327     merr = dmbp->hierarchy->get_coordinates(&connp[0], connp.size(), dmbp->hlevel, &pcoords[0]);  in DMCreateInterpolation_Moab()
332     std::vector<int> dofsp(connp.size()), dofsc(connc.size());  in DMCreateInterpolation_Moab()
334     PetscCall(DMMoabGetDofsBlocked(dmp, connp.size(), &connp[0], &dofsp[0]));  in DMCreateInterpolation_Moab()
349 …PetscCall(MatSetValues(*interpl, 1, &dofsc[tc], connp.size(), &dofsp[0], &values_phi_scalar[connp.…  in DMCreateInterpolation_Moab()
360       values_phi.resize(connp.size());  in DMCreateInterpolation_Moab()
364         for (unsigned tp = 0; tp < connp.size(); tp++) {  in DMCreateInterpolation_Moab()
375         for (unsigned tp = 0; tp < connp.size(); tp++) {  in DMCreateInterpolation_Moab()
376           if (values_phi[tp] > 1e11) values_phi[tp] = factor * 0.5 / connp.size();  in DMCreateInterpolation_Moab()
377           else values_phi[tp] = factor * values_phi[tp] * 0.5 / (connp.size() * normsum);  in DMCreateInterpolation_Moab()
380 …PetscCall(MatSetValues(*interpl, 1, &dofsc[tc], connp.size(), &dofsp[0], &values_phi_scalar2[0], A…  in DMCreateInterpolation_Moab()