aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim/extragalactic/clusters/solver.py
Commit message (Collapse)AuthorAgeFilesLines
* clusters/solver.py: Stop calc. at "tstop" by adding "tstep" parameterAaron LI2017-10-061-12/+19
|
* clusters/solver.py: Update comments/references/descriptions etc.Aaron LI2017-10-061-25/+31
|
* Improve logging messagesAaron LI2017-08-011-4/+4
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/solver.py: Use multiple data points to fit the power-lawAaron LI2017-07-251-17/+19
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters: Accept 1D numpy array and calculate values for all gamma'sAaron LI2017-07-231-4/+11
| | | | Signed-off-by: Aaron LI <aly@aaronly.me>
* clusters/halo.py: Constrain tau_acc to avoid zero diffusionAaron LI2017-07-231-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 upAaron LI2017-07-221-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 termAaron LI2017-06-251-5/+6
|
* solver.py: Add example for TDMAsolver()Aaron LI2017-06-251-0/+16
|
* clusters/solver.py: Split methods bound_w/Wplus/WminusAaron LI2017-06-251-14/+36
| | | | Also change an INFO message to be DEBUG level.
* clusters/solver.py: Fix "fix_boundary()" and remove a warningAaron LI2017-06-251-6/+5
|
* clusters/solver.py: Support escape termAaron LI2017-06-251-5/+12
|
* solver.py: Avoid possible overflow when w is too largeAaron LI2017-06-011-0/+10
|
* halo.py: Force a minimal value on acceleration coefficientAaron LI2017-06-011-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'sAaron LI2017-06-011-5/+15
|
* clusters/solver.py: Disable a RuntimeWarningAaron LI2017-06-011-2/+4
|
* Add clusters/solver.py to solve the Fokker-Planck equationAaron LI2017-06-011-0/+276
Adopt the finite difference scheme to solve the Fokker-Planck equation, following Park & Petrosian (1996, ApJS, 103, 255).