Commit 6adbe7cc authored by steffen.schotthoefer's avatar steffen.schotthoefer
Browse files

added MN linesource test case

parent 62955d62
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Linesource Benchmarking File MN %
% Author <Steffen Schotthöfer> %
% Date 01.07.2020 %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% ---- File specifications ----
OUTPUT_DIR = ../../result
OUTPUT_FILE = rtsn_test_linesource_MN
LOG_DIR = ../../result/logs
MESH_FILE = linesource.su2
% ---- Solver specifications ----
CFL_NUMBER = 0.5
TIME_FINAL = 0.3
PROBLEM = LINESOURCE
SOLVER = MN_SOLVER
MAX_MOMENT_SOLVER = 2
ENTROPY_FUNCTIONAL = QUADRATIC
ENTROPY_OPTIMIZER = NEWTON
%
% ---- Boundary Conditions ----
BC_DIRICHLET = ( void )
QUAD_TYPE = GAUSS_LEGENDRE_TENSORIZED
QUAD_ORDER = 2
This source diff could not be displayed because it is too large. You can view the blob instead.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Linesource Benchmarking File PN %
% Author <Steffen Schotthöfer> %
% Date 01.07.2020 %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% ---- File specifications ----
OUTPUT_DIR = ../../result
OUTPUT_FILE = rtsn_test_linesource_PN
......
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Linesource Benchmarking File MN %
% Author <Steffen Schotthöfer> %
% Date 01.07.2020 %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% ---- File specifications ----
OUTPUT_DIR = ../../result
OUTPUT_FILE = rtsn_test_linesource
......
......@@ -81,3 +81,21 @@ TEST_CASE( "linesource_PN", "[testcases]" ) {
REQUIRE( std::fabs( test[i] - reference[i] ) < eps );
}
}
TEST_CASE( "linesource_MN", "[testcases]" ) {
char config_file_name[MAX_STRING_SIZE] = "../tests/input/linesource_MN.cfg";
Config* config = new Config( config_file_name );
Solver* solver = Solver::Create( config );
solver->Solve();
solver->Save();
auto test = readVTKFile( "../result/rtsn_test_linesource_PN.vtk" );
auto reference = readVTKFile( "../tests/input/linesource_MN_reference.vtk" );
double eps = 1e-3;
REQUIRE( test.size() == reference.size() );
for( unsigned i = 0; i < test.size(); ++i ) {
REQUIRE( std::fabs( test[i] - reference[i] ) < eps );
}
}
Supports Markdown
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