<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><br><div>
Thanks. That worked out okay!<br></div><br><br><div><hr id="stopSpelling">Date: Wed, 7 Sep 2011 21:43:17 -0500<br>From: abhyshr@hawk.iit.edu<br>To: petsc-users@mcs.anl.gov<br>Subject: Re: [petsc-users] PETSc Error when writing a vector object to a binary file<br><br>Use VecView(), <a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Vec/VecView.html" target="_blank">http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Vec/VecView.html</a>, to write the vector to a binary file and then use PetscBinaryRead() in Matlab<br>
<div><br></div><div>ierr = PetscViewerBinaryOpen(PETSC_COMM_WORLD, filesol, FILE_MODE_WRITE, &binviewersol) ; CHKERRQ(ierr) ; </div><div>ierr = VecView(b,binviewersol);CHKERRQ(ierr);</div><div><br></div><div>In Matlab,</div>
<div>b = PetscBinaryRead(binviewersol,'complex', true);</div><div><br></div><div>See the function PetscBinaryRead() in $PETSC_DIR/bin/matlab</div><div><br><div class="ecxgmail_quote">On Wed, Sep 7, 2011 at 9:17 PM, amrit poudel <span dir="ltr"><<a href="mailto:amrit_pou@hotmail.com">amrit_pou@hotmail.com</a>></span> wrote:<br>
<blockquote class="ecxgmail_quote" style="border-left:1px #ccc solid;padding-left:1ex">
<div><div dir="ltr">
I am using the following sequence of commands to write a Vec b to a binary file and then eventually read it in MATLAB, but I get an error when I try writing this vector to a file in PETSc.<br><br>ierr = PetscViewerBinaryOpen(PETSC_COMM_WORLD, filesol, FILE_MODE_WRITE, &binviewersol) ; CHKERRQ(ierr) ; <br>
ierr = PetscViewerBinaryGetDescriptor(binviewersol, &fd); CHKERRQ(ierr) ;<br>ierr = PetscBinaryWrite(fd, b, 5, PETSC_SCALAR, PETSC_FALSE); CHKERRQ(ierr) ;<br><br><br><br>Here is the error :<br><br><br>[1]PETSC ERROR: Write to file failed!<br>
[1]PETSC ERROR: Error writing to file.!<br><br><br><br>NOTE : Petsc is configured with complex option and Vec b has complex entries. Length of Vec b is 5 .<br> </div></div>
</blockquote></div><br></div></div> </div></body>
</html>