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
|
/*=========================================================================
Program: Visualization Toolkit
Module: TestGenericCell.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 "vtkGenericCell.h"
#include "vtkMath.h"
int TestGenericCell(int , char *[])
{
int rval = 0;
vtkGenericCell *cell = vtkGenericCell::New();
for(int i=0; i<VTK_NUMBER_OF_CELL_TYPES;++i)
{
cell->SetCellType( i );
if( cell->RequiresInitialization() )
{
cell->Initialize();
}
cell->Print( cout );
int numPts = cell->GetNumberOfPoints();
int numEdges = cell->GetNumberOfEdges();
int numFaces = cell->GetNumberOfFaces();
double center[3];
int a = cell->GetParametricCenter(center);
(void)a;
double *pcoords = cell->GetParametricCoords();
if( cell->GetCellType() != VTK_EMPTY_CELL
&& cell->GetCellType() != VTK_POLY_VERTEX // FIXME
&& cell->GetCellType() != VTK_POLY_LINE // FIXME
&& cell->GetCellType() != VTK_TRIANGLE_STRIP // FIXME
&& cell->GetCellType() != VTK_POLYGON // FIXME
&& cell->GetCellType() != VTK_CONVEX_POINT_SET // FIXME
&& cell->GetCellType() != VTK_POLYHEDRON) // FIXME
{
double m[3] = {0., 0., 0.};
// We add all the points since
// Those on the corner points indeed define the parametric center
// The dof node (center mid points) by definition have the same parametric center
// and taking into account the center point only add a 0 vector to the sum
// therefore we do not need to differenciate corner from the rest in this sum:
for(int j=0; j<numPts; ++j)
{
double *point = pcoords + 3*j;
m[0] += point[0];
m[1] += point[1];
m[2] += point[2];
}
m[0] /= numPts;
m[1] /= numPts;
m[2] /= numPts;
if( fabs( center[0] - m[0] ) > 1e-6
|| fabs( center[1] - m[1] ) > 1e-6
|| fabs( center[2] - m[2] ) > 1e-6)
{
cerr << "Cell: " << i << endl;
cerr << "Center: " << center[0] << "," << center[1] << "," << center[2] << endl;
cerr << "M : " << m[0] << "," << m[1] << "," << m[2] << endl;
++rval;
}
}
int p = cell->IsPrimaryCell();
(void)p;
int cellDim = cell->GetCellDimension();
(void)cellDim;
int l = cell->IsLinear();
(void)l;
for(int e=0; e<numEdges; ++e)
{
vtkCell *c = cell->GetEdge(e);
c->Print( cout );
}
for(int f=0; f<numFaces; ++f)
{
vtkCell *c = cell->GetFace(f);
c->Print( cout );
}
if( cell->GetCellType() != i && cell->GetCellType() != VTK_EMPTY_CELL )
{
++rval;
}
}
cell->Delete();
return rval;
}
|
