aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim/extragalactic/clusters/solver.py
diff options
context:
space:
mode:
Diffstat (limited to 'fg21sim/extragalactic/clusters/solver.py')
-rw-r--r--fg21sim/extragalactic/clusters/solver.py7
1 files changed, 7 insertions, 0 deletions
diff --git a/fg21sim/extragalactic/clusters/solver.py b/fg21sim/extragalactic/clusters/solver.py
index 6c0e653..dea7f3f 100644
--- a/fg21sim/extragalactic/clusters/solver.py
+++ b/fg21sim/extragalactic/clusters/solver.py
@@ -102,6 +102,13 @@ class FokkerPlanckSolver:
by extrapolating an power law to avoid unphysical accumulations.
This fix is ignored if this parameter is not specified.
+ NOTE
+ ----
+ The diffusion coefficients (i.e., calculated by ``f_diffusion()``)
+ should be *positive* (i.e., C(x) > 0), otherwise unstable or wrong
+ results may occur, due to the current numerical scheme/algorithm
+ adopted.
+
References
----------
[1] Park & Petrosian 1996, ApJS, 103, 255