I get an stackoverflow error when I try to apply an adjoint of a blockbanded operator on a function. For instance
Op=Derivative(Chebyshev())⊗Derivative(Chebyshev())
f=Fun(sin)⊗Fun(sin)
Op*f # works fine
Op'*f # Stackoverflow
I believe this is due to blockbandwidths(P::AdjointOperator) not being implemented for adjoint operators (while bandwidths(P::AdjointOperator) is implemented in OperatorLayout.jl) and thus the operator is not considered blockbanded. Indeed
I get an stackoverflow error when I try to apply an adjoint of a blockbanded operator on a function. For instance
I believe this is due to
blockbandwidths(P::AdjointOperator)
not being implemented for adjoint operators (whilebandwidths(P::AdjointOperator)
is implemented in OperatorLayout.jl) and thus the operator is not considered blockbanded. IndeedThe following implementation did the job for me:
and
Op'*f
works.