<div dir="ltr">On Tue, Mar 5, 2013 at 2:30 AM, Dharmendar Reddy <span dir="ltr"><<a href="mailto:dharmareddy84@gmail.com" target="_blank">dharmareddy84@gmail.com</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">On a related note, How does one handle auxilary data on the mesh ? Say in a Poisson problem, I want to assign dielectric constant to each cell. Of course i can create a indexset or vec tor but, I want the aux data to be distributed along with the mesh. </blockquote>
<div><br></div><div style>1) I have created an F90 version of DMPlexCreateSection()</div><div style><br></div><div style>2) My plan was to have them assigned on the distributed mesh, since this is scalable. However,</div>
<div style> assigning them on the serial mesh will be straightforward. You will use DMPlexDistributeField().</div><div style> The only thing you will need is the PetscSF created during distribution, which I will preserve, since</div>
<div style> right now it is destroyed at the end.</div><div style><br></div><div style> Thanks,</div><div style><br></div><div style> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="HOEnZb"><div class="h5"><div class="gmail_quote">On Tue, Mar 5, 2013 at 1:18 AM, Dharmendar Reddy <span dir="ltr"><<a href="mailto:dharmareddy84@gmail.com" target="_blank">dharmareddy84@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hello,<br> Trying to understand the DMPlexCreateSection first. I have created a test case, can you please have a look at it? Actually the code has a link error on my system with petsc-dev. I used DMPlexeCreateBoxMesh to create the mesh. Looks like there is no Fortran interface to that function. If i comment out the dmplexcreateboxmesh. It does create and executable.<br>
<br><span style="color:rgb(51,102,255)">petscDMTest_v2.o: In function `MAIN__':<br>petscDMTest_v2.F90:(.text+0x9b): undefined reference to `dmplexcreateboxmesh_'<br>make: [testDMMeshv2] Error 1 (ignored)</span><br>
<br>Also, If include finclude/petsc.h90 instead of finclude/petsch.h i get compile error :<br><br><span style="color:rgb(51,51,255)"><span>petscDMTest_v2.F90(73): error #6691: A pointer dummy argument may only be argument a<br>
ssociated with a pointer. [NUMCOMP]<br> call DMPlexCreateSection(dmMesh, topologyDim, numField, numComp, numDof, &<br>----------------------------------------------------------^<br>petscDMTest_v2.F90(73): error #6691: A pointer dummy argument may only be argument a<br>
ssociated with a pointer. [NUMDOF]<br> call DMPlexCreateSection(dmMesh, topologyDim, numField, numComp, numDof, &<br>-------------------------------------------------------------------^<br>petscDMTest_v2.F90(74): error #6691: A pointer dummy argument may only be argument a<br>
ssociated with a pointer. [BCFIELD]<br> numBC, bcField, bcPointIS,section,ierr)<br>----------------------------------^<br>petscDMTest_v2.F90(74): error #6691: A pointer dummy argument may only be argument a<br>
ssociated with a pointer. [BCPOINTIS]<br> numBC, bcField, bcPointIS,section,ierr)</span></span><br>-------------------------------------------^<br><br>Please have a look at the attached test case and Makefile<div>
<div><br>
<br><div class="gmail_quote">On Mon, Mar 4, 2013 at 7:39 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</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"><div>On Mon, Mar 4, 2013 at 7:50 PM, Dharmendar Reddy <span dir="ltr"><<a href="mailto:dharmareddy84@gmail.com" target="_blank">dharmareddy84@gmail.com</a>></span> wrote:<br></div><div class="gmail_extra">
<div class="gmail_quote"><div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br> Is there a DMPlex based fortran example for solving a pde such as laplace ? Even a c example will do. I am not able to understand the terminology of chart, cone etc. Is there a paper or user guide to learn more about the definitions of terms such as chart cone stratum ? I remember reading a presentation file by mathew about sieve. Is there a newer tutorial with terminology mapped to the new DMPlex/DMMesh functions. <br>
</blockquote><div><br></div></div><div>There is a fully worked out Stokes example (SNES ex62) and a Finite Volume example (TS ex11),</div><div>and there are a few tests. You should not need much of that stuff for defining a problem. I am</div>
<div>writing up a paper about ex11 right now. There is a short manual section on unstructured grids as well.</div><div><br></div><div>I think the key thing to understand is how to define a data layout over the mesh using a PetscSection.</div>
<div>Once that is done, all the DM functions work as normal, and the last hard part is coding your residual</div><div>loop. I have support for FEM schemes, but its undocumented and slightly complex in order to allow</div>
<div>GPU assembly, so you may be more comfortable with your own element loop to start.</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Thanks<span><font color="#888888"><br>Reddy <br><br clear="all"><br>-- <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: <a href="tel:%2B1-512-350-9082" value="+15123509082" target="_blank">+1-512-350-9082</a><br>
United States of America.<br>Homepage: <a href="https://webspace.utexas.edu/~dpr342" target="_blank">https://webspace.utexas.edu/~dpr342</a><br>
</font></span></blockquote></div></div><span><font color="#888888"><br><br clear="all"><div><br></div>-- <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
</font></span></div></div>
</blockquote></div><br><br clear="all"><br>-- <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: <a href="tel:%2B1-512-350-9082" value="+15123509082" target="_blank">+1-512-350-9082</a><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></blockquote></div><br><br clear="all"><br>-- <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: <a href="tel:%2B1-512-350-9082" value="+15123509082" target="_blank">+1-512-350-9082</a><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></blockquote></div><br><br clear="all"><div><br></div>-- <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
</div></div>