[Scalapack] A new storage scheme to solve general five diagonal matrix
From: ehsan jamshidy
Date: Wed, 9 Feb 2011 06:50:17 +0330
I have developed a general five diagonal n*n matrix solver by
allocating six, one dimensional arrays (nx1).
This matrix is used excessively for solving pressure and velocity in CFD
(Computational Fluid Dynamics). A five diagonal matrix which is used in CFD
is a bounded matrix with a major diagonal plus two upper
diagonal, that there are some optional zeroes diagonals between them, and
the same for lower diagonals.
In traditional solvers like (DMCHB or DGELB subroutines) we hove to store
the zeroes between two upper and lower diagonals, but in the present method
we don't need to allocate these zeroes.
Have you seen yet a program by this features?