[petsc-dev] [Fortran] MatLoad() error: MatLoad is not supported for type: mpiaij!

[Barry Smith]

   Did you get the matrix in the writer code using a DA or some other code that puts the vector from the 2d or 3d problem in the natural ordering on the binary file.

   Barry

On Sep 1, 2010, at 9:55 AM, Leo van Kampenhout wrote:

> Hi all, 
> 
> I'm having trouble with MatLoad on Petsc-Dev. Not sure if this problem also exists in the main PETSc 3.1 since my code won't compile on that. 
> The problematic calls are as follows: 
> 
>       call PetscViewerBinaryOpen(PETSC_COMM_WORLD,file2,
>      &     FILE_MODE_READ,fd,ierr)
>       call DAGetMatrix(da,MATMPIAIJ,A,ierr)
>       call MatLoad(A,fd,ierr)
> 
> which supposedly will load a matrix A which is stored by the following commands:
> 
>       call PetscViewerBinaryOpen(PETSC_COMM_WORLD,file2,
>      &     FILE_MODE_WRITE,fd,ierr)
>       call MatView(A,fd,ierr)
> 
> by another program (on a single core). The MatLoad() does work perfectly on a single core, however two or more cores make the program crash with the following error (see below)
> The programcall is 
>  mpirun -n 2 ./main
> 
> Is something broken with Petsc or do I need to set some extra parameters? 
> Thanks in advance, 
> 
> Leo van Kampenhout
> 
> PS. I'm using a build from last friday, since when I updated petsc-dev today using Mercurial I couldnt get it to configure or compile :S (which is another subject)
> 
> 
> [0]PETSC ERROR: --------------------- Error Message ------------------------------------
> [0]PETSC ERROR: No support for this operation for this object type!
> [0]PETSC ERROR: MatLoad is not supported for type: mpiaij!
> [0]PETSC ERROR: ------------------------------------------------------------------------
> [0]PETSC ERROR: Petsc Development HG revision: 84fec3ca75412473b06251029b58f2a68f1334f3  HG Date: Fri Aug 27 07:06:08 2010 +0200
> [0]PETSC ERROR: See docs/changes/index.html for recent updates.
> [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
> [0]PETSC ERROR: See docs/index.html for manual pages.
> [0]PETSC ERROR: ------------------------------------------------------------------------
> [0]PETSC ERROR: ./main on a arch-linu named wingtip72 by csg4035 Wed Sep  1 15:52:11 2010
> [0]PETSC ERROR: Libraries linked from /net/users/csg/csg4035/install/petsc-dev/arch-linux-gnu-c-debug/lib
> [0]PETSC ERROR: Configure run at Fri Aug 27 08:37:32 2010
> [0]PETSC ERROR: Configure options 
> [0]PETSC ERROR: ------------------------------------------------------------------------
> [0]PETSC ERROR: MatLoad() line 843 in src/mat/interface/matrix.c
> 
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20100901/93150d91/attachment.html>