On Tue, 1 Jun 2010 00:11:41 +0800, "Gong Ding" <gdiso at ustc.edu> wrote: > Dear all, > I need to solve Poisson's equation many times with same Matrix A and > different RHS b. > The most efficient way is do pre-factorization I think. > > For serial problem, the PCLU seems work as I want to. > However, how to do for parallel situation? > Should I use an external LU solver such as superlu_dist?
Sure, or MUMPS or PaStiX. As the problem size gets bigger, eventually a multigrid will beat even the triangular solves. Jed
