Home
last modified time | relevance | path

Searched refs:beta (Results 1 – 25 of 171) sorted by relevance

1234567

/petsc/src/dm/dt/tests/
H A Dex1.c35 static PetscErrorCode CheckQuadrature_Basics(PetscInt npoints, PetscReal alpha, PetscReal beta, con… in CheckQuadrature_Basics() argument
41 … PetscInt_FMT ", x[i] = %g, x[i-1] = %g", npoints, (double)alpha, (double)beta, i, (double)x[i], (… in CheckQuadrature_Basics()
44 …ta = %g, i = %" PetscInt_FMT ", w[i] = %g", npoints, (double)alpha, (double)beta, i, (double)w[i]); in CheckQuadrature_Basics()
49 static PetscErrorCode CheckQuadrature(PetscInt npoints, PetscReal alpha, PetscReal beta, const Pets… in CheckQuadrature() argument
59 PetscCall(PetscDTJacobiEval(npoints, alpha, beta, x, 1, &i, Pi, NULL, NULL)); in CheckQuadrature()
64 PetscCall(PetscDTJacobiEval(npoints, alpha, beta, x, 1, &j, Pj, NULL, NULL)); in CheckQuadrature()
70 I_exact = PetscPowReal(2.0, alpha + beta + 1.) / (2. * i + alpha + beta + 1.); in CheckQuadrature()
72 …cExpReal(PetscLGamma(i + alpha + 1.) + PetscLGamma(i + beta + 1.) - (PetscLGamma(i + alpha + beta … in CheckQuadrature()
75 PetscInt ibeta = (PetscInt)beta; in CheckQuadrature()
77 …PetscCheck((PetscReal)ibeta == beta, PETSC_COMM_SELF, PETSC_ERR_SUP, "lgamma() - math routine is u… in CheckQuadrature()
[all …]
/petsc/src/snes/tutorials/
H A Dex18.c42 PetscReal beta, bm1, coef; /* nonlinear diffusivity parameterizations */ member
65 user.beta = 2.5; in main()
68 PetscCall(PetscOptionsGetReal(NULL, NULL, "-beta", &user.beta, NULL)); in main()
69 user.bm1 = user.beta - 1.0; in main()
70 user.coef = user.beta / 2.0; in main()
139 PetscScalar tleft, tright, beta; in FormFunction() local
154 beta = user->beta; in FormFunction()
173 dw = PetscPowScalar(aw, beta); in FormFunction()
178 de = PetscPowScalar(ae, beta); in FormFunction()
183 ds = PetscPowScalar(as, beta); in FormFunction()
[all …]
/petsc/src/snes/tests/
H A Dex20.c42 PetscReal beta, bm1, coef; /* nonlinear diffusivity parameterizations */ member
67 user.beta = 2.5; in main()
70 PetscCall(PetscOptionsGetReal(NULL, NULL, "-beta", &user.beta, NULL)); in main()
73 user.bm1 = user.beta - 1.0; in main()
74 user.coef = user.beta / 2.0; in main()
139 PetscScalar tleft, tright, beta, td, ad, dd, fd = 0.0, tu, au, du = 0.0, fu = 0.0; in FormFunction() local
156 beta = user->beta; in FormFunction()
176 dw = PetscPowScalar(aw, beta); in FormFunction()
181 de = PetscPowScalar(ae, beta); in FormFunction()
186 ds = PetscPowScalar(as, beta); in FormFunction()
[all …]
/petsc/src/mat/impls/baij/seq/baijmkl/
H A Dbaijmkl.h25 …kl_sparse_x_mv(operation, alpha, A, descr, x, beta, y) mkl_sparse_s_mv(operation, alpha, A, descr,… argument
27 …kl_sparse_x_mv(operation, alpha, A, descr, x, beta, y) mkl_sparse_d_mv(operation, alpha, A, descr,… argument
31 …kl_sparse_x_mv(operation, alpha, A, descr, x, beta, y) mkl_sparse_c_mv(operation, alpha, A, descr,… argument
33 …kl_sparse_x_mv(operation, alpha, A, descr, x, beta, y) mkl_sparse_z_mv(operation, alpha, A, descr,… argument
/petsc/src/ksp/ksp/impls/cg/
H A Dcg.c122 PetscScalar dpi = 0.0, a = 1.0, beta, betaold = 1.0, b = 0, *e = NULL, *d = NULL, dpiold; in KSPSolve_CG() local
180 PetscCall(VecXDot(Z, R, &beta)); /* beta <- z'*r */ in KSPSolve_CG()
181 KSPCheckDot(ksp, beta); in KSPSolve_CG()
182 dp = PetscSqrtReal(PetscAbsScalar(beta)); /* dp <- r'*z = r'*B*r = e'*A'*B*A*e */ in KSPSolve_CG()
215 PetscCall(VecXDot(Z, R, &beta)); /* beta <- z'*r */ in KSPSolve_CG()
216 KSPCheckDot(ksp, beta); in KSPSolve_CG()
222 if (beta == 0.0) { in KSPSolve_CG()
227 } else if ((i > 0) && (beta * betaold < 0.0)) { in KSPSolve_CG()
228 …to indefinite preconditioner, beta %g, betaold %g", (double)PetscRealPart(beta), (double)PetscReal… in KSPSolve_CG()
244 b = beta / betaold; in KSPSolve_CG()
[all …]
/petsc/src/mat/impls/aij/seq/aijmkl/
H A Daijmkl.h33 …, beta, y) mkl_scsrmv(transa, m, k, alpha, matdescra, (MatScalar *)val, (PetscInt *)indx, (PetscIn… argument
35 …, beta, y) mkl_dcsrmv(transa, m, k, alpha, matdescra, (MatScalar *)val, (PetscInt *)indx, (PetscIn… argument
39 …, beta, y) mkl_ccsrmv(transa, m, k, alpha, matdescra, (MatScalar *)val, (PetscInt *)indx, (PetscIn… argument
41 …, beta, y) mkl_zcsrmv(transa, m, k, alpha, matdescra, (MatScalar *)val, (PetscInt *)indx, (PetscIn… argument
/petsc/src/ksp/ksp/impls/symmlq/
H A Dsymmlq.c17 …PetscScalar alpha, beta, ibeta, betaold, beta1, ceta = 0, ceta_oold = 0.0, ceta_old = 0.0, ceta_ba… in KSPSolve_SYMMLQ() local
73 beta = dp; /* beta <- sqrt(r'*z) */ in KSPSolve_SYMMLQ()
74 beta1 = beta; in KSPSolve_SYMMLQ()
79 ibeta = 1.0 / beta; in KSPSolve_SYMMLQ()
104 PetscCall(VecScale(V, 1.0 / beta)); /* v <- ibeta*r; */ in KSPSolve_SYMMLQ()
106 PetscCall(VecScale(U, 1.0 / beta)); /* u <- ibeta*z; */ in KSPSolve_SYMMLQ()
126 PetscCall(VecAXPY(R, -beta, VOLD)); /* r <- r - beta * v_old; */ in KSPSolve_SYMMLQ()
127 PetscCall(VecAXPY(Z, -beta, UOLD)); /* z <- z - beta * u_old; */ in KSPSolve_SYMMLQ()
128 betaold = beta; /* beta_k */ in KSPSolve_SYMMLQ()
142 beta = PetscSqrtScalar(dp); /* beta = sqrt(dp); */ in KSPSolve_SYMMLQ()
[all …]
/petsc/src/tao/bound/impls/bncg/
H A Dbncg.c439 PetscReal gamma = 1.0, tau_k, beta; in TaoBNCGStepDirectionUpdate() local
521 beta = tau_k * gkp1_yk / dk_yk; in TaoBNCGStepDirectionUpdate()
522 PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient)); in TaoBNCGStepDirectionUpdate()
526 beta = gkp1_yk / dk_yk; in TaoBNCGStepDirectionUpdate()
527 PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work)); in TaoBNCGStepDirectionUpdate()
539 beta = tau_k * gnorm2 / gnorm2_old; in TaoBNCGStepDirectionUpdate()
540 PetscCall(VecAXPBY(tao->stepdirection, -tau_k, beta, tao->gradient)); in TaoBNCGStepDirectionUpdate()
545 beta = tmp / gnorm2_old; in TaoBNCGStepDirectionUpdate()
546 PetscCall(VecAXPBY(tao->stepdirection, -1.0, beta, cg->g_work)); in TaoBNCGStepDirectionUpdate()
556 beta = tau_k * gkp1_yk / gnorm2_old; in TaoBNCGStepDirectionUpdate()
[all …]
/petsc/lib/petsc/bin/
H A DtasClasses.py63 row[item] = '{:.3g}'.format(field.beta)
74 def __init__(self, fileName, fieldName, alpha=0, cRate=0, beta=0): argument
80 self.beta = beta
85 def setBeta(self, beta): argument
86 self.beta = beta
/petsc/src/ksp/ksp/impls/bcgs/pipebcgs/
H A Dpipebcgs.c18 PetscScalar rho, rhoold, alpha, beta, omega = 0.0, d1, d2, d3; in KSPSolve_PIPEBCGS() local
88 beta = 0.0; in KSPSolve_PIPEBCGS()
99 …PetscCall(VecAXPBYPCZ(P2, 1.0, -beta * omega, beta, R2, S2)); /* p2 <- beta * p2 + r2 - beta * ome… in KSPSolve_PIPEBCGS()
100 …PetscCall(VecAXPBYPCZ(S, 1.0, -beta * omega, beta, W, Z)); /* s <- beta * s + w - beta * ome… in KSPSolve_PIPEBCGS()
101 …PetscCall(VecAXPBYPCZ(S2, 1.0, -beta * omega, beta, W2, Z2)); /* s2 <- beta * s2 + w2 - beta * ome… in KSPSolve_PIPEBCGS()
102 …PetscCall(VecAXPBYPCZ(Z, 1.0, -beta * omega, beta, T, V)); /* z <- beta * z + t - beta * ome… in KSPSolve_PIPEBCGS()
160 …PetscCheck(d2 + beta * d1 - beta * omega * d3 != 0.0, PetscObjectComm((PetscObject)ksp), PETSC_ERR… in KSPSolve_PIPEBCGS()
162 beta = (rho / rhoold) * (alpha / omega); in KSPSolve_PIPEBCGS()
163 …alpha = rho / (d2 + beta * d1 - beta * omega * d3); /* alpha <- rho / (d2 + beta * d1 - beta * ome… in KSPSolve_PIPEBCGS()
/petsc/src/ts/tutorials/
H A Dex42.c17 PetscReal beta; member
30 PetscReal alpha, beta; in RHSFunction() local
40 beta = user->beta; in RHSFunction()
55 da = alpha * a * a / (1. + beta * h) + rho_a - mu_a * a; in RHSFunction()
79 PetscReal alpha, beta; in RHSJacobian() local
90 beta = user->beta; in RHSJacobian()
103 va[0] = 2 * alpha * a / (1. + beta * h) - mu_a; in RHSJacobian()
106 va[1] = -beta * alpha * a * a / ((1. + beta * h) * (1. + beta * h)); in RHSJacobian()
215 user.beta = 1.; in main()
226 PetscCall(PetscOptionsReal("-beta", "Inhibition factor", "ex42.c", user.beta, &user.beta, NULL)); in main()
[all …]
/petsc/src/ksp/ksp/impls/tcqmr/
H A Dtcqmr.c14 PetscScalar deltmp, rho, beta, eptmp, ta, s, c, tau_n, delta; in KSPSolve_TCQMR() local
70 beta = rho / (f * rhom1); in KSPSolve_TCQMR()
75 PetscCall(VecAXPY(utmp, -2.0 * beta, p)); in KSPSolve_TCQMR()
78 PetscCall(VecAXPY(up1, f * beta * beta, um1)); in KSPSolve_TCQMR()
83 PetscCall(VecAYPX(p, -beta, z)); /* p = f*(z-beta*p) */ in KSPSolve_TCQMR()
87 beta = beta / Gamma; in KSPSolve_TCQMR()
88 eptmp = beta; in KSPSolve_TCQMR()
91 PetscCall(VecAXPY(vp1, -beta, vm1)); in KSPSolve_TCQMR()
103 theta = sl1 * beta; in KSPSolve_TCQMR()
104 eptmp = -cl1 * beta; in KSPSolve_TCQMR()
/petsc/src/ksp/ksp/utils/lmvm/blas_cyclic/
H A Dblas_cyclic.c16 …Private(PetscInt m, PetscScalar alpha, const PetscScalar x[], PetscScalar beta, PetscScalar y[], P… in AXPBY_Private() argument
18 for (PetscInt i = 0; i < m; i++) y[i * y_stride] = alpha * x[i] + beta * y[i * y_stride]; in AXPBY_Private()
21 …est, PetscInt next, PetscScalar alpha, const PetscScalar x[], PetscScalar beta, PetscScalar y[], P… in AXPBYCylic_Private() argument
28 AXPBY_Private(m, alpha, x, beta, y, y_stride); in AXPBYCylic_Private()
30 AXPBY_Private(i_next - i_oldest, alpha, &x[i_oldest], beta, &y[i_oldest * y_stride], y_stride); in AXPBYCylic_Private()
32 AXPBY_Private(i_next, alpha, x, beta, y, y_stride); in AXPBYCylic_Private()
33 AXPBY_Private(m - i_oldest, alpha, &x[i_oldest], beta, &y[i_oldest * y_stride], y_stride); in AXPBYCylic_Private()
38 …e VecAXPBYCyclic(PetscInt oldest, PetscInt next, PetscScalar alpha, Vec x, PetscScalar beta, Vec y) in VecAXPBYCyclic() argument
61 if (m_local == m) PetscCall(AXPBYCyclic_CUPM_Private(m, oldest, next, alpha, x_, beta, y_, 1)); in VecAXPBYCyclic()
71 PetscCall(AXPBYCylic_Private(m, oldest, next, alpha, x_, beta, y_, 1)); in VecAXPBYCyclic()
[all …]
/petsc/src/ksp/ksp/utils/lmvm/blas_cyclic/cupm/
H A Dblas_cyclic_cupm.cxx30 …est, PetscInt next, PetscScalar alpha, const PetscScalar x[], PetscScalar beta, PetscScalar y[], P… in AXPBYCyclic_CUPM_Private() argument
41 …Petsc::device::cupm::DeviceType::CUDA>::axpby(dctx, m, oldest, next, alpha, x, beta, y, y_stride)); in AXPBYCyclic_CUPM_Private()
46 …:Petsc::device::cupm::DeviceType::HIP>::axpby(dctx, m, oldest, next, alpha, x, beta, y, y_stride)); in AXPBYCyclic_CUPM_Private()
55 …PetscScalar alpha, const PetscScalar A[], const PetscScalar x[], PetscScalar beta, PetscScalar y[]) in DMVCyclic_CUPM_Private() argument
66 …e::cupm::DeviceType::CUDA>::dmv(dctx, hermitian_transpose, m, oldest, next, alpha, A, x, beta, y)); in DMVCyclic_CUPM_Private()
71 …ce::cupm::DeviceType::HIP>::dmv(dctx, hermitian_transpose, m, oldest, next, alpha, A, x, beta, y)); in DMVCyclic_CUPM_Private()
130 …pha, const PetscScalar A[], PetscInt lda, const PetscScalar x[], PetscScalar beta, PetscScalar y[]) in HEMVCyclic_CUPM_Private() argument
141 …<::Petsc::device::cupm::DeviceType::CUDA>::hemv(dctx, m, oldest, next, alpha, A, lda, x, beta, y)); in HEMVCyclic_CUPM_Private()
146 …c<::Petsc::device::cupm::DeviceType::HIP>::hemv(dctx, m, oldest, next, alpha, A, lda, x, beta, y)); in HEMVCyclic_CUPM_Private()
155 …pha, const PetscScalar A[], PetscInt lda, const PetscScalar x[], PetscScalar beta, PetscScalar y[]) in GEMVCyclic_CUPM_Private() argument
[all …]
/petsc/src/ksp/ksp/impls/cg/cgne/
H A Dcgne.c44 PetscScalar dpi, a = 1.0, beta, betaold = 1.0, b = 0, *e = NULL, *d = NULL; in KSPSolve_CGNE() local
96 PetscCall(VecXDot(Z, R, &beta)); in KSPSolve_CGNE()
97 KSPCheckDot(ksp, beta); in KSPSolve_CGNE()
98 dp = PetscSqrtReal(PetscAbsScalar(beta)); in KSPSolve_CGNE()
109 PetscCall(VecXDot(Z, R, &beta)); /* beta <- r'z */ in KSPSolve_CGNE()
110 KSPCheckDot(ksp, beta); in KSPSolve_CGNE()
111 if (beta == 0.0) { in KSPSolve_CGNE()
116 } else if (beta < 0.0) { in KSPSolve_CGNE()
126 b = beta / betaold; in KSPSolve_CGNE()
133 betaold = beta; in KSPSolve_CGNE()
[all …]
/petsc/src/ksp/ksp/utils/lmvm/
H A Dlmbasis.c262 …sGEMV(LMBasis A, PetscInt oldest, PetscInt next, PetscScalar alpha, Vec x, PetscScalar beta, Vec y) in LMBasisGEMV() argument
279 if (beta != 1.0 && beta != 0.0) PetscCall(VecScale(y, beta)); in LMBasisGEMV()
282 if (beta == 0.0) PetscCall(MatMult(A->vecs, x_, y)); in LMBasisGEMV()
286 if (beta == 0.0) PetscCall(MatMultColumnRange(A->vecs, x_, y, oldest_idx, next_idx)); in LMBasisGEMV()
289 if (beta == 0.0) PetscCall(MatMultColumnRange(A->vecs, x_, y, 0, next_idx)); in LMBasisGEMV()
298 …GEMVH(LMBasis A, PetscInt oldest, PetscInt next, PetscScalar alpha, Vec x, PetscScalar beta, Vec y) in LMBasisGEMVH() argument
316 PetscCall(VecAXPBYCyclic(oldest, next, alpha, A->cached_product, beta, y)); in LMBasisGEMVH()
321 if (alpha != 1.0 || (beta != 1.0 && beta != 0.0)) PetscCall(LMBasisGetWorkRow(A, &y_)); in LMBasisGEMVH()
324 if (alpha == 1.0 && beta == 1.0) PetscCall(MatMultHermitianTransposeAdd(A->vecs, x, y_, y_)); in LMBasisGEMVH()
328 …if (alpha == 1.0 && beta == 1.0) PetscCall(MatMultHermitianTransposeAddColumnRange(A->vecs, x, y_,… in LMBasisGEMVH()
[all …]
/petsc/src/ksp/ksp/tests/output/
H A Dex5_mat_type-lmvmbfgs.out21 Rescale params: alpha=1., beta=0.5, rho=1.
38 Rescale params: alpha=1., beta=0.5, rho=1.
55 Rescale params: alpha=1., beta=0.5, rho=1.
72 Rescale params: alpha=1., beta=0.5, rho=1.
89 Rescale params: alpha=1., beta=0.5, rho=1.
106 Rescale params: alpha=1., beta=0.5, rho=1.
123 Rescale params: alpha=1., beta=0.5, rho=1.
140 Rescale params: alpha=1., beta=0.5, rho=1.
157 Rescale params: alpha=1., beta=0.5, rho=1.
174 Rescale params: alpha=1., beta=0.5, rho=1.
[all …]
H A Dex5_mat_type-lmvmsymbroyden.out22 Rescale params: alpha=1., beta=0.5, rho=1.
40 Rescale params: alpha=1., beta=0.5, rho=1.
58 Rescale params: alpha=1., beta=0.5, rho=1.
76 Rescale params: alpha=1., beta=0.5, rho=1.
94 Rescale params: alpha=1., beta=0.5, rho=1.
112 Rescale params: alpha=1., beta=0.5, rho=1.
130 Rescale params: alpha=1., beta=0.5, rho=1.
148 Rescale params: alpha=1., beta=0.5, rho=1.
166 Rescale params: alpha=1., beta=0.5, rho=1.
184 Rescale params: alpha=1., beta=0.5, rho=1.
[all …]
H A Dex5_mat_type-lmvmdfp.out21 Rescale params: alpha=1., beta=0.5, rho=1.
38 Rescale params: alpha=1., beta=0.5, rho=1.
55 Rescale params: alpha=1., beta=0.5, rho=1.
72 Rescale params: alpha=1., beta=0.5, rho=1.
89 Rescale params: alpha=1., beta=0.5, rho=1.
106 Rescale params: alpha=1., beta=0.5, rho=1.
123 Rescale params: alpha=1., beta=0.5, rho=1.
140 Rescale params: alpha=1., beta=0.5, rho=1.
157 Rescale params: alpha=1., beta=0.5, rho=1.
174 Rescale params: alpha=1., beta=0.5, rho=1.
[all …]
H A Dex5_mat_type-lmvmdiagbroyden.out21 Rescale params: alpha=1., beta=0.5, rho=1.
38 Rescale params: alpha=1., beta=0.5, rho=1.
55 Rescale params: alpha=1., beta=0.5, rho=1.
72 Rescale params: alpha=1., beta=0.5, rho=1.
89 Rescale params: alpha=1., beta=0.5, rho=1.
106 Rescale params: alpha=1., beta=0.5, rho=1.
123 Rescale params: alpha=1., beta=0.5, rho=1.
140 Rescale params: alpha=1., beta=0.5, rho=1.
157 Rescale params: alpha=1., beta=0.5, rho=1.
174 Rescale params: alpha=1., beta=0.5, rho=1.
[all …]
H A Dex5_mat_type-lmvmsymbadbroyden.out22 Rescale params: alpha=1., beta=0.5, rho=1.
40 Rescale params: alpha=1., beta=0.5, rho=1.
58 Rescale params: alpha=1., beta=0.5, rho=1.
76 Rescale params: alpha=1., beta=0.5, rho=1.
94 Rescale params: alpha=1., beta=0.5, rho=1.
112 Rescale params: alpha=1., beta=0.5, rho=1.
130 Rescale params: alpha=1., beta=0.5, rho=1.
148 Rescale params: alpha=1., beta=0.5, rho=1.
166 Rescale params: alpha=1., beta=0.5, rho=1.
184 Rescale params: alpha=1., beta=0.5, rho=1.
[all …]
/petsc/src/sys/utils/
H A Dmathfit.c23 PetscReal *X, *Y, beta[2]; in PetscLinearRegression() local
64 beta[i] = 0.0; in PetscLinearRegression()
65 for (PetscInt k = 0; k < n; ++k) beta[i] += Y[i * n + k] * y[k]; in PetscLinearRegression()
68 *intercept = beta[0]; in PetscLinearRegression()
69 *slope = beta[1]; in PetscLinearRegression()
/petsc/src/snes/impls/ngmres/
H A Dngmresfunc.c44 PetscScalar *beta = ngmres->beta; in SNESNGMRESFormCombinedSolution_Private() local
57 PetscCall(VecMDotBegin(Fdot[ivec], l, Fdot, beta)); in SNESNGMRESFormCombinedSolution_Private()
59 PetscCall(VecMDotEnd(Fdot[ivec], l, Fdot, beta)); in SNESNGMRESFormCombinedSolution_Private()
61 Q(i, ivec) = beta[i]; in SNESNGMRESFormCombinedSolution_Private()
62 Q(ivec, i) = beta[i]; in SNESNGMRESFormCombinedSolution_Private()
68 for (i = 0; i < l; i++) beta[i] = nu - xi[i]; in SNESNGMRESFormCombinedSolution_Private()
76 if (H(0, 0) != 0.) beta[0] = beta[0] / H(0, 0); in SNESNGMRESFormCombinedSolution_Private()
77 else beta[0] = 0.; in SNESNGMRESFormCombinedSolution_Private()
85 …_(&ngmres->m, &ngmres->n, &ngmres->nrhs, ngmres->h, &ngmres->lda, ngmres->beta, &ngmres->ldb, ngmr… in SNESNGMRESFormCombinedSolution_Private()
87 …_(&ngmres->m, &ngmres->n, &ngmres->nrhs, ngmres->h, &ngmres->lda, ngmres->beta, &ngmres->ldb, ngmr… in SNESNGMRESFormCombinedSolution_Private()
[all …]
/petsc/src/ksp/ksp/tests/
H A Dex38.c33 PetscReal h, gamma, beta; /* parameters */ in main() local
47 beta = 0.01; in main()
50 PetscCall(PetscOptionsGetReal(NULL, NULL, "-beta", &beta, NULL)); in main()
52 beta = beta / (h * h); in main()
91 co2 = beta * h * h; in main()
/petsc/src/ksp/ksp/impls/gmres/agmres/
H A Dagmresdeflation.c100 PetscScalar *beta = agmres->beta; in KSPAGMRESSchurForm() local
127 …hgeqz_("S", "I", "I", &KspSize, &ilo, &KspSize, A, &ldA, B, &ldB, wr, wi, beta, Q, &N, Z, &N, work… in KSPAGMRESSchurForm()
130 …PACKgges_("V", "V", "N", NULL, &KspSize, A, &ldA, B, &ldB, &sdim, wr, wi, beta, Q, &N, Z, &N, work… in KSPAGMRESSchurForm()
136 if (beta[i] != 0.0) { in KSPAGMRESSchurForm()
137 wr[i] /= beta[i]; in KSPAGMRESSchurForm()
138 wi[i] /= beta[i]; in KSPAGMRESSchurForm()
158 …tgsen_(&ijob, &wantQ, &wantZ, select, &KspSize, A, &ldA, B, &ldB, wr, wi, beta, Q, &N, Z, &N, &r, … in KSPAGMRESSchurForm()
187 PetscScalar alpha, beta; in KSPAGMRESComputeDeflationData() local
203 beta = 0.0; in KSPAGMRESComputeDeflationData()
207 …&bKspSize, &blC, &alpha, agmres->hes_origin, &blC, agmres->hes_origin, &blC, &beta, MatEigL, &bN)); in KSPAGMRESComputeDeflationData()

1234567