Commit 895cbd3a authored by Steffen Schotthöfer's avatar Steffen Schotthöfer
Browse files

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


Former-commit-id: 86298326
parent a1ba7cd2
......@@ -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 );
}
......
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
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
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 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 )
......@@ -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 ----
%
......
......@@ -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 ----
%
......
......@@ -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 ----
%
......
......@@ -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 ----
%
......
......@@ -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 ----
%
......
......@@ -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 ----
%
......
......@@ -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 ----
%
......
......@@ -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 ----
%
......
......@@ -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() );
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment