<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12pt"><div>Renumbering row will make an unknown not associate with its row (equation), will this affect the solver performance? For example, will the infill pattern change in ILU factorization?</div><div> </div><div>Thanks,</div><div>Qin </div><div><br></div> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 12pt;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 12pt;"> <div dir="ltr"> <div style="margin: 5px 0px; padding: 0px; border: 1px solid rgb(204, 204, 204); height: 0px; line-height: 0; font-size: 0px;" class="hr" contentEditable="false" readonly="true"></div> <font size="2" face="Arial"> <b><span style="font-weight: bold;">From:</span></b> Matthew
Knepley <knepley@gmail.com><br> <b><span style="font-weight: bold;">To:</span></b> Qin Lu <lu_qin_2000@yahoo.com> <br><b><span style="font-weight: bold;">Cc:</span></b> Barry Smith <bsmith@mcs.anl.gov>; petsc-users <petsc-users@mcs.anl.gov> <br> <b><span style="font-weight: bold;">Sent:</span></b> Wednesday, July 30, 2014 3:50 PM<br> <b><span style="font-weight: bold;">Subject:</span></b> Re: [petsc-users] Partition of parallel AIJ sparce matrix<br> </font> </div> <div class="y_msg_container"><br><div id="yiv7804996654"><div><div dir="ltr"><div class="yiv7804996654gmail_extra"><div class="yiv7804996654gmail_quote">On Wed, Jul 30, 2014 at 3:30 PM, Qin Lu <span dir="ltr"><<a href="mailto:lu_qin_2000@yahoo.com" shape="rect" rel="nofollow" target="_blank" ymailto="mailto:lu_qin_2000@yahoo.com">lu_qin_2000@yahoo.com</a>></span> wrote:<br clear="none">
<blockquote style="margin: 0px 0px 0px 0.8ex; padding-left: 1ex; border-left-color: rgb(204, 204, 204); border-left-width: 1px; border-left-style: solid;" class="yiv7804996654gmail_quote"><div style="color: rgb(0, 0, 0); font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);">
<div>In the context of domain decompostion, if the unknowns are ordered (to reduce the number of infills, for instance) in the way that each subdomain may not own consecutive unknown index, does this mean the partition of the domain will be different from the partition of the matrix? </div>
<div> </div><div>For example, if subdomain 1 (assigned to process 1) owns unknowns 1 and 3 (associated with equation 1 and 3), subdomain 2 (assigned to process 2) owns unknowns 2 and 4 (associated with equation 2 and 4) , how can I make each process own consecutive rows?</div>
</div></blockquote><div><br clear="none"></div><div>You renumber the rows once they are partitioned.</div><div><br clear="none"></div><div> Matt</div><div class="qtdSeparateBR"><br><br></div><div id="yiv7804996654yqtfd35508" class="yiv7804996654yqt4510240700"><div> </div></div><blockquote style="margin: 0px 0px 0px 0.8ex; padding-left: 1ex; border-left-color: rgb(204, 204, 204); border-left-width: 1px; border-left-style: solid;" class="yiv7804996654gmail_quote"><div id="yiv7804996654yqtfd81651" class="yiv7804996654yqt4510240700">
</div><div style="color: rgb(0, 0, 0); font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);"><div id="yiv7804996654yqtfd94662" class="yiv7804996654yqt4510240700"><div>Thanks,</div><div>Qin</div><div><br clear="none"></div> </div><div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 12pt;"><div id="yiv7804996654yqtfd41524" class="yiv7804996654yqt4510240700">
</div><div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 12pt;"><div id="yiv7804996654yqtfd84281" class="yiv7804996654yqt4510240700"> <div dir="ltr"> <div style="margin: 5px 0px; padding: 0px; border: 1px solid rgb(204, 204, 204); line-height: 0; font-size: 0px; min-height: 0px;">
</div> <font face="Arial"> <b><span style="font-weight: bold;">From:</span></b> Barry Smith <<a href="mailto:bsmith@mcs.anl.gov" shape="rect" rel="nofollow" target="_blank" ymailto="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>><br clear="none"> <b><span style="font-weight: bold;">To:</span></b> Qin Lu <<a href="mailto:lu_qin_2000@yahoo.com" shape="rect" rel="nofollow" target="_blank" ymailto="mailto:lu_qin_2000@yahoo.com">lu_qin_2000@yahoo.com</a>> <br clear="none">
<b><span style="font-weight: bold;">Cc:</span></b> petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" shape="rect" rel="nofollow" target="_blank" ymailto="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> <br clear="none"> <b><span style="font-weight: bold;">Sent:</span></b> Wednesday, July 30, 2014 2:49 PM<br clear="none">
<b><span style="font-weight: bold;">Subject:</span></b> Re: [petsc-users] Partition of parallel AIJ sparce matrix<br clear="none"> </font> </div> </div><div><div id="yiv7804996654yqtfd60774" class="yiv7804996654yqt4510240700"><br clear="none"><br clear="none">On Jul 30, 2014, at 11:08 AM, Qin Lu <<a href="mailto:lu_qin_2000@yahoo.com" shape="rect" rel="nofollow" target="_blank" ymailto="mailto:lu_qin_2000@yahoo.com">lu_qin_2000@yahoo.com</a>> wrote:<br clear="none">
<br clear="none">> Hello,<br clear="none">> <br clear="none">> Does a process have to own consecutive rows of the matrix? For example, suppose the global AIJ is 4x4, partitioned by 2 processes. Does process 1 have to own rows 1 and 2, process 2 own rows 3 and 4?<br clear="none">
<br clear="none"> Yes<br clear="none"><br clear="none">> Or process1 may own rows 1 and 3, and process 2 own row 2 and 4?<br clear="none"><br clear="none"> However, the numbering of degrees of freedom is arbitrary. Just renumber you degrees of freedom so the first set is on process 0, the next on process 1 etc.<br clear="none">
<br clear="none"> Barry</div><div><div id="yiv7804996654yqtfd03233" class="yiv7804996654yqt4510240700"><br clear="none"><br clear="none">> <br clear="none">> Thanks a lot for your help!<br clear="none">> <br clear="none">> Regards,<br clear="none">> Qin </div><br clear="none"></div>
<br clear="none"><br clear="none"></div> </div> </div> </div></blockquote></div><br clear="none"><br clear="all"><div><br clear="none"></div>-- <br clear="none">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br clear="none">
-- Norbert Wiener
</div></div></div></div><br><br></div> </div> </div> </div></body></html>