What are the properties of matrices being passed to kron_multiply?
i.e., what is the size and sparsity patterns of these matrices for the problem sizes of interest?
i.e., what is the size and sparsity patterns of these matrices for the problem sizes of interest?