<div dir="ltr"><a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatGetDiagonalBlock.html">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatGetDiagonalBlock.html</a><br></div><div class="gmail_extra">
<br><br><div class="gmail_quote">On Wed, Jan 16, 2013 at 3:51 PM, Fande Kong <span dir="ltr"><<a href="mailto:fd.kong@siat.ac.cn" target="_blank">fd.kong@siat.ac.cn</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Thanks, Jed,<div><br></div><div>Could I directly extract a local sub matrix from the parallel matrix? The local part should be a sequential matrix. <div><div class="h5"><br><br><div class="gmail_quote">On Wed, Jan 16, 2013 at 2:32 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">1. Partition the mesh<div>2. Create a sequential matrix representing local mesh adjacency</div><div>3. Use MATORDERINGRCM to get a new ordering for that sequential matrix</div>

<div>4. Broadcast the new global indices to neighbors</div>
</div><div><div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jan 16, 2013 at 2:26 PM, Fande Kong <span dir="ltr"><<a href="mailto:fd.kong@siat.ac.cn" target="_blank">fd.kong@siat.ac.cn</a>></span> wrote:<br>


<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Barry, Thank you very much.<div><br></div><div>What I want to do is to order the original mesh.  I have partitioned the mesh. But I want to reorder the vertices. Thus I create a matrix that represent the relationships of the vertices, and then I want to use the ordering methods in petsc to reorder the vertices by using matgetordering.</div>



<div><br></div><div>What should I do if I want to use the ordering methods (e.g. MATORDERINGRCM ) in petsc to reorder the mesh vertices?</div><div><div><div><br><div class="gmail_quote">On Wed, Jan 16, 2013 at 1:10 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br>



<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><br>
On Jan 16, 2013, at 1:53 PM, Fande Kong <<a href="mailto:fd.kong@siat.ac.cn" target="_blank">fd.kong@siat.ac.cn</a>> wrote:<br>
<br>
> Hi all,<br>
><br>
> I create a parallel mat with type MATMPIAIJ, and use the MatGetOrdering to reorder it. But it doesn't work. Who can tell me how to do that? I just want to reorder the local part of the matrix.<br>
<br>
</div>   We don't have code in place for this type of thing. Maybe if you explain why you would like it we could suggest alternatives.<br>
<br>
    Generally we recommend ordering the original mesh appropriately so that the resulting matrices have good orderings and don't recommend reordering the matrices directly (except for direct solvers).<br>
<span><font color="#888888"><br>
   Barry<br>
</font></span><div><div><br>
><br>
> Thanks,<br>
><br>
> --<br>
> Fande Kong<br>
> ShenZhen Institutes of Advanced Technology<br>
> Chinese Academy of Sciences<br>
<br>
<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div style="line-height:21px;font-size:14px;font-family:Verdana">Fande Kong</div><div style="line-height:21px;font-size:14px;font-family:Verdana">
ShenZhen Institutes of Advanced Technology</div><div style="line-height:21px;font-size:14px;font-family:Verdana">Chinese Academy of Sciences</div>
</div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div style="line-height:21px;font-size:14px;font-family:Verdana">Fande Kong</div><div style="line-height:21px;font-size:14px;font-family:Verdana">
ShenZhen Institutes of Advanced Technology</div><div style="line-height:21px;font-size:14px;font-family:Verdana">Chinese Academy of Sciences</div>
</div></div></div>
</blockquote></div><br></div>