<div dir="ltr">On Fri, Apr 5, 2013 at 2:33 PM, 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"><div dir="ltr"><div><div><div><div>Hello <br></div> Do i am currently cloned petsc using git. Looks like the function you mentioned is in knepley/pylith. How do i switch to that branch ?<br>
</div><div>I guess i need to pull the branch before switching ?<br>
</div>git branch knepley/pylith<br></div>and <br></div>reconfigure ?<br><div class="gmail_extra"><br></div><div class="gmail_extra">May not be of (immediate) concern but an observation about the LocatePoint_<Simplex> codes. <br>
<br>I see that test is done for every cell by calculating the ref cell mapped coordinates for a given point. It might help speed up things if there is a test to first check if the given point is with in the square/cube containing the simplex. <br>
</div></div></blockquote><div><br></div><div style>That does not change the order of the method. So far, this had not been worth optimizing, since its a very small fraction of time (< 1%).</div><div style>When it gets optimized, we will use a hierarchical division like a kd-tree.</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 dir="ltr">
<div class="gmail_extra"></div><div class="gmail_extra">The above test will make sure that the point is inside the cell or with in one of the adjoining cells.<br></div><div class="gmail_extra"><br></div><div class="gmail_extra">
I am assuming the DMPlex object has information about cells adjoining a given cell. Since the LocatePoint_<Simplex> calculates the ref cell mapped coordinates, I may be wrong here but i am thinking we can use that information to check outside which facet of the refcell the given point is and then try to locate the point in the cell sharing that facet. <br>
<br></div><div class="gmail_extra">Thanks<br>Reddy <br></div><div class="gmail_extra"><br></div><div class="gmail_extra"><div class="gmail_quote">On Thu, Apr 4, 2013 at 3:37 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 Fri, Apr 5, 2013 at 6:29 AM, 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"><div dir="ltr"><div><div><div><div><div>Hello,<br></div> Is there any DMPlex function which will tell me to which DMPlex point a given coordinate point belongs to ? <br>
</div>I am looking for cell id given a coordinates of a point.</div></div></div></div></blockquote><div><br></div></div><div>I have experimental code to do this, DMPlexLocatePoint(). There are tests, but they are not that stringent for</div>
<div>hexes. Let me know if this works for you.</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"><div dir="ltr">
<div><div><div>If not, I can think of writing of a test for point in cell<br></div></div>and loop over cells to find the cell id. Any suggestions for doing this ? <br><br><br></div><div>The physical problem i am looking at is a Poisson equation in a given domain. <br>
<br></div><div>- \Delta u(x) = rho(x) where rho is a sum of delta functions randomly located in the domain. <br></div><span><font color="#888888"><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: <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></font></span></div>
</blockquote></div></div><span><font color="#888888"><br><br clear="all"><span class="HOEnZb"><font color="#888888"><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></font></span></div></div><span class="HOEnZb"><font color="#888888">
</font></span></blockquote></div><span class="HOEnZb"><font color="#888888"><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></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>