Searched refs:lhh (Results 1 – 2 of 2) sorted by relevance
| /petsc/src/ksp/ksp/impls/gmres/ |
| H A D | borthog2.c | 39 PetscScalar *hh, *hes, *lhh; in KSPGMRESClassicalGramSchmidtOrthogonalization() local
46 lhh = gmres->orthogwork; in KSPGMRESClassicalGramSchmidtOrthogonalization()
62 PetscCall(VecMDot(VEC_VV(it + 1), it + 1, &(VEC_VV(0)), lhh)); /* <v,vnew> */ in KSPGMRESClassicalGramSchmidtOrthogonalization()
64 KSPCheckDot(ksp, lhh[j]); in KSPGMRESClassicalGramSchmidtOrthogonalization()
66 lhh[j] = -lhh[j]; in KSPGMRESClassicalGramSchmidtOrthogonalization()
73 PetscCall(VecMAXPY(VEC_VV(it + 1), it + 1, lhh, &VEC_VV(0))); in KSPGMRESClassicalGramSchmidtOrthogonalization()
76 hh[j] -= lhh[j]; /* hh += <v,vnew> */ in KSPGMRESClassicalGramSchmidtOrthogonalization()
77 hes[j] -= lhh[j]; /* hes += <v,vnew> */ in KSPGMRESClassicalGramSchmidtOrthogonalization()
86 for (j = 0; j <= it; j++) hnrm += PetscRealPart(lhh[j] * PetscConj(lhh[j])); in KSPGMRESClassicalGramSchmidtOrthogonalization()
99 PetscCall(VecMDot(VEC_VV(it + 1), it + 1, &(VEC_VV(0)), lhh)); /* <v,vnew> */ in KSPGMRESClassicalGramSchmidtOrthogonalization()
[all …]
|
| /petsc/src/ksp/ksp/impls/gmres/pipefgmres/ |
| H A D | pipefgmres.c | 52 PetscScalar *hh, *hes, *lhh, shift = pipefgmres->shift; in KSPPIPEFGMRESCycle() local
71 lhh = pipefgmres->orthogwork; in KSPPIPEFGMRESCycle()
140 PetscCall(VecMDotBegin(ZVEC(loc_it), loc_it + 2, redux, lhh)); in KSPPIPEFGMRESCycle()
164 PetscCall(VecMDotEnd(ZVEC(loc_it), loc_it + 2, redux, lhh)); in KSPPIPEFGMRESCycle()
165 tt = PetscRealPart(lhh[loc_it + 1]); in KSPPIPEFGMRESCycle()
173 hh[j] = lhh[j]; in KSPPIPEFGMRESCycle()
174 hes[j] = lhh[j]; in KSPPIPEFGMRESCycle()
176 hh[loc_it] = lhh[loc_it] + shift; in KSPPIPEFGMRESCycle()
177 hes[loc_it] = lhh[loc_it] + shift; in KSPPIPEFGMRESCycle()
180 for (j = 0; j <= loc_it; j++) lhh[j] = -lhh[j]; /* flip sign */ in KSPPIPEFGMRESCycle()
[all …]
|