Medical Imaging Interaction Toolkit  2016.11.0
Medical Imaging Interaction Toolkit
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
mitkVectorTest.cpp
Go to the documentation of this file.
1 /*===================================================================
2 
3 The Medical Imaging Interaction Toolkit (MITK)
4 
5 Copyright (c) German Cancer Research Center,
6 Division of Medical and Biological Informatics.
7 All rights reserved.
8 
9 This software is distributed WITHOUT ANY WARRANTY; without
10 even the implied warranty of MERCHANTABILITY or FITNESS FOR
11 A PARTICULAR PURPOSE.
12 
13 See LICENSE.txt or http://www.mitk.org for details.
14 
15 ===================================================================*/
16 
17 #include <mitkNumericTypes.h>
18 #include <mitkTestingMacros.h>
19 
20 #include <vnl/vnl_vector_fixed.txx>
21 
22 int mitkVectorTest(int /*argc*/, char * /*argv*/ [])
23 {
24  MITK_TEST_BEGIN("mitkVector");
25  // test itk vector equality methods
26  itk::Vector<mitk::ScalarType, 3> itkVector_1;
27  itkVector_1[0] = 4.6;
28  itkVector_1[1] = 9.76543;
29  itkVector_1[2] = 746.09;
30 
31  itk::Vector<mitk::ScalarType, 3> itkVector_2;
32  itk::Vector<mitk::ScalarType, 3> itkVector_3;
33  for (int i = 0; i < 3; i++)
34  {
35  itkVector_2[i] = itkVector_1[i] - mitk::eps * 1.1;
36  itkVector_3[i] = itkVector_1[i] - mitk::eps * 0.9;
37  }
38 
39  MITK_TEST_CONDITION(mitk::Equal(itkVector_1, itkVector_1),
40  "Test vector equality using the same vector with mitk::eps");
41  MITK_TEST_CONDITION(
42  !mitk::Equal(itkVector_1, itkVector_2),
43  "Test vector equality using different vectors with an element-wise difference greater than mitk::eps");
44  MITK_TEST_CONDITION(mitk::Equal(itkVector_1, itkVector_2, mitk::eps * 1.2),
45  "Vectors are equal for higher epsilon tolerance ( 1.2 * mitk::eps )");
46  MITK_TEST_CONDITION(
47  mitk::Equal(itkVector_1, itkVector_3),
48  "Test vector equality using different vectors with an element-wise difference less than mitk::eps");
49  // test itk point equality methods
50  itk::Point<mitk::ScalarType, 3> itkPoint_1;
51  itk::Point<mitk::ScalarType, 3> itkPoint_2;
52  itk::Point<mitk::ScalarType, 3> itkPoint_3;
53  for (int i = 0; i < 3; i++)
54  {
55  itkPoint_1[i] = itkVector_1[i];
56  itkPoint_2[i] = itkVector_2[i];
57  itkPoint_3[i] = itkVector_3[i];
58  }
59  MITK_TEST_CONDITION(mitk::Equal(itkPoint_1, itkPoint_1), "Test point equality using the same point with mitk::eps");
60  MITK_TEST_CONDITION(
61  !mitk::Equal(itkPoint_1, itkPoint_2),
62  "Test point equality using different points with an element-wise difference greater than mitk::eps");
63  MITK_TEST_CONDITION(mitk::Equal(itkPoint_1, itkPoint_2, mitk::eps * 1.2),
64  "Points are equal for higher epsilon tolerance ( 1.2 * mitk::eps )");
65  MITK_TEST_CONDITION(mitk::Equal(itkPoint_1, itkPoint_3),
66  "Test point equality using different points with an element-wise difference less than mitk::eps");
67  // test mitk vnl vector equality methods
68  mitk::VnlVector mitk_vnl_vector_1(3);
69  mitk::VnlVector mitk_vnl_vector_2(3);
70  mitk::VnlVector mitk_vnl_vector_3(3);
71  for (int i = 0; i < 3; i++)
72  {
73  mitk_vnl_vector_1.put(i, itkVector_1[i]);
74  mitk_vnl_vector_2.put(i, itkVector_2[i]);
75  mitk_vnl_vector_3.put(i, itkVector_1[i]);
76  }
77 
78  MITK_TEST_CONDITION(mitk::Equal(mitk_vnl_vector_1, mitk_vnl_vector_1),
79  "Test mitk vnl vector equality using the same mitk vnl vector with mitk::eps");
80  MITK_TEST_CONDITION(!mitk::Equal(mitk_vnl_vector_1, mitk_vnl_vector_2),
81  "Test mitk vnl vector equality using different mitk vnl vectors with an element-wise difference "
82  "greater than mitk::eps");
83  MITK_TEST_CONDITION(mitk::Equal(mitk_vnl_vector_1, mitk_vnl_vector_2, mitk::eps * 1.2),
84  "Vnl vectors are equal for higher epsilon tolerance ( 1.2 * mitk::eps )");
85  MITK_TEST_CONDITION(mitk::Equal(mitk_vnl_vector_1, mitk_vnl_vector_3),
86  "Test mitk vnl vector equality using "
87  "different mitk vnl vectors with an "
88  "element-wise difference less than mitk::eps");
89 
90  // test vnl_vector equality method
91  typedef mitk::ScalarType VnlValueType;
92  vnl_vector_fixed<VnlValueType, 7> vnlVector_1;
93  vnlVector_1[3] = 56.98;
94  vnlVector_1[4] = 22.32;
95  vnlVector_1[5] = 1.00;
96  vnlVector_1[6] = 746.09;
97  vnl_vector_fixed<VnlValueType, 7> vnlVector_2;
98  vnl_vector_fixed<VnlValueType, 7> vnlVector_3;
99  for (int i = 0; i < 7; i++)
100  {
101  if (i < 3)
102  {
103  vnlVector_1.put(i, itkVector_1[i]);
104  }
105 
106  vnlVector_2[i] = vnlVector_1[i] - mitk::eps * 1.1f;
107  vnlVector_3[i] = vnlVector_1[i] - mitk::eps * 0.9f;
108  }
109 
110  MITK_TEST_CONDITION((mitk::Equal<VnlValueType, 7>(vnlVector_1, vnlVector_1)), "vnl_fixed : v_1 == v_1 ");
111  // the v_2 is constructed so that the equality test fails for mitk::eps, the norm of the difference between the
112  // vectors is 7 * eps/6.9
113  MITK_TEST_CONDITION(!(mitk::Equal<VnlValueType, 7>(vnlVector_1, vnlVector_2)),
114  "vnl_fixed : v_1 != v_2 with mitk::eps ");
115  // increase the epsilon value used for testing equality - should now pass ( 1.2 * mitk::eps > 7 * mitk::eps/6.9 )
116  MITK_TEST_CONDITION((mitk::Equal<VnlValueType, 7>(vnlVector_1, vnlVector_2, mitk::eps * 1.2f)),
117  "vnl_fixed : v_1 == v_2 with eps = 1.2 * mitk::eps ");
118  MITK_TEST_CONDITION((mitk::Equal<VnlValueType, 7>(vnlVector_1, vnlVector_3, mitk::eps)),
119  "vnl_fixed : v_1 == v_3 with eps = 0.8 * mitk::eps ");
120  MITK_TEST_CONDITION(!(mitk::Equal<VnlValueType, 7>(vnlVector_1, vnlVector_3, mitk::eps * 0.8f)),
121  "vnl_fixed : v_1 != v_3 with eps = 0.8 * mitk::eps ");
122 
123  // test scalar equality method
124  mitk::ScalarType scalar1 = 0.5689;
125  mitk::ScalarType scalar2 = scalar1 + mitk::eps * 1.01;
126  mitk::ScalarType scalar3 = scalar1;
127  mitk::ScalarType scalar4 = scalar1 + mitk::eps * 0.95;
128  MITK_TEST_CONDITION(mitk::Equal(scalar1, scalar1), "Test scalar equality using the same scalar with mitk::eps");
129  MITK_TEST_CONDITION(!mitk::Equal(scalar1, scalar2),
130  "Test scalar equality using the different scalars with a difference greater than mitk::eps");
131  MITK_TEST_CONDITION(mitk::Equal(scalar1, scalar3),
132  "Test scalar equality using the different scalars with a difference equal to mitk::eps");
133  MITK_TEST_CONDITION(mitk::Equal(scalar1, scalar4),
134  "Test scalar equality using the different scalars with a difference less than mitk::eps");
135 
136  // test matrix equality methods
137  vnl_matrix_fixed<mitk::ScalarType, 3, 3> vnlMatrix3x3_1;
138  vnlMatrix3x3_1(0, 0) = 1.1;
139  vnlMatrix3x3_1(0, 1) = 0.4;
140  vnlMatrix3x3_1(0, 2) = 5.3;
141  vnlMatrix3x3_1(1, 0) = 2.7;
142  vnlMatrix3x3_1(1, 1) = 3578.56418;
143  vnlMatrix3x3_1(1, 2) = 123.56;
144  vnlMatrix3x3_1(2, 0) = 546.89;
145  vnlMatrix3x3_1(2, 1) = 0.0001;
146  vnlMatrix3x3_1(2, 2) = 1.0;
147  vnl_matrix_fixed<mitk::ScalarType, 3, 3> vnlMatrix3x3_2;
148  vnlMatrix3x3_2(0, 0) = 1.1000009;
149  vnlMatrix3x3_2(0, 1) = 0.4000009;
150  vnlMatrix3x3_2(0, 2) = 5.3000009;
151  vnlMatrix3x3_2(1, 0) = 2.7000009;
152  vnlMatrix3x3_2(1, 1) = 3578.5641809;
153  vnlMatrix3x3_2(1, 2) = 123.5600009;
154  vnlMatrix3x3_2(2, 0) = 546.8900009;
155  vnlMatrix3x3_2(2, 1) = 0.0001009;
156  vnlMatrix3x3_2(2, 2) = 1.0000009;
157 
158  mitk::ScalarType epsilon = 0.000001;
159  MITK_TEST_CONDITION(mitk::MatrixEqualElementWise(vnlMatrix3x3_1, vnlMatrix3x3_1, mitk::eps),
160  "Test for matrix equality with given epsilon=mitk::eps and exactly the same matrix elements");
161  MITK_TEST_CONDITION(!mitk::MatrixEqualElementWise(vnlMatrix3x3_1, vnlMatrix3x3_2, 0.0),
162  "Test for matrix equality with given epsilon=0.0 and slightly different matrix elements");
163  MITK_TEST_CONDITION(mitk::MatrixEqualElementWise(vnlMatrix3x3_1, vnlMatrix3x3_2, epsilon),
164  "Test for matrix equality with given epsilon and slightly different matrix elements");
165  MITK_TEST_CONDITION(!mitk::MatrixEqualRMS(vnlMatrix3x3_1, vnlMatrix3x3_2, 0.0),
166  "Test for matrix equality with given epsilon=0.0 and slightly different matrix elements");
167  MITK_TEST_CONDITION(mitk::MatrixEqualRMS(vnlMatrix3x3_1, vnlMatrix3x3_2, epsilon),
168  "Test for matrix equality with given epsilon and slightly different matrix elements");
169 
170  MITK_TEST_END();
171 }
mitk::MatrixEqualElementWise
bool MatrixEqualElementWise(const vnl_matrix_fixed< TCoordRep, NRows, NCols > &matrix1, const vnl_matrix_fixed< TCoordRep, NRows, NCols > &matrix2, mitk::ScalarType epsilon=mitk::eps)
Check for element-wise matrix equality with a user defined accuracy.
Definition: mitkMatrix.h:144
mitk::VnlVector
vnl_vector< ScalarType > VnlVector
Definition: mitkVector.h:138
mitk::ScalarType
double ScalarType
Definition: mitkNumericConstants.h:24
MITK_TEST_BEGIN
section GeneralTestsDeprecatedOldTestingStyle Deprecated macros All tests with MITK_TEST_BEGIN()
MITK_TEST_CONDITION
#define MITK_TEST_CONDITION(COND, MSG)
Definition: mitkTestingMacros.h:114
mitk::MatrixEqualRMS
bool MatrixEqualRMS(const vnl_matrix_fixed< TCoordRep, NRows, NCols > &matrix1, const vnl_matrix_fixed< TCoordRep, NRows, NCols > &matrix2, mitk::ScalarType epsilon=mitk::eps)
Check for matrix equality with a user defined accuracy. As an equality metric the root mean squared e...
Definition: mitkMatrix.h:100
mitk::Equal
MITKNEWMODULE_EXPORT bool Equal(mitk::ExampleDataStructure *leftHandSide, mitk::ExampleDataStructure *rightHandSide, mitk::ScalarType eps, bool verbose)
Returns true if the example data structures are considered equal.
Definition: mitkExampleDataStructure.cpp:60
mitkVectorTest
int mitkVectorTest(int, char *[])
Definition: mitkVectorTest.cpp:22
mitk::eps
MITKCORE_EXPORT const ScalarType eps
Definition: mitkNumericConstants.cpp:21
MITK_TEST_END
and MITK_TEST_END()