<div dir="ltr">Barry,<div>Thank you very much for your quick response.</div><div><br></div><div>The two links that you included clearly describe the formats, I have to write a reader for another code,</div><div>so if anyone has done this in a small code, any code snippet will be appreciated.</div><div>Cheers,</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Dec 29, 2020 at 1:15 PM Barry Smith <<a href="mailto:bsmith@petsc.dev">bsmith@petsc.dev</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;"><div><br></div>  <a href="https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Vec/VecLoad.html" target="_blank">https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Vec/VecLoad.html</a> and <a href="https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatLoad.html" target="_blank">https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatLoad.html</a> define the simple format. <div><br></div><div>  There are loaders for Matlab and Python distributed with PETSc so for those two languages you do not need to write your own custom loaders.</div><div><br></div><div>  Barry</div><div><br></div><div><br><div><br><blockquote type="cite"><div>On Dec 29, 2020, at 12:32 PM, C B <<a href="mailto:cebau.mail@gmail.com" target="_blank">cebau.mail@gmail.com</a>> wrote:</div><br><div><div dir="ltr">Hello PETSc community !<div><br></div><div>I would like to read in another code a matrix / rhs / sol  saved by PETSc.</div><div>I am dealing with large matrices and for precision I would like to use the binary format.</div><div><br></div><div>To save a matrix I am using the following lines (please let me know if I should make any changes), I am using an old version PETSc 2.1.5 for backward compatibility.</div><div></div><div><br></div><div> ierr = PetscViewerBinaryOpen(PETSC_COMM_WORLD,filename,PETSC_BINARY_CREATE,&petsc_view);CHKERRQ(ierr);<br> ierr = VecView(solut,petsc_view);CHKERRQ(ierr);<br>       ierr = PetscViewerDestroy(petsc_view);CHKERRQ(ierr);<br></div><div><br></div><div>And to save a vector:  </div><div><br></div><div> ierr = PetscViewerBinaryOpen(PETSC_COMM_WORLD, filename, PETSC_BINARY_CREATE, &petsc_view);CHKERRQ(ierr);<br> ierr = VecView(solut, petsc_view);CHKERRQ(ierr);<br>      ierr = PetscViewerDestroy(petsc_view);CHKERRQ(ierr);<br></div><div><br></div><div>The files are saved without any problems, and now I need to write the code to read the binary files.</div><div>Please, would anyone point me to the format of these binary files and/or a code snippet to take as an example ?  <br></div><div><br></div><div>Thank you very much in advance !</div><div>Cheers</div></div>
</div></blockquote></div><br></div></div></blockquote></div>