Andreas Varga

I solved the problem by simply making a copy of the input data. Regarding PR: I never prepared a PR and I am afraid I am not enough skilled to...

Let consider the inverse operator inv(L) of a Lyapunov operator L(X) =AX+XA'. This operator maps a symmetric matrix Y to a corresponding symmetric X which satisfies the Lyapunov equation AX+XA'+Y...

Yes, this could be an elegant solution to my problem. In the meantime, I think I found a less elegant way to manage the issue, by assuming that only the...

No matrices are involved, only vectors which are the results of product operations op*x and op'*x. So, it seems to be perfectly suited for LinearOperators. Dominique schrieb am Mi., 16....

Just an update: I formulated for my package an [issue](https://github.com/andreasvarga/MatrixEquations.jl/issues/2) which is related to the above discussion. Simply said, I was not able to ensure the condition `Matrix(T)' = Matrix(T')`...

Many thanks for your explanation. Abel Siqueira schrieb am Mo., 5. Okt. 2020, 17:00: > The issue here is that Rational does not provide enough information. It > should be...

In MatrixEquations.jl the operation exits with error for a Lyapunov operator of type Float32 and vector of type Float64. Abel Siqueira schrieb am Mo., 5. Okt. 2020, 17:07: > Regarding...

I updated my package [MatrixEquations.jl](https://github.com/andreasvarga/MatrixEquations.jl) to v2.0 and all tests involving LinearOperators failed! It is probably a trivial issue, but in this moment I would appreciate any hint on what...

Thanks. I wonder if this will be a simple exercise for me. Here is an example for the Lyapunov operator: ``` """ L = lyapop(A; disc = false, her =...

I am working with Julia 1.5. eltype(eltype(a)) does not work! Is this a new feature in Julia 1.6? john verzani schrieb am Sa., 6. Feb. 2021, 19:05: > I think...