Searched refs:Mk (Results 1 – 2 of 2) sorted by relevance
| /petsc/src/dm/dt/interface/ |
| H A D | dtaltv.c | 325 PetscInt i, j, Nk, Mk; in PetscDTAltVPullback() local
331 PetscCall(PetscDTBinomialInt(M, PetscAbsInt(k), &Mk)); in PetscDTAltVPullback()
339 for (j = 0; j < Mk; j++) sum += L[j * Nk + i] * w[j]; in PetscDTAltVPullback()
348 for (j = 0; j < Mk; j++) sum += L[(Mk - 1 - j) * Nk + (Nk - 1 - i)] * w[j] * mult[j]; in PetscDTAltVPullback()
360 …for (j = 0; j < Mk; j++) sum += (L[pairs[j][0] * N + pairs[i][0]] * L[pairs[j][1] * N + pairs[i][1… in PetscDTAltVPullback()
375 …for (j = 0; j < Mk; j++) sum += (L[pairs[offj + j][0] * N + pairs[offi + i][0]] * L[pairs[offj + j… in PetscDTAltVPullback()
392 PetscCall(PetscDTBinomialInt(M, PetscAbsInt(k), &Mk)); in PetscDTAltVPullback()
398 PetscCall(PetscMalloc1(Mk, &walloc)); in PetscDTAltVPullback()
406 for (i = 0; i < Mk; i++) { in PetscDTAltVPullback()
454 PetscInt Nk, Mk, Nf, i, j, l, p; in PetscDTAltVPullbackMatrix() local
[all …]
|
| /petsc/src/dm/dt/tests/ |
| H A D | ex7.c | 8 PetscInt Nk, Mk, i, j, l; in CheckPullback() local
18 PetscCall(PetscDTBinomialInt(M, k, &Mk)); in CheckPullback()
20 PetscCall(PetscMalloc1(Mk, &walloc)); in CheckPullback()
27 PetscCall(PetscMalloc2(Nk * Mk, &Lstar, Nk, &Lstarwcheck)); in CheckPullback()
52 for (j = 0; j < Mk; j++) sum += Lstar[i * Mk + j] * w[j]; in CheckPullback()
67 if (Nk * Mk > 0) PetscCall(PetscRealView(Nk * Mk, Lstar, viewer)); in CheckPullback()
354 PetscInt Mk, l; in main() local
356 PetscCall(PetscDTBinomialInt(M, k, &Mk)); in main()
357 PetscCall(PetscMalloc3(M * N, &L, Mk, &u, M * k, &x)); in main()
359 for (l = 0; l < Mk; l++) PetscCall(PetscRandomGetValueReal(rand, &u[l])); in main()
|