115 static PetscErrorCode CoordinateBasis(PetscInt dim, PetscInt npoints, const PetscReal points[], Pet…  in CoordinateBasis()  argument
118   PetscCall(PetscArrayzero(pScalar, (1 + dim) * Njet * npoints));  in CoordinateBasis()
123 …     for (PetscInt pt = 0; pt < npoints; pt++) pScalar[b * Njet * npoints + j * npoints + pt] = 1.;  in CoordinateBasis()
125 …for (PetscInt pt = 0; pt < npoints; pt++) pScalar[b * Njet * npoints + j * npoints + pt] = points[…  in CoordinateBasis()
128 …     for (PetscInt pt = 0; pt < npoints; pt++) pScalar[b * Njet * npoints + j * npoints + pt] = 1.;  in CoordinateBasis()
135 static PetscErrorCode PetscSpaceEvaluate_Polynomial(PetscSpace sp, PetscInt npoints, const PetscRea…  in PetscSpaceEvaluate_Polynomial()  argument
146     PetscCall(PetscSpaceEvaluate(sp, npoints, points, B, D, H));  in PetscSpaceEvaluate_Polynomial()
159   PetscCall(DMGetWorkArray(dm, Nb * Njet * npoints, MPIU_REAL, &pScalar));  in PetscSpaceEvaluate_Polynomial()
164     PetscCall(CoordinateBasis(dim, npoints, points, jet, Njet, pScalar));  in PetscSpaceEvaluate_Polynomial()
166     PetscCall(PetscDTPKDEvalJet(dim, npoints, points, sp->degree, jet, pScalar));  in PetscSpaceEvaluate_Polynomial()
172     PetscInt b_strr = Njet * npoints;  in PetscSpaceEvaluate_Polynomial()
175     PetscCall(PetscArrayzero(B, npoints * Nb));  in PetscSpaceEvaluate_Polynomial()
177 …for (PetscInt p = 0; p < npoints; p++) B[p * p_strl + b * b_strl] = pScalar[b * b_strr + p * p_str…  in PetscSpaceEvaluate_Polynomial()
185     PetscInt b_strr = Njet * npoints;  in PetscSpaceEvaluate_Polynomial()
186     PetscInt d_strr = npoints;  in PetscSpaceEvaluate_Polynomial()
189     PetscCall(PetscArrayzero(D, npoints * Nb * dim));  in PetscSpaceEvaluate_Polynomial()
192 …for (PetscInt p = 0; p < npoints; p++) D[p * p_strl + b * b_strl + d * d_strl] = pScalar[b * b_str…  in PetscSpaceEvaluate_Polynomial()
202     PetscInt b_strr = Njet * npoints;  in PetscSpaceEvaluate_Polynomial()
203     PetscInt j_strr = npoints;  in PetscSpaceEvaluate_Polynomial()
208     PetscCall(PetscArrayzero(H, npoints * Nb * dim * dim));  in PetscSpaceEvaluate_Polynomial()
218 …for (PetscInt p = 0; p < npoints; p++) H[p * p_strl + b * b_strl + d1 * d1_strl + d2 * d2_strl] = …  in PetscSpaceEvaluate_Polynomial()
224   PetscCall(DMRestoreWorkArray(dm, Nb * Njet * npoints, MPIU_REAL, &pScalar));  in PetscSpaceEvaluate_Polynomial()