Abstract type for statistics returned by a solver

Type for statistics returned by the majority of Krylov solvers, the attributes are:

• niter
• solved
• inconsistent
• residuals
• Aresiduals
• Acond
• status

Type for statistics returned by CG-LANCZOS, the attributes are:

• niter
• solved
• residuals
• indefinite
• Anorm
• Acond
• status

Type for statistics returned by CG-LANCZOS with shifts, the attributes are:

• niter
• solved
• residuals
• indefinite
• Anorm
• Acond
• status

Type for statistics returned by SYMMLQ, the attributes are:

• niter
• solved
• residuals
• residualscg
• errors
• errorscg
• Anorm
• Acond
• status

Type for statistics returned by adjoint systems solvers BiLQR and TriLQR, the attributes are:

• niter
• solved_primal
• solved_dual
• residuals_primal
• residuals_dual
• status

Type for statistics returned by the LNLQ method, the attributes are:

• niter
• solved
• residuals
• errorwithbnd
• errorbndx
• errorbndy
• status

Type for statistics returned by the LSLQ method, the attributes are:

• niter
• solved
• inconsistent
• residuals
• Aresiduals
• err_lbnds
• errorwithbnd
• errubndslq
• errubndscg
• status

Type for statistics returned by LSMR. The attributes are:

• niter
• solved
• inconsistent
• residuals
• Aresiduals
• Acond
• Anorm
• xNorm
• status

Abstract type for using Krylov solvers in-place

Type for storing the vectors required by the in-place version of MINRES.

The outer constructors

solver = MinresSolver(n, m, S; window :: Int=5)
solver = MinresSolver(A, b; window :: Int=5)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of CG.

The outer constructors

solver = CgSolver(n, m, S)
solver = CgSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of CR.

The outer constructors

solver = CrSolver(n, m, S)
solver = CrSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of SYMMLQ.

The outer constructors

solver = SymmlqSolver(n, m, S)
solver = SymmlqSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of CG-LANCZOS.

The outer constructors

solver = CgLanczosSolver(n, m, S)
solver = CgLanczosSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of CG-LANCZOS-SHIFT.

The outer constructors

solver = CgLanczosShiftSolver(n, m, nshifts, S)
solver = CgLanczosShiftSolver(A, b, nshifts)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of MINRES-QLP.

The outer constructors

solver = MinresQlpSolver(n, m, S)
solver = MinresQlpSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of DIOM.

The outer constructors

solver = DiomSolver(n, m, memory, S)
solver = DiomSolver(A, b, memory = 20)

may be used in order to create these vectors. memory is set to n if the value given is larger than n.

Type for storing the vectors required by the in-place version of FOM.

The outer constructors

solver = FomSolver(n, m, memory, S)
solver = FomSolver(A, b, memory = 20)

may be used in order to create these vectors. memory is set to n if the value given is larger than n.

Type for storing the vectors required by the in-place version of DQGMRES.

The outer constructors

solver = DqgmresSolver(n, m, memory, S)
solver = DqgmresSolver(A, b, memory = 20)

may be used in order to create these vectors. memory is set to n if the value given is larger than n.

Type for storing the vectors required by the in-place version of GMRES.

The outer constructors

solver = GmresSolver(n, m, memory, S)
solver = GmresSolver(A, b, memory = 20)

may be used in order to create these vectors. memory is set to n if the value given is larger than n.

Type for storing the vectors required by the in-place version of USYMLQ.

The outer constructors

solver = UsymlqSolver(n, m, S)
solver = UsymlqSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of USYMQR.

The outer constructors

solver = UsymqrSolver(n, m, S)
solver = UsymqrSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of TRICG.

The outer constructors

solver = TricgSolver(n, m, S)
solver = TricgSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of TRIMR.

The outer constructors

solver = TrimrSolver(n, m, S)
solver = TrimrSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of TRILQR.

The outer constructors

solver = TrilqrSolver(n, m, S)
solver = TrilqrSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of CGS.

The outer constructorss

solver = CgsSolver(n, m, S)
solver = CgsSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of BICGSTAB.

The outer constructors

solver = BicgstabSolver(n, m, S)
solver = BicgstabSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of BILQ.

The outer constructors

solver = BilqSolver(n, m, S)
solver = BilqSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of QMR.

The outer constructors

solver = QmrSolver(n, m, S)
solver = QmrSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of BILQR.

The outer constructors

solver = BilqrSolver(n, m, S)
solver = BilqrSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of CGLS.

The outer constructors

solver = CglsSolver(n, m, S)
solver = CglsSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of CRLS.

The outer constructors

solver = CrlsSolver(n, m, S)
solver = CrlsSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of CGNE.

The outer constructors

solver = CgneSolver(n, m, S)
solver = CgneSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of CRMR.

The outer constructors

solver = CrmrSolver(n, m, S)
solver = CrmrSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of LSLQ.

The outer constructors

solver = LslqSolver(n, m, S)
solver = LslqSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of LSQR.

The outer constructors

solver = LsqrSolver(n, m, S)
solver = LsqrSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of LSMR.

The outer constructors

solver = LsmrSolver(n, m, S)
solver = LsmrSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of LNLQ.

The outer constructors

solver = LnlqSolver(n, m, S)
solver = LnlqSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of CRAIG.

The outer constructors

solver = CraigSolver(n, m, S)
solver = CraigSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of CRAIGMR.

The outer constructors

solver = CraigmrSolver(n, m, S)
solver = CraigmrSolver(A, b)

may be used in order to create these vectors.

Type for storing the vectors required by the in-place version of GPMR.

The outer constructors

solver = GpmrSolver(n, m, memory, S)
solver = GpmrSolver(A, b, memory = 20)

may be used in order to create these vectors. memory is set to n + m if the value given is larger than n + m.

roots = roots_quadratic(q₂, q₁, q₀; nitref)

Find the real roots of the quadratic

q(x) = q₂ x² + q₁ x + q₀,

where q₂, q₁ and q₀ are real. Care is taken to avoid numerical cancellation. Optionally, nitref steps of iterative refinement may be performed to improve accuracy. By default, nitref=1.

(c, s, ρ) = sym_givens(a, b)

Numerically stable symmetric Givens reflection. Given a and b reals, return (c, s, ρ) such that

[ c  s ] [ a ] = [ ρ ]
[ s -c ] [ b ] = [ 0 ].

Numerically stable symmetric Givens reflection. Given a and b complexes, return (c, s, ρ) with c real and (s, ρ) complexes such that

[ c   s ] [ a ] = [ ρ ]
[ s̅  -c ] [ b ] = [ 0 ].

roots = to_boundary(x, d, radius; flip, xNorm2, dNorm2)

Given a trust-region radius radius, a vector x lying inside the trust-region and a direction d, return σ1 and σ2 such that

‖x + σi d‖ = radius, i = 1, 2

in the Euclidean norm. If known, ‖x‖² may be supplied in xNorm2.

If flip is set to true, σ1 and σ2 are computed such that

‖x - σi d‖ = radius, i = 1, 2.

s = vec2str(x; ndisp)

Display an array in the form

[ -3.0e-01 -5.1e-01  1.9e-01 ... -2.3e-01 -4.4e-01  2.4e-01 ]

with (ndisp - 1)/2 elements on each side.

S = ktypeof(v)

Return a dense storage type S based on the type of v.

v = kzeros(S, n)

Create an AbstractVector of storage type S of length n only composed of zero.

v = kones(S, n)

Create an AbstractVector of storage type S of length n only composed of one.