tridiagonal matrix
Peter Jay Salzman
p at belial.ucdavis.edu
Fri Jul 28 12:47:21 PDT 2000
dear all,
i'm just beginning to learn about MPI programming.
i have a series of linear equations which are represented by a tridiagonal
matrix. the order of the matrix is between 1000 and 10000. this is thfe
major bottleneck in my code, so i'd like to know if solving a tridiagonal
matrix is a parallelizable operation.
the only tridiagonal matrix algorithm i know of is the thomas algorithm.
the algorithm has a structure like:
get a result for a[i]
use a[i] to obtain a[i+1]
use a[i+1] to obtain a[i+2]
...
this looks decidedly non-parallelizable to me. i was wondering if anyone
knew of another algorithm which i could implement using MPI to get past this
bottleneck.
(drop in code for fortran or C++ would be fantabulous!).
sorry for typos, but i just moved, and dsl hasn't kicked in yet. i'm on a
painfully slow ppp connection.
much thanks!
pete
More information about the Beowulf
mailing list