Searched refs:aug1 (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/ksp/ksp/impls/gmres/dgmres/ |
| H A D | dgmres.c | 839 PetscInt aug1 = bmax + neig1; /* maximum size of the augmented invariant basis */ in KSPDGMRESImproveEig_DGMRES() local
855 PetscCall(PetscMalloc1(aug1 * aug1, &AUU)); in KSPDGMRESImproveEig_DGMRES()
856 PetscCall(PetscMalloc1(aug1 * aug1, &AUAU)); in KSPDGMRESImproveEig_DGMRES()
861 for (j = 0; j < r; j++) PetscCall(VecMDot(UU[j], r, MU, &AUU[j * aug1])); in KSPDGMRESImproveEig_DGMRES()
863 for (j = 0; j < neig; j++) PetscCall(VecMDot(XX[j], r, MU, &AUU[(r + j) * aug1])); in KSPDGMRESImproveEig_DGMRES()
865 for (j = 0; j < r; j++) PetscCall(VecMDot(UU[j], neig, MX, &AUU[j * aug1 + r])); in KSPDGMRESImproveEig_DGMRES()
867 for (j = 0; j < neig; j++) PetscCall(VecMDot(XX[j], neig, MX, &AUU[(r + j) * aug1 + r])); in KSPDGMRESImproveEig_DGMRES()
871 for (j = 0; j < r; j++) PetscCall(VecMDot(MU[j], r, MU, &AUAU[j * aug1])); in KSPDGMRESImproveEig_DGMRES()
873 for (j = 0; j < neig; j++) PetscCall(VecMDot(MX[j], r, MU, &AUAU[(r + j) * aug1])); in KSPDGMRESImproveEig_DGMRES()
875 for (j = 0; j < r; j++) PetscCall(VecMDot(MU[j], neig, MX, &AUAU[j * aug1 + r])); in KSPDGMRESImproveEig_DGMRES()
[all …]
|