607   PetscInt               i, j, l, i_start, j_start, l_start, m_f, n_f, p_f, Mx, My, Mz, dof;  in DMCreateInterpolation_DA_3D_Q0()  local
635   PetscCall(DMDAGetCorners(daf, &i_start, &j_start, &l_start, &m_f, &n_f, &p_f));  in DMCreateInterpolation_DA_3D_Q0()
661 …MatPreallocateBegin(PetscObjectComm((PetscObject)daf), m_f * n_f * p_f, col_scale * m_c * n_c * p_…  in DMCreateInterpolation_DA_3D_Q0()
662   for (l = l_start; l < l_start + p_f; l++) {  in DMCreateInterpolation_DA_3D_Q0()
697 …PetscCall(MatSetSizes(mat, m_f * n_f * p_f, col_scale * m_c * n_c * p_c, mx * my * mz, col_scale *…  in DMCreateInterpolation_DA_3D_Q0()
705   for (l = l_start; l < l_start + p_f; l++) {  in DMCreateInterpolation_DA_3D_Q0()
730   PetscCall(PetscLogFlops(13.0 * m_f * n_f * p_f));  in DMCreateInterpolation_DA_3D_Q0()
742   PetscInt               l_start, p_f, l_start_ghost, p_ghost, l_start_c, p_c;  in DMCreateInterpolation_DA_3D_Q1()  local
786   PetscCall(DMDAGetCorners(daf, &i_start, &j_start, &l_start, &m_f, &n_f, &p_f));  in DMCreateInterpolation_DA_3D_Q1()
797 …MatPreallocateBegin(PetscObjectComm((PetscObject)dac), m_f * n_f * p_f, m_c * n_c * p_c, dnz, onz);  in DMCreateInterpolation_DA_3D_Q1()
799   for (l = l_start; l < l_start + p_f; l++) {  in DMCreateInterpolation_DA_3D_Q1()
838   PetscCall(MatSetSizes(mat, m_f * n_f * p_f, m_c * n_c * p_c, mx * my * mz, Mx * My * Mz));  in DMCreateInterpolation_DA_3D_Q1()
847     for (l = l_start; l < l_start + p_f; l++) {  in DMCreateInterpolation_DA_3D_Q1()
928     for (l = l_start; l < l_start + p_f; l++) {  in DMCreateInterpolation_DA_3D_Q1()
1176   PetscInt               i, j, k, i_start, j_start, k_start, m_f, n_f, p_f, Mx, My, Mz;  in DMCreateInjection_DA_3D()  local
1217   PetscCall(DMDAGetCorners(daf, &i_start, &j_start, &k_start, &m_f, &n_f, &p_f));  in DMCreateInjection_DA_3D()
1229   PetscCall(PetscMalloc1(n_f * m_f * p_f, &cols));  in DMCreateInjection_DA_3D()
1339   PetscInt               i_start, j_start, l_start, m_f, n_f, p_f;  in DMDACreateAggregates()  local
1376   PetscCall(DMDAGetCorners(daf, &i_start, &j_start, &l_start, &m_f, &n_f, &p_f));  in DMDACreateAggregates()
1399 …PetscCall(MatCreateAIJ(PetscObjectComm((PetscObject)daf), m_c * n_c * p_c * dofc, m_f * n_f * p_f …  in DMDACreateAggregates()