51     PetscInt origDim, subDim, origNc, subNc, o, s;  in PetscSpaceView_Subspace()  local
56     PetscCall(PetscSpaceGetNumComponents(sp, &subNc));  in PetscSpaceView_Subspace()
79 …for (s = 0; s < subNc; s++) PetscCall(PetscViewerASCIIPrintf(viewer, " %g", (double)subsp->Ju[o * …  in PetscSpaceView_Subspace()
96   PetscInt             origDim, subDim, origNc, subNc, subNb, origNb, i, j, k, l, m, n, o;  in PetscSpaceEvaluate_Subspace()  local
103   PetscCall(PetscSpaceGetNumComponents(sp, &subNc));  in PetscSpaceEvaluate_Subspace()
131     for (i = 0; i < npoints * subNb * subNc * subDim; i++) D[i] = 0.0;  in PetscSpaceEvaluate_Subspace()
142           for (k = 0; k < subNc; k++) {  in PetscSpaceEvaluate_Subspace()
152           for (k = 0; k < subNc; k++) {  in PetscSpaceEvaluate_Subspace()
159           for (k = 0; k < subNc; k++) {  in PetscSpaceEvaluate_Subspace()
161 …for (m = 0; m < subDim * subDim; m++) H[((j * subNb + i) * subNc + k) * subDim * subDim + m] += su…  in PetscSpaceEvaluate_Subspace()
165           for (k = 0; k < PetscMin(subNc, origNc); k++) {  in PetscSpaceEvaluate_Subspace()
166 …for (l = 0; l < subDim * subDim; l++) H[((j * subNb + i) * subNc + k) * subDim * subDim + l] += ps…  in PetscSpaceEvaluate_Subspace()
171     PetscCall(DMRestoreWorkArray(sp->dm, subNc * origDim, MPIU_REAL, &psi));  in PetscSpaceEvaluate_Subspace()
180     for (i = 0; i < npoints * subNb * subNc * subDim; i++) D[i] = 0.0;  in PetscSpaceEvaluate_Subspace()
191           for (k = 0; k < subNc; k++) {  in PetscSpaceEvaluate_Subspace()
197           for (k = 0; k < subNc; k++) {  in PetscSpaceEvaluate_Subspace()
202           for (k = 0; k < subNc; k++) {  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()
208           for (k = 0; k < PetscMin(subNc, origNc); k++) {  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()
224         for (j = 0; j < subNc; j++) B[i * subNc + j] = subsp->u[j];  in PetscSpaceEvaluate_Subspace()
227       for (i = 0; i < npoints * subNb * subNc; i++) B[i] = 0.0;  in PetscSpaceEvaluate_Subspace()
238           for (k = 0; k < subNc; 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()
281   PetscInt             origDim, subDim, origNc, subNc, origNb, subNb, f, i, j, numPoints, offset;  in PetscSpaceSetUp_Subspace()  local
301   PetscCall(PetscDualSpaceGetNumComponents(dualSubspace, &subNc));  in PetscSpaceSetUp_Subspace()
337           for (j = 0; j < subNc; j++) allWeights[offset * origNc + i] += qw[j] * Ju[j * origNc + i];  in PetscSpaceSetUp_Subspace()
340         for (i = 0; i < PetscMin(subNc, origNc); i++) allWeights[offset * origNc + i] += qw[i];  in PetscSpaceSetUp_Subspace()
354         for (s = 0; s < subNc; s++) V[f * origNb + b] += qw[s] * u[s];  in PetscSpaceSetUp_Subspace()
450   PetscInt             origDim, subDim, origNc, subNc, subNb;  in PetscSpaceCreateSubspace()  local
467   PetscCall(PetscDualSpaceGetNumComponents(dualSubspace, &subNc));  in PetscSpaceCreateSubspace()
471   PetscCall(PetscSpaceSetNumComponents(*subspace, subNc));  in PetscSpaceCreateSubspace()
501       PetscCall(PetscMalloc1(subNc, &subsp->u_alloc));  in PetscSpaceCreateSubspace()
502       PetscCall(PetscArraycpy(subsp->u_alloc, u, subNc));  in PetscSpaceCreateSubspace()
506       PetscCall(PetscMalloc1(origNc * subNc, &subsp->Ju_alloc));  in PetscSpaceCreateSubspace()
507       PetscCall(PetscArraycpy(subsp->Ju_alloc, Ju, origNc * subNc));  in PetscSpaceCreateSubspace()