problems with MatLoad

[Andreas Grassl]

Matthew Knepley schrieb:
> On Wed, Apr 8, 2009 at 11:12 AM, Andreas Grassl
> <Andreas.Grassl at student.uibk.ac.at
> <mailto:Andreas.Grassl at student.uibk.ac.at>> wrote:
> 
>     Hello,
> 
>     I got some success on the localtoglobalmapping, but now I'm stuck
>     with writing
>     to/reading from files. In a sequential code I write out some
>     matrices with
> 
>     PetscViewerBinaryOpen(comms,matrixname,FILE_MODE_WRITE,&viewer);
>     for (k=0;k<np;k++){
>     MatView(AS[k],viewer);}
>     PetscViewerDestroy(viewer);
> 
>     and want to read them in in a parallel program, where each processor
>     should own
>     one matrix:
> 
>     ierr =
>     PetscViewerBinaryOpen(PETSC_COMM_WORLD,matrixname,FILE_MODE_READ,&viewer);CHKERRQ(ierr);
> 
> 
> The Viewer has COMM_WORLD, but you are reading a matrix with COMM_SELF.
> You need to create
> a separate viewer for each process to do what you want.
> 

Thank you for the fast answer. I resolved this issue now, but how could i gather
the Matrix from COMM_SELF to COMM_WORLD. I searched for functions doing such
matrix copying, but MatConvert and MatCopy act on the same communicator.

Thanks in advance

ando

-- 
 /"\                               Grassl Andreas
 \ /    ASCII Ribbon Campaign      Uni Innsbruck Institut f. Mathematik
  X      against HTML email        Technikerstr. 13 Zi 709
 / \                               +43 (0)512 507 6091