Commit 73005eb0 authored by Médéric Boquien's avatar Médéric Boquien

When calling np.sum(), the axis should be an integer, not a float, otherwise...

When calling np.sum(), the axis should be an integer, not a float, otherwise it triggers a very disgracious Numpy warning. Now it is gone.
parent a258d7c7
......@@ -17,6 +17,7 @@
- Remove the PSEUDO_D4000 filter which was incorrect. (Médéric Boquien)
- Indicate the correct transmission type for the PACS green and red filters. (Médéric Boquien)
- IRAS filters are defined in energy rather than photons. (Médéric Boquien)
- Remove a Numpy warning in the computation of the IGM absorption (Médéric Boquien)
### Optimised
- By default the MKL library created many threads for each for the parallel processes. Not only was this not necessary as a high-level parallelisation already exists, but it generated a strong oversubscription on the CPU and on the RAM. The slowdown was over a factor of ~2 in some cases. Now we mandate KML to use only 1 thread fo each process. (Médéric Boquien & Yannick Roehlly)
......
......@@ -119,7 +119,7 @@ def igm_transmission(wavelength, redshift):
tau_l_lls = np.zeros_like(wavelength)
tau_l_lls[w] = n0 * ((term1 - term2) * term3 - term4)
tau_taun = np.sum(tau_n[2:n_transitions_max, :], axis=0.)
tau_taun = np.sum(tau_n[2:n_transitions_max, :], axis=0)
lambda_min_igm = (1+redshift)*70.
w = np.where(wavelength < lambda_min_igm)
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment