[MOAB-dev] vertices of edges

[Grindeanu, Iulian R.]

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/080fd820/attachment.html>