96   PetscInt             origDim, subDim, origNc, subNc, subNb, origNb, i, j, k, l, m, n, o;  in PetscSpaceEvaluate_Subspace()  local
139 …for (l = 0; l < origNc * origDim * origDim; l++) phi[l] += inH[(j * origNb + k) * origNc * origDim…  in PetscSpaceEvaluate_Subspace()
143             for (l = 0; l < subDim; l++) {  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()
153             for (l = 0; l < PetscMin(subDim, origDim); l++) {  in PetscSpaceEvaluate_Subspace()
154 …; m < PetscMin(subDim, origDim); m++) psi[(k * subDim + l) * subDim + m] += phi[(k * origDim + l) …  in PetscSpaceEvaluate_Subspace()
160             for (l = 0; l < origNc; l++) {  in PetscSpaceEvaluate_Subspace()
161 …bNb + i) * subNc + k) * subDim * subDim + m] += subsp->Ju[k * origNc + l] * psi[l * subDim * subDi…  in PetscSpaceEvaluate_Subspace()
166 …for (l = 0; l < subDim * subDim; l++) H[((j * subNb + i) * subNc + k) * subDim * subDim + l] += ps…  in PetscSpaceEvaluate_Subspace()
188 …for (l = 0; l < origNc * origDim; l++) phi[l] += inD[(j * origNb + k) * origNc * origDim + l] * su…  in PetscSpaceEvaluate_Subspace()
192             for (l = 0; l < subDim; l++) {  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()
203             for (l = 0; l < origNc; 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()
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()