Home
last modified time | relevance | path

Searched refs:betas (Results 1 – 5 of 5) sorted by relevance

/petsc/src/ksp/ksp/impls/cheby/
H A Dbetas.c17 cheb->betas[0] = 1.12500000000000; in KSPChebyshevGetBetas_Private()
20 cheb->betas[0] = 1.02387287570313; in KSPChebyshevGetBetas_Private()
21 cheb->betas[1] = 1.26408905371085; in KSPChebyshevGetBetas_Private()
24 cheb->betas[0] = 1.00842544782028; in KSPChebyshevGetBetas_Private()
25 cheb->betas[1] = 1.08867839208730; in KSPChebyshevGetBetas_Private()
26 cheb->betas[2] = 1.33753125909618; in KSPChebyshevGetBetas_Private()
29 cheb->betas[0] = 1.00391310427285; in KSPChebyshevGetBetas_Private()
30 cheb->betas[1] = 1.04035811188593; in KSPChebyshevGetBetas_Private()
31 cheb->betas[2] = 1.14863498546254; in KSPChebyshevGetBetas_Private()
32 cheb->betas[3] = 1.38268869241000; in KSPChebyshevGetBetas_Private()
[all …]
H A Dchebyshevimpl.h13 …PetscReal *betas; /* store beta coefficients for 4th-kind Chebyshev smoothe… member
H A Dcheby.c566 PetscReal *betas = cheb->betas; in KSPSolve_Chebyshev_FourthKind() local
631 PetscCall(VecAXPBY(x, betas[i - 1], 1.0, d)); /* x = x + \beta_k d */ in KSPSolve_Chebyshev_FourthKind()
673 PetscCall(VecAXPBY(x, betas[ksp->max_it - 1], 1.0, d)); /* x = x + d */ in KSPSolve_Chebyshev_FourthKind()
768 PetscCall(PetscFree(cheb->betas)); in KSPSetUp_Chebyshev()
769 PetscCall(PetscMalloc1(ksp->max_it, &cheb->betas)); in KSPSetUp_Chebyshev()
774 for (PetscInt i = 0; i < ksp->max_it; i++) cheb->betas[i] = 1.0; in KSPSetUp_Chebyshev()
867 PetscCall(PetscFree(cheb->betas)); in KSPDestroy_Chebyshev()
/petsc/src/ksp/ksp/impls/fcg/pipefcg/
H A Dpipefcg.c84 PetscScalar alpha = 0.0, gamma, *betas, *dots; in KSPSolve_PIPEFCG_cycle() local
113 betas = dots; /* dots takes the result of all dot products of which the betas are a subset */ in KSPSolve_PIPEFCG_cycle()
219 …PetscCall(VecMXDotBegin(Z, j + 2, redux, betas)); /* Start split reduction… in KSPSolve_PIPEFCG_cycle()
225 …PetscCall(VecMXDotEnd(Z, j + 2, redux, betas)); /* Finish split reductio… in KSPSolve_PIPEFCG_cycle()
226 gamma = betas[j]; in KSPSolve_PIPEFCG_cycle()
227 delta = PetscRealPart(betas[j + 1]); in KSPSolve_PIPEFCG_cycle()
232 betas[j] /= -etas[kdx]; /* betak /= etak */ in KSPSolve_PIPEFCG_cycle()
233 *eta -= PetscAbsScalar(betas[j]) * PetscAbsScalar(betas[j]) * etas[kdx]; in KSPSolve_PIPEFCG_cycle()
250 … PetscCall(VecMAXPY(Pcurr, j, betas, pipefcg->Pold)); /* pi <- ui - sum_k beta_k p_k */ in KSPSolve_PIPEFCG_cycle()
251 … PetscCall(VecMAXPY(Scurr, j, betas, pipefcg->Sold)); /* si <- wi - sum_k beta_k s_k */ in KSPSolve_PIPEFCG_cycle()
[all …]
/petsc/src/ksp/ksp/impls/gcr/pipegcr/
H A Dpipegcr.c62 PetscScalar alpha = 0.0, gamma, *betas, *dots; in KSPSolve_PIPEGCR_cycle() local
86 betas = dots; /* dots takes the result of all dot products of which the betas are a subset */ in KSPSolve_PIPEGCR_cycle()
217 betas[j] /= -etas[kdx]; /* betak /= etak */ in KSPSolve_PIPEGCR_cycle()
218 *eta -= PetscAbsScalar(betas[j]) * PetscAbsScalar(betas[j]) * etas[kdx]; in KSPSolve_PIPEGCR_cycle()
238 PetscCall(VecMAXPY(t, j, betas, pipegcr->told)); /* ti <- n - sum_k beta_k t_k */ in KSPSolve_PIPEGCR_cycle()
240 PetscCall(VecMAXPY(p, j, betas, pipegcr->pold)); /* pi <- ui - sum_k beta_k p_k */ in KSPSolve_PIPEGCR_cycle()
241 PetscCall(VecMAXPY(s, j, betas, pipegcr->sold)); /* si <- wi - sum_k beta_k s_k */ in KSPSolve_PIPEGCR_cycle()
242 PetscCall(VecMAXPY(q, j, betas, pipegcr->qold)); /* qi <- m - sum_k beta_k q_k */ in KSPSolve_PIPEGCR_cycle()