Unary operations on adjoint fail about linearoperators.jl HOT 6 CLOSED

Thanks for the report @andreasvarga.

@dpo, define -A' as (-A)'?

from linearoperators.jl.

Yes

from linearoperators.jl.

Just a comment: Implementing -A' as (-A)' is OK. However, if you intend to implement 2A' as (2A)' , this may raise the question which operation should have priority: the multiplication with 2 or the adjoint. In the LinearAlgebra package of Julia, the multiplication defines the result, and therefore 2*Adjoint(a) does not lead to and Adjoint type, see below:

julia> 2*rand(2,2)'
2×2 Array{Float64,2}:
 1.25212  0.0511282
 1.77288  1.844

This approach has sense, because it is always possible to enforce the Adjoint type by computing (2*A)'.

from linearoperators.jl.

Do you mean that you're interested in cases where multiplication and transpose don't commute?

from linearoperators.jl.

from linearoperators.jl.

Because adjoint and transpose operators don't have explicit [ct]prod methods, I'm not sure how else it can be implemented.

from linearoperators.jl.