51     PetscInt origDim, subDim, origNc, subNc, o, s;  in PetscSpaceView_Subspace()  local
54     PetscCall(PetscSpaceGetNumComponents(subsp->origSpace, &origNc));  in PetscSpaceView_Subspace()
73 …for (o = 0; o < origNc; o++) PetscCall(PetscViewerASCIIPrintf(viewer, " %g\n", (double)subsp->u[o]…  in PetscSpaceView_Subspace()
77       for (o = 0; o < origNc; o++) {  in PetscSpaceView_Subspace()
96   PetscInt             origDim, subDim, origNc, subNc, subNb, origNb, i, j, k, l, m, n, o;  in PetscSpaceEvaluate_Subspace()  local
104   PetscCall(PetscSpaceGetNumComponents(origsp, &origNc));  in PetscSpaceEvaluate_Subspace()
122   if (B) PetscCall(DMGetWorkArray(sp->dm, npoints * origNb * origNc, MPIU_REAL, &inB));  in PetscSpaceEvaluate_Subspace()
123   if (D) PetscCall(DMGetWorkArray(sp->dm, npoints * origNb * origNc * origDim, MPIU_REAL, &inD));  in PetscSpaceEvaluate_Subspace()
124 …if (H) PetscCall(DMGetWorkArray(sp->dm, npoints * origNb * origNc * origDim * origDim, MPIU_REAL, …  in PetscSpaceEvaluate_Subspace()
129     PetscCall(DMGetWorkArray(sp->dm, origNc * origDim * origDim, MPIU_REAL, &phi));  in PetscSpaceEvaluate_Subspace()
130     PetscCall(DMGetWorkArray(sp->dm, origNc * subDim * subDim, MPIU_REAL, &psi));  in PetscSpaceEvaluate_Subspace()
136         for (k = 0; k < origNc * origDim; k++) phi[k] = 0.;  in PetscSpaceEvaluate_Subspace()
137         for (k = 0; k < origNc * subDim; k++) psi[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()
160             for (l = 0; l < origNc; l++) {  in PetscSpaceEvaluate_Subspace()
161 … H[((j * subNb + i) * subNc + k) * subDim * subDim + m] += subsp->Ju[k * origNc + l] * psi[l * sub…  in PetscSpaceEvaluate_Subspace()
165           for (k = 0; k < PetscMin(subNc, origNc); k++) {  in PetscSpaceEvaluate_Subspace()
172     PetscCall(DMRestoreWorkArray(sp->dm, origNc * origDim, MPIU_REAL, &phi));  in PetscSpaceEvaluate_Subspace()
173     PetscCall(DMRestoreWorkArray(sp->dm, npoints * origNb * origNc * origDim, MPIU_REAL, &inH));  in PetscSpaceEvaluate_Subspace()
178     PetscCall(DMGetWorkArray(sp->dm, origNc * origDim, MPIU_REAL, &phi));  in PetscSpaceEvaluate_Subspace()
179     PetscCall(DMGetWorkArray(sp->dm, origNc * subDim, MPIU_REAL, &psi));  in PetscSpaceEvaluate_Subspace()
185         for (k = 0; k < origNc * origDim; k++) phi[k] = 0.;  in PetscSpaceEvaluate_Subspace()
186         for (k = 0; k < origNc * subDim; k++) psi[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()
203             for (l = 0; l < origNc; l++) {  in PetscSpaceEvaluate_Subspace()
204 …Dim; m++) D[((j * subNb + i) * subNc + k) * subDim + m] += subsp->Ju[k * origNc + l] * psi[l * sub…  in PetscSpaceEvaluate_Subspace()
208           for (k = 0; k < PetscMin(subNc, origNc); k++) {  in PetscSpaceEvaluate_Subspace()
215     PetscCall(DMRestoreWorkArray(sp->dm, origNc * origDim, MPIU_REAL, &phi));  in PetscSpaceEvaluate_Subspace()
216     PetscCall(DMRestoreWorkArray(sp->dm, npoints * origNb * origNc * origDim, MPIU_REAL, &inD));  in PetscSpaceEvaluate_Subspace()
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()
247     PetscCall(DMRestoreWorkArray(sp->dm, npoints * origNb * origNc, MPIU_REAL, &inB));  in PetscSpaceEvaluate_Subspace()
281   PetscInt             origDim, subDim, origNc, subNc, origNb, subNb, f, i, j, numPoints, offset;  in PetscSpaceSetUp_Subspace()  local
295   PetscCall(PetscSpaceGetNumComponents(origSpace, &origNc));  in PetscSpaceSetUp_Subspace()
312 …alloc3(numPoints * origDim, &allPoints, numPoints * origNc, &allWeights, numPoints * origNb * orig…  in PetscSpaceSetUp_Subspace()
334       for (i = 0; i < origNc; i++) allWeights[origNc * offset + i] = 0.0;  in PetscSpaceSetUp_Subspace()
336         for (i = 0; i < origNc; i++) {  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()
361 …for (s = 0; s < origNc; s++) V[f * origNb + b] += B[(offset * origNb + b) * origNc + s] * allWeigh…  in PetscSpaceSetUp_Subspace()
450   PetscInt             origDim, subDim, origNc, subNc, subNb;  in PetscSpaceCreateSubspace()  local
462   PetscCall(PetscSpaceGetNumComponents(origSpace, &origNc));  in PetscSpaceCreateSubspace()
506       PetscCall(PetscMalloc1(origNc * subNc, &subsp->Ju_alloc));  in PetscSpaceCreateSubspace()
507       PetscCall(PetscArraycpy(subsp->Ju_alloc, Ju, origNc * subNc));  in PetscSpaceCreateSubspace()