space allocation for a parallel matrix

Satish Balay balay at mcs.anl.gov
Tue Feb 3 11:25:04 CST 2009


One issue that migh be cause for confusion is: for sparse mat-vec -
only a part of the vec from other procs is needed.

This part is communicated during Mat-Vec [using VecScatters]

So the 'whole 'vec' is not stored on each proc..

Satish

On Tue, 3 Feb 2009, Matthew Knepley wrote:

> On Tue, Feb 3, 2009 at 9:56 AM, Yixun Liu <enjoywm at cs.wm.edu> wrote:
> > It means that it is better to set n to be equal to the column to avoid
> > the communication. Right? If it does, I am confusing why we need this
> > parameter. Just for dealing with very large vector x?
> > Thanks.!
> 
> No. Please take a look at the book I recommended. If you do as you suggest,
> this is no longer parallel computing.
> 
>   Matt
> 
> 



More information about the petsc-users mailing list