I have implemented a basis transformation in C of the following form
B are matrices and
kron is the Kronecker product.
However, my naive code is not scalable, and I am looking for a library to achieve this. My starting point was to look in BLAS and LAPACK for Kronecker products and basis contractions (ie.
B*A*Transpose[B]) but from my reading of the documentation these are not in the standard build, is there an efficient way to achieve this using these packages (in particular the contraction)?