239   for (int tc = 0; tc < nlsizc; tc++) {  in DMCreateInterpolation_Moab()  local
241     nnz[tc] = std::min(nlsizp, nnz[tc]);  in DMCreateInterpolation_Moab()
242     onz[tc] = std::min(ngsizp - nlsizp, onz[tc]);  in DMCreateInterpolation_Moab()
244     PetscCall(PetscInfo(NULL, "  %d: NNZ = %d, ONZ = %d\n", tc, nnz[tc], onz[tc]));  in DMCreateInterpolation_Moab()
246     innz = (innz < nnz[tc] ? nnz[tc] : innz);  in DMCreateInterpolation_Moab()
247     ionz = (ionz < onz[tc] ? onz[tc] : ionz);  in DMCreateInterpolation_Moab()
294     for (unsigned tc = 0; tc < connc.size(); tc++) {  in DMCreateInterpolation_Moab()  local
295       const PetscInt offset = tc * 3;  in DMCreateInterpolation_Moab()
298 …, connp.size(), &pcoords[0], &ccoords[offset], &natparam[offset], &values_phi[connp.size() * tc]));  in DMCreateInterpolation_Moab()
347       for (unsigned tc = 0; tc < connc.size(); tc++) {  in DMCreateInterpolation_Moab()  local
349 …PetscCall(MatSetValues(*interpl, 1, &dofsc[tc], connp.size(), &dofsp[0], &values_phi_scalar[connp.…  in DMCreateInterpolation_Moab()
361       for (unsigned tc = 0; tc < connc.size(); tc++) {  in DMCreateInterpolation_Moab()  local
366 …gned k = 0; k < 3; k++) values_phi[tp] += std::pow(pcoords[tp * 3 + k] - ccoords[k + tc * 3], dim);  in DMCreateInterpolation_Moab()
380 …PetscCall(MatSetValues(*interpl, 1, &dofsc[tc], connp.size(), &dofsp[0], &values_phi_scalar2[0], A…  in DMCreateInterpolation_Moab()