<div dir="ltr">Thanks, that's what I needed!<div><br></div><div>Cheers,</div><div>Nico</div></div><br><div class="gmail_quote"><div dir="ltr">On Mon, Jan 4, 2016 at 3:23 PM Grindeanu, Iulian R. <<a href="mailto:iulian@mcs.anl.gov">iulian@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div>
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">by default, the last argument in get_adjacencies is intersection, you have to use union, something like this<br>
rval = mcomm_->get_moab()->get_adjacencies(edges, 0, false, verts, moab::Interface::UNION );<br>
<br>
or you can use connectivity, it is easier<br>
rval = mcomm_->get_moab()->get_connectivity( edges, verts);<br>
<br>
<div style="font-family:Times New Roman;color:#000000;font-size:16px">
<hr>
<div style="direction:ltr"><font size="2" color="#000000" face="Tahoma"><b>From:</b> <a href="mailto:moab-dev-bounces@mcs.anl.gov" target="_blank">moab-dev-bounces@mcs.anl.gov</a> [<a href="mailto:moab-dev-bounces@mcs.anl.gov" target="_blank">moab-dev-bounces@mcs.anl.gov</a>] on behalf of Nico Schlömer [<a href="mailto:nico.schloemer@gmail.com" target="_blank">nico.schloemer@gmail.com</a>]<br>
<b>Sent:</b> Sunday, January 03, 2016 5:32 PM<br>
<b>To:</b> <a href="mailto:moab-dev@mcs.anl.gov" target="_blank">moab-dev@mcs.anl.gov</a><br>
<b>Subject:</b> [MOAB-dev] vertices of edges<br>
</font><br>
</div></div></div></div><div><div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt"><div style="font-family:Times New Roman;color:#000000;font-size:16px">
<div></div>
<div>
<div dir="ltr">Given a bunch of edges, I'm trying to get the vertices on these edges. The output of
<div>```</div>
<div>
<div>  std::cout << "num boundary edges " << edges.size() << std::endl;</div>
<div>  moab::Range verts;</div>
<div>  rval = mcomm_->get_moab()->get_adjacencies(edges, 0, false, verts);</div>
<div>  if (rval != moab::MB_SUCCESS) {</div>
<div>    throw std::runtime_error("error in moab::get_adjacencies");</div>
<div>  }</div>
<div>  std::cout << "num boundary verts: " << verts.size() << std::endl;</div>
</div>
<div>```</div>
<div>is</div>
<div>```</div>
<div>
<div>num boundary edges 57</div>
<div>num boundary verts: 0</div>
</div>
<div>```</div>
<div>For the life of me I don't see what's wrong here. Any obvious mistakes?</div>
<div><br>
</div>
<div>Cheers,</div>
<div>Nico</div>
</div>
</div>
</div></div></div></blockquote></div>