From c010a5eecc06999b130d3ecc15f58b32153f2df1 Mon Sep 17 00:00:00 2001 From: Oliver Hahn Date: Wed, 15 Feb 2023 09:55:02 -0800 Subject: [PATCH] added output of 0 for delta_bc and theta_bc --- src/plugins/transfer_eisenstein.cc | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/src/plugins/transfer_eisenstein.cc b/src/plugins/transfer_eisenstein.cc index bec6cf1..6b37600 100644 --- a/src/plugins/transfer_eisenstein.cc +++ b/src/plugins/transfer_eisenstein.cc @@ -213,6 +213,7 @@ public: //! Computes the transfer function for k in Mpc/h by calling TFfit_onek inline double compute(double k, tf_type type) const { + if( type == theta_bc || type == delta_bc ) return 0.0; return etf_.at_k(k); } @@ -305,6 +306,7 @@ public: inline double compute(double k, tf_type type) const { + if( type == theta_bc || type == delta_bc ) return 0.0; return etf_.at_k(k) * pow(1.0 + pow(m_WDMalpha * k, 2.0 * wdmnu_), -5.0 / wdmnu_); } @@ -354,6 +356,8 @@ public: double kkfs2 = kkfs * kkfs; double kkd2 = (k / kd_) * (k / kd_); + if( type == theta_bc || type == delta_bc ) return 0.0; + // in principle the Green et al. (2004) works only up to k/k_fs < 1 // the fit crosses zero at (k/k_fs)**2 = 3/2, we just zero it there... if (kkfs2 < 1.5)