Has anyone developed some procedures for handling block matrix operations? I.e., blocks may be names of matrix stuctures (assigned or not) and operations on block matrices result in matrix multiplication and addition of the block rows and columns? This would require manipulating matrices of matrices (arrays of arrays) in the correct fashion.
The linalg[blockmatrix] function just helps build up regular matrices from smaller blocks but doesn’t preserve the block nature of the bigger matrix.
As a result of a request received by MUG I wrote some routines last summer. They can be found on my web page:
For those without net access I can email them. Comments and suggestions would be appreciated.