small indexing bug fix

Former-commit-id: 81071fa4

code/input/linesource_PN.cfg

@@ -40,7 +40,12 @@ RECONS_ORDER = 2
 % Dirichlet Boundary
 BC_DIRICHLET = ( void )
 %
+%
 % ----- Output ---- 
 %
-VOLUME_OUTPUT = (MINIMAL)
-
+VOLUME_OUTPUT = (MINIMAL, MOMENTS, DUAL_MOMENTS)
+VOLUME_OUTPUT_FREQUENCY = 1
+SCREEN_OUTPUT = (ITER, MASS, RMS_FLUX, VTK_OUTPUT, CSV_OUTPUT)
+SCREEN_OUTPUT_FREQUENCY = 1
+HISTORY_OUTPUT = (ITER, MASS, RMS_FLUX, VTK_OUTPUT, CSV_OUTPUT)
+HISTORY_OUTPUT_FREQUENCY = 1

code/input/pointSource_PNCSD.cfg

@@ -44,7 +44,7 @@ QUAD_ORDER = 6
 %
 % ----- Output ---- 
 %
-VOLUME_OUTPUT = (MINIMAL, MOMENTS, DUAL_MOMENTS)
+VOLUME_OUTPUT = (MINIMAL, MEDICAL)
 VOLUME_OUTPUT_FREQUENCY = 1
 SCREEN_OUTPUT = (ITER, MASS, RMS_FLUX, VTK_OUTPUT, CSV_OUTPUT)
 SCREEN_OUTPUT_FREQUENCY = 1

code/input/pointSource_sph_csd.cfg

-
-OUTPUT_DIR = ../result
-OUTPUT_FILE = IsotropicPointSource_sph_nq6
-LOG_DIR = ../result/logs
-MESH_FILE = meshes/phantom2D.su2
-PROBLEM = ISOTROPICPOINTSOURCE2D
-SOLVER = CSD_PN
-HYDROGEN_FILE = ENDL_H.txt
-OXYGEN_FILE = ENDL_O.txt
-KERNEL = ISOTROPIC
-CFL_NUMBER = 2.0
-TIME_FINAL = 1.0
-CLEAN_FLUX_MATRICES = NO
-
-BC_DIRICHLET = ( void )
-QUAD_TYPE = GAUSS_LEGENDRE_TENSORIZED
-QUAD_ORDER = 6

code/src/solvers/csdpnsolver.cpp

@@ -42,19 +42,21 @@ CSDPNSolver::CSDPNSolver( Config* settings ) : PNSolver( settings ) {
 
     Vector pos_beam = Vector{ 0.5, 0.5 };
     VectorVector IC( _nCells, Vector( _nSystem ) );
-    for( unsigned i = 0; i < _nCells; ++i ) {
-        double x            = _cellMidPoints[i][0];
-        double y            = _cellMidPoints[i][1];
+    for( unsigned idx_cell = 0; idx_cell < _nCells; ++idx_cell ) {
+        double x            = _cellMidPoints[idx_cell][0];
+        double y            = _cellMidPoints[idx_cell][1];
         const double stddev = .005;
         double f            = normpdf( x, pos_beam[0], stddev ) * normpdf( y, pos_beam[1], stddev );
-        for( unsigned j = 0; j < _nSystem; j++ ) {
-            IC[i][j] = f * StarMAPmoments[j];    // must be VectorVector
+        for( unsigned idx_sys = 0; idx_sys < _nSystem; idx_sys++ ) {
+            IC[idx_cell][idx_sys] = f * StarMAPmoments[idx_sys];    // must be VectorVector
         }
     }
 
+    printf( "%d", sigma_ref.rows() );
+
     Matrix sigma_t( _energies.size(), sigma_ref.rows() );
-    for( unsigned i = 0; i < _nSystem; ++i ) {
-        Vector xs_m = blaze::column( sigma_ref, i );
+    for( unsigned idx_degree = 0; idx_degree < _polyDegreeBasis; ++idx_degree ) {
+        Vector xs_m = blaze::column( sigma_ref, idx_degree );
         Interpolation interp( E_ref, xs_m );
         // blaze::column( sigma_t, i ) = interp( _energies );
     }