Concatenation of Kronecker Products

Open MKAbdElrahman opened this issue 4 years ago • 3 comments

I'm taking Kronecker products of sparse matrices. when I concatenate these matrices, full dense matrices are returned which is very unefficient.

Mar 08 '21 18:03 MKAbdElrahman

As in #76, I recommend doing "higher-order" manipulations like linear combinations and concatenation of Kronecker products with LinearMaps.jl. You are not restricted to use its own Kronecker product, but can safely use the types of Kronecker.jl.

Mar 08 '21 18:03 dkarrasch

This works, but doesn't mesh well with the iterativesolvers package. The performance is so slow even for small matrices. I used the kron function from LinearMaps.jl and performance is good.

Mar 08 '21 21:03 MKAbdElrahman

Strange, can you give a small working example?

Mar 09 '21 08:03 MichielStock