We propose a new way to treat the exponential/relative entropy cone using symmetric cone solvers. Our approach is based on highly accurate rational (Padé) approximations of the logarithm function. The key to this approach is that our rational approximations, by construction, inherit the (operator) concavity of the logarithm. Importantly, our method extends to the matrix logarithm and other derived functions such as the matrix relative entropy, giving new semidefinite optimization-based tools for convex optimization involving these functions. We include an implementation of our method for the MATLAB-based parser CVX. We compare our method to existing approximation schemes, and show that it can be much faster, especially for large problems. Preprint at https://arxiv.org/abs/1705.00812. Joint work with Hamza Fawzi (Cambridge) and James Saunderson (Monash).