diff options
| author | Aaron LI <aly@aaronly.me> | 2017-07-23 10:33:05 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2017-07-23 10:33:05 +0800 |
| commit | 7d1b88a07bee258ca3d93cbfc23d7b37f011065b (patch) | |
| tree | 5d791a82d483de1b8867f0ad54f2cfa77ff3d1fe /fg21sim/extragalactic/clusters/halo.py | |
| parent | 3bd9c730e20c9a70d17c1459ddf78ba71ee84f60 (diff) | |
| download | fg21sim-7d1b88a07bee258ca3d93cbfc23d7b37f011065b.tar.bz2 | |
clusters: Accept 1D numpy array and calculate values for all gamma's
Signed-off-by: Aaron LI <aly@aaronly.me>
Diffstat (limited to 'fg21sim/extragalactic/clusters/halo.py')
| -rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 33 |
1 files changed, 16 insertions, 17 deletions
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index 75789b6..ce4bcaf 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -182,7 +182,7 @@ class RadioHalo: Returns ------- - electron_spec : `~numpy.ndarray` + electron_spec : float 1D `~numpy.ndarray` The solved electron spectrum at ``zend``. Unit: [cm^-3] """ @@ -196,8 +196,7 @@ class RadioHalo: tstop = COSMO.age(zend) if n0_e is None: # Accumulated constantly injected electrons until ``tstart``. - n_inj = np.array([self.fp_injection(gm) - for gm in self.gamma]) + n_inj = self.fp_injection(self.gamma) n0_e = n_inj * tstart self.electron_spec = self.fpsolver.solve(u0=n0_e, tstart=tstart, @@ -218,8 +217,8 @@ class RadioHalo: Parameters ---------- - gamma : float - Lorentz factor of electrons + gamma : float, or float 1D `~numpy.ndarray` + Lorentz factors of electrons t : None Currently a constant injection rate is assumed, therefore this parameter is not used. Keep it for the consistency @@ -227,8 +226,8 @@ class RadioHalo: Returns ------- - Qe : float - Current electron injection rate at specified energy (gamma). + Qe : float, or float 1D `~numpy.ndarray` + Current electron injection rate at specified energies (gamma). Unit: [cm^-3 Gyr^-1] References @@ -251,16 +250,16 @@ class RadioHalo: Parameters ---------- - gamma : float - The Lorentz factor of electrons + gamma : float, or float 1D `~numpy.ndarray` + The Lorentz factors of electrons t : float Current (cosmic) time when solving the equation Unit: [Gyr] Returns ------- - diffusion : float - Diffusion coefficient + diffusion : float, or float 1D `~numpy.ndarray` + Diffusion coefficients Unit: [Gyr^-1] References @@ -283,8 +282,8 @@ class RadioHalo: Returns ------- - advection : float - Advection coefficient, describing the energy loss/gain rate. + advection : float, or float 1D `~numpy.ndarray` + Advection coefficients, describing the energy loss/gain rates. Unit: [Gyr^-1] """ advection = (abs(self._loss_ion(gamma, t)) + @@ -394,16 +393,16 @@ class RadioHalo: Parameters ---------- - gamma : float - The Lorentz factor of electrons + gamma : float, or float 1D `~numpy.ndarray` + The Lorentz factors of electrons t : float The cosmic time/age Unit: [Gyr] Returns ------- - loss : float - The energy loss rate + loss : float, or float 1D `~numpy.ndarray` + The energy loss rates Unit: [Gyr^-1] References |