Lines Matching refs:Nk

8   PetscInt         Nk, Mk, i, j, l;  in CheckPullback()  local
17 PetscCall(PetscDTBinomialInt(N, k, &Nk)); in CheckPullback()
26 PetscCall(PetscMalloc2(Nk, &Lstarw, M * k, &Lx)); in CheckPullback()
27 PetscCall(PetscMalloc2(Nk * Mk, &Lstar, Nk, &Lstarwcheck)); in CheckPullback()
33 PetscCall(PetscMalloc1(Nk, &sLsw)); in CheckPullback()
50 for (i = 0; i < Nk; i++) { in CheckPullback()
67 if (Nk * Mk > 0) PetscCall(PetscRealView(Nk * Mk, Lstar, viewer)); in CheckPullback()
72 if (Nk > 0) PetscCall(PetscRealView(Nk, Lstarw, viewer)); in CheckPullback()
133 PetscInt j, Nk, M; in main() local
137 PetscCall(PetscDTBinomialInt(N, k, &Nk)); in main()
140 …SCIIPrintf(viewer, "(%" PetscInt_FMT " choose %" PetscInt_FMT "): %" PetscInt_FMT "\n", N, k, Nk)); in main()
145 for (j = 0; j < Nk; j++) { in main()
168 PetscCall(PetscMalloc1(Nk, &w)); in main()
169 for (j = 0; j < Nk; j++) PetscCall(PetscRandomGetValueReal(rand, &w[j])); in main()
179 if (Nk) PetscCall(PetscRealView(Nk, w, viewer)); in main()
326 PetscCall(PetscMalloc2(Nk * Njk, &uWwmat, Njk, &uWwcheck)); in main()
331 if ((Nk * Njk) > 0) PetscCall(PetscRealView(Nk * Njk, uWwmat, viewer)); in main()
340 for (m = 0; m < Nk; m++) sum += uWwmat[l * Nk + m] * w[m]; in main()
382 …tscCall(PetscMalloc5(Nkm, &wIntv0, Nkm, &wIntv0check, Nk * Nkm, &intv0mat, Nk * Nkm, &matcheck, Nk in main()
389 for (l = 0; l < Nk * k; l++) { in main()
398 for (l = 0; l < Nkm * Nk; l++) matcheck[l] = 0.; in main()
399 for (l = 0; l < Nk * k; l++) { in main()
405 matcheck[row * Nk + col] = -v[-(x + 1)]; in main()
407 matcheck[row * Nk + col] = v[x]; in main()
412 for (l = 0; l < Nkm * Nk; l++) { in main()
424 for (m = 0; m < Nk; m++) sum += intv0mat[l * Nk + m] * w[m]; in main()
441 if (Nk * Nkm > 0) PetscCall(PetscRealView(Nk * Nkm, intv0mat, viewer)); in main()
456 PetscCall(PetscMalloc3(Nk, &u, Nk, &starw, Nk, &starstarw)); in main()
462 if (Nk) PetscCall(PetscRealView(Nk, starw, viewer)); in main()
467 if (Nk) PetscCall(PetscRealView(Nk, starstarw, viewer)); in main()
470 for (l = 0; l < Nk; l++) PetscCall(PetscRandomGetValueReal(rand, &u[l])); in main()
473 for (l = 0; l < Nk; l++) starwdotu += starw[l] * u[l]; in main()
479 for (l = 0; l < Nk; l++) { in main()