[petsc-dev] Vec distribution

[Matthew Knepley]

On Sat, Sep 14, 2013 at 5:55 PM, Jose David Bermeol <jbermeol at purdue.edu>wrote:

> Hi, right now I'm implementing the function MatSolve for MATMPIAIJ
> matrices, so I can assume that output vector x is an MPI vector. the
> problem is that my library returns the solution in an array to all my MPI
> process, so how do I know which elements of the solution should be safe in
> the local vector x.
>

I do not understand what you are asking. However, each Vec is divided
across processes by contiguous sets of rows. The rows
for a given process can be gotten from VecGetOwnershipRange().

   Matt


> Thanks
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130914/bda06f6b/attachment.html>