[petsc-users] Problem about the MatGetOrdering in Petsc

Matthew Knepley knepley at rice.edu
Fri Apr 22 10:23:29 CDT 2016


On Fri, Apr 22, 2016 at 12:56 AM, Wei Deng <deng106 at purdue.edu> wrote:

> Dear Prof. Knepley,
>
>
>
> I am a student from Prof. Sameh’s class, I am now trying to use RCM
> (Reverse Cuthill Mckee) to transform a sparse matrix into a banded one.
>
>
>
> In a test case which can be perfectly transformed into band in matlab, I
> can’t achieve it in Petsc using the same data. The code can work, but I
> don’t know why does this problem appear. Could you spend some to give me
> some advice about this?
>

This is indeed RCM. I have cleaned up the code and taken snapshots of the
matrices. I am not sure what you are using in Matlab, but here
we are just calling SPARSEPAK:


https://bitbucket.org/petsc/petsc/src/78026435dbba2cd6fc5f6af7e1f769fa9012b35f/src/mat/order/genrcm.c?at=master&fileviewer=file-view-default

  Thanks,

     Matt


> Really appreciate about that.
>
>
>
>
>
> Thanks a lot.
>
>
>
>
>
>
>
> Wei Deng
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160422/3050bb09/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: main.c
Type: text/x-csrc
Size: 2292 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160422/3050bb09/attachment-0001.c>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 2016-04-22 at 10.20.31 AM.png
Type: image/png
Size: 22348 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160422/3050bb09/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 2016-04-22 at 10.20.39 AM.png
Type: image/png
Size: 22501 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160422/3050bb09/attachment-0003.png>


More information about the petsc-users mailing list