[MOAB-dev] vertices of edges

Nico Schlömer nico.schloemer at gmail.com
Mon Jan 4 09:14:00 CST 2016


Thanks, that's what I needed!

Cheers,
Nico

On Mon, Jan 4, 2016 at 3:23 PM Grindeanu, Iulian R. <iulian at mcs.anl.gov>
wrote:

> by default, the last argument in get_adjacencies is intersection, you have
> to use union, something like this
> rval = mcomm_->get_moab()->get_adjacencies(edges, 0, false, verts,
> moab::Interface::UNION );
>
> or you can use connectivity, it is easier
> rval = mcomm_->get_moab()->get_connectivity( edges, verts);
>
> ------------------------------
> *From:* moab-dev-bounces at mcs.anl.gov [moab-dev-bounces at mcs.anl.gov] on
> behalf of Nico Schlömer [nico.schloemer at gmail.com]
> *Sent:* Sunday, January 03, 2016 5:32 PM
> *To:* moab-dev at mcs.anl.gov
> *Subject:* [MOAB-dev] vertices of edges
>
> Given a bunch of edges, I'm trying to get the vertices on these edges. The
> output of
> ```
>   std::cout << "num boundary edges " << edges.size() << std::endl;
>   moab::Range verts;
>   rval = mcomm_->get_moab()->get_adjacencies(edges, 0, false, verts);
>   if (rval != moab::MB_SUCCESS) {
>     throw std::runtime_error("error in moab::get_adjacencies");
>   }
>   std::cout << "num boundary verts: " << verts.size() << std::endl;
> ```
> is
> ```
> num boundary edges 57
> num boundary verts: 0
> ```
> For the life of me I don't see what's wrong here. Any obvious mistakes?
>
> Cheers,
> Nico
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/moab-dev/attachments/20160104/52dbe614/attachment.html>


More information about the moab-dev mailing list