Commit message (Collapse) | Author | Age | Files | Lines | |
---|---|---|---|---|---|
* | clusters: Accept 1D numpy array and calculate values for all gamma's | Aaron LI | 2017-07-23 | 1 | -4/+11 |
| | | | | Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/halo.py: Constrain tau_acc to avoid zero diffusion | Aaron LI | 2017-07-23 | 1 | -0/+7 |
| | | | | | | | | | | Zero or negative diffusion coefficient leads to unstable or wrong results due to numerical algorithm/scheme adopted to solve the Fokker-Planck equation. Also add a NOTE to the FokkerPlanckSolver class. Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | clusters/solver.py: Improve documents and clean up | Aaron LI | 2017-07-22 | 1 | -14/+49 |
| | | | | | | | | TODO: adaptively determine the proper time step instead of using a constant one. Signed-off-by: Aaron LI <aly@aaronly.me> | ||||
* | solver.py: Fix coefficients calculation w.r.t. escape term | Aaron LI | 2017-06-25 | 1 | -5/+6 |
| | |||||
* | solver.py: Add example for TDMAsolver() | Aaron LI | 2017-06-25 | 1 | -0/+16 |
| | |||||
* | clusters/solver.py: Split methods bound_w/Wplus/Wminus | Aaron LI | 2017-06-25 | 1 | -14/+36 |
| | | | | Also change an INFO message to be DEBUG level. | ||||
* | clusters/solver.py: Fix "fix_boundary()" and remove a warning | Aaron LI | 2017-06-25 | 1 | -6/+5 |
| | |||||
* | clusters/solver.py: Support escape term | Aaron LI | 2017-06-25 | 1 | -5/+12 |
| | |||||
* | solver.py: Avoid possible overflow when w is too large | Aaron LI | 2017-06-01 | 1 | -0/+10 |
| | |||||
* | halo.py: Force a minimal value on acceleration coefficient | Aaron LI | 2017-06-01 | 1 | -1/+1 |
| | | | | | | | To avoid the too small (or zero) values for the diffusion coefficient of the Fokker-Planck equation. Also change config "extragalactic/halo/pmax" from 1e4 to 1e5 | ||||
* | clusters/solver.py: Extrapolate the x grid to avoid NaN's | Aaron LI | 2017-06-01 | 1 | -5/+15 |
| | |||||
* | clusters/solver.py: Disable a RuntimeWarning | Aaron LI | 2017-06-01 | 1 | -2/+4 |
| | |||||
* | Add clusters/solver.py to solve the Fokker-Planck equation | Aaron LI | 2017-06-01 | 1 | -0/+276 |
Adopt the finite difference scheme to solve the Fokker-Planck equation, following Park & Petrosian (1996, ApJS, 103, 255). |