<div dir="ltr"><div><div>Hello,<br></div> I think i finally was able to go thorugh the process of solving a poisson problem using DMPlex object. Now i need to visualize the data. <br><br></div>I was looking at following lines from exmple 62:<br>
<br><pre width="80"><a name="line733">733: </a> <font color="#4169E1">if</font> (user.runType == RUN_FULL) {
<a name="line734">734: </a> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Viewer/PetscViewer.html#PetscViewer">PetscViewer</a> viewer;
<a name="line735">735: </a> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Vec/Vec.html#Vec">Vec</a> uLocal;
<a name="line737">737: </a> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Viewer/PetscViewerCreate.html#PetscViewerCreate">PetscViewerCreate</a>(<a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Sys/PETSC_COMM_WORLD.html#PETSC_COMM_WORLD">PETSC_COMM_WORLD</a>, &viewer);
<a name="line738">738: </a> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Viewer/PetscViewerSetType.html#PetscViewerSetType">PetscViewerSetType</a>(viewer, PETSCVIEWERVTK);
<a name="line739">739: </a> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Viewer/PetscViewerSetFormat.html#PetscViewerSetFormat">PetscViewerSetFormat</a>(viewer, PETSC_VIEWER_ASCII_VTK);
<a name="line740">740: </a> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Viewer/PetscViewerFileSetName.html#PetscViewerFileSetName">PetscViewerFileSetName</a>(viewer, <font color="#666666">"ex62_sol.vtk"</font>);
<a name="line742">742: </a> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/DM/DMGetLocalVector.html#DMGetLocalVector">DMGetLocalVector</a>(<a href="http://user.dm">user.dm</a>, &uLocal);
<a name="line743">743: </a> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/DM/DMGlobalToLocalBegin.html#DMGlobalToLocalBegin">DMGlobalToLocalBegin</a>(<a href="http://user.dm">user.dm</a>, u, <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Sys/INSERT_VALUES.html#INSERT_VALUES">INSERT_VALUES</a>, uLocal);
<a name="line744">744: </a> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/DM/DMGlobalToLocalEnd.html#DMGlobalToLocalEnd">DMGlobalToLocalEnd</a>(<a href="http://user.dm">user.dm</a>, u, <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Sys/INSERT_VALUES.html#INSERT_VALUES">INSERT_VALUES</a>, uLocal);
<a name="line746">746: </a> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Sys/PetscObjectReference.html#PetscObjectReference">PetscObjectReference</a>((<a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Sys/PetscObject.html#PetscObject">PetscObject</a>) <a href="http://user.dm">user.dm</a>); <font color="#B22222">/* Needed because viewer destroys the DM */</font>
<a name="line747">747: </a> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Sys/PetscObjectReference.html#PetscObjectReference">PetscObjectReference</a>((<a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Sys/PetscObject.html#PetscObject">PetscObject</a>) uLocal); <font color="#B22222">/* Needed because viewer destroys the <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Vec/Vec.html#Vec">Vec</a> */</font>
<a name="line748">748: </a> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Viewer/PetscViewerVTKAddField.html#PetscViewerVTKAddField">PetscViewerVTKAddField</a>(viewer, (<a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Sys/PetscObject.html#PetscObject">PetscObject</a>) <a href="http://user.dm">user.dm</a>, <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/DM/DMPlexVTKWriteAll.html#DMPlexVTKWriteAll">DMPlexVTKWriteAll</a>, PETSC_VTK_POINT_FIELD, (<a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Sys/PetscObject.html#PetscObject">PetscObject</a>) uLocal);
<a name="line749">749: </a> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/DM/DMRestoreLocalVector.html#DMRestoreLocalVector">DMRestoreLocalVector</a>(<a href="http://user.dm">user.dm</a>, &uLocal);
<a name="line750">750: </a> <a href="http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Viewer/PetscViewerDestroy.html#PetscViewerDestroy">PetscViewerDestroy</a>(&viewer);
<a name="line751">751: </a> }</pre>Can i do the same or some thing similar from fortran ? Some calls seem to have type casting. <br clear="all"><div><div><div><div><br></div><div>Othere wise if you can give me some idea on how to get the solution vector, including the boundary values.<br>
<br></div><div>As of now, i access the solution using solution vector passed to snes object but that does not give values at boundary nodes.<br><br></div><div>Thanks<br>Reddy<br></div><div>-- <br>-----------------------------------------------------<br>
Dharmendar Reddy Palle<br>Graduate Student<br>Microelectronics Research center,<br>University of Texas at Austin,<br>10100 Burnet Road, Bldg. 160<br>MER 2.608F, TX 78758-4445<br>e-mail: <a href="mailto:dharmareddy84@gmail.com" target="_blank">dharmareddy84@gmail.com</a><br>
Phone: +1-512-350-9082<br>United States of America.<br>Homepage: <a href="https://webspace.utexas.edu/~dpr342" target="_blank">https://webspace.utexas.edu/~dpr342</a><br>
</div></div></div></div></div>