96   PetscInt             origDim, subDim, origNc, subNc, subNb, origNb, i, j, k, l, m, n, o;  in PetscSpaceEvaluate_Subspace()  local
105   PetscCall(PetscSpaceGetDimension(sp, &subNb));  in PetscSpaceEvaluate_Subspace()
131     for (i = 0; i < npoints * subNb * subNc * subDim; i++) D[i] = 0.0;  in PetscSpaceEvaluate_Subspace()
132     for (i = 0; i < subNb; i++) {  in PetscSpaceEvaluate_Subspace()
161 …for (m = 0; m < subDim * subDim; m++) H[((j * subNb + i) * subNc + k) * subDim * subDim + m] += su…  in PetscSpaceEvaluate_Subspace()
166 …for (l = 0; l < subDim * subDim; l++) H[((j * subNb + i) * subNc + k) * subDim * subDim + l] += ps…  in PetscSpaceEvaluate_Subspace()
180     for (i = 0; i < npoints * subNb * subNc * subDim; i++) D[i] = 0.0;  in PetscSpaceEvaluate_Subspace()
181     for (i = 0; i < subNb; i++) {  in PetscSpaceEvaluate_Subspace()
204 …for (m = 0; m < subDim; m++) D[((j * subNb + i) * subNc + k) * subDim + m] += subsp->Ju[k * origNc…  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()
223       for (i = 0; i < npoints * subNb; i++) {  in PetscSpaceEvaluate_Subspace()
227       for (i = 0; i < npoints * subNb * subNc; i++) B[i] = 0.0;  in PetscSpaceEvaluate_Subspace()
229     for (i = 0; i < subNb; i++) {  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()
281   PetscInt             origDim, subDim, origNc, subNc, origNb, subNb, f, i, j, numPoints, offset;  in PetscSpaceSetUp_Subspace()  local
300   PetscCall(PetscDualSpaceGetDimension(dualSubspace, &subNb));  in PetscSpaceSetUp_Subspace()
303   for (f = 0, numPoints = 0; f < subNb; f++) {  in PetscSpaceSetUp_Subspace()
311   PetscCall(PetscMalloc1(subNb * origNb, &V));  in PetscSpaceSetUp_Subspace()
313   for (f = 0, offset = 0; f < subNb; f++) {  in PetscSpaceSetUp_Subspace()
345   for (f = 0, offset = 0; f < subNb; f++) {  in PetscSpaceSetUp_Subspace()
366   for (f = 0; f < subNb; f++) {  in PetscSpaceSetUp_Subspace()
374     for (j = f + 1; j < subNb; j++) {  in PetscSpaceSetUp_Subspace()
450   PetscInt             origDim, subDim, origNc, subNc, subNb;  in PetscSpaceCreateSubspace()  local
466   PetscCall(PetscDualSpaceGetDimension(dualSubspace, &subNb));  in PetscSpaceCreateSubspace()
475   subsp->Nb = subNb;  in PetscSpaceCreateSubspace()