[petsc-users] Outputting matrix for viewing in matlab
Patrick Sanan
patrick.sanan at gmail.com
Fri Nov 29 02:16:19 CST 2019
> Am 29.11.2019 um 09:14 schrieb Patrick Sanan <patrick.sanan at gmail.com>:
>
> PETSc has its own binary format, which is not the same as MATLAB's.
>
> However, PETSc includes some MATLAB/Octave scripts which will load these binary files.
>
> See $PETSC_DIR/share/matlab/PetscBinaryRead.m - there are some examples in the comments at the top of that file.
Correction: $PETSC_DIR/share/petsc/matlab/PetscBinaryRead.m
>
>
> Note that you will probably want to add $PETSC_DIR/share/matlab to your MATLAB path so that you can run the script. This is what I have for Octave, but I'm not sure if it this, precisely, works in MATLAB:
>
> $ cat ~/.octaverc
> PETSC_DIR=getenv('PETSC_DIR');
> if length(PETSC_DIR)==0
> PETSC_DIR='~/code/petsc'
> end
> addpath([PETSC_DIR,'/share/petsc/matlab'])
>
> (As an aside, note that there are also scripts included to load PETSc binary files to use with numpy/scipy in Python, e.g. $PETSC_DIR/lib/petsc/bin/PetscBinaryIO.py)
>
>> Am 29.11.2019 um 04:07 schrieb baikadi pranay <pranayreddy865 at gmail.com <mailto:pranayreddy865 at gmail.com>>:
>>
>> Hello PETSc users,
>>
>> I have a sparse matrix built and I want to output the matrix for viewing in matlab. However i'm having difficulty outputting the matrix. I am writing my program in Fortran90 and I've included the following lines to output the matrix.
>>
>> call PetscViewerBinaryOpen(PETSC_COMM_SELF,'matrix',FILE_MODE_WRITE,view,ierr)
>> call PetscViewerBinaryGetDescriptor(view,fd,ierr)
>> call PetscBinaryWrite(fd,ham,1,PETSC_SCALAR,PETSC_FALSE,ierr)
>>
>> These lines do create a matrix but matlab says its not a binary file. Could you please provide me some inputs on where I'm going wrong and how to proceed with this problem. I can provide any further information that you might need to help me solve this problem.
>>
>>
>> Thank you.
>>
>> Sincerely,
>> Pranay.
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20191129/472d9367/attachment-0001.html>
More information about the petsc-users
mailing list