On Thu, Aug 23, 2012 at 12:19 PM, Chris Eldred <span dir="ltr"><<a href="mailto:chris.eldred@gmail.com" target="_blank">chris.eldred@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I am working with a 2D unstructured mesh using sieve and I wanted to<br>
get some advice on how to determine adjacency relations using the<br>
topological operators provided by Sieve (cone, support, closure,<br>
star). I have figured out the following:<br>
<br>
Vertices for a given edge: cone(edge)<br>
Cells for a given edge: support(edge)<br>
<br>
Cells for a given vertex: star(vertex) \ support(vertex)<br>
Edges for a given vertex: support(vertex)<br>
<br>
Edges for a given cell: cone(cell)<br>
Vertices for a given cell: closure(cell) \ cone(cell)<br>
<br>
The set of edges associated with the cells for a given edge: U<br>
cone(support(edge)) - ie the union of the cones of the support of the<br>
edge<br>
<br>
Are there other topological operators that more naturally express<br>
these relations (especially the ones involving unions or complements)?</blockquote><div><br></div><div>The only other builtin operations are meet and join. You are correct that you</div><div>want to use union there. We could add routines like</div>
<div><br></div><div>  DMComplexGetConeUnion(dm, [p], &unionSize, &unionArray);</div><div><br></div><div>but it was not clear that would do anything other than bloat the interface. I</div><div>think the best way to procede is to write these convenience routines using</div>
<div>the lower level primitives, and move them into PETSc if they turn out to be</div><div>useful. I can help you write them.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<span class="HOEnZb"><font color="#888888"><br>
--<br>
Chris Eldred<br>
DOE Computational Science Graduate Fellow<br>
Graduate Student, Atmospheric Science, Colorado State University<br>
B.S. Applied Computational Physics, Carnegie Mellon University, 2009<br>
<a href="mailto:chris.eldred@gmail.com">chris.eldred@gmail.com</a><br>
</font></span></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<br>