I am trying to do a Quasi-Newton method with a basic cholesky adaptation, to guarantee B is positive definite. The code goes as follows: function choleksyadaption!(A,β = 10e-3, max_iter = 1000) @info "Adapting Cholesky Factorization to make Hessian positive definite." if (isposdef(A)) …

You can use Hermitian to tell the program that a matrix is supposed to be Hermitian. Then it will only look at half the matrix (by default, the upper triangle) and assume that the opposite triangle is conjugate-equal. julia> using LinearAlgebra julia> Hermitian([1 2; 3 4], :U) # :U is the default …

Struggling with isposdef and floating point errors

General Usage

stevengj March 13, 2024, 7:25pm 3

Topic		Replies	Views
Cholesky factorization for slightly non-Hermitian matrices Numerics linearalgebra	18	4126	January 28, 2019
`isposdef` inconsistent results, numerically unstable? Numerics bug	54	4554	February 21, 2018
Has anyone noticed numerical instability in Base.LinAlg.cholfact or Base.LinAlg.ishermitian? Statistics question , linearalgebra	3	899	April 4, 2018
How to force a matrix to be Hermitian? General Usage linearalgebra	7	4625	December 10, 2018
Is this a bug with Cholesky? General Usage linearalgebra	6	4784	June 1, 2019

Struggling with isposdef and floating point errors

Related topics