11 #ifndef PhaseSpaceDensity_H
12 #define PhaseSpaceDensity_H
17 #include "../Grid/Grid.h"
18 #include "../Matrix/Matrix.h"
19 #include "../DiffusionCoefficient/DiffusionCoefficient.h"
20 #include "../Parameters/Parameters.h"
21 #include "MatrixSolver.h"
46 void Load_initial_f(
GridElement &L,
GridElement &pc,
GridElement &alpha,
double tau,
double Kp,
double min_psd = 1.e-99,
string J_L7_function =
"J_L7",
double fb_out = 1,
double fb_in = 0);
48 void Load_initial_f_two_zone(
GridElement &L,
GridElement &pc,
GridElement &alpha,
double tau,
double tauLpp,
double Kp,
double min_psd = 1.e-99,
string J_L7_function =
"J_L7",
double fb_out = 1,
double fb_in = 0);
56 void DiffusionMixTermExplicit(
double dt,
double Lpp,
63 string pc_lowerBoundaryCondition_calculationType,
64 string pc_upperBoundaryCondition_calculationType,
65 string alpha_lowerBoundaryCondition_calculationType,
66 string alpha_upperBoundaryCondition_calculationType);
69 void Diffusion_alpha(
double dt,
double Lpp,
74 string alpha_lowerBoundaryCondition_calculationType,
75 string alpha_upperBoundaryCondition_calculationType);
77 void Diffusion_pc(
double dt,
double Lpp,
82 string pc_lowerBoundaryCondition_calculationType,
83 string pc_upperBoundaryCondition_calculationType);
85 void Diffusion_L(
double dt,
double Lpp,
90 string lowerBoundaryCondition_calculationType,
91 string upperBoundaryCondition_calculationType,
92 double tau,
double tauLpp);
95 void Diffusion_pc_alpha(
double dt,
double Lpp,
104 string pc_lowerBoundaryCondition_calculationType,
105 string pc_upperBoundaryCondition_calculationType,
106 string alpha_lowerBoundaryCondition_calculationType,
107 string alpha_upperBoundaryCondition_calculationType);
109 void Diffusion_pc_alpha_KC(
double dt,
double Lpp,
118 string pc_lowerBoundaryCondition_calculationType,
119 string pc_upperBoundaryCondition_calculationType,
120 string alpha_lowerBoundaryCondition_calculationType,
121 string alpha_upperBoundaryCondition_calculationType);
130 void SourcesAndLosses(
135 double Lpp,
double tau,
double tauLpp,
double Kp);
138 void Output_without_grid(
double time);
139 ofstream output_without_grid_file;
Interpolation parameters structure
Definition: Parameters.h:238
Array of values of coordinate axes.
Definition: Grid.h:28
Holds diffusion coefficient matrix and routines to load and calculate it.
Definition: DiffusionCoefficient.h:23
This matrix calculates the diagonal values and index given parameters for x, y, and z...
Definition: Matrix.h:310
Main parameters structure that holds smaller structures for individual parameters.
Definition: Parameters.h:87
Phase Space Density class.
Definition: PSD.h:48
Holds upper and lower boundary conditions.
Definition: BoundaryConditions.h:26
Computational grid composed of 3 different GridElement.
Definition: Grid.h:53
PSD parameters structure.
Definition: Parameters.h:178