<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'><div dir='ltr'>
<br>Hello everyone,<br><br>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?<br>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 :(<br><br> 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 :<br><br>PetscViewer binviewer ;<br>Mat C ;<br>char fileC[PETSC_MAX_PATH_LEN] ;<br><br>ierr = PetscOptionsGetString(PETSC_NULL, "-fC", fileC, PETSC_MAX_PATH_LEN-1, PETSC_NULL) ; CHKERRQ(ierr) ;<br><br>ierr = PetscViewerBinaryOpen(PETSC_COMM_WORLD, fileC, FILE_MODE_READ, &binviewer) ; CHKERRQ(ierr) ;<br><br>ierr = MatLoad(binviewer, MATMPIAIJ, &C) ; CHKERRQ(ierr) ;<br><br><br>I must add that I compiled PETSc with complex numbers. <br><br>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.<br><br>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. <br><br> I will truly appreciate any help or a pointer in this regard.<br> <br>Thank you.<br><br>Amrit <br><br> </div></body>
</html>