1099   PetscInt   degidx, kidx, d, pt;  in PetscDTPKDEvalJet()  local
1115   for (degidx = 0; degidx < Ndeg; degidx++) {  in PetscDTPKDEvalJet()
1124     PetscCall(PetscDTIndexToGradedOrder(dim, degidx, degtup));  in PetscDTPKDEvalJet()
1128       scales[degidx] = initscale;  in PetscDTPKDEvalJet()
1131         scales[degidx] /= norm;  in PetscDTPKDEvalJet()
1133       for (i = 0; i < npoints; i++) p[degidx * Nk * npoints + i] = 1.;  in PetscDTPKDEvalJet()
1134       for (i = 0; i < (Nk - 1) * npoints; i++) p[(degidx * Nk + 1) * npoints + i] = 0.;  in PetscDTPKDEvalJet()
1150     scales[degidx] = initscale;  in PetscDTPKDEvalJet()
1157       for (f = 0; f < degsum; f++) scales[degidx] *= 2.;  in PetscDTPKDEvalJet()
1159       scales[degidx] /= enorm;  in PetscDTPKDEvalJet()
1180 …p[(degidx * Nk + kidx) * npoints + pt] = (cnm1 * thetanm1 + cnm1x * thetanm1x) * p[(m1idx * Nk + k…  in PetscDTPKDEvalJet()
1181 …if (m2idx >= 0) p[(degidx * Nk + kidx) * npoints + pt] -= cnm2 * thetanm2 * p[(m2idx * Nk + kidx) …  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()
1222   for (degidx = 0; degidx < Ndeg; degidx++) {  in PetscDTPKDEvalJet()
1223     PetscReal scale = scales[degidx];  in PetscDTPKDEvalJet()
1226     for (i = 0; i < Nk * npoints; i++) p[degidx * Nk * npoints + i] *= scale;  in PetscDTPKDEvalJet()