Searched refs:Lstarw (Results 1 – 2 of 2) sorted by relevance
| /petsc/src/dm/dt/interface/ |
| H A D | dtaltv.c | 323 …back(PetscInt N, PetscInt M, const PetscReal *L, PetscInt k, const PetscReal *w, PetscReal *Lstarw) in PetscDTAltVPullback() argument
334 Lstarw[0] = w[0]; in PetscDTAltVPullback()
340 Lstarw[i] = sum; in PetscDTAltVPullback()
349 Lstarw[i] = mult[i] * sum; in PetscDTAltVPullback()
361 Lstarw[i] = sum; in PetscDTAltVPullback()
376 Lstarw[i] = sum; in PetscDTAltVPullback()
381 for (i = 0; i < Nk; i++) Lstarw[i] = detL * w[i]; in PetscDTAltVPullback()
405 for (i = 0; i < Nk; i++) Lstarw[i] = 0.; in PetscDTAltVPullback()
417 Lstarw[j] += prod; in PetscDTAltVPullback()
425 PetscCall(PetscDTAltVStar(N, N - k, -1, Lstarw, sLsw)); in PetscDTAltVPullback()
[all …]
|
| /petsc/src/dm/dt/tests/ |
| H A D | ex7.c | 9 PetscReal *Lstarw, *Lx, *Lstar, *Lstarwcheck, wLx, Lstarwx; in CheckPullback() local
26 PetscCall(PetscMalloc2(Nk, &Lstarw, M * k, &Lx)); in CheckPullback()
28 PetscCall(PetscDTAltVPullback(N, M, L, negative ? -k : k, w, Lstarw)); in CheckPullback()
34 PetscCall(PetscDTAltVStar(N, N - k, 1, Lstarw, sLsw)); in CheckPullback()
38 PetscCall(PetscDTAltVApply(N, k, Lstarw, x, &Lstarwx)); in CheckPullback()
54 diffMat += PetscSqr(PetscAbsReal(Lstarwcheck[i] - Lstarw[i])); in CheckPullback()
55 normMat += PetscSqr(Lstarwcheck[i]) + PetscSqr(Lstarw[i]); in CheckPullback()
72 if (Nk > 0) PetscCall(PetscRealView(Nk, Lstarw, viewer)); in CheckPullback()
80 PetscCall(PetscFree2(Lstarw, Lx)); in CheckPullback()
|