matrices


Ginodesilva
05-26-2005, 10:07 AM
Hi everybody,

ik have a problem and am looking for a nice solution (sequential or parallel doesn't matter)

I'm trying to write a recursive subroutine in which a lower-triangle matrix (diagonal is 1, above is 0, below doesn't matter) is divided into 4 equal-sized matrices untill the submatrix has a size of 2x2. In this case I can calculate the inverse of the 4x4 matrix easily by products of the 2x2 matrices.
I'm trying to achieve the following:
http://www.famhornes.nl/example.jpg

now, I'm puzzeled how to divide the matrix into 4 submatrices and back into 1 easy.

Does anybody have a clue?

Any help is appreciated.

Grtz, Gino