Home
last modified time | relevance | path

Searched refs:VEC_VV (Results 1 – 15 of 15) sorted by relevance

/petsc/src/ksp/ksp/impls/gmres/
H A Dborthog2.c62 PetscCall(VecMDot(VEC_VV(it + 1), it + 1, &(VEC_VV(0)), lhh)); /* <v,vnew> */ in KSPGMRESClassicalGramSchmidtOrthogonalization()
73 PetscCall(VecMAXPY(VEC_VV(it + 1), it + 1, lhh, &VEC_VV(0))); in KSPGMRESClassicalGramSchmidtOrthogonalization()
89 PetscCall(VecNorm(VEC_VV(it + 1), NORM_2, &wnrm)); in KSPGMRESClassicalGramSchmidtOrthogonalization()
99 PetscCall(VecMDot(VEC_VV(it + 1), it + 1, &(VEC_VV(0)), lhh)); /* <v,vnew> */ in KSPGMRESClassicalGramSchmidtOrthogonalization()
105 PetscCall(VecMAXPY(VEC_VV(it + 1), it + 1, lhh, &VEC_VV(0))); in KSPGMRESClassicalGramSchmidtOrthogonalization()
H A Dborthog.c43 PetscCall(VecDot(VEC_VV(it + 1), VEC_VV(j), hh)); in KSPGMRESModifiedGramSchmidtOrthogonalization()
48 PetscCall(VecAXPY(VEC_VV(it + 1), -(*hh++), VEC_VV(j))); in KSPGMRESModifiedGramSchmidtOrthogonalization()
H A Dgmres.c97 PetscCall(VecNormalize(VEC_VV(0), &res)); in KSPGMRESCycle()
133 PetscCall(KSP_PCApplyBAorAB(ksp, VEC_VV(it), VEC_VV(1 + it), VEC_TEMP_MATOP)); in KSPGMRESCycle()
140 PetscCall(VecNormalize(VEC_VV(it + 1), &tt)); in KSPGMRESCycle()
214 …PetscCall(KSPInitialResidual(ksp, ksp->vec_sol, VEC_TEMP, VEC_TEMP_MATOP, VEC_VV(0), ksp->vec_rhs)… in KSPSolve_GMRES()
223 PetscCall(VecDuplicateVecs(VEC_VV(0), N, &gmres->vecb)); in KSPSolve_GMRES()
226 for (i = 0; i < gmres->max_k + 1; i++) PetscCall(VecCopy(VEC_VV(i), gmres->vecb[i])); in KSPSolve_GMRES()
336 PetscCall(VecMAXPBY(VEC_TEMP, it + 1, nrs, 0, &VEC_VV(0))); in KSPGMRESBuildSoln()
515 x = VEC_VV(gmres->it + 1); in KSPGMRESMonitorKrylov()
H A Dgmresimpl.h88 #define VEC_VV(i) gmres->vecs[VEC_OFFSET + i] macro
H A Dgmreig.c258 …; j++) PetscCall(VecMAXPBY(S[j], bn, &SR[j * bn], 0, gmres->fullcycle ? gmres->vecb : &VEC_VV(0))); in KSPComputeRitz_GMRES()
/petsc/src/ksp/ksp/impls/gmres/pgmres/
H A Dpgmres.c26 PetscCall(VecNormalize(VEC_VV(0), &res_norm)); in KSPPGMRESCycle()
50 …Zcur = VEC_VV(it); /* Zcur is not yet orthogonal, but the VecMDot to orthogonalize it has bee… in KSPPGMRESCycle()
51 …Znext = VEC_VV(it + 1); /* This iteration will compute Znext, update with a deferred correction on… in KSPPGMRESCycle()
59 PetscCall(VecNormEnd(VEC_VV(it - 1), NORM_2, &newnorm)); in KSPPGMRESCycle()
63 PetscCall(VecMDotEnd(Zcur, it, &(VEC_VV(0)), HH(0, it - 1))); in KSPPGMRESCycle()
68 PetscCall(VecScale(VEC_VV(it - 1), 1. / *HH(it - 1, it - 2))); in KSPPGMRESCycle()
130 PetscCall(VecMAXPY(Znext, it + 1, work, &VEC_VV(0))); in KSPPGMRESCycle()
135 PetscCall(VecMAXPY(Zcur, it, work, &VEC_VV(0))); in KSPPGMRESCycle()
138 PetscCall(VecNormBegin(VEC_VV(it), NORM_2, &newnorm)); in KSPPGMRESCycle()
142 PetscCall(VecMDotBegin(Znext, it + 1, &VEC_VV(0), HH(0, it))); in KSPPGMRESCycle()
[all …]
H A Dpgmresimpl.h22 #define VEC_VV(i) pgmres->vecs[VEC_OFFSET + i] /* use to access othog basis vectors */ macro
/petsc/src/ksp/ksp/impls/gmres/lgmres/
H A Dlgmres.c123 PetscCall(VecNorm(VEC_VV(0), NORM_2, &res_norm)); in KSPLGMRESCycle()
140 PetscCall(VecScale(VEC_VV(0), tmp)); in KSPLGMRESCycle()
170 PetscCall(KSP_PCApplyBAorAB(ksp, VEC_VV(loc_it), VEC_VV(1 + loc_it), VEC_TEMP_MATOP)); in KSPLGMRESCycle()
180 PetscCall(VecCopy(A_AUGVEC(spot), VEC_VV(1 + loc_it))); in KSPLGMRESCycle()
189 PetscCall(VecNorm(VEC_VV(loc_it + 1), NORM_2, &tt)); in KSPLGMRESCycle()
199 PetscCall(VecScale(VEC_VV(loc_it + 1), tmp)); /* scale new direction by its norm */ in KSPLGMRESCycle()
289 PetscCall(VecMAXPBY(VEC_TEMP, it_total + 1, avec, 0, &VEC_VV(0))); /* answer is in VEC_TEMP */ in KSPLGMRESCycle()
324 …PetscCall(KSPInitialResidual(ksp, ksp->vec_sol, VEC_TEMP, VEC_TEMP_MATOP, VEC_VV(0), ksp->vec_rhs)… in KSPSolve_LGMRES()
410 PetscCall(VecMAXPBY(VEC_TEMP, it + 1, nrs, 0, &VEC_VV(0))); in KSPLGMRESBuildSoln()
413 PetscCall(VecMAXPBY(VEC_TEMP, it_arnoldi, nrs, 0, &VEC_VV(0))); in KSPLGMRESBuildSoln()
H A Dlgmresimpl.h49 #define VEC_VV(i) lgmres->vecs[VEC_OFFSET + i] /* use to access othog basis vectors */ macro
/petsc/src/ksp/ksp/impls/gmres/fgmres/
H A Dfgmres.c52 PetscCall(VecWAXPY(VEC_VV(0), -1.0, VEC_TEMP, ksp->vec_rhs)); in KSPFGMRESResidual()
78 PetscCall(VecNorm(VEC_VV(0), NORM_2, &res_norm)); in KSPFGMRESCycle()
97 PetscCall(VecScale(VEC_VV(0), 1.0 / res_norm)); in KSPFGMRESCycle()
123 PetscCall(KSP_PCApply(ksp, VEC_VV(loc_it), PREVEC(loc_it))); in KSPFGMRESCycle()
127 PetscCall(KSP_MatMult(ksp, Amat, PREVEC(loc_it), VEC_VV(1 + loc_it))); in KSPFGMRESCycle()
134 PetscCall(VecNorm(VEC_VV(loc_it + 1), NORM_2, &tt)); in KSPFGMRESCycle()
146 PetscCall(VecScale(VEC_VV(loc_it + 1), 1.0 / tt)); in KSPFGMRESCycle()
227 PetscCall(VecCopy(ksp->vec_rhs, VEC_VV(0))); in KSPSolve_FGMRES()
230 PetscCall(VecFlag(VEC_VV(0), ksp->reason == KSP_DIVERGED_PC_FAILED)); in KSPSolve_FGMRES()
H A Dfgmresimpl.h37 #define VEC_VV(i) \ macro
/petsc/src/ksp/ksp/impls/gmres/pipefgmres/
H A Dpipefgmres.c84 PetscCall(VecNorm(VEC_VV(0), NORM_2, &res_norm)); in KSPPIPEFGMRESCycle()
105 PetscCall(VecScale(VEC_VV(0), 1.0 / res_norm)); in KSPPIPEFGMRESCycle()
108 PetscCall(KSP_PCApply(ksp, VEC_VV(loc_it), PREVEC(loc_it))); in KSPPIPEFGMRESCycle()
111 PetscCall(VecAXPY(ZVEC(loc_it), -shift, VEC_VV(loc_it))); /* Note shift */ in KSPPIPEFGMRESCycle()
136 for (i = 0; i < loc_it + 1; i++) redux[i] = VEC_VV(i); in KSPPIPEFGMRESCycle()
208 PetscCall(VecCopy(ZVEC(loc_it), VEC_VV(loc_it + 1))); in KSPPIPEFGMRESCycle()
209 PetscCall(VecMAXPY(VEC_VV(loc_it + 1), loc_it + 1, lhh, &VEC_VV(0))); in KSPPIPEFGMRESCycle()
232 PetscCall(VecScale(VEC_VV(loc_it + 1), 1.0 / tt)); in KSPPIPEFGMRESCycle()
321 …PetscCall(KSPInitialResidual(ksp, ksp->vec_sol, VEC_TEMP, VEC_TEMP_MATOP, VEC_VV(0), ksp->vec_rhs)… in KSPSolve_PIPEFGMRES()
H A Dpipefgmresimpl.h41 #define VEC_VV(i) pipefgmres->vecs[VEC_OFFSET + i] /* use to access othog basis vectors. Note the o… macro
/petsc/src/ksp/ksp/impls/gmres/dgmres/
H A Ddgmres.c115 PetscCall(VecNormalize(VEC_VV(0), &res_norm)); in KSPDGMRESCycle()
148 PetscCall(KSP_PCApplyBAorAB(ksp, VEC_VV(it), VEC_TEMP, VEC_TEMP_MATOP)); in KSPDGMRESCycle()
150 PetscCall(KSPDGMRESApplyDeflation(ksp, VEC_TEMP, VEC_VV(1 + it))); in KSPDGMRESCycle()
152 PetscCall(KSPDGMRESApplyDeflation(ksp, VEC_VV(it), VEC_TEMP)); in KSPDGMRESCycle()
153 PetscCall(KSP_PCApplyBAorAB(ksp, VEC_TEMP, VEC_VV(1 + it), VEC_TEMP_MATOP)); in KSPDGMRESCycle()
156 PetscCall(KSP_PCApplyBAorAB(ksp, VEC_VV(it), VEC_VV(1 + it), VEC_TEMP_MATOP)); in KSPDGMRESCycle()
163 PetscCall(VecNormalize(VEC_VV(it + 1), &tt)); in KSPDGMRESCycle()
238 …PetscCall(KSPInitialResidual(ksp, ksp->vec_sol, VEC_TEMP, VEC_TEMP_MATOP, VEC_VV(0), ksp->vec_rhs)… in KSPSolve_DGMRES()
242 PetscCall(KSPDGMRESApplyDeflation(ksp, VEC_VV(0), VEC_TEMP)); in KSPSolve_DGMRES()
243 PetscCall(VecCopy(VEC_TEMP, VEC_VV(0))); in KSPSolve_DGMRES()
[all …]
H A Ddgmresimpl.h63 #define VEC_VV(i) dgmres->vecs[VEC_OFFSET + i] macro