Does anyone know an algorithm for computing kronecker products of two matrices? It's probably not that hard, but I feel like my head is about to explode ATM, so if you can help me out that'd be cool. I want to implement this in fortran... I'll give you an example; Say I want compute the kronecker product between two 2x2 matrices (http://en.wikipedia.org/wiki/Tensor_product#Kronecker_product_of_two_matrices). I know that in the end the result of the product is a tensor, so the straightforward way to do this would be to create a rank 4 tensor in the program, but that's unwieldy... What I want is to generate a matrix (4x4), starting from these 2x2 matrices.