[MOAB-dev] Setting connectivity for single vertex

[Grindeanu, Iulian R.]

Hello,
Moab does not have elements with a single vertex. Is this what you are trying to do?
A vertex (node) is represented with an entity handle of type MBVERTEX, and you can get coordinates for it

I am not sure what are you trying to do.
Best Regards,
Iulian

________________________________
From: moab-dev-bounces at mcs.anl.gov [moab-dev-bounces at mcs.anl.gov] on behalf of CORTIS M. [michael.cortis at durham.ac.uk]
Sent: Thursday, January 29, 2015 10:18 AM
To: moab-dev at mcs.anl.gov
Cc: Tautges, Timothy J.
Subject: Re: [MOAB-dev] Setting connectivity for single vertex

Thanks Tim,

I am still having a problem, may be the guys at Argonne can help with this.

I've tried a simpler one and trying to crete the vertex from point.

  double points[9]={0.5,0,0};
  EntityHandle vertex;
  rval=moab2.create_vertex(points,vertex);CHKERR_PETSC(rval);
  EntityHandle dummy;
  rval = moab2.create_element(MBVERTEX,&vertex,1,dummy); CHKERR_PETSC(rval);

But I am getting Error code 2, MB_TYPE_OUT_OF_RANGE !
Is it possible to use MBVERTEX here?

Thanks,

Michael

On 29 Jan 2015, at 15:04, Tim Tautges <tautges at mcs.anl.gov<mailto:tautges at mcs.anl.gov>> wrote:

[cc'ng moab-dev list, since I've left Argonne.]

dummy is 0, which is an invalid handle, but should be the handle(s) of the entities whose connectivity you're setting; I suspect what you want there is *vertices.begin().  And you can't set connectivity until you make an element, so you need to call create_element, probably with the number of verts per element as 1 and the connectivity just the single vertex defining the point element.

- tim

On 01/29/2015 03:34 AM, CORTIS M. wrote:
Dear Tim,

I've got the following problem:

I want to create vertex entities, ie create the points and vertex entities using those points.
So I've create a set of vertices from an array of coordinates.
And If you try to set the connectivity for every point you get Error code 4 (MB_ENTITY_NOT_FOUND).

This is my code:

double points[6]={1,0,0, 0,1,0, 0,0,1};
Range vertices;
rval = moab.create_vertices(points,2,vertices); CHKERR_PETSC(rval);
EntityHandle dummy=0;
rval = moab.set_connectivity(dummy,vertexlist,1); CHKERR_PETSC(rval); //Error code 4

Any suggestions how to overcome this problem.

Kind regards,

Michael







-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/moab-dev/attachments/20150129/5f5d937c/attachment.html>