1 // SMESH SMESH : implementaion of SMESH idl descriptions
3 // Copyright (C) 2003 OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
4 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
6 // This library is free software; you can redistribute it and/or
7 // modify it under the terms of the GNU Lesser General Public
8 // License as published by the Free Software Foundation; either
9 // version 2.1 of the License.
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 // Lesser General Public License for more details.
16 // You should have received a copy of the GNU Lesser General Public
17 // License along with this library; if not, write to the Free Software
18 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
20 // See http://www.opencascade.org/SALOME/ or email : webmaster.salome@opencascade.org
24 // File : SMESH_Gen.cxx
25 // Author : Paul RASCLE, EDF
29 #include "SMESH_Gen.hxx"
30 #include "SMESH_subMesh.hxx"
31 #include "SMDS_MeshElement.hxx"
32 #include "SMDS_MeshNode.hxx"
33 #include "SMESHDriver.h"
36 #include <BRep_Tool.hxx>
38 #include "utilities.h"
41 //=============================================================================
43 * default constructor:
45 //=============================================================================
47 SMESH_Gen::SMESH_Gen()
49 MESSAGE("SMESH_Gen::SMESH_Gen");
51 _hypothesisFactory.SetGen(this);
54 //=============================================================================
58 //=============================================================================
60 SMESH_Gen::~SMESH_Gen()
62 MESSAGE("SMESH_Gen::~SMESH_Gen");
65 //=============================================================================
69 //=============================================================================
71 SMESH_Hypothesis *SMESH_Gen::CreateHypothesis(const char *anHyp, int studyId)
72 throw(SALOME_Exception)
75 MESSAGE("CreateHypothesis("<<anHyp<<","<<studyId<<")");
76 // Get studyContext, create it if it does'nt exist, with a SMESHDS_Document
78 StudyContextStruct *myStudyContext = GetStudyContext(studyId);
80 // create a new hypothesis object, store its ref. in studyContext
82 SMESH_Hypothesis *myHypothesis = _hypothesisFactory.Create(anHyp, studyId);
83 int hypId = myHypothesis->GetID();
84 myStudyContext->mapHypothesis[hypId] = myHypothesis;
88 // store hypothesis in SMESHDS document
90 myStudyContext->myDocument->AddHypothesis(myHypothesis);
94 //=============================================================================
98 //=============================================================================
100 SMESH_Mesh *SMESH_Gen::Init(int studyId, const TopoDS_Shape & aShape, int meshID)
101 throw(SALOME_Exception)
103 MESSAGE("SMESH_Gen::Init");
104 // if (aShape.ShapeType() == TopAbs_COMPOUND)
106 // INFOS("Mesh Compound not yet implemented!");
107 // throw(SALOME_Exception(LOCALIZED("Mesh Compound not yet implemented!")));
110 // Get studyContext, create it if it does'nt exist, with a SMESHDS_Document
112 StudyContextStruct *myStudyContext = GetStudyContext(studyId);
114 // create a new SMESH_mesh object
118 else if(_localId<=meshID)
121 SMESH_Mesh *mesh = new SMESH_Mesh(meshID,
124 myStudyContext->myDocument);
125 myStudyContext->mapMesh[meshID] = mesh;
127 // associate a TopoDS_Shape to the mesh
129 mesh->ShapeToMesh(aShape);
133 //=============================================================================
135 * @TODO Doing a full update after computation is not optimal when doing a local
138 //=============================================================================
140 bool SMESH_Gen::Compute(SMESH_Mesh & aMesh, const TopoDS_Shape & aShape)
141 throw(SALOME_Exception)
143 MESSAGE("SMESH_Gen::Compute");
144 // bool isDone = false;
146 Algo : s'appuie ou non sur une geometrie
148 Vertex : rien �faire (range le point)
149 Edge, Wire, collection d'edge et wire : 1D
150 Face, Shell, collection de Face et Shells : 2D
151 Solid, Collection de Solid : 3D
153 // *** corriger commentaires
154 // check hypothesis associated to the mesh :
155 // - only one algo : type compatible with the type of the shape
156 // - hypothesis = compatible with algo
157 // - check if hypothesis are applicable to this algo
158 // - check contradictions within hypothesis
159 // (test if enough hypothesis is done further)
163 SMESH_subMesh *sm = aMesh.GetSubMesh(aShape);
165 SMESH_subMesh *smToCompute = sm->GetFirstToCompute();
168 TopoDS_Shape subShape = smToCompute->GetSubShape();
169 int dim = GetShapeDim(subShape);
173 bool ret1 = smToCompute->ComputeStateEngine(SMESH_subMesh::COMPUTE);
179 ASSERT(smToCompute->_vertexSet == false);
180 TopoDS_Vertex V1 = TopoDS::Vertex(subShape);
181 gp_Pnt P1 = BRep_Tool::Pnt(V1);
182 SMESHDS_Mesh * meshDS = aMesh.GetMeshDS();
183 //MESSAGE("point "<<nodeId<<" "<<P1.X()<<" "<<P1.Y()<<" "<<P1.Z());
184 SMDS_MeshNode * node = meshDS->AddNode(P1.X(), P1.Y(), P1.Z());
185 meshDS->SetNodeOnVertex(node, V1);
186 smToCompute->GetSubMeshDS();
187 smToCompute->_vertexSet = true;
188 smToCompute->ComputeStateEngine(SMESH_subMesh::COMPUTE);
190 smToCompute = sm->GetFirstToCompute();
193 aMesh.GetMeshDS()->logFullUpdate();
198 //=============================================================================
202 //=============================================================================
204 SMESH_Algo *SMESH_Gen::GetAlgo(SMESH_Mesh & aMesh, const TopoDS_Shape & aShape)
206 //MESSAGE("SMESH_Gen::GetAlgo");
208 const SMESHDS_Hypothesis *theHyp = NULL;
209 SMESH_Algo *algo = NULL;
210 const SMESHDS_Mesh * meshDS = aMesh.GetMeshDS();
215 // try shape first, then main shape
217 TopoDS_Shape mainShape = meshDS->ShapeToMesh();
218 const TopoDS_Shape *shapeToTry[2] = { &aShape, &mainShape };
220 for (int iShape = 0; iShape < 2; iShape++)
222 TopoDS_Shape tryShape = (*shapeToTry[iShape]);
224 const list<const SMESHDS_Hypothesis*>& listHyp =
225 meshDS->GetHypothesis(tryShape);
226 list<const SMESHDS_Hypothesis*>::const_iterator it=listHyp.begin();
229 int shapeDim = GetShapeDim(aShape);
230 int typeOfShape = aShape.ShapeType();
232 while (it!=listHyp.end())
234 const SMESHDS_Hypothesis *anHyp = *it;
235 hypType = anHyp->GetType();
237 if (hypType > SMESHDS_Hypothesis::PARAM_ALGO)
241 case SMESHDS_Hypothesis::ALGO_1D:
244 case SMESHDS_Hypothesis::ALGO_2D:
247 case SMESHDS_Hypothesis::ALGO_3D:
256 //SCRUTE(typeOfShape);
257 if (shapeDim == algoDim) // count only algos of shape dim.
258 { // discard algos for subshapes
259 hypId = anHyp->GetID(); // (of lower dim.)
260 ASSERT(_mapAlgo.find(hypId) != _mapAlgo.end());
261 SMESH_Algo *anAlgo = _mapAlgo[hypId];
262 //SCRUTE(anAlgo->GetShapeType());
263 //if (anAlgo->GetShapeType() == typeOfShape)
264 if ((anAlgo->GetShapeType()) & (1 << typeOfShape))
265 { // only specific TopoDS_Shape
271 if (nb_algo > 1) return NULL; // more than one algo
274 if (nb_algo == 1) // one algo found : OK
275 break; // do not try a parent shape
279 return NULL; // no algo found
281 hypType = theHyp->GetType();
282 hypId = theHyp->GetID();
284 ASSERT(_mapAlgo.find(hypId) != _mapAlgo.end());
285 algo = _mapAlgo[hypId];
286 //MESSAGE("Algo found " << algo->GetName() << " Id " << hypId);
290 //=============================================================================
294 //=============================================================================
296 StudyContextStruct *SMESH_Gen::GetStudyContext(int studyId)
298 // Get studyContext, create it if it does'nt exist, with a SMESHDS_Document
300 if (_mapStudyContext.find(studyId) == _mapStudyContext.end())
302 _mapStudyContext[studyId] = new StudyContextStruct;
303 _mapStudyContext[studyId]->myDocument = new SMESHDS_Document(studyId);
305 StudyContextStruct *myStudyContext = _mapStudyContext[studyId];
306 // ASSERT(_mapStudyContext.find(studyId) != _mapStudyContext.end());
307 return myStudyContext;
310 //=============================================================================
314 //=============================================================================
316 void SMESH_Gen::Save(int studyId, const char *aUrlOfFile)
320 //=============================================================================
324 //=============================================================================
326 void SMESH_Gen::Load(int studyId, const char *aUrlOfFile)
330 //=============================================================================
334 //=============================================================================
336 void SMESH_Gen::Close(int studyId)
340 //=============================================================================
344 //=============================================================================
346 const char *SMESH_Gen::ComponentDataType()
350 //=============================================================================
354 //=============================================================================
356 const char *SMESH_Gen::IORToLocalPersistentID(const char *IORString,
361 //=============================================================================
365 //=============================================================================
367 const char *SMESH_Gen::LocalPersistentIDToIOR(const char *aLocalPersistentID)
371 //=============================================================================
375 //=============================================================================
377 int SMESH_Gen::GetShapeDim(const TopoDS_Shape & aShape)
379 int shapeDim = -1; // Shape dimension: 0D, 1D, 2D, 3D
380 int type = aShape.ShapeType();
383 // case TopAbs_COMPOUND:
387 case TopAbs_COMPOUND:
388 case TopAbs_COMPSOLID:
395 // case TopAbs_SHELL:
418 * Import a mesh from a file
419 * @param fileName file name to be imported
420 * @param fileType Currently it could be either "DAT", "UNV" or "MED".
423 SMESH_Mesh * SMESH_Gen::Import(int studyId, const char * fileName,
424 const char * fileType)
426 MESSAGE("SMESH_Gen::Import("<<studyId<<","<<fileName<<","<<fileType<<")");
428 // Get studyContext, create it if it does'nt exist, with a SMESHDS_Document
429 StudyContextStruct *myStudyContext = GetStudyContext(studyId);
431 // will be used with document
432 /*Document_Reader * reader = SMESHDriver::GetDocumentReader(string(fileType));
433 reader->SetDocument(myStudyContext->myDocument);
434 reader->SetFile(string(fileName));
436 // currently we only read one mesh from a file (limitation on MED files).
438 // create a new SMESH_mesh object
439 SMESH_Mesh *mesh = new SMESH_Mesh(_localId++, studyId, this,
440 myStudyContext->myDocument);
441 myStudyContext->mapMesh[_localId] = mesh;
443 Mesh_Reader * reader = SMESHDriver::GetMeshReader(string(fileType));
444 reader->SetMesh(mesh->GetMeshDS());
445 reader->SetFile(string(fileName));
448 mesh->GetMeshDS()->logFullUpdate();