Speed up cold_gas and hsml by strategic use of numexpr, being
careful with casting: numexpr will cast float_array*double_scalar to
double, unlike numpy, which will leave it as float.