The additive 2nd term inside trace_inv has incompatible dimensions. For instance H' is 4 by 2 , and is multiplied by Rc, which is 4 by 4. So MATLAB (not even CVX) produces an error message, as it should, just based on double precision data. So I don’t know how you even got anything to run. Perhaps the inputs you ran are different than what’s in your post?
Also, your output shows CVX’s (unreliable) Successive Approximation method was used. You should use Mosek as solver if available, in which case Mosek’s native exponential cone capability will be used. if Mosek is not available to you, 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 .