The coordinates are stored according to the PETSc numbering. You get<br>them using<br><br> <a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/DA/DAGetCoordinates.html">http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/DA/DAGetCoordinates.html</a><br>
<br> Matt<br><br><div class="gmail_quote">On Fri, Feb 20, 2009 at 12:25 PM, Yixun Liu <span dir="ltr"><<a href="mailto:enjoywm@cs.wm.edu">enjoywm@cs.wm.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi,<br>
In ex2, the following codes Change the element vertex information to the<br>
PETSc numbering.<br>
<br>
705:AOCreateBasic <<a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-3.0.0/docs/manualpages/AO/AOCreateBasic.html#AOCreateBasic" target="_blank">http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-3.0.0/docs/manualpages/AO/AOCreateBasic.html#AOCreateBasic</a>>(PETSC_COMM_WORLD <<a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-3.0.0/docs/manualpages/Sys/PETSC_COMM_WORLD.html#PETSC_COMM_WORLD" target="_blank">http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-3.0.0/docs/manualpages/Sys/PETSC_COMM_WORLD.html#PETSC_COMM_WORLD</a>>,gdata->nlocal,gdata->localvert,PETSC_NULL <<a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-3.0.0/docs/manualpages/Sys/PETSC_NULL.html#PETSC_NULL" target="_blank">http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-3.0.0/docs/manualpages/Sys/PETSC_NULL.html#PETSC_NULL</a>>,&ao);<br>
<<a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-3.0.0/docs/manualpages/AO/AOCreateBasic.html#AOCreateBasic" target="_blank">http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-3.0.0/docs/manualpages/AO/AOCreateBasic.html#AOCreateBasic</a>>710:AOApplicationToPetsc <<a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-3.0.0/docs/manualpages/AO/AOApplicationToPetsc.html#AOApplicationToPetsc" target="_blank">http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-3.0.0/docs/manualpages/AO/AOApplicationToPetsc.html#AOApplicationToPetsc</a>>(ao,3*gdata->mlocal_ele,gdata->ele);<br>
<br>
However, I cannot find the codes which are used to create an array or vector, which stores the coordinates according to the PETSc numbering. Without this array how we get the coordinate of the vertex of an element?<br>
<br>
Thanks.<br>
<br>
<br>
</blockquote></div><br><br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener<br>