[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