Commit 81071fa4 authored by Steffen Schotthöfer's avatar Steffen Schotthöfer
Browse files

small indexing bug fix

parent 263a242f
Pipeline #142966 passed with stage
in 20 minutes and 38 seconds
......@@ -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
......@@ -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
......
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
......@@ -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 );
}
......
Markdown is supported
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