small changes

README.md

@@ -98,7 +98,7 @@ Append `HINTS VTK_INSTALL_DIR` to the `find_package( VTK ... )` line in the CMak
 find_package( VTK REQUIRED COMPONENTS vtkIOGeometry vtkFiltersCore HINTS ~/VTK-install )
 ```
 
-Compile it
+Compile it (make sure you use GCC as compiler: module load compiler/gnu/10.2) 
 ```bash
 module load devel/cmake/3.16
 module load compiler/gnu/9.2

code/input/exampleDataGen.cfg

@@ -7,9 +7,9 @@
 %
 % ---- Global settings ----
 DATA_GENERATOR_MODE = YES
-TRAINING_SET_SIZE = 500
-MAX_VALUE_FIRST_MOMENT = 100
-BOUNDARY_DISTANCE_REALIZABLE_SET = 0.1
+TRAINING_SET_SIZE = 100000
+MAX_VALUE_FIRST_MOMENT = 10
+BOUNDARY_DISTANCE_REALIZABLE_SET = 0.01
 %
 % ---- File specifications ----
 % 
@@ -23,10 +23,10 @@ LOG_DIR = ../result/logs
 % --- Spherical Basis ----
 %
 % Choice of basis functions
-SPHERICAL_BASIS = SPHERICAL_MONOMIALS
+SPHERICAL_BASIS = SPHERICAL_HARMONICS
 %
 % Maximal Moment degree
-MAX_MOMENT_SOLVER = 1
+MAX_MOMENT_SOLVER = 0
 %
 % --- Entropy settings ----
 ENTROPY_FUNCTIONAL = MAXWELL_BOLTZMANN

code/input/exampleSN.cfg

@@ -10,28 +10,25 @@
 % Output directory
 OUTPUT_DIR = ../result
 % Output file
-OUTPUT_FILE = example_IC3_lCFL
+OUTPUT_FILE = example_SN
 % Log directory
 LOG_DIR = ../result/logs
 % Mesh File
 MESH_FILE = linesource.su2
 %
-PROBLEM = ELECTRONRT
+PROBLEM = LINESOURCE
 %
 % ---- Solver specifications ----
 %
 % Solver type
-SOLVER = CSD_SN_SOLVER
-% Slowing down approximation
-CONTINUOUS_SLOWING_DOWN = YES
+SOLVER = SN_SOLVER
 % CFL number
-CFL_NUMBER = 0.9
+CFL_NUMBER = 0.7
 % Final time for simulation
 TIME_FINAL = 0.5
 % Reconstruction order
 RECONS_ORDER = 2
 %
-CLEAN_FLUX_MATRICES = NO
 
 % ---- Boundary Conditions ----
 % Example: BC_DIRICLET = (dummyMarker1, dummyMarker2)
@@ -40,11 +37,11 @@ BC_DIRICHLET = ( void )
 %
 %% Quadrature Specifications
 % Quadrature Type
-QUAD_TYPE = LEBEDEV
+QUAD_TYPE = GAUSS_LEGENDRE_TENSORIZED
 % Quadrature Order
 QUAD_ORDER = 8
 %
 % ----- Output ---- 
 %
 VOLUME_OUTPUT = (ANALYTIC,  MINIMAL)
-OUTPUT_FREQUENCY = 1
+VOLUME_OUTPUT_FREQUENCY = 1

code/src/optimizers/newtonoptimizer.cpp

@@ -89,7 +89,7 @@ void NewtonOptimizer::SolveMultiCell( VectorVector& lambda, VectorVector& sol, c
 
 void NewtonOptimizer::Solve( Vector& lambda, Vector& sol, const VectorVector& moments, unsigned idx_cell ) {
 
-    /* solve the problem argmin ( <eta(alpha*m)>-alpha*u))
+    /* solve the problem argmin ( <eta_*(alpha*m)>-alpha*u))
      * where alpha = Lagrange multiplier
      *           m = moment basis
      *           u = current "moment solution"

code/src/toolboxes/datagenerator.cpp

@@ -158,7 +158,7 @@ void nnDataGenerator::SampleSolutionU() {
             }
         }
     }
-    else {
+    if( _LMaxDegree > 1 ) {
         ErrorMessages::Error( "Sampling for order higher than 1 is not yet supported", CURRENT_FUNCTION );
     }
 }