Searched refs:bell (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/ksp/ksp/impls/bcgsl/ |
| H A D | bcgsl.c | 164 PetscBLASInt ione = 1, bell; in KSPSolve_BCGSL() local
165 PetscCall(PetscBLASIntCast(bcgsl->ell, &bell)); in KSPSolve_BCGSL()
170 …APACKgesvd", LAPACKgesvd_("A", "A", &bell, &bell, &MZa[1 + ldMZ], &ldMZ, bcgsl->s, bcgsl->u, &bell… in KSPSolve_BCGSL()
172 …APACKgesvd", LAPACKgesvd_("A", "A", &bell, &bell, &MZa[1 + ldMZ], &ldMZ, bcgsl->s, bcgsl->u, &bell… in KSPSolve_BCGSL()
180 for (i = 1; i < bell; i++) { in KSPSolve_BCGSL()
184 pinv_tol = bell * max_s * PETSC_MACHINE_EPSILON; in KSPSolve_BCGSL()
185 PetscCall(PetscArrayzero(&AY0c[1], bell)); in KSPSolve_BCGSL()
186 for (i = 0; i < bell; i++) { in KSPSolve_BCGSL()
189 for (j = 0; j < bell; j++) utb += MZb[1 + j] * bcgsl->u[i * bell + j]; in KSPSolve_BCGSL()
191 for (j = 0; j < bell; j++) AY0c[1 + j] += utb / bcgsl->s[i] * bcgsl->v[j * bell + i]; in KSPSolve_BCGSL()
[all …]
|