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

use PN flux

parent 9916779a
Pipeline #143090 passed with stage
in 20 minutes and 28 seconds
......@@ -32,7 +32,7 @@ class CSDPNSolver : public PNSolver
void IterPreprocessing( unsigned /*idx_iter*/ ) override;
void IterPostprocessing( unsigned /*idx_iter*/ ) override;
void FluxUpdate() override;
// void FluxUpdate() override;
void FVMUpdate( unsigned idx_energy ) override;
void PrepareVolumeOutput() override;
......
......@@ -24,7 +24,7 @@ KERNEL = ISOTROPIC
% Solver type
SOLVER = CSD_PN
% CFL number
CFL_NUMBER = 0.7
CFL_NUMBER = 0.6
% Final time for simulation
TIME_FINAL = 0.1
% Maximal Moment degree
......
......@@ -11,6 +11,7 @@
#include "toolboxes/textprocessingtoolbox.h"
// externals
#include "spdlog/spdlog.h"
#include <iostream>
#include <mpi.h>
double normpdf( double x, double mu, double sigma ) {
......@@ -104,11 +105,10 @@ void CSDPNSolver::IterPostprocessing( unsigned idx_iter ) {
else {
_dose[j] += _dE * _fluxNew[j] * _s[_nEnergies - n - 1] / _density[j];
}
_solverOutput[j] = _fluxNew[j]; // Carefull here
}
}
void CSDPNSolver::FluxUpdate() {
/*void CSDPNSolver::FluxUpdate() {
// _mesh->ReconstructSlopesU( _nSystem, _solDx, _solDy, _sol );
#pragma omp parallel for
......@@ -146,7 +146,7 @@ void CSDPNSolver::FluxUpdate() {
}
}
}
*/
void CSDPNSolver::FVMUpdate( unsigned idx_energy ) {
// loop over all spatial cells
#pragma omp parallel for
......
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