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
|
/*=========================================================================
Program: Visualization Toolkit
Module: TestMappedGridDeepCopy
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.
=========================================================================*/
/*----------------------------------------------------------------------------
This test was written by Menno Deij - van Rijswijk (MARIN).
----------------------------------------------------------------------------*/
#include "vtkCell.h" // for cell types
#include "vtkCellIterator.h"
#include "vtkDataArray.h"
#include "vtkDebugLeaks.h"
#include "vtkIdList.h"
#include "vtkIdTypeArray.h"
#include "vtkInformation.h"
#include "vtkMappedUnstructuredGridGenerator.h"
#include "vtkNew.h"
#include "vtkPoints.h"
#include "vtkTestUtilities.h"
#include "vtkUnstructuredGrid.h"
#include "vtkXMLUnstructuredGridReader.h"
#include "vtkXMLUnstructuredGridWriter.h"
#include <algorithm>
#include <fstream>
int TestMappedGridDeepCopy(int vtkNotUsed(argc), char*[] vtkNotUsed(argv))
{
vtkUnstructuredGrid* original;
vtkMappedUnstructuredGridGenerator::GenerateUnstructuredGrid(&original);
vtkUnstructuredGridBase* mg;
vtkMappedUnstructuredGridGenerator::GenerateMappedUnstructuredGrid(&mg);
// copy the mapped grid into a normal unstructured grid.
// copying will proceed via the super class
// vtkUnstructuredGridBase::DeepCopy function
// implementation that uses a cell iterator. This will
// invoke to InsertNextCell function with face list
// for *all* cells (even if they are not VTK_POLYHEDRON).
// In the old implementation this gave copy errors. The fix
// proposed together with this test addresses that issue.
vtkNew<vtkUnstructuredGrid> copy;
copy->Allocate(mg->GetNumberOfCells());
copy->DeepCopy(mg);
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->GetFaces();
vtkIdList* cFaces = cIt->GetFaces();
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();
mg->Delete();
return EXIT_SUCCESS;
}
|
