Integration of OCCT 6.5.0 from SVN
[occt.git] / src / StlMesh / StlMesh.cxx
1 // File:        StlMesh.cxx
2 // Created:     Fri Jun 21 15:48:55 1996
3 // Author:      Bruno TACCI
4 //              <bti@sgi44>
5
6 #include <StlMesh.ixx>
7 #include <StlMesh_SequenceOfMeshTriangle.hxx>
8 #include <StlMesh_MeshTriangle.hxx>
9 #include <TColgp_SequenceOfXYZ.hxx>
10
11 //=======================================================================
12 //function : Merge
13 //purpose  : 
14 //=======================================================================
15
16 Handle(StlMesh_Mesh) StlMesh::Merge(const Handle(StlMesh_Mesh)& mesh1, const Handle(StlMesh_Mesh)& mesh2)
17 {
18   Handle(StlMesh_Mesh) mergeMesh = new StlMesh_Mesh;
19   StlMesh_SequenceOfMeshTriangle aSeqOfTriangle;
20   TColgp_SequenceOfXYZ aSeqOfVertex;
21   Standard_Real xn,yn,zn;
22   Standard_Integer v1,v2,v3;
23
24   // Chargement de mesh1 dans mergeMesh
25   // Boucle sur les domaines puis sur les triangles
26
27   Standard_Integer idom;
28   for (idom = 1; idom <= mesh1->NbDomains(); idom++) {
29     aSeqOfTriangle = mesh1->Triangles(idom);
30     aSeqOfVertex = mesh1->Vertices(idom);
31     mergeMesh->AddDomain(mesh1->Deflection(idom));
32     
33     for (Standard_Integer itri = 1; itri <= mesh1->NbTriangles(idom); itri++) {
34       const Handle(StlMesh_MeshTriangle) aTrian = aSeqOfTriangle.Value(itri);
35       aTrian->GetVertexAndOrientation(v1,v2,v3,xn,yn,zn);
36       mergeMesh->AddTriangle(v1,v2,v3,xn,yn,zn);
37     }
38     
39     for (Standard_Integer iver = 1; iver <= mesh1->NbVertices(idom); iver++) {
40       mergeMesh->AddVertex(aSeqOfVertex.Value(iver).X(),
41                            aSeqOfVertex.Value(iver).Y(),
42                            aSeqOfVertex.Value(iver).Z());
43     }
44     
45   }
46   // Idem avec mesh2
47   
48   for (idom = 1; idom <= mesh2->NbDomains(); idom++) {
49     aSeqOfTriangle = mesh2->Triangles(idom);
50     aSeqOfVertex = mesh2->Vertices(idom);
51     mergeMesh->AddDomain(mesh2->Deflection(idom));
52     
53     for (Standard_Integer itri = 1; itri <= mesh2->NbTriangles(idom); itri++) {
54       const Handle(StlMesh_MeshTriangle) aTrian = aSeqOfTriangle.Value(itri);
55       aTrian->GetVertexAndOrientation(v1,v2,v3,xn,yn,zn);
56       mergeMesh->AddTriangle(v1,v2,v3,xn,yn,zn);
57     }
58     
59     for (Standard_Integer iver = 1; iver <= mesh2->NbVertices(idom); iver++) {
60       mergeMesh->AddVertex(aSeqOfVertex.Value(iver).X(),
61                            aSeqOfVertex.Value(iver).Y(),
62                            aSeqOfVertex.Value(iver).Z());
63     }
64   }
65   return mergeMesh;
66 }