<html><head><meta http-equiv="Content-Type" content="text/html charset=windows-1252"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Thanks for your quick reply.<br><div apple-content-edited="true">
<div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br></div></div><div><blockquote type="cite"><blockquote type="cite">I'm using Petsc interface to ParMetis, as detailed in section 3.5 of the user's manual.<br><br>I have a problem with the creation of the AO mapps using the sentence <br><br>AOCreateBasicIS(isg,PETSC NULL,&ao);<br></blockquote><br>This line in the user's manual is incorrect. See Section 2.3.1 for a<br>proper explanation.</blockquote><div><br></div><div>OK, good to know ;).</div><br><blockquote type="cite"><blockquote type="cite">where isg contains for each local node the new global (petsc) number of that node.<br>The problem is that the resulting AO map is inverted, and hence the results obtained with <br>AOApplicationToPetsc and AOPetscToApplication are not correct. <br></blockquote><br>You have the arguments backward.<br></blockquote><div><br></div><div><br></div><div>Yes, I agree, and this explains why the ApplicationToPetsc and PetscToApplication maps are inverted. </div><div>The problem is that you can neither call AOCreateBasicIS(PETSC_NULL,isg,&ao);</div><div>since it seems that PETSC_NULL is only allowed in the second argument.</div><div><br></div><blockquote type="cite"><a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/AO/AOCreateBasic.html">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/AO/AOCreateBasic.html</a><br><br>It's not clear from your description above whether you have an<br>"application ordering" or just a local to global map.<br></blockquote></div><br><div>I have an application ordering and a petsc ordering, both global, in the spirit of Section 2.3.1. </div><div>Here, the application ordering is the natural ordering (0, 1, 2, … ), resulting from a simple distribution </div><div>of the nodes across processors, before partition with ParMetis.</div><div><br></div><div>Best,</div><div>Miguel</div><div><br></div><div><br></div><div><br></div><div>--</div><div><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; font-size: 12px; border-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px; "><span class="Apple-style-span" style="font-family: 'Trebuchet MS', Tahoma, Arial, Helvetica, sans-serif; font-size: medium; ">INRIA Paris - Rocquencourt <br>Building 16, office 11</span></span></span></div><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; font-size: 12px; border-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px; "><span class="Apple-style-span" style="font-family: 'Trebuchet MS', Tahoma, Arial, Helvetica, sans-serif; font-size: medium; ">78153 Le Chesnay Cedex (France)<br>phone: +33 (0)1 3963 5470<br>fax: +33 (0)1 3963 5882<br><a href="http://www-rocq.inria.fr/who/Miguel.Fernandez" style="text-decoration: none; ">www-rocq.inria.fr/who/Miguel.Fernandez</a></span></span></span></div></div></body></html>