Andreas Varga

`u` is not modified, but certainly `a`: ``` julia> using BenchmarkTools julia> using LinearAlgebra julia> a = rand(300,300); u = triu(rand(300,300)); julia> A = copy(a); U = copy(u); julia> @benchmark...

Sorry if I am wrong, but here is what I am expecting. For a matrix based linear map I have: ``` julia> B = LinearMap(rand(2,2)) 2×2 LinearMaps.WrappedMap{Float64} of 2×2 Matrix{Float64}...

Actually I already used [https://diffeq.sciml.ai/stable/features/diffeq_operator/ ](https://diffeq.sciml.ai/stable/features/diffeq_operator/ ) to implement the function for a particular integration method (Magnus' mthod). I think it will be possible to implement function matrices, as for...

My understanding is that` \lstset` is a convenient way to modify some predefined language settings for generating a code listing in that language. So, this feature would be desirable to...

Thanks for your time and I fully agree with your points. However, for my purposes the output of the above code is not satisfactory and therefore I need more flexibility...

And even ``` julia> Ty = Double64; c = rand(Ty,1) ; d =rand(Ty,1); c[1]*d[1]-d[1]*c[1] 1.5407439555097887e-33 ``` Thus, the multiplication is not commutative?

Could you provide a link to `orthogonal_transformation_to`? What mathematical condition guarantees the existence of an orthogonal solution?

The fact that A and B are commuting (i.e., A*B = B*A) does guarantee the existence of an orthogonal solution? I was not able to find such a result in...

Could you provide an example for using this function?

It is a good idea to add some solvers for large Sylvester and Lyapunov equations. The provided implementation has in my opinion some issues, as for example, (1) the lack...