You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
julia> adjoint(opEye()) * randn(10)
ERROR: type opEye has no field nrow
Stacktrace:
[1] getproperty(x::opEye, f::Symbol)
@ Base ./Base.jl:37
[2] size(op::opEye)
@ LinearOperators ~/.julia/packages/LinearOperators/WS4fM/src/abstract.jl:194
[3] size(A::AdjointLinearOperator{Any, opEye})
@ LinearOperators ~/.julia/packages/LinearOperators/WS4fM/src/adjtrans.jl:62
[4] *(op::AdjointLinearOperator{Any, opEye}, v::Vector{Float64})
@ LinearOperators ~/.julia/packages/LinearOperators/WS4fM/src/operations.jl:41
[5] top-level scope
@ REPL[3]:1
I think the problem is that *(A::opEye, x) has a special implementation, but for the adjoint, it falls back to the generic implementation of *(A::AdjointLinearOperator, x) which assumes a fully-described LinearOperator.
The bug should be fixed (or rather avoided) when adding the following code
adjoint(A::opEye) = A
transpose(A::opEye) = A
conj(A::opEye) = A
which will ensure that the *(A::opEye, x) implementations are used. But I am sure there are other solutions.
Happy to create a PR -- LMK. Thanks.
-ja
The text was updated successfully, but these errors were encountered:
Hi everyone,
I observed the following bug:
I think the problem is that
*(A::opEye, x)
has a special implementation, but for the adjoint, it falls back to the generic implementation of*(A::AdjointLinearOperator, x)
which assumes a fully-described LinearOperator.The bug should be fixed (or rather avoided) when adding the following code
which will ensure that the
*(A::opEye, x)
implementations are used. But I am sure there are other solutions.Happy to create a PR -- LMK. Thanks.
-ja
The text was updated successfully, but these errors were encountered: