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

low cfl test


Former-commit-id: 0758f531
parent a090dc2b
Pipeline #147662 failed with stage
......@@ -8,7 +8,7 @@
% ---- File specifications ----
%
OUTPUT_DIR = ../result
OUTPUT_FILE = IsotropicPointSource_sph_nq6
OUTPUT_FILE = volOut_07_moments
LOG_DIR = ../result/logs
MESH_FILE = meshes/phantom2D.su2
%
......@@ -24,7 +24,7 @@ KERNEL = ISOTROPIC
% Solver type
SOLVER = CSD_PN
% CFL number
CFL_NUMBER = 0.6
CFL_NUMBER = 0.006
% Final time for simulation
TIME_FINAL = 0.1
% Maximal Moment degree
......@@ -45,7 +45,7 @@ QUAD_ORDER = 6
%
% ----- Output ----
%
VOLUME_OUTPUT = (MINIMAL, MEDICAL)
VOLUME_OUTPUT = (MINIMAL, MEDICAL, MOMENTS)
VOLUME_OUTPUT_FREQUENCY = 1
SCREEN_OUTPUT = (ITER, MASS, RMS_FLUX, VTK_OUTPUT, CSV_OUTPUT)
SCREEN_OUTPUT_FREQUENCY = 1
......
......@@ -47,11 +47,15 @@ CSDPNSolver::CSDPNSolver( Config* settings ) : PNSolver( settings ) {
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 idx_sys = 0; idx_sys < _nSystem; idx_sys++ ) {
_sol[idx_cell][idx_sys] = f * StarMAPmoments[idx_sys]; // must be VectorVector
//_sol[idx_cell][idx_sys] = 0;
_sol[idx_cell][0] = f * StarMAPmoments[0];
for( unsigned idx_sys = 1; idx_sys < _nSystem; idx_sys++ ) {
_sol[idx_cell][idx_sys] = f * 0; // must be VectorVector
//_sol[idx_cell][idx_sys] = 0;
}
}
_solNew = _sol;
_sigmaS = VectorVector( _nEnergies, Vector( _polyDegreeBasis ) );
......
......@@ -340,7 +340,7 @@ void PNSolver::ComputeFluxComponents() {
// std::cout << "Spectral Radius Z " << blaze::max( blaze::abs( eigenValues ) ) << "\n";
//_combinedSpectralRadius = blaze::max( blaze::abs( eigenValues + eigenValuesX + eigenValuesY ) );
// std::cout << "Spectral Radius combined " << _combinedSpectralRadius << "\n";
std::cout << "Spectral Radius combined " << blaze::max( blaze::abs( eigenValues + eigenValuesX + eigenValuesY ) ) << "\n";
}
void PNSolver::ComputeScatterMatrix() {
......
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