Home
last modified time | relevance | path

Searched refs:alpha (Results 1 – 25 of 440) sorted by relevance

12345678910>>...18

/petsc/src/binding/petsc4py/demo/legacy/taosolve/
H A Drosenbrock.py28 def __init__(self, n=2, alpha=99.0): argument
30 self.alpha = float(alpha)
34 alpha = self.alpha
40 ff += alpha*t1*t1 + t2*t2;
45 alpha = self.alpha
51 G[2*i] = -4*alpha*t1*x[2*i] - 2*t2;
52 G[2*i+1] = 2*alpha*t1;
56 alpha = self.alpha
63 ff += alpha*t1*t1 + t2*t2;
64 G[2*i] = -4*alpha*t1*x[2*i] - 2*t2;
[all …]
/petsc/src/tao/unconstrained/tutorials/
H A Drosenbrock3.c20 PetscReal alpha; /* condition parameter */ member
48 user.alpha = 99.0; in main()
52 PetscCall(PetscOptionsGetReal(NULL, NULL, "-alpha", &user.alpha, &flg)); in main()
137 PetscReal ff = 0, t1, t2, alpha = user->alpha; in FormFunctionGradient() local
151 ff += PetscSqr(1 - x[i]) + alpha * t1 * t1; in FormFunctionGradient()
152 g[i] += -2 * (1 - x[i]) + 2 * alpha * t1 * (-2 * x[i]); in FormFunctionGradient()
153 g[i + 1] = 2 * alpha * t1; in FormFunctionGradient()
159 ff += alpha * t1 * t1 + t2 * t2; in FormFunctionGradient()
160 g[2 * i] = -4 * alpha * t1 * x[2 * i] - 2.0 * t2; in FormFunctionGradient()
161 g[2 * i + 1] = 2 * alpha * t1; in FormFunctionGradient()
[all …]
H A Drosenbrock2.c20 PetscReal alpha; /* condition parameter */ member
48 user.alpha = 99.0; in main()
52 PetscCall(PetscOptionsGetReal(NULL, NULL, "-alpha", &user.alpha, &flg)); in main()
132 PetscReal ff = 0, t1, t2, alpha = user->alpha; in FormFunctionGradient() local
146 ff += PetscSqr(1 - x[i]) + alpha * t1 * t1; in FormFunctionGradient()
147 g[i] += -2 * (1 - x[i]) + 2 * alpha * t1 * (-2 * x[i]); in FormFunctionGradient()
148 g[i + 1] = 2 * alpha * t1; in FormFunctionGradient()
154 ff += alpha * t1 * t1 + t2 * t2; in FormFunctionGradient()
155 g[2 * i] = -4 * alpha * t1 * x[2 * i] - 2.0 * t2; in FormFunctionGradient()
156 g[2 * i + 1] = 2 * alpha * t1; in FormFunctionGradient()
[all …]
H A Drosenbrock1.c20 PetscReal alpha; /* condition parameter */ member
51 user.alpha = 99.0; in main()
55 PetscCall(PetscOptionsGetReal(NULL, NULL, "-alpha", &user.alpha, &flg)); in main()
141 PetscReal ff = 0, t1, t2, alpha = user->alpha; in FormFunctionGradient() local
155 ff += PetscSqr(1 - x[i]) + alpha * t1 * t1; in FormFunctionGradient()
156 g[i] += -2 * (1 - x[i]) + 2 * alpha * t1 * (-2 * x[i]); in FormFunctionGradient()
157 g[i + 1] = 2 * alpha * t1; in FormFunctionGradient()
163 ff += alpha * t1 * t1 + t2 * t2; in FormFunctionGradient()
164 g[2 * i] = -4 * alpha * t1 * x[2 * i] - 2.0 * t2; in FormFunctionGradient()
165 g[2 * i + 1] = 2 * alpha * t1; in FormFunctionGradient()
[all …]
H A Drosenbrock1f.F9042 PetscReal alpha
44 common/params/alpha, n
58 ff = ff + alpha*t1*t1 + t2*t2
59 g_v(1 + 2*i) = -4*alpha*t1*x_v(1 + 2*i) - 2.0*t2
60 g_v(1 + 2*i + 1) = 2.0*alpha*t1
110 PetscReal alpha
112 common/params/alpha, n
129 v(1, 1) = 2.0*alpha
130 v(0, 0) = -4.0*alpha*(x_v(1 + 2*i + 1) - 3*x_v(1 + 2*i)*x_v(1 + 2*i)) + 2
131 v(1, 0) = -4.0*alpha*x_v(1 + 2*i)
[all …]
/petsc/src/dm/impls/swarm/tests/
H A Dex8.c79 PetscReal alpha, gmin, gmax; in TestDistribution() local
109 PetscCall(PetscProbComputeKSStatisticWeighted(locx, locw, PetscCDFConstant1D, &alpha)); in TestDistribution()
110 …alpha < confidenceLevel) PetscCall(PetscPrintf(comm, "The KS test for X rejects the null hypothesi… in TestDistribution()
111 …or X accepts the null hypothesis at level %.2g (%.2g)\n", (double)confidenceLevel, (double)alpha)); in TestDistribution()
112 PetscCall(PetscProbComputeKSStatisticWeighted(locv, locw, PetscCDFGaussian1D, &alpha)); in TestDistribution()
113 …alpha < confidenceLevel) PetscCall(PetscPrintf(comm, "The KS test for V rejects the null hypothesi… in TestDistribution()
114 …or V accepts the null hypothesis at level %.2g (%.2g)\n", (double)confidenceLevel, (double)alpha)); in TestDistribution()
115 PetscCall(PetscProbComputeKSStatisticMagnitude(locv, cdf, &alpha)); in TestDistribution()
116 …alpha < confidenceLevel) PetscCall(PetscPrintf(comm, "The KS test for |V| rejects the null hypothe… in TestDistribution()
117 … |V| accepts the null hypothesis at level %.2g (%.2g)\n", (double)confidenceLevel, (double)alpha)); in TestDistribution()
/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 …ta = %g, i = %" PetscInt_FMT ", x[i] = %g, x[i-1] = %g", npoints, (double)alpha, (double)beta, i, … in CheckQuadrature_Basics()
44 …lpha = %g, beta = %g, i = %" PetscInt_FMT ", w[i] = %g", npoints, (double)alpha, (double)beta, 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 …I_exact *= PetscExpReal(PetscLGamma(i + alpha + 1.) + PetscLGamma(i + beta + 1.) - (PetscLGamma(i … in CheckQuadrature()
78 for (k = 0; k < ibeta; k++) I_exact *= (i + 1. + k) / (i + alpha + 1. + k); in CheckQuadrature()
82 PetscCall(PetscDTJacobiNorm(alpha, beta, i, &norm)); in CheckQuadrature()
[all …]
/petsc/src/mat/impls/baij/seq/baijmkl/
H A Dbaijmkl.h25 …#define mkl_sparse_x_mv(operation, alpha, A, descr, x, beta, y) mkl_sparse_s_mv(operation, alpha, … argument
27 …#define mkl_sparse_x_mv(operation, alpha, A, descr, x, beta, y) mkl_sparse_d_mv(operation, alpha, … argument
31 …#define mkl_sparse_x_mv(operation, alpha, A, descr, x, beta, y) mkl_sparse_c_mv(operation, alpha, … argument
33 …#define mkl_sparse_x_mv(operation, alpha, A, descr, x, beta, y) mkl_sparse_z_mv(operation, alpha, … argument
/petsc/src/mat/tests/output/
H A Dex2_24.out17 MatAXPY: B = B + alpha * A
25 MatAYPX: B = alpha*B + A
33 MatAXPY: C = C + alpha * A, C=A, SAME_NONZERO_PATTERN
42 MatAXPY: (C^T)^T = (C^T)^T + alpha * A, C=A, SAME_NONZERO_PATTERN
50 MatAXPY: C = C + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN
58 MatAXPY: (C^T)^T = (C^T)^T + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN
75 MatAXPY: (C^H)^H = (C^H)^H + alpha * A, C=A, SAME_NONZERO_PATTERN
83 MatAXPY: C = C + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN
91 MatAXPY: (C^H)^H = (C^H)^H + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN
123 MatAXPY: B = B + alpha * A, SUBSET_NONZERO_PATTERN
H A Dex2_21.out17 MatAXPY: B = B + alpha * A
25 MatAYPX: B = alpha*B + A
33 MatAXPY: C = C + alpha * A, C=A, SAME_NONZERO_PATTERN
42 MatAXPY: (C^T)^T = (C^T)^T + alpha * A, C=A, SAME_NONZERO_PATTERN
50 MatAXPY: C = C + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN
58 MatAXPY: (C^T)^T = (C^T)^T + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN
75 MatAXPY: (C^H)^H = (C^H)^H + alpha * A, C=A, SAME_NONZERO_PATTERN
83 MatAXPY: C = C + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN
91 MatAXPY: (C^H)^H = (C^H)^H + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN
123 MatAXPY: B = B + alpha * A, SUBSET_NONZERO_PATTERN
H A Dex2_22.out17 MatAXPY: B = B + alpha * A
25 MatAYPX: B = alpha*B + A
33 MatAXPY: C = C + alpha * A, C=A, SAME_NONZERO_PATTERN
42 MatAXPY: (C^T)^T = (C^T)^T + alpha * A, C=A, SAME_NONZERO_PATTERN
50 MatAXPY: C = C + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN
58 MatAXPY: (C^T)^T = (C^T)^T + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN
75 MatAXPY: (C^H)^H = (C^H)^H + alpha * A, C=A, SAME_NONZERO_PATTERN
83 MatAXPY: C = C + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN
91 MatAXPY: (C^H)^H = (C^H)^H + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN
123 MatAXPY: B = B + alpha * A, SUBSET_NONZERO_PATTERN
H A Dex2_23.out17 MatAXPY: B = B + alpha * A
25 MatAYPX: B = alpha*B + A
33 MatAXPY: C = C + alpha * A, C=A, SAME_NONZERO_PATTERN
42 MatAXPY: (C^T)^T = (C^T)^T + alpha * A, C=A, SAME_NONZERO_PATTERN
50 MatAXPY: C = C + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN
58 MatAXPY: (C^T)^T = (C^T)^T + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN
75 MatAXPY: (C^H)^H = (C^H)^H + alpha * A, C=A, SAME_NONZERO_PATTERN
83 MatAXPY: C = C + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN
91 MatAXPY: (C^H)^H = (C^H)^H + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN
123 MatAXPY: B = B + alpha * A, SUBSET_NONZERO_PATTERN
H A Dex2_12_B.out15 MatAXPY: C = C + alpha * A, C=A, SAME_NONZERO_PATTERN
24 MatAXPY: (C^T)^T = (C^T)^T + alpha * A, C=A, SAME_NONZERO_PATTERN
32 MatAXPY: C = C + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN
40 MatAXPY: (C^T)^T = (C^T)^T + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN
57 MatAXPY: (C^H)^H = (C^H)^H + alpha * A, C=A, SAME_NONZERO_PATTERN
65 MatAXPY: C = C + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN
73 MatAXPY: (C^H)^H = (C^H)^H + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN
105 MatAXPY: B = B + alpha * A, SUBSET_NONZERO_PATTERN
H A Dex2_4.out19 MatAXPY: C = C + alpha * A, C=A, SAME_NONZERO_PATTERN
28 MatAXPY: (C^T)^T = (C^T)^T + alpha * A, C=A, SAME_NONZERO_PATTERN
36 MatAXPY: C = C + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN
44 MatAXPY: (C^T)^T = (C^T)^T + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN
61 MatAXPY: (C^H)^H = (C^H)^H + alpha * A, C=A, SAME_NONZERO_PATTERN
69 MatAXPY: C = C + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN
77 MatAXPY: (C^H)^H = (C^H)^H + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN
109 MatAXPY: B = B + alpha * A, SUBSET_NONZERO_PATTERN
H A Dex2_12_A.out19 MatAXPY: C = C + alpha * A, C=A, SAME_NONZERO_PATTERN
28 MatAXPY: (C^T)^T = (C^T)^T + alpha * A, C=A, SAME_NONZERO_PATTERN
36 MatAXPY: C = C + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN
44 MatAXPY: (C^T)^T = (C^T)^T + alpha * (A^T)^T, C=A, SAME_NONZERO_PATTERN
61 MatAXPY: (C^H)^H = (C^H)^H + alpha * A, C=A, SAME_NONZERO_PATTERN
69 MatAXPY: C = C + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN
77 MatAXPY: (C^H)^H = (C^H)^H + alpha * (A^H)^H, C=A, SAME_NONZERO_PATTERN
109 MatAXPY: B = B + alpha * A, SUBSET_NONZERO_PATTERN
/petsc/src/binding/petsc4py/demo/legacy/bratu2d/
H A Dbratu2d.py12 def __init__(self, nx, ny, alpha, impl='python'): argument
15 self.alpha = alpha
29 alpha = self.alpha
33 self.compute(alpha, x, f)
41 alpha = OptDB.getReal('alpha', 6.8) variable
46 appc = Bratu2D(nx, ny, alpha, impl)
/petsc/src/vec/vec/impls/seq/
H A Dbvec1.c51 PetscErrorCode VecScale_Seq(Vec xin, PetscScalar alpha) in VecScale_Seq() argument
54 if (alpha == (PetscScalar)0.0) { in VecScale_Seq()
55 PetscCall(VecSet_Seq(xin, alpha)); in VecScale_Seq()
56 } else if (alpha != (PetscScalar)1.0) { in VecScale_Seq()
64 PetscCallBLAS("BLASscal", BLASscal_(&bn, &alpha, xarray, &one)); in VecScale_Seq()
70 PetscErrorCode VecAXPY_Seq(Vec yin, PetscScalar alpha, Vec xin) in VecAXPY_Seq() argument
74 if (alpha != (PetscScalar)0.0) { in VecAXPY_Seq()
84 PetscCallBLAS("BLASaxpy", BLASaxpy_(&bn, &alpha, xarray, &one, yarray, &one)); in VecAXPY_Seq()
120 PetscErrorCode VecAXPBYPCZ_Seq(Vec zin, PetscScalar alpha, PetscScalar beta, PetscScalar gamma, Vec… in VecAXPBYPCZ_Seq() argument
131 if (alpha == (PetscScalar)1.0) { in VecAXPBYPCZ_Seq()
[all …]
/petsc/src/ts/tutorials/
H A Dex53.c49 PetscScalar alpha; /* Biot effective stress coefficient */ member
151 const PetscReal alpha = PetscRealPart(constants[2]); in f0_quadratic_linear_u() local
153 const PetscReal K_d = K_u - alpha * alpha * M; in f0_quadratic_linear_u()
157 for (d = 0; d < dim - 1; ++d) f0[d] -= 2.0 * G - alpha * t; in f0_quadratic_linear_u()
158 f0[dim - 1] -= 2.0 * lambda + 4.0 * G - alpha * t; in f0_quadratic_linear_u()
163 const PetscReal alpha = PetscRealPart(constants[2]); in f0_quadratic_linear_p() local
169 f0[0] += u_t ? alpha * u_t[uOff[1]] : 0.0; in f0_quadratic_linear_p()
233 const PetscReal alpha = PetscRealPart(constants[2]); in f0_quadratic_trig_u() local
235 const PetscReal K_d = K_u - alpha * alpha * M; in f0_quadratic_trig_u()
239 for (d = 0; d < dim - 1; ++d) f0[d] -= 2.0 * G - alpha * PetscCosReal(t); in f0_quadratic_trig_u()
[all …]
/petsc/src/ksp/ksp/impls/cg/
H A Dcgls.c32 PetscReal alpha, gamma, oldgamma; in KSPSolve_CGLS() local
63 PetscCall(VecNorm(q, NORM_2, &alpha)); in KSPSolve_CGLS()
64 KSPCheckNorm(ksp, alpha); in KSPSolve_CGLS()
65 alpha = alpha * alpha; /* alpha = norm2(q)^2 */ in KSPSolve_CGLS()
66 alpha = gamma / alpha; /* alpha = gamma / alpha */ in KSPSolve_CGLS()
67 PetscCall(VecAXPY(x, alpha, p)); /* x += alpha * p */ in KSPSolve_CGLS()
68 PetscCall(VecAXPY(r, -alpha, q)); /* r -= alpha * q */ in KSPSolve_CGLS()
/petsc/src/ksp/ksp/utils/lmvm/blas_cyclic/
H A Dblas_cyclic.c16 static inline void AXPBY_Private(PetscInt m, PetscScalar alpha, const PetscScalar x[], PetscScalar … 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 …XPBYCylic_Private(PetscInt m, PetscInt oldest, PetscInt next, PetscScalar alpha, const PetscScalar… 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 PETSC_INTERN PetscErrorCode VecAXPBYCyclic(PetscInt oldest, PetscInt next, PetscScalar alpha, Vec x… 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/impls/cg/stcg/
H A Dstcg.c19 PetscReal alpha, beta, kappa, rz, rzm1; in KSPCGSolve_STCG()
88 alpha = 1.0; in KSPCGSolve_STCG()
91 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_STCG()
95 PetscCall(VecAXPY(d, alpha, r)); /* d = d + alpha r */ in KSPCGSolve_STCG()
123 alpha = 1.0; in KSPCGSolve_STCG()
126 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_STCG()
130 PetscCall(VecAXPY(d, alpha, r)); /* d = d + alpha r */ in KSPCGSolve_STCG()
201 alpha = 1.0; in KSPCGSolve_STCG()
204 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_STCG()
208 PetscCall(VecAXPY(d, alpha, r)); /* d = d + alpha r */ in KSPCGSolve_STCG()
[all …]
/petsc/src/ksp/ksp/impls/cg/nash/
H A Dnash.c20 PetscReal alpha, beta, kappa, rz, rzm1; in KSPCGSolve_NASH()
91 alpha = 1.0; in KSPCGSolve_NASH()
94 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_NASH()
98 PetscCall(VecAXPY(d, alpha, r)); /* d = d + alpha r */ in KSPCGSolve_NASH()
126 alpha = 1.0; in KSPCGSolve_NASH()
129 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_NASH()
133 PetscCall(VecAXPY(d, alpha, r)); /* d = d + alpha r */ in KSPCGSolve_NASH()
204 alpha = 1.0; in KSPCGSolve_NASH()
207 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_NASH()
211 PetscCall(VecAXPY(d, alpha, r)); /* d = d + alpha r */ in KSPCGSolve_NASH()
[all …]
/petsc/src/mat/tests/
H A Dex2.c5 static PetscErrorCode TransposeAXPY(Mat C, PetscScalar alpha, Mat mat, PetscErrorCode (*f)(Mat, Mat… in TransposeAXPY() argument
20 PetscCall(MatAXPY(E, alpha, mat, SAME_NONZERO_PATTERN)); in TransposeAXPY()
38 PetscCall(MatAXPY(C, alpha, E, SAME_NONZERO_PATTERN)); in TransposeAXPY()
53 PetscCall(MatAXPY(E, alpha, G, SAME_NONZERO_PATTERN)); in TransposeAXPY()
69 PetscCall(MatAXPY(F, alpha, mat, SAME_NONZERO_PATTERN)); in TransposeAXPY()
84 PetscScalar v, alpha; in main() local
147 alpha = 1.0; in main()
148 PetscCall(PetscOptionsGetScalar(NULL, NULL, "-alpha", &alpha, NULL)); in main()
150 PetscCall(MatAXPY(tmat, alpha, mat, DIFFERENT_NONZERO_PATTERN)); in main()
154 PetscCall(MatAYPX(tmat, alpha, mat, DIFFERENT_NONZERO_PATTERN)); in main()
[all …]
/petsc/src/mat/impls/aij/seq/aijmkl/
H A Daijmkl.h33 …#define mkl_xcsrmv(transa, m, k, alpha, matdescra, val, indx, pntrb, pntre, x, beta, y) mkl_scsrmv… argument
35 …#define mkl_xcsrmv(transa, m, k, alpha, matdescra, val, indx, pntrb, pntre, x, beta, y) mkl_dcsrmv… argument
39 …#define mkl_xcsrmv(transa, m, k, alpha, matdescra, val, indx, pntrb, pntre, x, beta, y) mkl_ccsrmv… argument
41 …#define mkl_xcsrmv(transa, m, k, alpha, matdescra, val, indx, pntrb, pntre, x, beta, y) mkl_zcsrmv… argument
/petsc/src/ksp/ksp/impls/lsqr/
H A Dlsqr.c72 PetscReal beta, alpha, rnorm; in KSPSolve_LSQR() local
125 PetscCall(VecNorm(V, NORM_2, &alpha)); in KSPSolve_LSQR()
130 PetscCall(VecDotRealPart(V, Z, &alpha)); in KSPSolve_LSQR()
131 if (alpha <= 0.0) { in KSPSolve_LSQR()
133 PetscCall(PetscInfo(ksp, "Diverging due to breakdown alpha (%g) <= 0\n", (double)alpha)); in KSPSolve_LSQR()
134 …m((PetscObject)ksp), PETSC_ERR_NOT_CONVERGED, "KSPSolve breakdown alpha (%g) <= 0", (double)alpha); in KSPSolve_LSQR()
137 alpha = PetscSqrtReal(alpha); in KSPSolve_LSQR()
138 PetscCall(VecScale(Z, 1.0 / alpha)); in KSPSolve_LSQR()
140 PetscCall(VecScale(V, 1.0 / alpha)); in KSPSolve_LSQR()
151 lsqr->arnorm = alpha * beta; in KSPSolve_LSQR()
[all …]

12345678910>>...18