127     PetscReal *phi, *psi;  in PetscSpaceEvaluate_Subspace()  local
129     PetscCall(DMGetWorkArray(sp->dm, origNc * origDim * origDim, MPIU_REAL, &phi));  in PetscSpaceEvaluate_Subspace()
136         for (k = 0; k < origNc * origDim; k++) phi[k] = 0.;  in PetscSpaceEvaluate_Subspace()
139 …for (l = 0; l < origNc * origDim * origDim; l++) phi[l] += inH[(j * origNb + k) * origNc * origDim…  in PetscSpaceEvaluate_Subspace()
146 …) * subDim + n] += subsp->Jx[m * subDim + l] * subsp->Jx[o * subDim + n] * phi[(k * origDim + m) *…  in PetscSpaceEvaluate_Subspace()
154 … m < PetscMin(subDim, origDim); m++) psi[(k * subDim + l) * subDim + m] += phi[(k * origDim + l) *…  in PetscSpaceEvaluate_Subspace()
172     PetscCall(DMRestoreWorkArray(sp->dm, origNc * origDim, MPIU_REAL, &phi));  in PetscSpaceEvaluate_Subspace()
176     PetscReal *phi, *psi;  in PetscSpaceEvaluate_Subspace()  local
178     PetscCall(DMGetWorkArray(sp->dm, origNc * origDim, MPIU_REAL, &phi));  in PetscSpaceEvaluate_Subspace()
185         for (k = 0; k < origNc * origDim; k++) phi[k] = 0.;  in PetscSpaceEvaluate_Subspace()
188 …for (l = 0; l < origNc * origDim; l++) phi[l] += inD[(j * origNb + k) * origNc * origDim + l] * su…  in PetscSpaceEvaluate_Subspace()
193 …for (m = 0; m < origDim; m++) psi[k * subDim + l] += subsp->Jx[m * subDim + l] * phi[k * origDim +…  in PetscSpaceEvaluate_Subspace()
198 …       for (l = 0; l < PetscMin(subDim, origDim); l++) psi[k * subDim + l] += phi[k * origDim + l];  in PetscSpaceEvaluate_Subspace()
215     PetscCall(DMRestoreWorkArray(sp->dm, origNc * origDim, MPIU_REAL, &phi));  in PetscSpaceEvaluate_Subspace()
219     PetscReal *phi;  in PetscSpaceEvaluate_Subspace()  local
221     PetscCall(DMGetWorkArray(sp->dm, origNc, MPIU_REAL, &phi));  in PetscSpaceEvaluate_Subspace()
233         for (k = 0; k < origNc; k++) phi[k] = 0.;  in PetscSpaceEvaluate_Subspace()
235           for (l = 0; l < origNc; l++) phi[l] += inB[(j * origNb + k) * origNc + l] * subq[k];  in PetscSpaceEvaluate_Subspace()
239 … for (l = 0; l < origNc; l++) B[(j * subNb + i) * subNc + k] += subsp->Ju[k * origNc + l] * phi[l];  in PetscSpaceEvaluate_Subspace()
242           for (k = 0; k < PetscMin(subNc, origNc); k++) B[(j * subNb + i) * subNc + k] += phi[k];  in PetscSpaceEvaluate_Subspace()
246     PetscCall(DMRestoreWorkArray(sp->dm, origNc, MPIU_REAL, &phi));  in PetscSpaceEvaluate_Subspace()