Dominique
Dominique
> LDLFactorizations is not used in LinearOperators.jl That's easy to change.
@geoffroyleconte Should be all fixed. Do you mind taking another quick look?
It's expected that FreeBSD/nightly fail at the moment.
Thank you @Joshua-Wolff and @geoffroyleconte !
Is every update accepted? https://github.com/JuliaSmoothOptimizers/LinearOperators.jl/blob/main/src/lsr1.jl#L150
Thank you for this analysis. > I made several tests on my SR1 implementation using collinear pairs s,y, and the resulting matrices satisfied the secant equation and were well-defined. How...
In fact, I don't remember where `scaling_condition` came from. Maybe it was an experiment and should be removed. Clearly, it prevents updates when s and y are collinear.
If I understand your question correctly, symmetry of the solution of your Lyapunov equation doesn't really depend on the operator, but rather on how you model the equation. Couldn't you...
So it's the default implementation of `Matrix(op::AbstractLinearOperator)` that is causing you trouble. An easy solution might be to define your own operator type, say `LyapunovOperator
> I think I found a less elegant way to manage the issue, by assuming that only the upper triangular part of X is packed in the vector x Yes,...