General Question about Linear Solvers

Matthew Knepley knepley at
Tue Mar 6 15:44:42 CST 2007

In general, iterative information cannot be reused, or it is of little use.
Until you get truly large, I would use sparse direct packages like
MUMPS or SuperLU_dist through PETSc. They are simple (hopefully)
with configure.


On 3/6/07, Manav Bhatia <manav at> wrote:
> Hi
>     I have an application where a linear system of equations has to be
> solved for many right hand sides. I am curious to know about how the
> different solvers compare with each other for such a case?
> `   So far, I have only used the LU decompositionm, where, once
> decomposed,
> the factorization can be use each time. I have been solving small systems
> so far, hence, cost of factorization has not been my concern. However, I
> will be moving towards much larger problems, requiring parallel
> processing.
>     I would greatly appreciate if someone could comment on this.
> Thanks,
> Manav Bhatia
> PhD Candidate
> Department of Aeronautics and Astronautics
> Guggenheim Building
> University of Washington
> Seattle, WA 98195-2400
> Ph (O) (206) 685 8063
> e-mail: manav at

One trouble is that despite this system, anyone who reads journals widely
and critically is forced to realize that there are scarcely any bars to eventual
publication. There seems to be no study too fragmented, no hypothesis too
trivial, no literature citation too biased or too egotistical, no design too
warped, no methodology too bungled, no presentation of results too
inaccurate, too obscure, and too contradictory, no analysis too self-serving,
no argument too circular, no conclusions too trifling or too unjustified, and
no grammar and syntax too offensive for a paper to end up in print. --
Drummond Rennie

More information about the petsc-users mailing list