LinearOperators.jl
LinearOperators.jl copied to clipboard
JuliaSmoothOptimizers
Linear Operators for Julia
Hello all, I did test the LSR1 operators for myself, and it seems that there is a problem with the LSR1Operator of dimension 1. In my work, I generate automatically...
Proposition to fix #225 fix LSR1Operator when s and y are collinear by cancelling data.scaling
Hello! Context: - I have to update several `LBFGSOperator` with differents `(s,y)` for a partitioned update [PLBFGS](https://github.com/paraynaud/PartitionedStructures.jl/blob/main/src/algorithms/PartitionedLOQuasiNewton.jl). - I cannot garantee that each `(s,y)` verify the cruvature condition. I wonder...
What is the best way to broadcast `LinearOperators`? For instance, ``` Jx = jac_op(nlp, x) # where nlp is an AbstractNLPModel, and x a vector ``` and then ``` Jx...
I would be nice to add in-place solves to LinearOperators in the linalg.jl file.
`*` between operators and Numbers is assumed to be a commutative operation, but this could not be the expected behaviour for some specific Numbers (quaternions? ).
The following code illustrates this issue: ``` M = ones(Rational,3,3) v = ones(Int,3) M*v # OK, the result is a rational vector LinearOperator(M)*v # failure ``` The execution of this...
I implemented several linear operators related to solve control problems. An example is the Lyapunov operator L: X -> AX+XA', where A is a squre matrix. This operator acts usually...