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
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
|
// SPDX-FileCopyrightText: Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
// SPDX-License-Identifier: BSD-3-Clause
/*----------------------------------------------------------------------------
This was taken from TestMappedGridDeepCopy and modified by Bengt Rosenberger.
----------------------------------------------------------------------------*/
#include "vtkCell.h" // for cell types
#include "vtkCellData.h"
#include "vtkCellIterator.h"
#include "vtkDataArray.h"
#include "vtkDebugLeaks.h"
#include "vtkDoubleArray.h"
#include "vtkIdList.h"
#include "vtkIdTypeArray.h"
#include "vtkInformation.h"
#include "vtkIntArray.h"
#include "vtkMappedUnstructuredGridGenerator.h"
#include "vtkNew.h"
#include "vtkPointData.h"
#include "vtkPoints.h"
#include "vtkTestUtilities.h"
#include "vtkUnstructuredGrid.h"
#include "vtkXMLUnstructuredGridReader.h"
#include "vtkXMLUnstructuredGridWriter.h"
#include <algorithm>
#include <fstream>
int TestMappedGridShallowCopy(int vtkNotUsed(argc), char*[] vtkNotUsed(argv))
{
vtkUnstructuredGridBase* original;
vtkMappedUnstructuredGridGenerator::GenerateMappedUnstructuredGrid(&original);
// This simulates the executive calling vtkUnstructuredGrid::Initialize
// when preparing an algorithm output
vtkUnstructuredGrid* copy = vtkUnstructuredGrid::New();
copy->Initialize();
// Make sure shallow copy succeeds after call to Initialize
copy->ShallowCopy(original);
// Compare number of points
if (copy->GetNumberOfPoints() != original->GetNumberOfPoints())
{
cerr << "Number of points do not match" << endl;
return EXIT_FAILURE;
}
// Compare number of cells
if (copy->GetNumberOfCells() != original->GetNumberOfCells())
{
cerr << "Number of cells do not match" << endl;
return EXIT_FAILURE;
}
// Comparison code below taken from TestMappedGridDeepCopy.cxx
vtkCellIterator* oIt = original->NewCellIterator();
vtkCellIterator* cIt = copy->NewCellIterator();
vtkNew<vtkGenericCell> orig, copied;
for (cIt->InitTraversal(), oIt->InitTraversal();
!cIt->IsDoneWithTraversal() && !oIt->IsDoneWithTraversal();
cIt->GoToNextCell(), oIt->GoToNextCell())
{
oIt->GetCell(orig.GetPointer());
cIt->GetCell(copied.GetPointer());
if (cIt->GetCellType() != oIt->GetCellType())
{
cerr << "Cell types do not match" << endl;
return EXIT_FAILURE;
}
if (cIt->GetCellType() == VTK_POLYHEDRON)
{
vtkIdList* oFaces = oIt->GetSerializedCellFaces();
vtkIdList* cFaces = cIt->GetSerializedCellFaces();
if (cFaces->GetNumberOfIds() != oFaces->GetNumberOfIds())
{
cerr << "Face id list length does not match" << endl;
cerr << "Original: ";
for (vtkIdType i = 0; i < oFaces->GetNumberOfIds(); ++i)
{
cerr << oFaces->GetId(i) << " ";
}
cerr << endl;
cerr << "Copied: ";
for (vtkIdType i = 0; i < cFaces->GetNumberOfIds(); ++i)
cerr << cFaces->GetId(i) << " ";
cerr << endl;
return EXIT_FAILURE;
}
for (vtkIdType i = 0; i < cFaces->GetNumberOfIds(); ++i)
{
vtkIdType c = cFaces->GetId(i);
vtkIdType o = oFaces->GetId(i);
if (c != o)
{
cerr << "Face id list content does not match at" << i << endl;
return EXIT_FAILURE;
}
}
}
}
oIt->Delete();
cIt->Delete();
original->Delete();
copy->Delete();
return EXIT_SUCCESS;
}
|
