127     PetscReal *phi, *psi;  in PetscSpaceEvaluate_Subspace()  local
130     PetscCall(DMGetWorkArray(sp->dm, origNc * subDim * subDim, MPIU_REAL, &psi));  in PetscSpaceEvaluate_Subspace()
137         for (k = 0; k < origNc * subDim; k++) psi[k] = 0.;  in PetscSpaceEvaluate_Subspace()
146 …for (o = 0; o < origDim; o++) psi[(k * subDim + l) * subDim + n] += subsp->Jx[m * subDim + l] * su…  in PetscSpaceEvaluate_Subspace()
154 …for (m = 0; m < PetscMin(subDim, origDim); m++) psi[(k * subDim + l) * subDim + m] += phi[(k * ori…  in PetscSpaceEvaluate_Subspace()
161 …Nb + i) * subNc + k) * subDim * subDim + m] += subsp->Ju[k * origNc + l] * psi[l * subDim * subDim…  in PetscSpaceEvaluate_Subspace()
166 …m * subDim; l++) H[((j * subNb + i) * subNc + k) * subDim * subDim + l] += psi[k * subDim * subDim…  in PetscSpaceEvaluate_Subspace()
171     PetscCall(DMRestoreWorkArray(sp->dm, subNc * origDim, MPIU_REAL, &psi));  in PetscSpaceEvaluate_Subspace()
176     PetscReal *phi, *psi;  in PetscSpaceEvaluate_Subspace()  local
179     PetscCall(DMGetWorkArray(sp->dm, origNc * subDim, MPIU_REAL, &psi));  in PetscSpaceEvaluate_Subspace()
186         for (k = 0; k < origNc * subDim; k++) psi[k] = 0.;  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()
204 …) D[((j * subNb + i) * subNc + k) * subDim + m] += subsp->Ju[k * origNc + l] * psi[l * subDim + m];  in PetscSpaceEvaluate_Subspace()
209 … for (l = 0; l < subDim; l++) D[((j * subNb + i) * subNc + k) * subDim + l] += psi[k * subDim + l];  in PetscSpaceEvaluate_Subspace()
214     PetscCall(DMRestoreWorkArray(sp->dm, subNc * origDim, MPIU_REAL, &psi));  in PetscSpaceEvaluate_Subspace()