51     PetscInt origDim, subDim, origNc, subNc, o, s;  in PetscSpaceView_Subspace()  local
55     PetscCall(PetscSpaceGetNumVariables(sp, &subDim));  in PetscSpaceView_Subspace()
64         PetscCall(PetscViewerASCIIPrintf(viewer, " %g", (double)subsp->Jx[o * subDim + 0]));  in PetscSpaceView_Subspace()
66 …for (s = 1; s < subDim; s++) PetscCall(PetscViewerASCIIPrintf(viewer, " %g", (double)subsp->Jx[o *…  in PetscSpaceView_Subspace()
96   PetscInt             origDim, subDim, origNc, subNc, subNb, origNb, i, j, k, l, m, n, o;  in PetscSpaceEvaluate_Subspace()  local
101   PetscCall(PetscSpaceGetNumVariables(sp, &subDim));  in PetscSpaceEvaluate_Subspace()
116 …for (k = 0; k < subDim; k++) inpoints[i * origDim + j] += subsp->Jx[j * subDim + k] * points[i * s…  in PetscSpaceEvaluate_Subspace()
119 …for (j = 0; j < PetscMin(subDim, origDim); j++) inpoints[i * origDim + j] += points[i * subDim + j…  in PetscSpaceEvaluate_Subspace()
130     PetscCall(DMGetWorkArray(sp->dm, origNc * subDim * subDim, MPIU_REAL, &psi));  in PetscSpaceEvaluate_Subspace()
131     for (i = 0; i < npoints * subNb * subNc * subDim; i++) D[i] = 0.0;  in PetscSpaceEvaluate_Subspace()
137         for (k = 0; k < origNc * subDim; k++) psi[k] = 0.;  in PetscSpaceEvaluate_Subspace()
143             for (l = 0; l < subDim; l++) {  in PetscSpaceEvaluate_Subspace()
145                 for (n = 0; n < subDim; n++) {  in PetscSpaceEvaluate_Subspace()
146 …or (o = 0; o < origDim; o++) psi[(k * subDim + l) * subDim + n] += subsp->Jx[m * subDim + l] * sub…  in PetscSpaceEvaluate_Subspace()
153             for (l = 0; l < PetscMin(subDim, origDim); l++) {  in PetscSpaceEvaluate_Subspace()
154 …for (m = 0; m < PetscMin(subDim, origDim); m++) psi[(k * subDim + l) * subDim + m] += phi[(k * ori…  in PetscSpaceEvaluate_Subspace()
161 …m = 0; m < subDim * subDim; m++) H[((j * subNb + i) * subNc + k) * subDim * subDim + m] += subsp->…  in PetscSpaceEvaluate_Subspace()
166 …for (l = 0; l < subDim * subDim; l++) H[((j * subNb + i) * subNc + k) * subDim * subDim + l] += ps…  in PetscSpaceEvaluate_Subspace()
179     PetscCall(DMGetWorkArray(sp->dm, origNc * subDim, MPIU_REAL, &psi));  in PetscSpaceEvaluate_Subspace()
180     for (i = 0; i < npoints * subNb * subNc * subDim; i++) D[i] = 0.0;  in PetscSpaceEvaluate_Subspace()
186         for (k = 0; k < origNc * subDim; k++) psi[k] = 0.;  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()
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()
281   PetscInt             origDim, subDim, origNc, subNc, origNb, subNb, f, i, j, numPoints, offset;  in PetscSpaceSetUp_Subspace()  local
298   PetscCall(DMGetDimension(dm, &subDim));  in PetscSpaceSetUp_Subspace()
329 …        for (j = 0; j < subDim; j++) allPoints[origDim * offset + i] += Jx[i * subDim + j] * qp[j];  in PetscSpaceSetUp_Subspace()
332         for (i = 0; i < PetscMin(subDim, origDim); i++) allPoints[origDim * offset + i] += qp[i];  in PetscSpaceSetUp_Subspace()
393       PetscInt subDim, origDim, i, j;  in PetscSpacePolynomialGetTensor_Subspace()  local
397       PetscCall(PetscSpaceGetNumVariables(sp, &subDim));  in PetscSpacePolynomialGetTensor_Subspace()
402         for (j = 0; j < subDim; j++) nnz += (subsp->Jx[i * subDim + j] != 0.);  in PetscSpacePolynomialGetTensor_Subspace()
405       for (j = 0; j < subDim; j++) {  in PetscSpacePolynomialGetTensor_Subspace()
408         for (i = 0; i < origDim; i++) nnz += (subsp->Jx[i * subDim + j] != 0.);  in PetscSpacePolynomialGetTensor_Subspace()
450   PetscInt             origDim, subDim, origNc, subNc, subNb;  in PetscSpaceCreateSubspace()  local
465   PetscCall(DMGetDimension(dm, &subDim));  in PetscSpaceCreateSubspace()
470   PetscCall(PetscSpaceSetNumVariables(*subspace, subDim));  in PetscSpaceCreateSubspace()
496       PetscCall(PetscMalloc1(origDim * subDim, &subsp->Jx_alloc));  in PetscSpaceCreateSubspace()
497       PetscCall(PetscArraycpy(subsp->Jx_alloc, Jx, origDim * subDim));  in PetscSpaceCreateSubspace()