I have a matrix A of size n^2 by n^2, and I wanted to know if for a given accuracy (or a number r) there is a way to express A as the sum of Bi kron Ci for i=1...R
where Bi, Ci are n by n? i.e.
There are some toolboxes for cases where A is a vector, but I do not want to reshape A as a vector since it grows badly with n. So I want to deal with A directly.
How can I decompose this matrix?

v = M(:).