-// Copyright (C) 2007-2014 CEA/DEN, EDF R&D
+// Copyright (C) 2007-2015 CEA/DEN, EDF R&D
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
#include <SMESH_MesherHelper.hxx>
#include <StdMeshers_FaceSide.hxx>
#include <StdMeshers_ViscousLayers2D.hxx>
+#include <SMESH_File.hxx>
#include <utilities.h>
*/
//=============================================================================
-inline std::string to_string(double d)
+inline std::string val_to_string(double d)
{
std::ostringstream o;
o << d;
return o.str();
}
-inline std::string to_string_rel(double d)
+inline std::string val_to_string_rel(double d)
{
std::ostringstream o;
o << d;
return o.str();
}
-inline std::string to_string(int i)
+inline std::string val_to_string(int i)
{
std::ostringstream o;
o << i;
return o.str();
}
-inline std::string to_string_rel(int i)
+inline std::string val_to_string_rel(int i)
{
std::ostringstream o;
o << i;
}
double _smp_phy_size;
-// #if BLSURF_VERSION_LONG >= "3.1.1"
-// // sizemap_t *geo_sizemap_e, *geo_sizemap_f;
-// sizemap_t *iso_sizemap_p, *iso_sizemap_e, *iso_sizemap_f;
-// // sizemap_t *clean_geo_sizemap_e, *clean_geo_sizemap_f;
-// sizemap_t *clean_iso_sizemap_p, *clean_iso_sizemap_e, *clean_iso_sizemap_f;
-// #endif
status_t size_on_surface(integer face_id, real *uv, real *size, void *user_data);
status_t size_on_edge(integer edge_id, real t, real *size, void *user_data);
status_t size_on_vertex(integer vertex_id, real *size, void *user_data);
return face_ids;
}
+BLSURFPlugin_BLSURF::TListOfIDs _getSubShapeIDsInMainShape(SMESH_Mesh* theMesh,
+ TopoDS_Shape theSubShape,
+ TopAbs_ShapeEnum theShapeType)
+{
+ BLSURFPlugin_BLSURF::TListOfIDs face_ids;
+
+ for (TopExp_Explorer face_iter(theSubShape,theShapeType);face_iter.More();face_iter.Next())
+ {
+ int face_id = theMesh->GetMeshDS()->ShapeToIndex(face_iter.Current());
+ if (face_id == 0)
+ throw SALOME_Exception ( SMESH_Comment("Sub_shape not found in main_shape"));
+ face_ids.push_back(face_id);
+ }
+
+ return face_ids;
+}
+
void BLSURFPlugin_BLSURF::addCoordsFromVertices(const std::vector<std::string> &theVerticesEntries, std::vector<double> &theVerticesCoords)
{
for (std::vector<std::string>::const_iterator it = theVerticesEntries.begin(); it != theVerticesEntries.end(); it++)
/////////////////////////////////////////////////////////
void BLSURFPlugin_BLSURF::createVerticesPeriodicity(TopoDS_Shape theGeomShape, BLSURFPlugin_Hypothesis::TEntry theEdge1, BLSURFPlugin_Hypothesis::TEntry theVertex1,
- BLSURFPlugin_Hypothesis::TEntry theEdge2, BLSURFPlugin_Hypothesis::TEntry theVertex2)
+ BLSURFPlugin_Hypothesis::TEntry theEdge2, BLSURFPlugin_Hypothesis::TEntry theVertex2)
{
MESSAGE("BLSURFPlugin_BLSURF::createVerticesPeriodicity");
/////////////////////////////////////////////////////////
-void BLSURFPlugin_BLSURF::SetParameters(
-// #if BLSURF_VERSION_LONG >= "3.1.1"
-// cad_t * c,
-// #endif
- const BLSURFPlugin_Hypothesis* hyp,
+void BLSURFPlugin_BLSURF::SetParameters(const BLSURFPlugin_Hypothesis* hyp,
cadsurf_session_t * css,
precad_session_t * pcs,
const TopoDS_Shape& theGeomShape,
- bool * use_precad
- )
+ bool * use_precad
+ )
{
// rnc : Bug 1457
// Clear map so that it is not stored in the algorithm with old enforced vertices in it
EnfVertexCoords2EnfVertexList.clear();
-
- double diagonal = SMESH_Mesh::GetShapeDiagonalSize( theGeomShape );
- double bbSegmentation = _gen->GetBoundaryBoxSegmentation();
- int _physicalMesh = BLSURFPlugin_Hypothesis::GetDefaultPhysicalMesh();
- int _geometricMesh = BLSURFPlugin_Hypothesis::GetDefaultGeometricMesh();
- double _phySize = BLSURFPlugin_Hypothesis::GetDefaultPhySize(diagonal, bbSegmentation);
- bool _phySizeRel = BLSURFPlugin_Hypothesis::GetDefaultPhySizeRel();
- double _minSize = BLSURFPlugin_Hypothesis::GetDefaultMinSize(diagonal);
- bool _minSizeRel = BLSURFPlugin_Hypothesis::GetDefaultMinSizeRel();
- double _maxSize = BLSURFPlugin_Hypothesis::GetDefaultMaxSize(diagonal);
- bool _maxSizeRel = BLSURFPlugin_Hypothesis::GetDefaultMaxSizeRel();
- double _gradation = BLSURFPlugin_Hypothesis::GetDefaultGradation();
- bool _quadAllowed = BLSURFPlugin_Hypothesis::GetDefaultQuadAllowed();
- double _angleMesh = BLSURFPlugin_Hypothesis::GetDefaultAngleMesh();
- double _chordalError = BLSURFPlugin_Hypothesis::GetDefaultChordalError(diagonal);
- bool _anisotropic = BLSURFPlugin_Hypothesis::GetDefaultAnisotropic();
- double _anisotropicRatio = BLSURFPlugin_Hypothesis::GetDefaultAnisotropicRatio();
- bool _removeTinyEdges = BLSURFPlugin_Hypothesis::GetDefaultRemoveTinyEdges();
- double _tinyEdgeLength = BLSURFPlugin_Hypothesis::GetDefaultTinyEdgeLength(diagonal);
- bool _badElementRemoval = BLSURFPlugin_Hypothesis::GetDefaultBadElementRemoval();
- double _badElementAspectRatio = BLSURFPlugin_Hypothesis::GetDefaultBadElementAspectRatio();
- bool _optimizeMesh = BLSURFPlugin_Hypothesis::GetDefaultOptimizeMesh();
- bool _quadraticMesh = BLSURFPlugin_Hypothesis::GetDefaultQuadraticMesh();
- int _verb = BLSURFPlugin_Hypothesis::GetDefaultVerbosity();
- int _topology = BLSURFPlugin_Hypothesis::GetDefaultTopology();
+
+ double diagonal = SMESH_Mesh::GetShapeDiagonalSize( theGeomShape );
+ double bbSegmentation = _gen->GetBoundaryBoxSegmentation();
+ int _physicalMesh = BLSURFPlugin_Hypothesis::GetDefaultPhysicalMesh();
+ int _geometricMesh = BLSURFPlugin_Hypothesis::GetDefaultGeometricMesh();
+ double _phySize = BLSURFPlugin_Hypothesis::GetDefaultPhySize(diagonal, bbSegmentation);
+ bool _phySizeRel = BLSURFPlugin_Hypothesis::GetDefaultPhySizeRel();
+ double _minSize = BLSURFPlugin_Hypothesis::GetDefaultMinSize(diagonal);
+ bool _minSizeRel = BLSURFPlugin_Hypothesis::GetDefaultMinSizeRel();
+ double _maxSize = BLSURFPlugin_Hypothesis::GetDefaultMaxSize(diagonal);
+ bool _maxSizeRel = BLSURFPlugin_Hypothesis::GetDefaultMaxSizeRel();
+ double _gradation = BLSURFPlugin_Hypothesis::GetDefaultGradation();
+ bool _quadAllowed = BLSURFPlugin_Hypothesis::GetDefaultQuadAllowed();
+ double _angleMesh = BLSURFPlugin_Hypothesis::GetDefaultAngleMesh();
+ double _chordalError = BLSURFPlugin_Hypothesis::GetDefaultChordalError(diagonal);
+ bool _anisotropic = BLSURFPlugin_Hypothesis::GetDefaultAnisotropic();
+ double _anisotropicRatio = BLSURFPlugin_Hypothesis::GetDefaultAnisotropicRatio();
+ bool _removeTinyEdges = BLSURFPlugin_Hypothesis::GetDefaultRemoveTinyEdges();
+ double _tinyEdgeLength = BLSURFPlugin_Hypothesis::GetDefaultTinyEdgeLength(diagonal);
+ bool _badElementRemoval = BLSURFPlugin_Hypothesis::GetDefaultBadElementRemoval();
+ double _badElementAspectRatio = BLSURFPlugin_Hypothesis::GetDefaultBadElementAspectRatio();
+ bool _optimizeMesh = BLSURFPlugin_Hypothesis::GetDefaultOptimizeMesh();
+ bool _quadraticMesh = BLSURFPlugin_Hypothesis::GetDefaultQuadraticMesh();
+ int _verb = BLSURFPlugin_Hypothesis::GetDefaultVerbosity();
+ int _topology = BLSURFPlugin_Hypothesis::GetDefaultTopology();
// PreCAD
- int _precadMergeEdges = BLSURFPlugin_Hypothesis::GetDefaultPreCADMergeEdges();
- int _precadProcess3DTopology = BLSURFPlugin_Hypothesis::GetDefaultPreCADProcess3DTopology();
- int _precadDiscardInput = BLSURFPlugin_Hypothesis::GetDefaultPreCADDiscardInput();
+ int _precadMergeEdges = BLSURFPlugin_Hypothesis::GetDefaultPreCADMergeEdges();
+ int _precadProcess3DTopology = BLSURFPlugin_Hypothesis::GetDefaultPreCADProcess3DTopology();
+ int _precadDiscardInput = BLSURFPlugin_Hypothesis::GetDefaultPreCADDiscardInput();
if (hyp) {
}
if (hyp->GetGradation() > 0)
_gradation = hyp->GetGradation();
- _quadAllowed = hyp->GetQuadAllowed();
- if (hyp->GetAngleMesh() > 0)
- _angleMesh = hyp->GetAngleMesh();
- if (hyp->GetChordalError() > 0)
- _chordalError = hyp->GetChordalError();
- _anisotropic = hyp->GetAnisotropic();
- if (hyp->GetAnisotropicRatio() >= 0)
- _anisotropicRatio = hyp->GetAnisotropicRatio();
- _removeTinyEdges = hyp->GetRemoveTinyEdges();
- if (hyp->GetTinyEdgeLength() > 0)
- _tinyEdgeLength = hyp->GetTinyEdgeLength();
- _badElementRemoval = hyp->GetBadElementRemoval();
- if (hyp->GetBadElementAspectRatio() >= 0)
- _badElementAspectRatio = hyp->GetBadElementAspectRatio();
- _optimizeMesh = hyp->GetOptimizeMesh();
- _quadraticMesh = hyp->GetQuadraticMesh();
+ _quadAllowed = hyp->GetQuadAllowed();
+ if (hyp->GetAngleMesh() > 0)
+ _angleMesh = hyp->GetAngleMesh();
+ if (hyp->GetChordalError() > 0)
+ _chordalError = hyp->GetChordalError();
+ _anisotropic = hyp->GetAnisotropic();
+ if (hyp->GetAnisotropicRatio() >= 0)
+ _anisotropicRatio = hyp->GetAnisotropicRatio();
+ _removeTinyEdges = hyp->GetRemoveTinyEdges();
+ if (hyp->GetTinyEdgeLength() > 0)
+ _tinyEdgeLength = hyp->GetTinyEdgeLength();
+ _badElementRemoval = hyp->GetBadElementRemoval();
+ if (hyp->GetBadElementAspectRatio() >= 0)
+ _badElementAspectRatio = hyp->GetBadElementAspectRatio();
+ _optimizeMesh = hyp->GetOptimizeMesh();
+ _quadraticMesh = hyp->GetQuadraticMesh();
_verb = hyp->GetVerbosity();
- _topology = (int) hyp->GetTopology();
- // PreCAD
- _precadMergeEdges = hyp->GetPreCADMergeEdges();
- _precadProcess3DTopology = hyp->GetPreCADProcess3DTopology();
- _precadDiscardInput = hyp->GetPreCADDiscardInput();
+ _topology = (int) hyp->GetTopology();
+ // PreCAD
+ _precadMergeEdges = hyp->GetPreCADMergeEdges();
+ _precadProcess3DTopology = hyp->GetPreCADProcess3DTopology();
+ _precadDiscardInput = hyp->GetPreCADDiscardInput();
- const BLSURFPlugin_Hypothesis::TOptionValues & opts = hyp->GetOptionValues();
+ const BLSURFPlugin_Hypothesis::TOptionValues& opts = hyp->GetOptionValues();
BLSURFPlugin_Hypothesis::TOptionValues::const_iterator opIt;
for ( opIt = opts.begin(); opIt != opts.end(); ++opIt )
if ( !opIt->second.empty() ) {
MESSAGE("cadsurf_set_param(): " << opIt->first << " = " << opIt->second);
set_param(css, opIt->first.c_str(), opIt->second.c_str());
}
-
- const BLSURFPlugin_Hypothesis::TOptionValues & preCADopts = hyp->GetPreCADOptionValues();
+
+ const BLSURFPlugin_Hypothesis::TOptionValues& custom_opts = hyp->GetCustomOptionValues();
+ for ( opIt = custom_opts.begin(); opIt != custom_opts.end(); ++opIt )
+ if ( !opIt->second.empty() ) {
+ MESSAGE("cadsurf_set_param(): " << opIt->first << " = " << opIt->second);
+ set_param(css, opIt->first.c_str(), opIt->second.c_str());
+ }
+
+ const BLSURFPlugin_Hypothesis::TOptionValues& preCADopts = hyp->GetPreCADOptionValues();
for ( opIt = preCADopts.begin(); opIt != preCADopts.end(); ++opIt )
if ( !opIt->second.empty() ) {
*use_precad = true;
MESSAGE("precad_set_param(): " << opIt->first << " = " << opIt->second);
precad_set_param(pcs, opIt->first.c_str(), opIt->second.c_str());
}
+
+ const BLSURFPlugin_Hypothesis::TOptionValues& custom_preCADopts = hyp->GetCustomPreCADOptionValues();
+ for ( opIt = custom_preCADopts.begin(); opIt != custom_preCADopts.end(); ++opIt )
+ if ( !opIt->second.empty() ) {
+ *use_precad = true;
+ MESSAGE("precad_set_param(): " << opIt->first << " = " << opIt->second);
+ precad_set_param(pcs, opIt->first.c_str(), opIt->second.c_str());
+ }
}
// else {
// //0020968: EDF1545 SMESH: Problem in the creation of a mesh group on geometry
// PreCAD
if (_topology == BLSURFPlugin_Hypothesis::PreCAD) {
*use_precad = true;
- precad_set_param(pcs, "verbose", to_string(_verb).c_str());
+ precad_set_param(pcs, "verbose", val_to_string(_verb).c_str());
precad_set_param(pcs, "merge_edges", _precadMergeEdges ? "1" : "0");
precad_set_param(pcs, "process_3d_topology", _precadProcess3DTopology ? "1" : "0");
precad_set_param(pcs, "discard_input_topology", _precadDiscardInput ? "1" : "0");
{
case BLSURFPlugin_Hypothesis::PhysicalGlobalSize:
set_param(css, "physical_size_mode", "global");
- set_param(css, "global_physical_size", _phySizeRel ? to_string_rel(_phySize).c_str() : to_string(_phySize).c_str());
+ set_param(css, "global_physical_size", _phySizeRel ? val_to_string_rel(_phySize).c_str() : val_to_string(_phySize).c_str());
break;
case BLSURFPlugin_Hypothesis::PhysicalLocalSize:
set_param(css, "physical_size_mode", "local");
- set_param(css, "global_physical_size", _phySizeRel ? to_string_rel(_phySize).c_str() : to_string(_phySize).c_str());
+ set_param(css, "global_physical_size", _phySizeRel ? val_to_string_rel(_phySize).c_str() : val_to_string(_phySize).c_str());
useGradation = true;
break;
default:
{
case BLSURFPlugin_Hypothesis::GeometricalGlobalSize:
set_param(css, "geometric_size_mode", "global");
- set_param(css, "geometric_approximation", to_string(_angleMesh).c_str());
- set_param(css, "chordal_error", to_string(_chordalError).c_str());
+ set_param(css, "geometric_approximation", val_to_string(_angleMesh).c_str());
+ set_param(css, "chordal_error", val_to_string(_chordalError).c_str());
useGradation = true;
break;
case BLSURFPlugin_Hypothesis::GeometricalLocalSize:
set_param(css, "geometric_size_mode", "local");
- set_param(css, "geometric_approximation", to_string(_angleMesh).c_str());
- set_param(css, "chordal_error", to_string(_chordalError).c_str());
+ set_param(css, "geometric_approximation", val_to_string(_angleMesh).c_str());
+ set_param(css, "chordal_error", val_to_string(_chordalError).c_str());
useGradation = true;
break;
default:
double mins = _minSizeRel ? _minSize * diagonal : _minSize;
// - min size should not be greater than user size
if ( _phySize < mins )
- set_param(css, "min_size", _phySizeRel ? to_string_rel(_phySize).c_str() : to_string(_phySize).c_str());
+ set_param(css, "min_size", _phySizeRel ? val_to_string_rel(_phySize).c_str() : val_to_string(_phySize).c_str());
else
- set_param(css, "min_size", _minSizeRel ? to_string_rel(_minSize).c_str() : to_string(_minSize).c_str());
+ set_param(css, "min_size", _minSizeRel ? val_to_string_rel(_minSize).c_str() : val_to_string(_minSize).c_str());
// - compute absolute max size
double maxs = _maxSizeRel ? _maxSize * diagonal : _maxSize;
// - max size should not be less than user size
if ( _phySize > maxs )
- set_param(css, "max_size", _phySizeRel ? to_string_rel(_phySize).c_str() : to_string(_phySize).c_str());
+ set_param(css, "max_size", _phySizeRel ? val_to_string_rel(_phySize).c_str() : val_to_string(_phySize).c_str());
else
- set_param(css, "max_size", _maxSizeRel ? to_string_rel(_maxSize).c_str() : to_string(_maxSize).c_str());
+ set_param(css, "max_size", _maxSizeRel ? val_to_string_rel(_maxSize).c_str() : val_to_string(_maxSize).c_str());
}
else {
// user size is not explicitly specified
// - if minsize is not explicitly specified, we pass default value computed automatically, in this case "relative" flag is ignored
- set_param(css, "min_size", _minSizeRel ? to_string_rel(_minSize).c_str() : to_string(_minSize).c_str());
+ set_param(css, "min_size", _minSizeRel ? val_to_string_rel(_minSize).c_str() : val_to_string(_minSize).c_str());
// - if maxsize is not explicitly specified, we pass default value computed automatically, in this case "relative" flag is ignored
- set_param(css, "max_size", _maxSizeRel ? to_string_rel(_maxSize).c_str() : to_string(_maxSize).c_str());
+ set_param(css, "max_size", _maxSizeRel ? val_to_string_rel(_maxSize).c_str() : val_to_string(_maxSize).c_str());
}
// anisotropic and quadrangle mesh requires disabling gradation
if ( _anisotropic && _quadAllowed )
useGradation = false; // limitation of V1.3
if ( useGradation )
- set_param(css, "gradation", to_string(_gradation).c_str());
+ set_param(css, "gradation", val_to_string(_gradation).c_str());
set_param(css, "element_generation", _quadAllowed ? "quad_dominant" : "triangle");
set_param(css, "metric", _anisotropic ? "anisotropic" : "isotropic");
if ( _anisotropic )
- set_param(css, "anisotropic_ratio", to_string(_anisotropicRatio).c_str());
+ set_param(css, "anisotropic_ratio", val_to_string(_anisotropicRatio).c_str());
set_param(css, "remove_tiny_edges", _removeTinyEdges ? "1" : "0");
if ( _removeTinyEdges )
- set_param(css, "tiny_edge_length", to_string(_tinyEdgeLength).c_str());
+ set_param(css, "tiny_edge_length", val_to_string(_tinyEdgeLength).c_str());
set_param(css, "force_bad_surface_element_removal", _badElementRemoval ? "1" : "0");
if ( _badElementRemoval )
- set_param(css, "bad_surface_element_aspect_ratio", to_string(_badElementAspectRatio).c_str());
+ set_param(css, "bad_surface_element_aspect_ratio", val_to_string(_badElementAspectRatio).c_str());
set_param(css, "optimisation", _optimizeMesh ? "yes" : "no");
set_param(css, "element_order", _quadraticMesh ? "quadratic" : "linear");
- set_param(css, "verbose", to_string(_verb).c_str());
+ set_param(css, "verbose", val_to_string(_verb).c_str());
_smp_phy_size = _phySizeRel ? _phySize*diagonal : _phySize;
if ( _verb > 0 )
{
MESSAGE("USING PRECAD FOR PERIODICITY")
*use_precad = true;
- precad_set_param(pcs, "verbose", to_string(_verb).c_str());
+ precad_set_param(pcs, "verbose", val_to_string(_verb).c_str());
}
MESSAGE("SetParameters facesPeriodicityVector");
/* create a distene context (generic object) */
status_t status = STATUS_ERROR;
+ myMesh = &aMesh;
SMESHDS_Mesh* meshDS = aMesh.GetMeshDS();
SMESH_MesherHelper helper( aMesh );
// do not call helper.IsQuadraticSubMesh() because sub-meshes
TSubMeshSet edgeSubmeshes;
TSubMeshSet& mergeSubmeshes = edgeSubmeshes;
- TopTools_IndexedMapOfShape fmap;
- TopTools_IndexedMapOfShape emap;
- TopTools_IndexedMapOfShape pmap;
+ TopTools_IndexedMapOfShape pmap, emap, fmap;
// Issue 0019864. On DebianSarge, FE signals do not obey to OSD::SetSignal(false)
#ifndef WIN32
MESSAGE("BEGIN SetParameters");
bool use_precad = false;
- SetParameters(
- // #if BLSURF_VERSION_LONG >= "3.1.1"
- // c,
- // #endif
- _hypothesis, css, pcs, aShape, &use_precad);
+ SetParameters(_hypothesis, css, pcs, aShape, &use_precad);
MESSAGE("END SetParameters");
MESSAGE("_preCadFacesIDsPeriodicityVector.size() = " << _preCadFacesIDsPeriodicityVector.size());
vector<Handle(Geom2d_Curve)> curves;
vector<Handle(Geom_Surface)> surfaces;
- fmap.Clear();
emap.Clear();
pmap.Clear();
FaceId2PythonSmp.clear();
if (f.Orientation() != TopAbs_FORWARD && f.Orientation() != TopAbs_REVERSED )
f.Orientation(TopAbs_FORWARD);
- if (fmap.FindIndex(f) > 0)
- continue;
iface = fmap.Add(f);
// std::string aFileName = "fmap_face_";
-// aFileName.append(to_string(iface));
+// aFileName.append(val_to_string(iface));
// aFileName.append(".brep");
// BRepTools::Write(f,aFileName.c_str());
cad_face_t *fce = cad_face_new(c, iface, surf_fun, surfaces.back());
/* by default a face has no tag (color).
- The following call sets it to the same value as the face_id : */
- cad_face_set_tag(fce, iface);
+ The following call sets it to the same value as the Geom module ID : */
+ const int faceTag = meshDS->ShapeToIndex(f);
+ cad_face_set_tag(fce, faceTag);
/* Set face orientation (optional if you want a well oriented output mesh)*/
if(f.Orientation() != TopAbs_FORWARD)
else
cad_face_set_orientation(fce, CAD_ORIENTATION_FORWARD);
- if (HasSizeMapOnFace && !use_precad)
+ if (HasSizeMapOnFace /*&& !use_precad*/) //22903: use_precad seems not to interfere
{
// -----------------
// Classic size map
if (FaceId2SizeMap.find(faceKey)!=FaceId2SizeMap.end()) {
- MESSAGE("A size map is defined on face :"<<faceKey)
- theSizeMapStr = FaceId2SizeMap[faceKey];
+ MESSAGE("A size map is defined on face :"<<faceKey);
+ theSizeMapStr = FaceId2SizeMap[faceKey];
// check if function ends with "return"
if (theSizeMapStr.find(bad_end) == (theSizeMapStr.size()-bad_end.size()-1))
continue;
ic = emap.Add(e);
// std::string aFileName = "fmap_edge_";
-// aFileName.append(to_string(ic));
+// aFileName.append(val_to_string(ic));
// aFileName.append(".brep");
// BRepTools::Write(e,aFileName.c_str());
}
// std::string aFileName = "fmap_vertex_";
-// aFileName.append(to_string(*ip));
+// aFileName.append(val_to_string(*ip));
// aFileName.append(".brep");
// BRepTools::Write(v,aFileName.c_str());
GMFFileName = _hypothesis->GetGMFFile();
if (GMFFileName != "") {
// bool GMFFileMode = _hypothesis->GetGMFFileMode();
- bool asciiFound = (GMFFileName.find(".mesh",GMFFileName.length()-5) != std::string::npos);
+ bool asciiFound = (GMFFileName.find(".mesh", GMFFileName.length()-5) != std::string::npos);
bool binaryFound = (GMFFileName.find(".meshb",GMFFileName.length()-6) != std::string::npos);
if (!asciiFound && !binaryFound)
GMFFileName.append(".mesh");
mesh_write_mesh(msh, GMFFileName.c_str());
+ //cout << GMFFileName << " created = " << SMESH_File(GMFFileName).exists() << endl;
}
/* retrieve mesh data (see meshgems/mesh.h) */
mesh_get_triangle_extra_vertices(msh, it, &type, evtri);
mesh_get_triangle_tag(msh, it, &tag);
if (tags[vtx[0]]) {
- meshDS->SetNodeOnFace(nodes[vtx[0]], TopoDS::Face(fmap(tag)));
+ meshDS->SetNodeOnFace(nodes[vtx[0]], tag);
tags[vtx[0]] = false;
};
if (tags[vtx[1]]) {
- meshDS->SetNodeOnFace(nodes[vtx[1]], TopoDS::Face(fmap(tag)));
+ meshDS->SetNodeOnFace(nodes[vtx[1]], tag);
tags[vtx[1]] = false;
};
if (tags[vtx[2]]) {
- meshDS->SetNodeOnFace(nodes[vtx[2]], TopoDS::Face(fmap(tag)));
+ meshDS->SetNodeOnFace(nodes[vtx[2]], tag);
tags[vtx[2]] = false;
};
if (type == MESHGEMS_MESH_ELEMENT_TYPE_TRIA6) {
// QUADRATIC TRIANGLE
if (tags[evtri[0]]) {
- meshDS->SetNodeOnFace(nodes[evtri[0]], TopoDS::Face(fmap(tag)));
+ meshDS->SetNodeOnFace(nodes[evtri[0]], tag);
tags[evtri[0]] = false;
}
if (tags[evtri[1]]) {
- meshDS->SetNodeOnFace(nodes[evtri[1]], TopoDS::Face(fmap(tag)));
+ meshDS->SetNodeOnFace(nodes[evtri[1]], tag);
tags[evtri[1]] = false;
}
if (tags[evtri[2]]) {
- meshDS->SetNodeOnFace(nodes[evtri[2]], TopoDS::Face(fmap(tag)));
+ meshDS->SetNodeOnFace(nodes[evtri[2]], tag);
tags[evtri[2]] = false;
}
tri = meshDS->AddFace(nodes[vtx[0]], nodes[vtx[1]], nodes[vtx[2]],
else {
tri = helper.AddFace(nodes[vtx[0]], nodes[vtx[1]], nodes[vtx[2]]);
}
- meshDS->SetMeshElementOnShape(tri, TopoDS::Face(fmap(tag)));
+ meshDS->SetMeshElementOnShape(tri, tag);
}
/* enumerate quadrangles */
mesh_get_quadrangle_extra_vertices(msh, it, &type, evquad);
mesh_get_quadrangle_tag(msh, it, &tag);
if (tags[vtx[0]]) {
- meshDS->SetNodeOnFace(nodes[vtx[0]], TopoDS::Face(fmap(tag)));
+ meshDS->SetNodeOnFace(nodes[vtx[0]], tag);
tags[vtx[0]] = false;
};
if (tags[vtx[1]]) {
- meshDS->SetNodeOnFace(nodes[vtx[1]], TopoDS::Face(fmap(tag)));
+ meshDS->SetNodeOnFace(nodes[vtx[1]], tag);
tags[vtx[1]] = false;
};
if (tags[vtx[2]]) {
- meshDS->SetNodeOnFace(nodes[vtx[2]], TopoDS::Face(fmap(tag)));
+ meshDS->SetNodeOnFace(nodes[vtx[2]], tag);
tags[vtx[2]] = false;
};
if (tags[vtx[3]]) {
- meshDS->SetNodeOnFace(nodes[vtx[3]], TopoDS::Face(fmap(tag)));
+ meshDS->SetNodeOnFace(nodes[vtx[3]], tag);
tags[vtx[3]] = false;
};
if (type == MESHGEMS_MESH_ELEMENT_TYPE_QUAD9) {
// QUADRATIC QUADRANGLE
std::cout << "This is a quadratic quadrangle" << std::endl;
if (tags[evquad[0]]) {
- meshDS->SetNodeOnFace(nodes[evquad[0]], TopoDS::Face(fmap(tag)));
+ meshDS->SetNodeOnFace(nodes[evquad[0]], tag);
tags[evquad[0]] = false;
}
if (tags[evquad[1]]) {
- meshDS->SetNodeOnFace(nodes[evquad[1]], TopoDS::Face(fmap(tag)));
+ meshDS->SetNodeOnFace(nodes[evquad[1]], tag);
tags[evquad[1]] = false;
}
if (tags[evquad[2]]) {
- meshDS->SetNodeOnFace(nodes[evquad[2]], TopoDS::Face(fmap(tag)));
+ meshDS->SetNodeOnFace(nodes[evquad[2]], tag);
tags[evquad[2]] = false;
}
if (tags[evquad[3]]) {
- meshDS->SetNodeOnFace(nodes[evquad[3]], TopoDS::Face(fmap(tag)));
+ meshDS->SetNodeOnFace(nodes[evquad[3]], tag);
tags[evquad[3]] = false;
}
if (tags[evquad[4]]) {
- meshDS->SetNodeOnFace(nodes[evquad[4]], TopoDS::Face(fmap(tag)));
+ meshDS->SetNodeOnFace(nodes[evquad[4]], tag);
tags[evquad[4]] = false;
}
quad = meshDS->AddFace(nodes[vtx[0]], nodes[vtx[1]], nodes[vtx[2]], nodes[vtx[3]],
else {
quad = helper.AddFace(nodes[vtx[0]], nodes[vtx[1]], nodes[vtx[2]], nodes[vtx[3]]);
}
- meshDS->SetMeshElementOnShape(quad, TopoDS::Face(fmap(tag)));
+ meshDS->SetMeshElementOnShape(quad, tag);
}
/* release the mesh object, the rest is released by cleaner */
_comment = "No mesh elements assigned to a face";
}
bool badFaceFound = false;
- for ( int i = 1; i <= fmap.Extent(); ++i )
+ for (TopExp_Explorer face_iter(aShape,TopAbs_FACE);face_iter.More();face_iter.Next())
{
- SMESH_subMesh* sm = aMesh.GetSubMesh( fmap(i) );
+ TopoDS_Face f = TopoDS::Face(face_iter.Current());
+ SMESH_subMesh* sm = aMesh.GetSubMesh( f );
if ( !sm->GetSubMeshDS() || sm->GetSubMeshDS()->NbElements() == 0 )
{
sm->GetComputeError().reset( new SMESH_ComputeError( err, _comment, this ));
