115 … dim, PetscInt npoints, const PetscReal points[], PetscInt jet, PetscInt Njet, PetscReal pScalar[])  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()
141   PetscReal       *pScalar;  in PetscSpaceEvaluate_Polynomial()  local
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()
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()
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()
218 …npoints; p++) H[p * p_strl + b * b_strl + d1 * d1_strl + d2 * d2_strl] = pScalar[b * b_strr + j * …  in PetscSpaceEvaluate_Polynomial()
224   PetscCall(DMRestoreWorkArray(dm, Nb * Njet * npoints, MPIU_REAL, &pScalar));  in PetscSpaceEvaluate_Polynomial()