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

removed bug where energies are overwritten by csdsolvers

parent c68eeaed
Pipeline #142880 passed with stage
in 20 minutes and 27 seconds
......@@ -9,8 +9,7 @@ class CSDSNSolver : public SNSolver
std::vector<double> _dose; /*!< @brief TODO */
// Physics acess
Vector _energies; /*!< @brief energy levels for CSD, lenght = _nEnergies */
Vector _angle; /*!< @brief angles for SN */
Vector _angle; /*!< @brief angles for SN */
std::vector<Matrix> _sigmaSE; /*!< @brief scattering cross section for all energies*/
Vector _sigmaTE; /*!< @brief total cross section for all energies*/
......
......@@ -11,8 +11,7 @@ class CSDSNSolverFP : public SNSolver
std::vector<double> _dose; /*!< @brief TODO */
// Physics acess
Vector _energies; /*!< @brief energy levels for CSD, length = _nEnergies */
Vector _angle; /*!< @brief angles for SN */
Vector _angle; /*!< @brief angles for SN */
std::vector<Matrix> _sigmaSE; /*!< @brief scattering cross section for all energies*/
Vector _sigmaTE; /*!< @brief total cross section for all energies*/
......
......@@ -11,8 +11,7 @@ class CSDSNSolverNoTrafo : public SNSolver
std::vector<double> _dose; /*!< @brief TODO */
// Physics acess
Vector _energies; /*!< @brief energy levels for CSD, lenght = _nEnergies */
Vector _angle; /*!< @brief angles for SN */
Vector _angle; /*!< @brief angles for SN */
std::vector<Matrix> _sigmaSE; /*!< @brief scattering cross section for all energies*/
Vector _sigmaTE; /*!< @brief total cross section for all energies*/
......
......@@ -11,8 +11,7 @@ class CSDSolverTrafoFP : public SNSolver
std::vector<double> _dose; /*!< @brief TODO */
// Physics acess
Vector _energies; /*!< @brief energy levels for CSD, lenght = _nEnergies */
Vector _angle; /*!< @brief angles for SN */
Vector _angle; /*!< @brief angles for SN */
std::vector<Matrix> _sigmaSE; /*!< @brief scattering cross section for all energies*/
Vector _sigmaTE; /*!< @brief total cross section for all energies*/
......
......@@ -11,8 +11,7 @@ class CSDSolverTrafoFP2D : public SNSolver
std::vector<double> _dose; /*!< @brief TODO */
// Physics acess
Vector _energies; /*!< @brief energy levels for CSD, lenght = _nEnergies */
Vector _angle; /*!< @brief angles for SN */
Vector _angle; /*!< @brief angles for SN */
std::vector<Matrix> _sigmaSE; /*!< @brief scattering cross section for all energies*/
Vector _sigmaTE; /*!< @brief total cross section for all energies*/
......@@ -20,9 +19,7 @@ class CSDSolverTrafoFP2D : public SNSolver
Matrix _L; /*!< @brief Laplace Beltrami Matrix */
Matrix _IL; /*!< @brief Laplace Beltrami Matrix */
VectorVector _quadPoints;
VectorVector _quadPointsSphere;
Vector _weights;
Vector _mu;
Vector _phi;
Vector _wp;
......
......@@ -22,8 +22,7 @@ class CSDSolverTrafoFPSH2D : public SNSolver
std::vector<double> _dose; /*!< @brief TODO */
// Physics acess
Vector _energies; /*!< @brief energy levels for CSD, lenght = _nEnergies */
Vector _angle; /*!< @brief angles for SN */
Vector _angle; /*!< @brief angles for SN */
std::vector<Matrix> _sigmaSE; /*!< @brief scattering cross section for all energies*/
Vector _sigmaTE; /*!< @brief total cross section for all energies*/
......@@ -35,9 +34,7 @@ class CSDSolverTrafoFPSH2D : public SNSolver
Matrix _S;
Matrix _M;
VectorVector _quadPoints;
VectorVector _quadPointsSphere;
Vector _weights;
Vector _mu;
Vector _phi;
Vector _wp;
......
......@@ -71,14 +71,14 @@ void CSDPNSolver::IterPreprocessing( unsigned /*idx_iter*/ ) {
// TODO
}
void CSDPNSolver::IterPostprocessing( unsigned /*idx_iter*/ ) {
void CSDPNSolver::IterPostprocessing( unsigned idx_iter ) {
// --- Update Solution ---
_sol = _solNew;
// --- Compute Flux for solution and Screen Output ---
ComputeRadFlux();
unsigned n = idx_pseudotime;
unsigned n = idx_iter;
// -- Compute Dose
for( unsigned j = 0; j < _nCells; ++j ) {
if( n > 0 ) {
......
......@@ -119,7 +119,6 @@ SolverBase* SolverBase::Create( Config* settings ) {
void SolverBase::Solve() {
// --- Preprocessing ---
PrepareVolumeOutput();
DrawPreSolverOutput();
......
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