[petsc-users] Loading Binary file problem

[Jin, Shuangshuang]

Hello everyone,

     I created a 4*3 matrix A in matlab, saved it to binary file "myA.dat" using PetscBinaryWrite('myA.dat',sparse(A)). Everything seems fine to me because I can read the binary file 'myA.dat' by PetscBinaryRead('myA.dat')correctly.

A = [1     2     0;
     2     0     5;
     0     3     0;
     1     4     1]
PetscBinaryWrite('myA.dat',sparse(A))
PetscBinaryRead('myA.dat')

However, when I tried to load the binary file "myA.dat" in PETSc using the example code "ex12.c" under "petsc-3.3-p0/src/mat/examples/tutorials", I got error message below:

[d3m956 at olympus tutorials]$ ./ex12 -f0 myA.dat
[0]PETSC ERROR: --------------------- Error Message ------------------------------------
[0]PETSC ERROR: Read from file failed!
[0]PETSC ERROR: Read past end of file!
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 0, Tue Jun  5 14:20:42 CDT 2012
[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: ./ex12 on a arch-linu named olympus.local by d3m956 Tue Jun 26 13:51:28 2012
[0]PETSC ERROR: Libraries linked from /pic/projects/mca/ss/PETSC/petsc-3.3-p0/arch-linux2-c-debug/lib
[0]PETSC ERROR: Configure run at Thu Jun 14 17:00:19 2012
[0]PETSC ERROR: Configure options --with-cc=gcc --with-fc=gfortran --download-f-blas-lapack --download-mpich
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: PetscBinaryRead() line 271 in src/sys/fileio/sysio.c
[0]PETSC ERROR: PetscBinarySynchronizedRead() line 558 in src/sys/fileio/sysio.c
[0]PETSC ERROR: PetscViewerBinaryRead() line 764 in src/sys/viewer/impls/binary/binv.c
[0]PETSC ERROR: PetscViewerBinaryReadVecHeader_Private() line 24 in src/vec/vec/utils/vecio.c
[0]PETSC ERROR: VecLoad_Binary() line 97 in src/vec/vec/utils/vecio.c
[0]PETSC ERROR: VecLoad_Default() line 348 in src/vec/vec/utils/vecio.c
[0]PETSC ERROR: VecLoad() line 1111 in src/vec/vec/interface/vector.c
[0]PETSC ERROR: main() line 92 in src/mat/examples/tutorials/ex12.c

It said "Read from file failed! Read past end of file!" Anything wrong? I was using the example code without changing anything. Please help me to figure out the problem.

Thanks,
Shuangshuang


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120626/b52eba6e/attachment.html>