Does there exist functionality analogous to those listed in https://julialinearalgebra.github.io/LinearMaps.jl/stable/types/#BlockMap-and-BlockDiagonalMap?
In particular, lazy hcat/vcat/hvcat/diagonal op formation, etc.
We should also revisit this line of code:
|
LinearAlgebra.Diagonal(L::MatrixOperator) = MatrixOperator(Diagonal(L.A)) |