I am trying to solve systems of linear equations Au=b where the coefficient matrix A possesses a block tridiagonal structure with each of the block matrices is block tridiagonal as well. Anyone who can help me to solve this system of linear equations. is there any available algorithm for the inversion of block-block tridiagonal matrix?

Don't know if there's a special algorithm for the inner set of blocks, but the 'standard' block triagonal algorithm should work.

Are you sure you need the inverse? Factoring the matrix into upper and lower triangular parts is much more efficient if you really just need the vector "u".

Beyond this answer, some help from a linear algebra expert is in order.

Good luck!

Hi,

The discretization of a Poisson Equation results in these Block-Block Tridiagonal Matrices. There are lots of open source packages which you can use to solve these type of equations. LAPACK, SuperLU, PETSc.. Actually you will not need to construct the inverse to solve these type of problems.. Direct methods such as LU and iterative methods such as GMRES/Conjugate Gradient should solve your problem.

-Dominic

what about TDMA....???