First working beyond-limber. Still need a couple of cosmetics

marcobonici · Apr 18, 2021 · 0f5e90b · 0f5e90b
1 parent c3447b1
commit 0f5e90b
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 6 deletions.
diff --git a/src/PowerSpectrum.jl b/src/PowerSpectrum.jl
@@ -73,8 +73,10 @@ function InterpolateAndEvaluatePowerSpectrum(CosmologicalGrid::CosmologicalGrid,
             PowerSpectrum.InterpolatedPowerSpectrumBeyondLimber[iidx, lidx, :] =
             10 .^(InterpPmm.(log10.(
             CosmologicalGrid.KBeyondLimberArray[lidx, :]),
-            (ConvolvedDensity.ZBinArray[iidx+1]+
-            ConvolvedDensity.ZBinArray[iidx])/2))
+            #(ConvolvedDensity.ZBinArray[iidx+1]+
+            #ConvolvedDensity.ZBinArray[iidx])/2)
+            0.001
+            ))
         end
     end
 end

diff --git a/src/TransferFunctions.jl b/src/TransferFunctions.jl
@@ -3,13 +3,14 @@ function EvaluateTransferFunction(CosmologicalGrid::CosmologicalGridStruct,
     ConvolvedDensity::ConvolvedDensityStruct,
     κTransferFunction::κTransferFunctionStruct,
     PowerSpectrum::PowerSpectrumStruct)
-    GrowthFactorZ = Dierckx.Spline1D(CosmologicalGrid.ZArray,
-    PowerSpectrum.GrowthFactor)
+    #GrowthFactorZ = Dierckx.Spline1D(CosmologicalGrid.ZArray,
+    #PowerSpectrum.GrowthFactor)
     for iidx in 1:length(ConvolvedDensity.ZBinArray)-1
         FFTLog = FFTLogStruct(XArray = BackgroundQuantities.rZArray, FXArray =
         κTransferFunction.LensingSourceFunction.SourceFunctionArray[iidx, :] .*
-        PowerSpectrum.GrowthFactor ./ GrowthFactorZ((
-        ConvolvedDensity.ZBinArray[iidx+1]+ConvolvedDensity.ZBinArray[iidx])/2))
+        PowerSpectrum.GrowthFactor)# ./ GrowthFactorZ((
+        #ConvolvedDensity.ZBinArray[iidx+1]+ConvolvedDensity.ZBinArray[iidx])/2))
+        #TODO Understand why the elimination of this term looks to be the fix
         Kl, κTransferFunction.TransferFunctionArray[iidx, :, :] =
         CosmoCentral.EvaluateFFTLog(FFTLog, CosmologicalGrid.MultipolesArray)
         κTransferFunction.TransferFunctionArray[iidx, :, :] =