[petsc-users] Write complex matrix/vector from Matlab

[amrit poudel]

Hello everyone,

Does anyone here know how one can write a complex or real, matrix or vector (matrix could be sparse or full) in MATLAB  to a binary file  such that it can be read correctly by PETSc?
I tried using  PetscBinaryWrite.m that came with PETSc ( /bin/matlab) to write a simple real (or complex), full (or sparse) matrix in binary file named 'matrix.dat' and tried reading it in PETSc, but it absolutely does not work :(

 I did PetscBinaryWrite('matrix.dat', C) to write a matrix C in an empty file named 'matrix.dat' in MATLAB. Then I tried reading 'matrix.dat' file in PETSc using the following command :

PetscViewer        binviewer ;
Mat                    C ;
char            fileC[PETSC_MAX_PATH_LEN] ;

ierr = PetscOptionsGetString(PETSC_NULL, "-fC", fileC, PETSC_MAX_PATH_LEN-1, PETSC_NULL) ; CHKERRQ(ierr) ;

ierr = PetscViewerBinaryOpen(PETSC_COMM_WORLD, fileC, FILE_MODE_READ, &binviewer) ; CHKERRQ(ierr) ;

ierr = MatLoad(binviewer, MATMPIAIJ, &C) ; CHKERRQ(ierr) ;


I must add that I compiled PETSc with complex numbers.  

I plan on assembling my matrix (very large square sparse matrix, size greater than 1e+8) in MATLAB, which does not take too long to assemble in MATLAB at all. Then I will write the matrix to a binary file and read that in PETSc and solve a very large linear system of equations in parallel in PETSc using MPI. Hence, I really need to be able to write the matrix in MATLAB to a binary file  such that PETSc can correctly read the matrix.

I am really surprised to find out that PETSc cannot read/write a matrix or vector, full or sparse, real or complex in MATRIX MARKET FORMAT. I have gone though the manual, googled around quite a bit, but I just have not found anything useful. 

 I will truly appreciate any help or a pointer in this regard.
 
Thank you.

Amrit 

 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110905/3aa0f1e9/attachment.htm>