0024428: Implementation of LGPL license
[occt.git] / src / StlMesh / StlMesh.cxx
1 // Created on: 1996-06-21
2 // Created by: Bruno TACCI
3 // Copyright (c) 1996-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
5 //
6 // This file is part of Open CASCADE Technology software library.
7 //
8 // This library is free software; you can redistribute it and / or modify it
9 // under the terms of the GNU Lesser General Public version 2.1 as published
10 // by the Free Software Foundation, with special exception defined in the file
11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
12 // distribution for complete text of the license and disclaimer of any warranty.
13 //
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
16
17 #include <StlMesh.ixx>
18 #include <StlMesh_SequenceOfMeshTriangle.hxx>
19 #include <StlMesh_MeshTriangle.hxx>
20 #include <TColgp_SequenceOfXYZ.hxx>
21
22 //=======================================================================
23 //function : Merge
24 //purpose  : 
25 //=======================================================================
26
27 Handle(StlMesh_Mesh) StlMesh::Merge(const Handle(StlMesh_Mesh)& mesh1, const Handle(StlMesh_Mesh)& mesh2)
28 {
29   Handle(StlMesh_Mesh) mergeMesh = new StlMesh_Mesh;
30   StlMesh_SequenceOfMeshTriangle aSeqOfTriangle;
31   TColgp_SequenceOfXYZ aSeqOfVertex;
32   Standard_Real xn,yn,zn;
33   Standard_Integer v1,v2,v3;
34
35   // Chargement de mesh1 dans mergeMesh
36   // Boucle sur les domaines puis sur les triangles
37
38   Standard_Integer idom;
39   for (idom = 1; idom <= mesh1->NbDomains(); idom++) {
40     aSeqOfTriangle = mesh1->Triangles(idom);
41     aSeqOfVertex = mesh1->Vertices(idom);
42     mergeMesh->AddDomain(mesh1->Deflection(idom));
43     
44     for (Standard_Integer itri = 1; itri <= mesh1->NbTriangles(idom); itri++) {
45       const Handle(StlMesh_MeshTriangle) aTrian = aSeqOfTriangle.Value(itri);
46       aTrian->GetVertexAndOrientation(v1,v2,v3,xn,yn,zn);
47       mergeMesh->AddTriangle(v1,v2,v3,xn,yn,zn);
48     }
49     
50     for (Standard_Integer iver = 1; iver <= mesh1->NbVertices(idom); iver++) {
51       mergeMesh->AddVertex(aSeqOfVertex.Value(iver).X(),
52                            aSeqOfVertex.Value(iver).Y(),
53                            aSeqOfVertex.Value(iver).Z());
54     }
55     
56   }
57   // Idem avec mesh2
58   
59   for (idom = 1; idom <= mesh2->NbDomains(); idom++) {
60     aSeqOfTriangle = mesh2->Triangles(idom);
61     aSeqOfVertex = mesh2->Vertices(idom);
62     mergeMesh->AddDomain(mesh2->Deflection(idom));
63     
64     for (Standard_Integer itri = 1; itri <= mesh2->NbTriangles(idom); itri++) {
65       const Handle(StlMesh_MeshTriangle) aTrian = aSeqOfTriangle.Value(itri);
66       aTrian->GetVertexAndOrientation(v1,v2,v3,xn,yn,zn);
67       mergeMesh->AddTriangle(v1,v2,v3,xn,yn,zn);
68     }
69     
70     for (Standard_Integer iver = 1; iver <= mesh2->NbVertices(idom); iver++) {
71       mergeMesh->AddVertex(aSeqOfVertex.Value(iver).X(),
72                            aSeqOfVertex.Value(iver).Y(),
73                            aSeqOfVertex.Value(iver).Z());
74     }
75   }
76   return mergeMesh;
77 }