408     PetscInt     li, nxi, lj, neta;  in DMCreateInterpolation_DA_2D_Q1()  local
416 …  for (lj = 0; lj < neta; lj++) eta[lj] = -1.0 + (PetscScalar)lj * (2.0 / (PetscScalar)(neta - 1));  in DMCreateInterpolation_DA_2D_Q1()
430         lj = j - ratioj * (j / ratioj);  in DMCreateInterpolation_DA_2D_Q1()
431         if (j == my - 1) lj = neta - 1;  in DMCreateInterpolation_DA_2D_Q1()
438           if ((lj == 0) || (lj == neta - 1)) {  in DMCreateInterpolation_DA_2D_Q1()
455         Ni[0] = 0.25 * (1.0 - xi[li]) * (1.0 - eta[lj]);  in DMCreateInterpolation_DA_2D_Q1()
456         Ni[1] = 0.25 * (1.0 + xi[li]) * (1.0 - eta[lj]);  in DMCreateInterpolation_DA_2D_Q1()
457         Ni[2] = 0.25 * (1.0 - xi[li]) * (1.0 + eta[lj]);  in DMCreateInterpolation_DA_2D_Q1()
458         Ni[3] = 0.25 * (1.0 + xi[li]) * (1.0 + eta[lj]);  in DMCreateInterpolation_DA_2D_Q1()
914     PetscInt     li, nxi, lj, neta, lk, nzeta, n;  in DMCreateInterpolation_DA_3D_Q1()  local
925 …  for (lj = 0; lj < neta; lj++) eta[lj] = -1.0 + (PetscScalar)lj * (2.0 / (PetscScalar)(neta - 1));  in DMCreateInterpolation_DA_3D_Q1()
941           lj = j - ratioj * (j / ratioj);  in DMCreateInterpolation_DA_3D_Q1()
942           if (j == my - 1) lj = neta - 1;  in DMCreateInterpolation_DA_3D_Q1()
951             if ((lj == 0) || (lj == neta - 1)) {  in DMCreateInterpolation_DA_3D_Q1()
976           Ni[0] = 0.125 * (1.0 - xi[li]) * (1.0 - eta[lj]) * (1.0 - zeta[lk]);  in DMCreateInterpolation_DA_3D_Q1()
977           Ni[1] = 0.125 * (1.0 + xi[li]) * (1.0 - eta[lj]) * (1.0 - zeta[lk]);  in DMCreateInterpolation_DA_3D_Q1()
978           Ni[2] = 0.125 * (1.0 - xi[li]) * (1.0 + eta[lj]) * (1.0 - zeta[lk]);  in DMCreateInterpolation_DA_3D_Q1()
979           Ni[3] = 0.125 * (1.0 + xi[li]) * (1.0 + eta[lj]) * (1.0 - zeta[lk]);  in DMCreateInterpolation_DA_3D_Q1()
981           Ni[4] = 0.125 * (1.0 - xi[li]) * (1.0 - eta[lj]) * (1.0 + zeta[lk]);  in DMCreateInterpolation_DA_3D_Q1()
982           Ni[5] = 0.125 * (1.0 + xi[li]) * (1.0 - eta[lj]) * (1.0 + zeta[lk]);  in DMCreateInterpolation_DA_3D_Q1()
983           Ni[6] = 0.125 * (1.0 - xi[li]) * (1.0 + eta[lj]) * (1.0 + zeta[lk]);  in DMCreateInterpolation_DA_3D_Q1()
984           Ni[7] = 0.125 * (1.0 + xi[li]) * (1.0 + eta[lj]) * (1.0 + zeta[lk]);  in DMCreateInterpolation_DA_3D_Q1()