<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Apr 22, 2016 at 12:56 AM, Wei Deng <span dir="ltr"><<a href="mailto:deng106@purdue.edu" target="_blank">deng106@purdue.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div lang="EN-US" link="#0563C1" vlink="#954F72"><p class="MsoNormal"><span style="font-size:13pt;font-family:'Helvetica Neue',serif">Dear Prof. Knepley,<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:13pt;font-family:'Helvetica Neue',serif"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:13pt;font-family:'Helvetica Neue',serif">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.<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:13pt;font-family:'Helvetica Neue',serif"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:13pt;font-family:'Helvetica Neue',serif">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?</span></p></div></blockquote><div><br></div><div>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</div><div>we are just calling SPARSEPAK:</div><div><br></div><div>  <a href="https://bitbucket.org/petsc/petsc/src/78026435dbba2cd6fc5f6af7e1f769fa9012b35f/src/mat/order/genrcm.c?at=master&fileviewer=file-view-default">https://bitbucket.org/petsc/petsc/src/78026435dbba2cd6fc5f6af7e1f769fa9012b35f/src/mat/order/genrcm.c?at=master&fileviewer=file-view-default</a></div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> <span style="font-family:'Helvetica Neue',serif;font-size:13pt"> </span></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div lang="EN-US" link="#0563C1" vlink="#954F72"><p class="MsoNormal"><span style="font-size:13pt;font-family:'Helvetica Neue',serif">Really appreciate about that.<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:13pt;font-family:'Helvetica Neue',serif"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:13pt;font-family:'Helvetica Neue',serif"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:13pt;font-family:'Helvetica Neue',serif">Thanks a lot.<span class=""><font color="#888888"><u></u><u></u></font></span></span></p><span class=""><font color="#888888"><p class="MsoNormal"><span style="font-size:13pt;font-family:'Helvetica Neue',serif"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:13pt;font-family:'Helvetica Neue',serif"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:13pt;font-family:'Helvetica Neue',serif"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:13pt;font-family:'Helvetica Neue',serif">Wei Deng<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif"><u></u> <u></u></span></p></font></span></div></blockquote></div><br></div></div>