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
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
|
// SPDX-FileCopyrightText: Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
// SPDX-License-Identifier: BSD-3-Clause
#include "vtkSetGet.h"
#include "vtkMathUtilities.h"
#include "vtkVector.h"
#include "vtkVectorOperators.h"
//------------------------------------------------------------------------------
int TestVector(int, char*[])
{
// Store up any errors, return non-zero if something fails.
int retVal = 0;
// Test out the general vector data types, give nice API and great memory use
vtkVector2i vec2i;
int arr2i[2] = { 0, 0 };
vec2i.Set(arr2i[0], arr2i[1]);
if (sizeof(vec2i) != sizeof(arr2i))
{
// The two should be the same size and memory layout - error out if not
cerr << "vtkVector2i should be the same size as int[2]." << endl
<< "sizeof(vec2i) = " << sizeof(vec2i) << endl
<< "sizeof(int[2]) = " << sizeof(arr2i) << endl;
++retVal;
}
vtkVector<float, 3> vector3f;
if (vector3f.GetSize() != 3)
{
cerr << "Incorrect size of vector3f, should be 3, but is " << vector3f.GetSize() << endl;
++retVal;
}
// Test out vtkVector3i and ensure the various access methods are the same
vtkVector3i vec3i(0, 6, 9);
if (vec3i.GetX() != vec3i[0] || vec3i.GetX() != 0)
{
cerr << "vec3i.GetX() should equal vec3i.GetData()[0] which should equal 0."
<< "\nvec3i.GetX() = " << vec3i.GetX() << endl
<< "vec3i[0] = " << vec3i[0] << endl;
++retVal;
}
if (vec3i.GetY() != vec3i[1] || vec3i.GetY() != 6)
{
cerr << "vec3i.GetY() should equal vec3i.GetData()[1] which should equal 6."
<< "\nvec3i.GetY() = " << vec3i.GetY() << endl
<< "vec3i[1] = " << vec3i[1] << endl;
++retVal;
}
if (vec3i.GetZ() != vec3i[2] || vec3i.GetZ() != 9)
{
cerr << "vec3i.GetZ() should equal vec3i.GetData()[2] which should equal 9."
<< "\nvec3i.GetZ() = " << vec3i.GetZ() << endl
<< "vec3i[2] = " << vec3i[2] << endl;
++retVal;
}
// Assign the data to an int array and ensure the two ways of referencing are
// the same.
int* intPtr = vec3i.GetData();
for (int i = 0; i < 3; ++i)
{
if (vec3i[i] != intPtr[i] || vec3i(i) != vec3i[i])
{
cerr << "Error: vec3i[i] != intPtr[i]" << endl
<< "vec3i[i] = " << vec3i[i] << endl
<< "intPtr[i] = " << intPtr[i] << endl;
++retVal;
}
}
// Test out casting...
vtkVector<float, 3> castVec = vec3i.Cast<float>();
vtkVector3d castVecd(castVec.Cast<double>().GetData());
if (castVecd[0] < -0.0000001 || castVecd[0] > 0.0000001)
{
// Then the number did not make it through within reasonable precision.
cerr << "Error: castVecd value incorrect. Should be ~0.0 for component 1."
<< "\ncastVecd[0] = " << castVecd[0] << endl;
++retVal;
}
cout << "castVecd[0] = " << castVecd[0] << endl;
// Test the normalize and normalized functions.
vtkVector3d normy(1, 2, 3);
vtkVector3d normed = normy.Normalized();
double dotted = normy.Dot(normed);
if (!vtkMathUtilities::FuzzyCompare(dotted, 3.74166, 0.0001))
{
cerr << "The dot product of " << normy << " and " << normed << " was " << dotted
<< ", expected 3.74166." << endl;
++retVal;
}
if (!normed.Compare(vtkVector3d(0.267261, 0.534522, 0.801784), 0.0001))
{
cerr << "Error vtkVector3d::Normalized() failed: " << normed << endl;
++retVal;
}
normy.Normalize();
if (!normy.Compare(normed, 0.0001))
{
cerr << "Error vtkVector3d::Normalize() failed: " << normy << endl;
}
if (!vtkMathUtilities::FuzzyCompare(normy.Norm(), 1.0, 0.0001))
{
cerr << "Normalized length should always be ~= 1.0, value is " << normy.Norm() << endl;
++retVal;
}
if (!vtkMathUtilities::FuzzyCompare(dotted, 3.74166, 0.0001))
{
cerr << "The dot product of " << normy.Norm() << endl;
++retVal;
}
if (!vtkMathUtilities::FuzzyCompare(normy.Dot(normed), 1.0, 0.0001))
{
cerr << "The dot product of " << normy << " and " << normed << " was " << normy.Dot(normed)
<< ", expected 1.0." << endl;
++retVal;
}
// Some cross product stuff now...
if (!normy.Cross(normed).Compare(vtkVector3d(0, 0, 0), 0.0001))
{
cerr << normy << " cross " << normed << " expected to be 0, got " << normy.Cross(normed)
<< endl;
++retVal;
}
if (!normy.Cross(vtkVector3d(0, 1, 0)).Compare(vtkVector3d(-0.801784, 0, 0.267261), 0.0001))
{
cerr << normy << " cross (0, 1, 0) expected to be (-0.801784, 0, 0.267261), got "
<< normy.Cross(normed) << endl;
++retVal;
}
return retVal;
}
|
