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	4226	January 28, 2019
`isposdef` inconsistent results, numerically unstable? Numerics bug	54	4816	February 21, 2018
Isposdef and eigvals do not agree Statistics linearalgebra , distributions , potential-bug	9	246	August 16, 2024
A * B * A' is not hermitian, even when B is General Usage linearalgebra	8	1120	November 1, 2021
How to force a matrix to be Hermitian? General Usage linearalgebra	7	4701	December 10, 2018

Struggling with isposdef and floating point errors

Related topics