Compute partial derivative of real parts to log10(rho0)

def dre_dlog10rho0(self, pars):
        """Compute partial derivative of real parts to log10(rho0)
        """

        # first call the linear response to set the parameters
        linear_response = self.dre_drho0(pars)
        result = np.log(10) * self.rho0 * linear_response
        return result