!----------------------------------------------------------------------------- ! ! Compute sigma^2(R) from the linear matter power spectrum, and ! fit it by Chebyshev polynomials (single precision) ! ! << sample code >> ! ! USE linearpk ! USE sigma ! REAL :: chebev ! REAL :: lnsigma2,lnRh,Rh,dlnsigma2dlnRh ! character(len=128) :: filename ! integer :: n ! filename='wmap5baosn_max_likelihood_matterpower.dat' ! n=896 ! # of lines in the file ! CALL get_linearpk(filename,n) ! CALL compute_sigma2 ! Rh=8. ! h^-1 Mpc ! lnRh = alog(Rh) ! if(lnRh>lnR2.or.lnRh