JuliaSmoothOptimizers / JuliaSmoothOptimizers/LDLFactorizations.jl

possible method error upper triangle

Open
#76 1 comment 0 reactions 0 assignees View on GitHub
documentation
Dominant language
Julia
Stars
36
Forks
13
PR merge metrics
No merged PRs in 30d

Description

```
n, m = size(A)
A = sprand(Float64, 10, 5, 0.25)
λ = 1.0e-3

K1 = triu([sparse(1.0I,n,n) A; A' -λ^2 * I])
S1 = ldl(Symmetric(K1, :U))

K2 = [sparse(1.0I,n,n) A; A' -λ^2 * I]
S2 = ldl(Symmetric(K2, :U))

b = ones(m+n)

norm(S1 \ b - S2 \ b) # not zero
```

Hi! I think this may be a method error. Do we still need L521 - L577 ?

Contributor guide

Open the contributing guide

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.