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