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
Dealing with non-positive definite matrices in Cholesky decomposition: R vs. Julia General Usage linearalgebra	4	1336	March 14, 2022
`isposdef` inconsistent results, numerically unstable? Numerics bug	54	4378	February 21, 2018
Is this a bug with Cholesky? General Usage linearalgebra	6	4757	June 1, 2019
A * B * A' is not hermitian, even when B is General Usage linearalgebra	8	1047	November 1, 2021
Isposdef() true with small negative eigenvalues? General Usage linearalgebra	3	212	March 12, 2024

Struggling with isposdef and floating point errors

Related topics