space allocation for a parallel matrix

Yixun Liu enjoywm at cs.wm.edu
Tue Feb 3 09:25:39 CST 2009


Hi,
I am confusing to the allocation for a parallel matrix in PETSC manual,

"If PETSC DECIDE is not used for the arguments m and n, then the user
must ensure that they are chosen
to be compatible with the vectors. To do this, one first considers the
matrix-vector product y = Ax. The m
that is used in the matrix creation routine MatCreateMPIAIJ() must match
the local size used in the vector
creation routine VecCreateMPI() for y. Likewise, the n used must match
that used as the local size in
VecCreateMPI() for x."

Since the process holds all the columns, why does it only hold part of
vector x? Is it convenient for the matrix vector multiplication?

Thanks.

Yixun


More information about the petsc-users mailing list