aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim/extragalactic/clusters/solver.py
Commit message (Collapse)AuthorAgeFilesLines
* clusters/solver.py: Add FokkerPlanckTests with 3 casesAaron LI2017-10-121-9/+158
| | | | | | | | The FokkerPlanckSolver is validated with all the 3 test cases! References: * Park & Petrosian 1996, ApJS, 103, 255 * Donnert & Brunetti 2014, MNRAS, 443, 3564
* clusters/solver.py: significantly improve fix_boundary()Aaron LI2017-10-121-10/+30
| | | | Also suggest that ``buffer_np`` be specified to 5%-10% of ``x_np``.
* clusters/solver.py: Small tweaksAaron LI2017-10-071-8/+6
|
* clusters/solver.py: Adjust parameter tc,uc ordersAaron LI2017-10-061-10/+7
|
* 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).