> The following code uses the for-loop method I mentioned earlier ... it will save you memory space, but at the cost of slower execution than with the use of the 'kron' function.
Yes, using multiple loops was actually my starting point, but then I removed the loops in an attempt to speed up execution. This worked pretty well for modestly-sized matrices (allowing execution to finish within a day or two). However, I've hit a wall with this slightly larger data size in that I now run out of memory.
I suppose a brute-force approach would be to get my hands on a 64-bit architecture. Or I could try keeping "some" loops and not others... but it's difficult to identify a good trade-off.