linesource corrected

code/include/problems/problembase.h

@@ -49,7 +49,9 @@ class ProblemBase
      *        in vector energy
      * @param energy is the energy the cross section is queried for
      */
-    virtual VectorVector GetScatteringXSE( const Vector& energies, const Vector& angles ) { return VectorVector( 1, Vector( 1 ) ); }
+    virtual VectorVector GetScatteringXSE( const Vector& energies, const Vector& angles ) {
+        return VectorVector( energies.size(), Vector( angles.size() ) );
+    }
 
     /**
      * @brief GetExternalSource gives back vector of vectors of source terms for each

code/input/example.cfg

@@ -37,4 +37,4 @@ BC_DIRICHLET = ( void )
 % Quadrature Type
 QUAD_TYPE = LEBEDEV
 % Quadrature Order
-QUAD_ORDER = 47
+QUAD_ORDER = 4

code/input/example_csd.cfg

 OUTPUT_DIR = ../result
 OUTPUT_FILE = example_csd
 LOG_DIR = ../result/logs
-MESH_FILE = linesource_pseudo_1d.su2
-PROBLEM = WATERPHANTOM
+MESH_FILE = linesource_pseudo_1D.su2
+PROBLEM = LINESOURCE_PSEUDO_1D
 SOLVER = CSD_SN_SOLVER
-CONTINUOUS_SLOWING_DOWN = YES
+CONTINUOUS_SLOWING_DOWN = NO
 HYDROGEN_FILE = ENDL_H.txt
 OXYGEN_FILE = ENDL_O.txt
 KERNEL = ISOTROPIC_1D
-CFL_NUMBER = 0.0001
-TIME_FINAL = 0.00001
+CFL_NUMBER = 0.1
+TIME_FINAL = 1.0
 CLEAN_FLUX_MATRICES = NO
 
 BC_DIRICHLET = ( void )

code/src/problems/linesource.cpp

@@ -43,7 +43,7 @@ VectorVector LineSource_SN_Pseudo1D::SetupIC() {
     double t      = 3.2e-4;    // pseudo time for gaussian smoothing
     for( unsigned j = 0; j < cellMids.size(); ++j ) {
         double x = cellMids[j][0];
-        psi[j]   = 1.0 / ( 4.0 * M_PI * t ) * std::exp( -( x * x ) / ( 4 * t ) );
+        psi[j]   = 1.0 / ( 2.0 * M_PI * t ) * std::exp( -( x * x ) / ( 2.0 * t ) );
     }
     return psi;
 }