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
6 // This file is part of Open CASCADE Technology software library.
8 // This library is free software; you can redistribute it and/or modify it under
9 // the terms of the GNU Lesser General Public License 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.
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
18 #include <StlMesh.hxx>
19 #include <StlMesh_Mesh.hxx>
20 #include <StlMesh_MeshTriangle.hxx>
21 #include <StlMesh_SequenceOfMeshTriangle.hxx>
22 #include <TColgp_SequenceOfXYZ.hxx>
24 //=======================================================================
27 //=======================================================================
28 Handle(StlMesh_Mesh) StlMesh::Merge(const Handle(StlMesh_Mesh)& mesh1, const Handle(StlMesh_Mesh)& mesh2)
30 Handle(StlMesh_Mesh) mergeMesh = new StlMesh_Mesh;
31 StlMesh_SequenceOfMeshTriangle aSeqOfTriangle;
32 TColgp_SequenceOfXYZ aSeqOfVertex;
33 Standard_Real xn,yn,zn;
34 Standard_Integer v1,v2,v3;
36 // Chargement de mesh1 dans mergeMesh
37 // Boucle sur les domaines puis sur les triangles
39 Standard_Integer idom;
40 for (idom = 1; idom <= mesh1->NbDomains(); idom++) {
41 aSeqOfTriangle = mesh1->Triangles(idom);
42 aSeqOfVertex = mesh1->Vertices(idom);
43 mergeMesh->AddDomain(mesh1->Deflection(idom));
45 for (Standard_Integer itri = 1; itri <= mesh1->NbTriangles(idom); itri++) {
46 const Handle(StlMesh_MeshTriangle) aTrian = aSeqOfTriangle.Value(itri);
47 aTrian->GetVertexAndOrientation(v1,v2,v3,xn,yn,zn);
48 mergeMesh->AddTriangle(v1,v2,v3,xn,yn,zn);
51 for (Standard_Integer iver = 1; iver <= mesh1->NbVertices(idom); iver++) {
52 mergeMesh->AddVertex(aSeqOfVertex.Value(iver).X(),
53 aSeqOfVertex.Value(iver).Y(),
54 aSeqOfVertex.Value(iver).Z());
60 for (idom = 1; idom <= mesh2->NbDomains(); idom++) {
61 aSeqOfTriangle = mesh2->Triangles(idom);
62 aSeqOfVertex = mesh2->Vertices(idom);
63 mergeMesh->AddDomain(mesh2->Deflection(idom));
65 for (Standard_Integer itri = 1; itri <= mesh2->NbTriangles(idom); itri++) {
66 const Handle(StlMesh_MeshTriangle) aTrian = aSeqOfTriangle.Value(itri);
67 aTrian->GetVertexAndOrientation(v1,v2,v3,xn,yn,zn);
68 mergeMesh->AddTriangle(v1,v2,v3,xn,yn,zn);
71 for (Standard_Integer iver = 1; iver <= mesh2->NbVertices(idom); iver++) {
72 mergeMesh->AddVertex(aSeqOfVertex.Value(iver).X(),
73 aSeqOfVertex.Value(iver).Y(),
74 aSeqOfVertex.Value(iver).Z());