changed mesh file directory, added csd waterphantom test. waterphantom test...

changed mesh file directory, added csd waterphantom test. waterphantom test has error with absolute paths in icru class

code/src/solvers/mnsolver.cpp

@@ -109,8 +109,6 @@ Vector MNSolver::ConstructFlux( unsigned idx_cell ) {
             entropyFlux += _g->Flux( _quadPoints[idx_quad], entropyL, entropyR, _normals[idx_cell][idx_neigh] );
         }
         flux += _moments[idx_quad] * ( _weights[idx_quad] * entropyFlux );
-
-        // ------- Relizablity Reconstruction Step ----
     }
     return flux;
 }
@@ -128,10 +126,9 @@ void MNSolver::ComputeRealizableSolution( unsigned idx_cell ) {
 void MNSolver::IterPreprocessing() {
 
     // ------- Reconstruction Step ----------------
-
     _optimizer->SolveMultiCell( _alpha, _sol, _moments );
 
-    // ------- Relizablity Reconstruction Step ----
+    // ------- Relizablity Preservation Step ----
     for( unsigned idx_cell = 0; idx_cell < _nCells; idx_cell++ ) {
         ComputeRealizableSolution( idx_cell );
     }

code/tests/input/validation_tests/CSD_SN_FP_solver/stopping_power.txt

+1.000E-02 2.256E+01
+1.250E-02 1.898E+01
+1.500E-02 1.647E+01
+1.750E-02 1.461E+01
+2.000E-02 1.318E+01
+2.500E-02 1.110E+01
+3.000E-02 9.657E+00
+3.500E-02 8.596E+00
+4.000E-02 7.781E+00
+4.500E-02 7.134E+00
+5.000E-02 6.607E+00
+5.500E-02 6.170E+00
+6.000E-02 5.801E+00
+7.000E-02 5.211E+00
+8.000E-02 4.761E+00
+9.000E-02 4.407E+00
+1.000E-01 4.119E+00
+1.250E-01 3.596E+00
+1.500E-01 3.242E+00
+1.750E-01 2.988E+00
+2.000E-01 2.798E+00
+2.500E-01 2.533E+00
+3.000E-01 2.360E+00
+3.500E-01 2.241E+00
+4.000E-01 2.154E+00
+4.500E-01 2.090E+00
+5.000E-01 2.041E+00
+5.500E-01 2.003E+00
+6.000E-01 1.972E+00
+7.000E-01 1.926E+00
+8.000E-01 1.896E+00
+9.000E-01 1.876E+00
+1.000E+00 1.862E+00
+1.250E+00 1.845E+00
+1.500E+00 1.841E+00
+1.750E+00 1.844E+00
+2.000E+00 1.850E+00
+2.500E+00 1.868E+00
+3.000E+00 1.889E+00
+3.500E+00 1.910E+00
+4.000E+00 1.931E+00
+4.500E+00 1.951E+00
+5.000E+00 1.971E+00
+5.500E+00 1.991E+00
+6.000E+00 2.010E+00
+7.000E+00 2.047E+00
+8.000E+00 2.082E+00
+9.000E+00 2.116E+00
+1.000E+01 2.149E+00
+1.250E+01 2.230E+00
+1.500E+01 2.306E+00
+1.750E+01 2.381E+00
+2.000E+01 2.454E+00
+2.500E+01 2.598E+00
+3.000E+01 2.738E+00
+3.500E+01 2.876E+00
+4.000E+01 3.013E+00
+4.500E+01 3.150E+00
+5.000E+01 3.286E+00
+5.500E+01 3.421E+00
+6.000E+01 3.556E+00
+7.000E+01 3.827E+00
+8.000E+01 4.096E+00
+9.000E+01 4.366E+00
+1.000E+02 4.636E+00
+1.250E+02 5.311E+00
+1.500E+02 5.987E+00
+1.750E+02 6.663E+00
+2.000E+02 7.341E+00
+2.500E+02 8.698E+00
+3.000E+02 1.006E+01
+3.500E+02 1.142E+01
+4.000E+02 1.278E+01
+4.500E+02 1.414E+01
+5.000E+02 1.551E+01
+5.500E+02 1.688E+01
+6.000E+02 1.824E+01
+7.000E+02 2.098E+01
+8.000E+02 2.371E+01
+9.000E+02 2.645E+01
+1.000E+03 2.919E+01

code/tests/input/validation_tests/CSD_SN_FP_solver/waterphantom_1D.cfg

+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%  Waterphantom Benchmarking File CSD  %
+%  Author <P. Stammer, J. Kusch>       %
+%  Date   01.12.2020                   %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%
+% ---- File specifications ----
+%
+OUTPUT_DIR = ../../../result
+OUTPUT_FILE = rtsn_test_waterphantom_1D_CSD_FP
+LOG_DIR = ../../../result/logs
+MESH_FILE = ../../mesh_files/waterphantom1D.su2
+%
+% ---- Problem specifications ----
+%
+PROBLEM = WATERPHANTOM
+HYDROGEN_FILE = ENDL_H.txt
+OXYGEN_FILE = ENDL_O.txt
+STOPPING_POWER_FILE = stopping_power.txt
+%
+% ---- Solver specifications ----
+%
+SOLVER = CSD_SN_FOKKERPLANCK_TRAFO_SOLVER
+KERNEL = ISOTROPIC_1D
+CFL_NUMBER = 0.008
+TIME_FINAL = 1.0
+%
+% ---- Quadrature ----
+%
+QUAD_TYPE = GAUSS_LEGENDRE_1D
+QUAD_ORDER = 26
+%
+BC_DIRICHLET = ( dirichlet )
+BC_NEUMANN = ( wall_low, wall_up )

code/tests/input/validation_tests/MN_solver/checkerboard_MN.cfg

@@ -9,7 +9,7 @@
 OUTPUT_DIR = ../../../result
 OUTPUT_FILE = rtsn_test_checkerboard_MN
 LOG_DIR = ../../../result/logs
-MESH_FILE = ../meshes/checkerboard.su2
+MESH_FILE = ../../mesh_files/checkerboard.su2
 %
 % ---- Problem specifications ----
 %

code/tests/input/validation_tests/MN_solver/linesource_MN_MB.cfg

@@ -9,7 +9,7 @@
 OUTPUT_DIR = ../../../result
 OUTPUT_FILE = rtsn_test_linesource_MN_MB
 LOG_DIR = ../../../result/logs
-MESH_FILE = ../meshes/linesource.su2
+MESH_FILE = ../../mesh_files/linesource.su2
 %
 % ---- Problem specifications ----
 %

code/tests/input/validation_tests/MN_solver/linesource_MN_Quad.cfg

@@ -9,7 +9,7 @@
 OUTPUT_DIR = ../../../result
 OUTPUT_FILE = rtsn_test_linesource_MN_Quad
 LOG_DIR = ../../../result/logs
-MESH_FILE = ../meshes/linesource.su2
+MESH_FILE = ../../mesh_files/linesource.su2
 %
 % ---- Problem specifications ----
 %

code/tests/input/validation_tests/PN_solver/checkerboard_PN.cfg

@@ -9,7 +9,7 @@
 OUTPUT_DIR = ../../../result
 OUTPUT_FILE = rtsn_test_checkerboard_PN
 LOG_DIR = ../../../result/logs
-MESH_FILE = ../meshes/checkerboard.su2
+MESH_FILE = ../../mesh_files/checkerboard.su2
 %
 % ---- Problem specifications ----
 %

code/tests/input/validation_tests/PN_solver/linesource_PN.cfg

@@ -8,7 +8,7 @@
 OUTPUT_DIR = ../../../result
 OUTPUT_FILE = rtsn_test_linesource_PN
 LOG_DIR = ../../../result/logs
-MESH_FILE = ../meshes/linesource.su2
+MESH_FILE = ../../mesh_files/linesource.su2
 %
 % ---- Problem specifications ----
 %

code/tests/input/validation_tests/SN_solver/checkerboard_SN.cfg

@@ -9,7 +9,7 @@
 OUTPUT_DIR = ../../../result
 OUTPUT_FILE = rtsn_test_checkerboard_SN
 LOG_DIR = ../../../result/logs
-MESH_FILE = ../meshes/checkerboard.su2
+MESH_FILE = ../../mesh_files/checkerboard.su2
 %
 % ---- Problem specifications ----
 %

code/tests/input/validation_tests/SN_solver/linesource_SN.cfg

@@ -8,7 +8,7 @@
 OUTPUT_DIR = ../../../result
 OUTPUT_FILE = rtsn_test_linesource_SN
 LOG_DIR = ../../../result/logs
-MESH_FILE = ../meshes/linesource.su2
+MESH_FILE = ../../mesh_files/linesource.su2
 %
 % ---- Problem specifications ----
 %

code/tests/input/validation_tests/output/validate_logger.cfg

@@ -9,7 +9,7 @@ OUTPUT_DIR = ../../../result
 OUTPUT_FILE = validate_screen_output
 LOG_DIR = ../../../result/logs
 LOG_FILE = validate_logger_output
-MESH_FILE = ../meshes/linesource.su2
+MESH_FILE = ../../mesh_files/linesource.su2
 %
 % ---- Problem specifications ----
 %

code/tests/test_cases.cpp

@@ -190,18 +190,22 @@ TEST_CASE( "MN_SOLVER", "[validation_tests]" ) {
 }
 
 TEST_CASE( "CSD_SN_FP_SOLVER", "[validation_tests]" ) {
-    std::string sn_fileDir = "input/validation_tests/CSD_SN_FP_solver/";
+    std::string csd_sn_fileDir = "input/validation_tests/CSD_SN_FP_solver/";
     SECTION( "waterphantom 1D" ) {
-        std::string config_file_name = std::string( TESTS_PATH ) + sn_fileDir + "waterphantom_1D.cfg";
+        std::cout << "here\n";
+
+        std::string config_file_name = std::string( TESTS_PATH ) + csd_sn_fileDir + "waterphantom_1D.cfg";
+        std::cout << "here\n";
 
         Config* config = new Config( config_file_name );
         Solver* solver = Solver::Create( config );
-
+        std::cout << "here\n";
         solver->Solve();
+        std::cout << "here\n";
         solver->PrintVolumeOutput();
-
+        std::cout << "here\n";
         auto test      = readVTKFile( std::string( TESTS_PATH ) + "result/rtsn_test_waterphantom_1D_CSD_FP_1.vtk" );
-        auto reference = readVTKFile( std::string( TESTS_PATH ) + sn_fileDir + "waterphantom_1D_CSD_SN_FP_reference.vtk" );
+        auto reference = readVTKFile( std::string( TESTS_PATH ) + csd_sn_fileDir + "waterphantom_1D_CSD_SN_FP_reference.vtk" );
 
         double eps = 1e-3;
         REQUIRE( test.size() == reference.size() );