Searched refs:Lstar (Results 1 – 2 of 2) sorted by relevance
| /petsc/src/dm/dt/interface/ |
| H A D | dtaltv.c | 452 …PetscDTAltVPullbackMatrix(PetscInt N, PetscInt M, const PetscReal *L, PetscInt k, PetscReal *Lstar) in PetscDTAltVPullbackMatrix() argument
469 Lstar[0] = 1.; in PetscDTAltVPullbackMatrix()
472 for (j = 0; j < Mk; j++) Lstar[i * Mk + j] = L[j * Nk + i]; in PetscDTAltVPullbackMatrix()
476 …for (j = 0; j < Mk; j++) Lstar[i * Mk + j] = L[(Mk - 1 - j) * Nk + (Nk - 1 - i)] * mult[i] * mult[… in PetscDTAltVPullbackMatrix()
486 …for (j = 0; j < Mk; j++) Lstar[i * Mk + j] = L[pairs[j][0] * N + pairs[i][0]] * L[pairs[j][1] * N … in PetscDTAltVPullbackMatrix()
498 …for (j = 0; j < Mk; j++) Lstar[i * Mk + j] = L[pairs[offj + j][0] * N + pairs[offi + i][0]] * L[pa… in PetscDTAltVPullbackMatrix()
503 for (i = 0; i < Nk; i++) Lstar[i] = detL; in PetscDTAltVPullbackMatrix()
514 for (i = 0; i < Nk * Mk; i++) Lstar[i] = 0.; in PetscDTAltVPullbackMatrix()
536 Lstar[jidx * Mk + iidx] += prod; in PetscDTAltVPullbackMatrix()
|
| /petsc/src/dm/dt/tests/ |
| H A D | ex7.c | 9 PetscReal *Lstarw, *Lx, *Lstar, *Lstarwcheck, wLx, Lstarwx; in CheckPullback() local
27 PetscCall(PetscMalloc2(Nk * Mk, &Lstar, Nk, &Lstarwcheck)); in CheckPullback()
29 PetscCall(PetscDTAltVPullbackMatrix(N, M, L, negative ? -k : k, Lstar)); 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()
79 PetscCall(PetscFree2(Lstar, Lstarwcheck)); in CheckPullback()
|