497   PetscInt         dim, Nc, imageNc, formSize, Ncopies, imageFormSize, Npoints, pt, i, j, c;  in PetscQuadraturePushForward()  local
518   for (pt = 0; pt < Npoints; pt++) {  in PetscQuadraturePushForward()
519     const PetscReal *point      = PetscSafePointerPlusOffset(points, pt * dim);  in PetscQuadraturePushForward()
520     PetscReal       *imagePoint = &imagePoints[pt * imageDim];  in PetscQuadraturePushForward()
529       const PetscReal *form      = &weights[pt * Nc + c * formSize];  in PetscQuadraturePushForward()
530       PetscReal       *imageForm = &imageWeights[pt * imageNc + c * imageFormSize];  in PetscQuadraturePushForward()
1099   PetscInt   degidx, kidx, d, pt;  in PetscDTPKDEvalJet()  local
1163     for (pt = 0; pt < npoints; pt++) {  in PetscDTPKDEvalJet()
1167       thetanm1x = dim - (d + 1) + 2. * points[pt * dim + d];  in PetscDTPKDEvalJet()
1168       for (e = d + 1; e < dim; e++) thetanm1x += points[pt * dim + e];  in PetscDTPKDEvalJet()
1171       for (e = d + 1; e < dim; e++) thetanm1 -= points[pt * dim + e];  in PetscDTPKDEvalJet()
1180 …p[(degidx * Nk + kidx) * npoints + pt] = (cnm1 * thetanm1 + cnm1x * thetanm1x) * p[(m1idx * Nk + k…  in PetscDTPKDEvalJet()
1181 …2idx >= 0) p[(degidx * Nk + kidx) * npoints + pt] -= cnm2 * thetanm2 * p[(m2idx * Nk + kidx) * npo…  in PetscDTPKDEvalJet()
1196 …p[(degidx * Nk + kidx) * npoints + pt] += mplty * 0.5 * (cnm1x - cnm1) * p[(m1idx * Nk + km1idx) *…  in PetscDTPKDEvalJet()
1198 …p[(degidx * Nk + kidx) * npoints + pt] += mplty * cnm2 * thetanm1 * p[(m2idx * Nk + km1idx) * npoi…  in PetscDTPKDEvalJet()
1210 …p[(degidx * Nk + kidx) * npoints + pt] -= 0.5 * factor * cnm2 * p[(m2idx * Nk + km2idx) * npoints …  in PetscDTPKDEvalJet()
1215 … p[(degidx * Nk + kidx) * npoints + pt] += mplty * cnm1x * p[(m1idx * Nk + km1idx) * npoints + pt];  in PetscDTPKDEvalJet()
1347 …for (PetscInt pt = 0; pt < npoints; pt++) p_jet[pt] += scale * h_jet[pt] * (points[pt * dim + v] -…  in PetscDTPTrimmedEvalJet_Internal()  local
1355             for (PetscInt pt = 0; pt < npoints; pt++) p_jet[pt] += scale * mult * h_jet[pt];  in PetscDTPTrimmedEvalJet_Internal()  local
2038     for (PetscInt pt = 0, l = 0; l < totprev; l++) {  in PetscDTStroudConicalQuadrature()  local
2040         for (PetscInt m = 0; m < totrem; m++, pt++) {  in PetscDTStroudConicalQuadrature()
2041 …for (PetscInt k = 0; k < i; k++) x[pt * dim + k] = (x[pt * dim + k] + 1.) * (1. - p1[j]) * 0.5 - 1…  in PetscDTStroudConicalQuadrature()
2042           x[pt * dim + i] = p1[j];  in PetscDTStroudConicalQuadrature()
2043           for (PetscInt c = 0; c < Nc; c++) w[pt * Nc + c] *= mul * w1[j];  in PetscDTStroudConicalQuadrature()