VERB4D
|
Calculating the diffusion in 2D using implicit implementation of diagonal terms, explicit implementation of mixed terms. More...
Functions | |
bool | Diffusion_2D_ADI1 (Matrix2D< double > &psd, Matrix2D< double > x, Matrix2D< double > y, int x_size, int y_size, Matrix1D< double > x_LBC, Matrix1D< double > x_UBC, Matrix1D< double > y_LBC, Matrix1D< double > y_UBC, string x_LBC_type, string x_UBC_type, string y_LBC_type, string y_UBC_type, Matrix2D< double > Dxx, Matrix2D< double > Dyy, Matrix2D< double > Dxy, Matrix2D< double > Dyx, Matrix2D< double > G, Matrix2D< double > Sources, Matrix2D< double > Losses, double dt) |
Calculating the diffusion in 2D using implicit implementation of diagonal terms, explicit implementation of mixed terms.
Created on: May 28, 2011 Author: dimath
bool Diffusion_2D_ADI1 | ( | Matrix2D< double > & | psd, |
Matrix2D< double > | x, | ||
Matrix2D< double > | y, | ||
int | x_size, | ||
int | y_size, | ||
Matrix1D< double > | x_LBC, | ||
Matrix1D< double > | x_UBC, | ||
Matrix1D< double > | y_LBC, | ||
Matrix1D< double > | y_UBC, | ||
string | x_LBC_type, | ||
string | x_UBC_type, | ||
string | y_LBC_type, | ||
string | y_UBC_type, | ||
Matrix2D< double > | Dxx, | ||
Matrix2D< double > | Dyy, | ||
Matrix2D< double > | Dxy, | ||
Matrix2D< double > | Dyx, | ||
Matrix2D< double > | G, | ||
Matrix2D< double > | Sources, | ||
Matrix2D< double > | Losses, | ||
double | dt | ||
) |
Create model matrices and solve the system, Implicit implementation of diagonal terms, explicit implementation of mixed terms
Method:
For the x direction and then for the y direction
psd | - phase space density |
x | - one dimensional slice |
y | - one dimensional slice |
x_size | - size of x slice |
y_size | - size of y slice |
x_LBC | - lower boundary condition for param x |
x_UBC | - upper boundary condition for param x |
y_LBC | - lower boundary condition for param y |
y_UBC | - upper boundary condition for param y |
x_LBC_type | - type of lower boundary condition for param x |
x_UBC_type | - type of upper boundary condition for param x |
y_LBC_type | - type of lower boundary condition for param y |
y_UBC_type | - type of upper boundary condition for param y |
Dxx | - 2D Diffusion matrix |
Dyy | - 2D Diffusion matrix |
Dxy | - 2D Diffusion matrix |
Dyx | - 2D Diffusion matrix |
G | - 2D used for Jacobian to normalize matrix |
Sources | - matrix used for Sources |
Losses | - Matrix used for Losses (loss cone) |
dt | - change in time of single time step |
< Rearranged PSD into one vector of unknown variables