http://camb.readthedocs.org/en/latest/CAMBdemo.html
there is a function in
Code: Select all
camb.model
Code: Select all
set_for_lmax(lmax, max_eta_k=None, lens_potential_accuracy=0, lens_margin=150, k_eta_fac=2.5, lens_k_eta_reference=18000.0)[source]
Code: Select all
import camb
pars.set_for_lmax(512, lens_potential_accuracy=0)
#get dictionary of CAMB power spectra
powers = results.get_cmb_power_spectra(pars)
totCL=powers['total']
print totCL.shape
(551, 4)
That's unexpected. Shouldn't the length of the array for power spectra by 513, as that's the lmax value?
What's going on here? Was that in the original CAMB Fortran code?