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 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.
  13 //
  14 // Alternatively, this file may be used under the terms of Open CASCADE
  15 // commercial license or contractual agreement.
  16
  17
  18 #include <StlMesh.hxx>
  19 #include <StlMesh_Mesh.hxx>
  20 #include <StlMesh_MeshTriangle.hxx>
  21 #include <StlMesh_SequenceOfMeshTriangle.hxx>
  22 #include <TColgp_SequenceOfXYZ.hxx>
  23
  24 //=======================================================================
  25 //function : Merge
  26 //purpose  :
  27 //=======================================================================
  28 Handle(StlMesh_Mesh) StlMesh::Merge(const Handle(StlMesh_Mesh)& mesh1, const Handle(StlMesh_Mesh)& mesh2)
  29 {
  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;
  35
  36   // Chargement de mesh1 dans mergeMesh
  37   // Boucle sur les domaines puis sur les triangles
  38
  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));
  44
  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);
  49     }
  50
  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());
  55     }
  56
  57   }
  58   // Idem avec mesh2
  59
  60   for (idom = 1; idom <= mesh2->NbDomains(); idom++) {
  61     aSeqOfTriangle = mesh2->Triangles(idom);
  62     aSeqOfVertex = mesh2->Vertices(idom);
  63     mergeMesh->AddDomain(mesh2->Deflection(idom));
  64
  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);
  69     }
  70
  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());
  75     }
  76   }
  77   return mergeMesh;
  78 }