In this work, we adopt the fitting formula proposed in [1]}, [2]} that approximates the full transfer function as the sum of the
baryon and cold DM contribution on all scales
\(T(\mathrm {k})=\mathrm {\frac{\Omega _{b,0}}{\Omega _{M,0}}}T_{b}(\mathrm {k})+\mathrm {\frac{\Omega _{DM,0}}{\Omega _{M,0}}}T_{DM}(\mathrm {k}).\) 
