From bce68a6e007ddd39f0c87833c9e2de9b2cae029f Mon Sep 17 00:00:00 2001 From: jaimerzp Date: Fri, 15 Nov 2024 17:07:56 +0000 Subject: [PATCH] extrapolate bug --- src/tracers.jl | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/tracers.jl b/src/tracers.jl index abcb98e..8fcf01a 100755 --- a/src/tracers.jl +++ b/src/tracers.jl @@ -90,7 +90,7 @@ WeakLensingTracer(cosmo::Cosmology, z_n, nz; if IA_params != [0.0, 0.0] hz = Hmpc(cosmo, z_w) - As = get_IA(cosmo, z_w,IA_params) + As = get_IA(cosmo, z_w, IA_params) corr = @. As * (nz_w * hz / nz_norm) w_arr = @. w_arr - corr end @@ -160,13 +160,13 @@ function nz_interpolate(z, nz, res; mode="linear") if mode!="none" if mode=="linear" nz_int = linear_interpolation(z, nz; - extrapolation_bc=Line()) + extrapolation_bc=0.0) end if mode=="cubic" dz = mean(z[2:end] - z[1:end-1]) z_range = z[1]:dz:z[end] nz_int = cubic_spline_interpolation(z_range, nz; - extrapolation_bc=Line()) + extrapolation_bc=0.0) end zz_range = range(0.00001, stop=z[end], length=res) nzz = nz_int(zz_range)