Searched refs:subsetw (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/dm/dt/interface/ |
| H A D | dtaltv.c | 386 PetscInt *subsetw, *subsetv; in PetscDTAltVPullback() local
404 PetscCall(PetscMalloc5(k, &subsetw, k, &subsetv, k, &perm, N * k, &Lw, k * k, &Lwv)); in PetscDTAltVPullback()
407 PetscCall(PetscDTEnumSubset(M, k, i, subsetw)); in PetscDTAltVPullback()
416 for (l = 0; l < k; l++) prod *= L[subsetw[perm[l]] * N + subsetv[l]]; in PetscDTAltVPullback()
429 PetscCall(PetscFree5(subsetw, subsetv, perm, Lw, Lwv)); in PetscDTAltVPullback()
456 PetscInt *subsetw, *subsetv; in PetscDTAltVPullbackMatrix() local
513 PetscCall(PetscMalloc5(M, &subsetw, N, &subsetv, k, &perm, N * k, &Lw, k * k, &Lwv)); in PetscDTAltVPullbackMatrix()
519 PetscCall(PetscDTEnumSplit(M, k, i, subsetw, &iOdd)); in PetscDTAltVPullbackMatrix()
535 for (l = 0; l < k; l++) prod *= L[subsetw[perm[l]] * N + subsetv[l]]; in PetscDTAltVPullbackMatrix()
540 PetscCall(PetscFree5(subsetw, subsetv, perm, Lw, Lwv)); in PetscDTAltVPullbackMatrix()
|