Yes, you can have 0 length IS arguments. Matt
On Wed, Dec 31, 2008 at 10:21 AM, Yujie <recrusader at gmail.com> wrote: > Hi, PETSc Developers > > In parallel mode, I have 3 MPIDense matrices, A1, A2, A3. Now, I need to > exact submatrices B1, B2, B3 in parallel mode from A1, A2, A3 respectively. > I know I should use MatGetSubMatrix(). My problem is in > > MatGetSubMatrix(Mat mat,IS isrow,IS iscol,PetscInt csize,MatReuse cll,Mat > *newmat) > > , the parameter "isrow" means "rows this processor should obtain". if Bi in > some nodes (cpus) of the cluster (Ai have rows in them) don't have rows, is > it work? > > After getting Bi, I will create MPIDense-based C and combine them into C > like > > B1 > C= B2 > B3 > MatGetRow(); > MatSetValues() should work, right? > > Could you give me some comments about these operations? thanks a lot. Happy > new year! > > Regards, > Yujie > -- 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-users/attachments/20081231/53fe1174/attachment.htm>
