I tried to modify camb with non-zero dark matter sound speed. for that I was looking for fderivs subroutine in equation.f90 I could not find it ..
does anyone know where I can put non-tribial dark matter sound speed.
say c_s^2 = .01 ( constant) and then c_s^2 = g(k) , function of k