Lines Matching refs:Nk
1100 PetscInt Nk, Ndeg;
1106 PetscCall(PetscDTBinomialInt(dim + k, k, &Nk));
1133 for (i = 0; i < npoints; i++) p[degidx * Nk * npoints + i] = 1.;
1134 for (i = 0; i < (Nk - 1) * npoints; i++) p[(degidx * Nk + 1) * npoints + i] = 0.;
1175 for (kidx = 0; kidx < Nk; kidx++) {
1180 p[(degidx * Nk + kidx) * npoints + pt] = (cnm1 * thetanm1 + cnm1x * thetanm1x) * p[(m1idx * Nk + kidx) * npoints + pt];
1181 if (m2idx >= 0) p[(degidx * Nk + kidx) * npoints + pt] -= cnm2 * thetanm2 * p[(m2idx * Nk + kidx) * npoints + pt];
1196 p[(degidx * Nk + kidx) * npoints + pt] += mplty * 0.5 * (cnm1x - cnm1) * p[(m1idx * Nk + km1idx) * npoints + pt];
1198 p[(degidx * Nk + kidx) * npoints + pt] += mplty * cnm2 * thetanm1 * p[(m2idx * Nk + km1idx) * npoints + pt];
1210 p[(degidx * Nk + kidx) * npoints + pt] -= 0.5 * factor * cnm2 * p[(m2idx * Nk + km2idx) * npoints + pt];
1215 p[(degidx * Nk + kidx) * npoints + pt] += mplty * cnm1x * p[(m1idx * Nk + km1idx) * npoints + pt];
1226 for (i = 0; i < Nk * npoints; i++) p[degidx * Nk * npoints + i] *= scale;
1283 PetscInt Nk;
1284 PetscCall(PetscDTBinomialInt(dim + jetDegree, dim, &Nk));
1285 PetscCall(PetscArrayzero(p, Nbpt * Nf * Nk * npoints));
1290 PetscCall(PetscMalloc1(Nbpm1 * Nk * npoints, &p_scalar));
1296 const PetscReal *src = &p_scalar[i * Nk * npoints];
1298 PetscReal *dest = &p[(total++ * Nf + f) * Nk * npoints];
1299 PetscCall(PetscArraycpy(dest, src, Nk * npoints));
1321 const PetscReal *h_scalar = &p_scalar[(Nbpm1 - Nh) * Nk * npoints];
1323 const PetscReal *h_s = &h_scalar[b * Nk * npoints];
1331 PetscReal *p_f = &p[total++ * Nf * Nk * npoints];
1342 PetscReal *p_i = &p_f[i * Nk * npoints];
1343 for (PetscInt jet = 0; jet < Nk; jet++) {
1353 if (l >= Nk) continue;