96   PetscInt             origDim, subDim, origNc, subNc, subNb, origNb, i, j, k, l, m, n, o;  in PetscSpaceEvaluate_Subspace()  local
106   PetscCall(PetscSpaceGetDimension(origsp, &origNb));  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()
133       const PetscReal *subq = &subsp->Q[i * origNb];  in PetscSpaceEvaluate_Subspace()
138         for (k = 0; k < origNb; k++) {  in PetscSpaceEvaluate_Subspace()
139 …for (l = 0; l < origNc * origDim * origDim; l++) phi[l] += inH[(j * origNb + k) * origNc * origDim…  in PetscSpaceEvaluate_Subspace()
173     PetscCall(DMRestoreWorkArray(sp->dm, npoints * origNb * origNc * origDim, MPIU_REAL, &inH));  in PetscSpaceEvaluate_Subspace()
182       const PetscReal *subq = &subsp->Q[i * origNb];  in PetscSpaceEvaluate_Subspace()
187         for (k = 0; k < origNb; k++) {  in PetscSpaceEvaluate_Subspace()
188 …for (l = 0; l < origNc * origDim; l++) phi[l] += inD[(j * origNb + k) * origNc * origDim + l] * su…  in PetscSpaceEvaluate_Subspace()
216     PetscCall(DMRestoreWorkArray(sp->dm, npoints * origNb * origNc * origDim, MPIU_REAL, &inD));  in PetscSpaceEvaluate_Subspace()
230       const PetscReal *subq = &subsp->Q[i * origNb];  in PetscSpaceEvaluate_Subspace()
234         for (k = 0; k < origNb; k++) {  in PetscSpaceEvaluate_Subspace()
235           for (l = 0; l < origNc; l++) phi[l] += inB[(j * origNb + k) * origNc + l] * subq[k];  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
299   PetscCall(PetscSpaceGetDimension(origSpace, &origNb));  in PetscSpaceSetUp_Subspace()
311   PetscCall(PetscMalloc1(subNb * origNb, &V));  in PetscSpaceSetUp_Subspace()
312 …umPoints * origDim, &allPoints, numPoints * origNc, &allWeights, numPoints * origNb * origNc, &B));  in PetscSpaceSetUp_Subspace()
353       for (b = 0; b < origNb; b++) {  in PetscSpaceSetUp_Subspace()
354         for (s = 0; s < subNc; s++) V[f * origNb + b] += qw[s] * u[s];  in PetscSpaceSetUp_Subspace()
357       for (b = 0; b < origNb; b++) V[f * origNb + b] = 0.0;  in PetscSpaceSetUp_Subspace()
360       for (b = 0; b < origNb; b++) {  in PetscSpaceSetUp_Subspace()
361 …for (s = 0; s < origNc; s++) V[f * origNb + b] += B[(offset * origNb + b) * origNc + s] * allWeigh…  in PetscSpaceSetUp_Subspace()
369     for (i = 0; i < origNb; i++) rho += PetscSqr(V[f * origNb + i]);  in PetscSpaceSetUp_Subspace()
373     for (i = 0; i < origNb; i++) V[f * origNb + i] *= scal;  in PetscSpaceSetUp_Subspace()
375       for (i = 0, scal = 0.; i < origNb; i++) scal += V[f * origNb + i] * V[j * origNb + i];  in PetscSpaceSetUp_Subspace()
376       for (i = 0; i < origNb; i++) V[j * origNb + i] -= V[f * origNb + i] * scal;  in PetscSpaceSetUp_Subspace()