Commit a7458217 authored by pia.stammer's avatar pia.stammer
Browse files

Commented csdsnsolverfp.h and csdsolvertrafofp.h

parent d36b1405
......@@ -12,7 +12,7 @@ class CSDSNSolverFP : public SNSolver
std::vector<double> _dose; /*! @brief: TODO */
// Physics acess
Vector _energies; /*! @brief: energy levels for CSD, lenght = _nEnergies */
Vector _energies; /*! @brief: energy levels for CSD, length = _nEnergies */
Vector _angle; /*! @brief: angles for SN */
std::vector<Matrix> _sigmaSE; /*! @brief scattering cross section for all energies*/
......@@ -21,17 +21,15 @@ class CSDSNSolverFP : public SNSolver
Matrix _L; /*! @brief Laplace Beltrami Matrix */
Matrix _IL; /*! @brief Laplace Beltrami Matrix */
double _alpha;
double _beta;
double _alpha; /*! @brief Coefficient of GFP operators (see Olbrant 2010, eq. (8)*/
double _beta; /*! @brief Coefficient of GFP operators (see Olbrant 2010, eq. (8)*/
Vector _xi1;
Vector _xi2;
Matrix _xi;
Matrix _xi; /*! @brief matrix of transport coefficients */
bool _RT;
bool _RT; /*! @brief radiotherapy application (on/off), if true use crosssections + stopping powers from database */
double _energyMin;
double _energyMax;
double _energyMin; /*! @brief minimal energy in energy grid*/
double _energyMax; /*! @brief maximal energy in energy grid*/
public:
/**
......
......@@ -21,20 +21,18 @@ class CSDSolverTrafoFP : public SNSolver
Matrix _L; /*! @brief Laplace Beltrami Matrix */
Matrix _IL; /*! @brief Laplace Beltrami Matrix */
double _alpha;
double _alpha2;
double _beta;
double _alpha; /*! @brief Coefficient of GFP operators (see Olbrant 2010, Appendix B)*/
double _alpha2; /*! @brief Coefficient of GFP operators (see Olbrant 2010, Appendix B)*/
double _beta; /*! @brief Coefficient of GFP operators (see Olbrant 2010, Appendix B)*/
Vector _xi1;
Vector _xi2;
Matrix _xi;
Matrix _xi; /*! @brief matrix of transport coefficients */
unsigned _FPMethod;
unsigned _FPMethod; /*! @brief Encodes different ways of computing coefficients alpha, alpha2 & beta, _FPMethod == 1, 2 ,3 stand for methods with increasing accuracy (see Olbrant 2010, Appendix B)*/
bool _RT;
bool _RT; /*! @brief radiotherapy application (on/off), if true use crosssections + stopping powers from database */
double _energyMin;
double _energyMax;
double _energyMin; /*! @brief minimal energy in energy grid*/
double _energyMax; /*! @brief maximal energy in energy grid*/
void GenerateEnergyGrid( bool refinement );
......
......@@ -17,9 +17,9 @@ CSDSNSolverFP::CSDSNSolverFP( Config* settings ) : SNSolver( settings ) {
// Set angle and energies
_angle = Vector( _settings->GetNQuadPoints(), 0.0 ); // my
_energies = Vector( _nEnergies, 0.0 ); // equidistant
_energyMin = 1e-4;
//_energyMax = 10.0;
_energyMax = 5.0;
_energyMin = 1e-4;
//_energyMax = 10.0;
_energyMax = 5.0;
// write equidistant energy grid
_dE = ComputeTimeStep( settings->GetCFL() );
......@@ -62,8 +62,6 @@ CSDSNSolverFP::CSDSNSolverFP( Config* settings ) : SNSolver( settings ) {
double g = 0.8;
// determine momente of Heney-Greenstein
_xi1 = Vector(_nEnergies, 1.0 - g); // paper Olbrant, Frank (11)
_xi2 = Vector(_nEnergies, 4.0 / 3.0 - 2.0 * g + 2.0 / 3.0 * g * g);
_xi = Matrix(3,_nEnergies);
_s = Vector( _nEnergies, 1.0 );
......
......@@ -54,7 +54,7 @@ CSDSolverTrafoFP::CSDSolverTrafoFP( Config* settings ) : SNSolver( settings ) {
// Heney-Greenstein parameter
double g = 0.8;
// determine momente of Heney-Greenstein
// determine transport coefficients of Heney-Greenstein
_xi1 = Vector( _nEnergies, 1.0 - g ); // paper Olbrant, Frank (11)
_xi2 = Vector( _nEnergies, 4.0 / 3.0 - 2.0 * g + 2.0 / 3.0 * g * g );
_xi = Matrix( 4, _nEnergies );
......
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