2525 FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
2626 DEALINGS IN THE SOFTWARE.
2727*/
28-
28+
2929#include < gtest/gtest.h>
3030#include < Core/Math/MiscMath.h>
3131#include < Core/Datatypes/Legacy/Field/VField.h>
@@ -49,46 +49,46 @@ namespace
4949{
5050 FieldHandle CreateGetMeshSolution ()
5151 {
52- return loadFieldFromFile (TestResources::rootDir () / " generateroistatistics/mesh_sol.fld" );
52+ return loadFieldFromFile (TestResources::rootDir () / " Fields/ generateroistatistics/mesh_sol.fld" );
5353 }
54-
54+
5555 FieldHandle CreateCoordinatesCubesShiftedBy10 ()
5656 {
57- return loadFieldFromFile (TestResources::rootDir () / " generateroistatistics/coordinates_cubes_shifted_by_10.fld" );
57+ return loadFieldFromFile (TestResources::rootDir () / " Fields/ generateroistatistics/coordinates_cubes_shifted_by_10.fld" );
5858 }
59-
59+
6060 FieldHandle CreateCoordinatesCubes ()
6161 {
62- return loadFieldFromFile (TestResources::rootDir () / " generateroistatistics/coordinates_cubes.fld" );
63- }
64-
62+ return loadFieldFromFile (TestResources::rootDir () / " Fields/ generateroistatistics/coordinates_cubes.fld" );
63+ }
64+
6565 FieldHandle CreateAtlas ()
6666 {
67- return loadFieldFromFile (TestResources::rootDir () / " generateroistatistics/atlas.fld" );
68- }
69-
70- DenseMatrixHandle expected_result_no_roi_spec ()
67+ return loadFieldFromFile (TestResources::rootDir () / " Fields/ generateroistatistics/atlas.fld" );
68+ }
69+
70+ DenseMatrixHandle expected_result_no_roi_spec ()
7171 {
7272 DenseMatrixHandle expected_result (boost::make_shared<DenseMatrix>(8 , 5 ));
7373
74- (*expected_result)(0 , 0 ) = 0.5398800608953688 ; (*expected_result)(1 , 0 ) = 0.7546866819823609 ; (*expected_result)(2 , 0 ) = 0.3784870349158202 ; (*expected_result)(3 , 0 ) = 0.5215146855270543 ; (*expected_result)(4 , 0 ) = 0.6204598403245096 ; (*expected_result)(5 , 0 ) = 0.4103867315361441 ; (*expected_result)(6 , 0 ) = 0.6164145705872455 ; (*expected_result)(7 , 0 ) = 0.3773992798882833 ;
75- (*expected_result)(0 , 1 ) = 0.2023192248404992 ; (*expected_result)(1 , 1 ) = std::numeric_limits<double >::quiet_NaN (); (*expected_result)(2 , 1 ) = 0.270036948169627 ; (*expected_result)(3 , 1 ) = 0.2819928747450936 ; (*expected_result)(4 , 1 ) = 0.2464919130115238 ; (*expected_result)(5 , 1 ) = 0.348358459345021 ; (*expected_result)(6 , 1 ) = 0.3382031144990396 ; (*expected_result)(7 , 1 ) = 0.1699586184078796 ;
76- (*expected_result)(0 , 2 ) = 0.1868726045543786 ; (*expected_result)(1 , 2 ) = 0.7546866819823609 ; (*expected_result)(2 , 2 ) = 0.1189976815583766 ; (*expected_result)(3 , 2 ) = 0.223811939491137 ; (*expected_result)(4 , 2 ) = 0.2550951154592691 ; (*expected_result)(5 , 2 ) = 0.1386244428286791 ; (*expected_result)(6 , 2 ) = 0.2435249687249893 ; (*expected_result)(7 , 2 ) = 0.1965952504312082 ;
77- (*expected_result)(0 , 3 ) = 0.7951999011370632 ; (*expected_result)(1 , 3 ) = 0.7546866819823609 ; (*expected_result)(2 , 3 ) = 0.6797026768536748 ; (*expected_result)(3 , 3 ) = 0.9597439585160811 ; (*expected_result)(4 , 3 ) = 0.8909032525357985 ; (*expected_result)(5 , 3 ) = 0.9592914252054443 ; (*expected_result)(6 , 3 ) = 0.9292636231872278 ; (*expected_result)(7 , 3 ) = 0.6160446761466392 ;
74+ (*expected_result)(0 , 0 ) = 0.5398800608953688 ; (*expected_result)(1 , 0 ) = 0.7546866819823609 ; (*expected_result)(2 , 0 ) = 0.3784870349158202 ; (*expected_result)(3 , 0 ) = 0.5215146855270543 ; (*expected_result)(4 , 0 ) = 0.6204598403245096 ; (*expected_result)(5 , 0 ) = 0.4103867315361441 ; (*expected_result)(6 , 0 ) = 0.6164145705872455 ; (*expected_result)(7 , 0 ) = 0.3773992798882833 ;
75+ (*expected_result)(0 , 1 ) = 0.2023192248404992 ; (*expected_result)(1 , 1 ) = std::numeric_limits<double >::quiet_NaN (); (*expected_result)(2 , 1 ) = 0.270036948169627 ; (*expected_result)(3 , 1 ) = 0.2819928747450936 ; (*expected_result)(4 , 1 ) = 0.2464919130115238 ; (*expected_result)(5 , 1 ) = 0.348358459345021 ; (*expected_result)(6 , 1 ) = 0.3382031144990396 ; (*expected_result)(7 , 1 ) = 0.1699586184078796 ;
76+ (*expected_result)(0 , 2 ) = 0.1868726045543786 ; (*expected_result)(1 , 2 ) = 0.7546866819823609 ; (*expected_result)(2 , 2 ) = 0.1189976815583766 ; (*expected_result)(3 , 2 ) = 0.223811939491137 ; (*expected_result)(4 , 2 ) = 0.2550951154592691 ; (*expected_result)(5 , 2 ) = 0.1386244428286791 ; (*expected_result)(6 , 2 ) = 0.2435249687249893 ; (*expected_result)(7 , 2 ) = 0.1965952504312082 ;
77+ (*expected_result)(0 , 3 ) = 0.7951999011370632 ; (*expected_result)(1 , 3 ) = 0.7546866819823609 ; (*expected_result)(2 , 3 ) = 0.6797026768536748 ; (*expected_result)(3 , 3 ) = 0.9597439585160811 ; (*expected_result)(4 , 3 ) = 0.8909032525357985 ; (*expected_result)(5 , 3 ) = 0.9592914252054443 ; (*expected_result)(6 , 3 ) = 0.9292636231872278 ; (*expected_result)(7 , 3 ) = 0.6160446761466392 ;
7878 (*expected_result)(0 , 4 ) = 9 ; (*expected_result)(1 , 4 ) = 1 ; (*expected_result)(2 , 4 ) = 5 ; (*expected_result)(3 , 4 ) = 5 ; (*expected_result)(4 , 4 ) = 5 ; (*expected_result)(5 , 4 ) = 5 ; (*expected_result)(6 , 4 ) = 5 ; (*expected_result)(7 , 4 ) = 5 ; (*expected_result)(8 , 4 ) = 5 ;
79-
79+
8080 return expected_result;
8181 }
82-
83- DenseMatrixHandle expected_result_roi_spec ()
82+
83+ DenseMatrixHandle expected_result_roi_spec ()
8484 {
8585 DenseMatrixHandle expected_result (boost::make_shared<DenseMatrix>(1 , 5 ));
8686 (*expected_result)(0 , 0 ) = 0.5046221600763444 ;
8787 (*expected_result)(0 , 1 ) = 0.2271478773489659 ;
8888 (*expected_result)(0 , 2 ) = 0.1626117351946306 ;
8989 (*expected_result)(0 , 3 ) = 0.7951999011370632 ;
9090 (*expected_result)(0 , 4 ) = 14 ;
91-
91+
9292 return expected_result;
9393 }
9494}
@@ -97,22 +97,22 @@ TEST(GenerateROIStatisticsAlgorithm, ValidInput_TestAllROI)
9797{
9898 GenerateROIStatisticsAlgorithm algo;
9999 FieldHandle atlas = CreateAtlas ();
100- FieldHandle createcoordcubes = CreateCoordinatesCubes ();
101- FieldHandle createmeshsolution = CreateGetMeshSolution ();
100+ FieldHandle createcoordcubes = CreateCoordinatesCubes ();
101+ FieldHandle createmeshsolution = CreateGetMeshSolution ();
102102 const std::string& atlasMeshLabelsStr = std::string (" Region1;Region2;Region3;Region4;Region5;Region6;Region7;Region8" );
103103 DenseMatrixHandle specROI (new DenseMatrix (5 , 1 ));
104104 (*specROI)(0 , 0 ) = 0 ;
105105 (*specROI)(1 , 0 ) = 0 ;
106106 (*specROI)(2 , 0 ) = 0 ;
107- (*specROI)(3 , 0 ) = 0 ;
108- (*specROI)(4 , 0 ) = 0 ;
107+ (*specROI)(3 , 0 ) = 0 ;
108+ (*specROI)(4 , 0 ) = 0 ;
109109 auto output = algo.run (createmeshsolution,atlas,createcoordcubes,atlasMeshLabelsStr, specROI);
110110 DenseMatrixHandle outputMatrix = output.get <0 >();
111-
111+
112112 DenseMatrixHandle expected_result = expected_result_no_roi_spec ();
113113 EXPECT_EQ (outputMatrix->rows (),expected_result->rows ());
114- EXPECT_EQ (outputMatrix->cols (),expected_result->cols ());
115-
114+ EXPECT_EQ (outputMatrix->cols (),expected_result->cols ());
115+
116116 for (int i = 0 ; i < outputMatrix->rows (); i++)
117117 for (int j = 0 ; j < outputMatrix->cols (); j++)
118118 if (i==1 && j==1 )
@@ -121,62 +121,62 @@ TEST(GenerateROIStatisticsAlgorithm, ValidInput_TestAllROI)
121121 {
122122 FAIL () << " ERROR: Test on NaN failed!" << std::endl;
123123 }
124- }
124+ }
125125 else
126126 {
127127 EXPECT_NEAR ((*outputMatrix)(i, j), (*expected_result)(i,j), 1e-10 );
128128 }
129-
129+
130130}
131131
132- TEST (GenerateROIStatisticsAlgorithm, ValidInput_TestSpecROI)
132+ TEST (GenerateROIStatisticsAlgorithm, ValidInput_TestSpecROI)
133133{
134134 GenerateROIStatisticsAlgorithm algo;
135135 FieldHandle atlas = CreateAtlas ();
136- FieldHandle createcoordcubes = CreateCoordinatesCubes ();
137- FieldHandle createmeshsolution = CreateGetMeshSolution ();
136+ FieldHandle createcoordcubes = CreateCoordinatesCubes ();
137+ FieldHandle createmeshsolution = CreateGetMeshSolution ();
138138 const std::string& atlasMeshLabelsStr = std::string (" Region1;Region2;Region3;Region4;Region5;Region6;Region7;Region8" );
139139 DenseMatrixHandle specROI (new DenseMatrix (5 , 1 ));
140140 (*specROI)(0 , 0 ) = 0.5 ;
141141 (*specROI)(1 , 0 ) = 0.5 ;
142142 (*specROI)(2 , 0 ) = 0.5 ;
143- (*specROI)(3 , 0 ) = 0 ;
144- (*specROI)(4 , 0 ) = 1 ;
143+ (*specROI)(3 , 0 ) = 0 ;
144+ (*specROI)(4 , 0 ) = 1 ;
145145 auto output = algo.run (createmeshsolution,atlas,createcoordcubes,atlasMeshLabelsStr, specROI);
146146 DenseMatrixHandle outputMatrix = output.get <0 >();
147-
147+
148148 DenseMatrixHandle expected_result = expected_result_roi_spec ();
149149 EXPECT_EQ (outputMatrix->rows (),expected_result->rows ());
150- EXPECT_EQ (outputMatrix->cols (),expected_result->cols ());
151-
150+ EXPECT_EQ (outputMatrix->cols (),expected_result->cols ());
151+
152152 for (int i = 0 ; i < outputMatrix->rows (); i++)
153153 for (int j = 0 ; j < outputMatrix->cols (); j++)
154154 EXPECT_NEAR ((*outputMatrix)(i, j), (*expected_result)(i,j), 1e-10 );
155-
155+
156156}
157157
158- TEST (GenerateROIStatisticsAlgorithm, ValidInput_TestCoordinateSpace)
158+ TEST (GenerateROIStatisticsAlgorithm, ValidInput_TestCoordinateSpace)
159159{
160160 GenerateROIStatisticsAlgorithm algo;
161- FieldHandle atlas = CreateAtlas ();
161+ FieldHandle atlas = CreateAtlas ();
162162 FieldHandle createcoordcubesshiftby10 = CreateCoordinatesCubesShiftedBy10 ();
163- FieldHandle createmeshsolution = CreateGetMeshSolution ();
163+ FieldHandle createmeshsolution = CreateGetMeshSolution ();
164164 const std::string& atlasMeshLabelsStr = std::string (" Region1;Region2;Region3;Region4;Region5;Region6;Region7;Region8" );
165165 DenseMatrixHandle specROI (new DenseMatrix (5 , 1 ));
166166 (*specROI)(0 , 0 ) = 10.5 ;
167167 (*specROI)(1 , 0 ) = 10.5 ;
168168 (*specROI)(2 , 0 ) = 10.5 ;
169- (*specROI)(3 , 0 ) = 0 ;
170- (*specROI)(4 , 0 ) = 1.0 ;
169+ (*specROI)(3 , 0 ) = 0 ;
170+ (*specROI)(4 , 0 ) = 1.0 ;
171171 auto output = algo.run (createmeshsolution,atlas,createcoordcubesshiftby10,atlasMeshLabelsStr, specROI);
172172 DenseMatrixHandle outputMatrix = output.get <0 >();
173-
173+
174174 DenseMatrixHandle expected_result = expected_result_roi_spec ();
175175 EXPECT_EQ (outputMatrix->rows (),expected_result->rows ());
176- EXPECT_EQ (outputMatrix->cols (),expected_result->cols ());
177-
176+ EXPECT_EQ (outputMatrix->cols (),expected_result->cols ());
177+
178178 for (int i = 0 ; i < outputMatrix->rows (); i++)
179179 for (int j = 0 ; j < outputMatrix->cols (); j++)
180180 EXPECT_NEAR ((*outputMatrix)(i, j), (*expected_result)(i,j), 1e-10 );
181-
181+
182182}
0 commit comments