Skip to content
Snippets Groups Projects
Commit e3eb32a1 authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

*** empty log message ***

parent 365395dc
No related branches found
No related tags found
No related merge requests found
...@@ -599,20 +599,21 @@ static void applyOCCMeshConstraints(GModel *m, const void *constraints) ...@@ -599,20 +599,21 @@ static void applyOCCMeshConstraints(GModel *m, const void *constraints)
// set the mesh as immutable // set the mesh as immutable
ge->meshAttributes.Method == MESH_NONE; ge->meshAttributes.Method == MESH_NONE;
// set the correct tags on the boundary vertices // set the correct tags on the boundary vertices
Msg::Debug("Applying mesh contraints on edge %d first node %d", int numbeg = nodeNum.Value(1);
ge->tag(), nodeNum.Value(1)); int numend = nodeNum.Value(n);
ge->getBeginVertex()->mesh_vertices[0]->setNum(nodeNum.Value(1)); Msg::Debug("Applying mesh contraints on edge %d: beg=%d end=%d",
Msg::Debug("Applying mesh contraints on edge %d last node %d", ge->tag(), numbeg, numend);
ge->tag(), nodeNum.Value(n)); ge->getBeginVertex()->mesh_vertices[0]->setNum(numbeg);
ge->getEndVertex()->mesh_vertices[0]->setNum(nodeNum.Value(n)); ge->getEndVertex()->mesh_vertices[0]->setNum(numend);
// set the mesh on the edge // set the mesh on the edge
for(int i = 2; i < n; i++){ for(int i = 2; i < n; i++){
Msg::Debug("Applying mesh contraints on edge %d node %d", int num = nodeNum.Value(i);
ge->tag(), nodeNum.Value(i));
double u = nodePar.Value(i); double u = nodePar.Value(i);
GPoint p = ge->point(u); GPoint p = ge->point(u);
Msg::Debug("... adding vertex on edge %d: num=%d u=%g xyz=(%g,%g,%g)",
ge->tag(), num, u, p.x(), p.y(), p.z());
MEdgeVertex *v = new MEdgeVertex(p.x(), p.y(), p.z(), ge, u); MEdgeVertex *v = new MEdgeVertex(p.x(), p.y(), p.z(), ge, u);
v->setNum(nodeNum.Value(i)); v->setNum(num);
ge->mesh_vertices.push_back(v); ge->mesh_vertices.push_back(v);
} }
for(unsigned int i = 0; i < ge->mesh_vertices.size() + 1; i++){ for(unsigned int i = 0; i < ge->mesh_vertices.size() + 1; i++){
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment