update default cfg files, prepare new solver output

Former-commit-id: 8015e5ea

code/input/exampleMN.cfg

+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%  Linesource Benchmarking File PN     %
+%  Author <Steffen Schotthöfer>		   %
+%  Date   01.07.2020                   %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%
+%
+% ---- File specifications ----
+% 
+% Output directory
+OUTPUT_DIR = ../result
+% Output file
+OUTPUT_FILE = exampleMN_ML
+% Log directory
+LOG_DIR = ../result/logs
+% Mesh File
+MESH_FILE = linesource.su2
+%MESH_FILE = linesource_debug.su2
+%
+PROBLEM = LINESOURCE
+%
+% ---- Solver specifications ----
+%
+% Solver type
+SOLVER = MN_SOLVER
+% CFL number
+CFL_NUMBER = 0.7
+% Final time for simulation
+TIME_FINAL = 0.3
+% Maximal Moment degree
+MAX_MOMENT_SOLVER = 1
+%
+%% Entropy settings
+ENTROPY_FUNCTIONAL = MAXWELL_BOLZMANN
+ENTROPY_OPTIMIZER = NEWTON
+%
+% ----- Newton Solver Specifications ----
+%
+NEWTON_FAST_MODE = NO
+NEWTON_ITER = 1000
+NEWTON_EPSILON = 0.01
+NEWTON_STEP_SIZE = 0.7
+NEWTON_LINE_SEARCH_ITER = 1000
+%
+% ---- Boundary Conditions ----
+% Example: BC_DIRICLET = (dummyMarker1, dummyMarker2)
+% Dirichlet Boundary
+BC_DIRICHLET = ( void )
+%
+% ----- Quadrature ----
+% 
+% Quadrature Rule
+%QUAD_TYPE = MONTE_CARLO
+QUAD_TYPE = GAUSS_LEGENDRE_TENSORIZED
+%
+% Quadrature Order
+QUAD_ORDER = 12

code/input/examplePN.cfg

+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%  Linesource Benchmarking File PN     %
+%  Author <Steffen Schotthöfer>		   %
+%  Date   01.07.2020                   %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%
+%
+% ---- File specifications ----
+% 
+% Output directory
+OUTPUT_DIR = ../result
+% Output file
+OUTPUT_FILE = examplePN
+% Log directory
+LOG_DIR = ../result/logs
+% Mesh File
+MESH_FILE = linesource.su2
+%MESH_FILE = linesource_debug.su2
+
+%
+PROBLEM = LINESOURCE
+%
+% ---- Solver specifications ----
+%
+% Solver type
+SOLVER = PN_SOLVER
+% CFL number
+CFL_NUMBER = 0.6
+% Final time for simulation
+TIME_FINAL = 0.3
+% Maximal Moment degree
+MAX_MOMENT_SOLVER = 2
+
+% ---- Boundary Conditions ----
+% Example: BC_DIRICLET = (dummyMarker1, dummyMarker2)
+% Dirichlet Boundary
+BC_DIRICHLET = ( void )
+%
+

code/input/exampleSN.cfg

+% This is a comment
+%%%%%%%%%%%%%%%%%%%%%%
+% Example config file
+%
+% ---- File specifications ----
+% 
+% Output directory
+OUTPUT_DIR = ../result
+% Output file
+OUTPUT_FILE = exampleSN
+% Log directory
+LOG_DIR = ../result/logs
+% Mesh File
+MESH_FILE = linesource.su2
+%
+PROBLEM = LINESOURCE
+%
+% ---- Solver specifications ----
+%
+% Solver type
+SOLVER = SN_SOLVER
+% CFL number
+CFL_NUMBER = 0.4
+% Final time for simulation
+TIME_FINAL = 0.3
+
+% ---- Boundary Conditions ----
+% Example: BC_DIRICLET = (dummyMarker1, dummyMarker2)
+% Dirichlet Boundary
+BC_DIRICHLET = ( void )
+%
+%% Quadrature Specifications
+% Quadrature Type
+QUAD_TYPE = LEBEDEV
+% Quadrature Order
+QUAD_ORDER = 5

code/src/solvers/mnsolver.cpp

@@ -223,9 +223,16 @@ void MNSolver::Save() const {
 }
 
 void MNSolver::Save( int currEnergy ) const {
-    std::vector<std::string> fieldNames{ "flux" };
+    std::vector<std::string> fieldNames = { "flux" };
+
+    // Multifield ooutput
+    // std::vector<std::string> fieldNames( _nTotalEntries, "flux" );
+    // for( unsigned idx_sys = 1; idx_sys < _nTotalEntries; idx_sys++ ) {
+    //    fieldNames[idx_sys] = std::string( "flux_moment_" + std::to_string( idx_sys ) );
+    //}
+
     std::vector<std::vector<double>> scalarField( 1, _solverOutput );
-    std::vector<std::vector<std::vector<double>>> results{ scalarField };
+    // std::vector<std::vector<std::vector<double>>> results{ _outputFields };
     ExportVTK( _settings->GetOutputFile() + "_" + std::to_string( currEnergy ), results, fieldNames, _mesh );
 }