Standard_Real y0 = xyzPoint.Y();
Standard_Real z0 = xyzPoint.Z();
MESSAGE("Projected Vertex: " << x0 << ", " << y0 << ", " << z0);
+ MESSAGE("Parametric coordinates: " << u0 << ", " << v0 );
coords.push_back(u0);
coords.push_back(v0);
coords.push_back(x0);
MESSAGE("Face indice: " << iface);
MESSAGE("Adding enforced vertices");
evl = evmIt->second;
- MESSAGE("Number of vertices to add: "<< evl.size())
+ MESSAGE("Number of vertices to add: "<< evl.size());
std::set< std::vector<double> >::const_iterator evlIt = evl.begin();
for (; evlIt != evl.end(); ++evlIt) {
// ev = *evlIt;
mesh_get_vertex_coordinates(msh, iv, xyz);
mesh_get_vertex_tag(msh, iv, &tag);
// Issue 0020656. Use vertex coordinates
- if ( tag ) {
- gp_Pnt p = BRep_Tool::Pnt( TopoDS::Vertex(pmap(tag)));
- xyz[0] = p.X(); xyz[1] = p.Y(); xyz[2] = p.Z();
- }
+ // GDD: commented because enforced vertices
+ // are badly taken into account
+ //if ( tag ) {
+ // gp_Pnt p = BRep_Tool::Pnt( TopoDS::Vertex(pmap(tag)));
+ // xyz[0] = p.X(); xyz[1] = p.Y(); xyz[2] = p.Z();
+ //}
nodes[iv] = meshDS->AddNode(xyz[0], xyz[1], xyz[2]);
// internal point are tagged to zero
if(tag){