If you have Mosek 9.x available, i recommend you use that (latest version of Mosek 9.2.x) as solver. If not, install CVXQUAD’s exponential.m replacement, and follow the instructions at CVXQUAD: How to use CVXQUAD's Pade Approximant instead of CVX's unreliable Successive Approximation for GP mode, log, exp, entr, rel_entr, kl_div, log_det, det_rootn, exponential cone. CVXQUAD's Quantum (Matrix) Entropy & Matrix Log related functions .
After doing that, if the result is still infeasible, follow the guidance at https://yalmip.github.io/debugginginfeasible/ , all of which also applies to CVX, excpet for section 1.