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
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
|
/*=========================================================================
Program: Visualization Toolkit
Module: TestPolyhedron3.cxx
Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
All rights reserved.
See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notice for more information.
=========================================================================*/
#include "vtkPlane.h"
#include "vtkPolyData.h"
#include "vtkPolyhedron.h"
#include "vtkUnstructuredGrid.h"
#include "vtkClipDataSet.h"
#include "vtkNew.h"
#include "vtkTestUtilities.h"
#include "vtkUnstructuredGridReader.h"
#include "vtkXMLUnstructuredGridWriter.h"
const char inputDataStream[] =
"# vtk DataFile Version 3.0\n"
"vtk output\n"
"ASCII\n"
"DATASET UNSTRUCTURED_GRID\n"
"POINTS 8 float\n"
"1337.72 1586.34 914.4 1337.72 1586.34 1371.6 1261.68 1606.71 914.4 \n"
"1261.68 1606.71 1371.6 1337.72 1484.47 914.4 1337.72 1484.47 1371.6 \n"
"1261.68 1464.1 914.4 1261.68 1464.1 1371.6 \n"
"CELLS 1 32\n"
"31 6 4 4 6 2 0 4 1 3 7 5 4 0 2 3 1 4 2 6 7 3 4 6 4 5 7 4 4 0 1 5 \n"
"CELL_TYPES 1\n"
"42\n";
// Test of contour/clip of vtkPolyhedron. uses input from
// https://gitlab.kitware.com/vtk/vtk/-/issues/15026
int TestPolyhedron3(int argc, char* argv[])
{
(void)argc;
(void)argv;
vtkNew<vtkUnstructuredGridReader> reader;
reader->SetInputString(inputDataStream);
reader->ReadFromInputStringOn();
vtkNew<vtkPlane> plane;
plane->SetNormal(0.847934330264784, 0.530022019598814, -0.00916680417631942);
plane->SetOrigin(1254.0760499239, 1489.93486006017, 1143.9780493697);
vtkNew<vtkClipDataSet> clip;
clip->SetInputConnection(reader->GetOutputPort());
clip->SetClipFunction(plane);
clip->Update();
vtkUnstructuredGrid* result = clip->GetOutput(0);
if (!result)
return 1;
if (result->GetNumberOfCells() != 1)
{
std::cout << "Expected 1 but found " << result->GetNumberOfCells()
<< " cells in intersected polyhedron" << std::endl;
return EXIT_FAILURE;
}
vtkCell* clipped = result->GetCell(0);
if (!clipped)
return 1;
if (clipped->GetNumberOfFaces() != 7)
{
std::cout << "Expected 7 but found " << clipped->GetNumberOfFaces()
<< " faces on in intersected polyhedron" << std::endl;
return EXIT_FAILURE;
}
return EXIT_SUCCESS;
}
|
