I am currently working on optimizing an objective function which contains quantum relative entropy (sometimes referred to matrix relative entropy) and is described in https://en.wikipedia.org/wiki/Quantum_relative_entropy.

The above may be written in MATLAB code as `trace(rho*logm(rho) - rho*logm(zrho))`

where `rho`

and `zrho`

are density matrices.

Quantum relative entropy is known to be convex, but it seems that CVX is currently unable to support such a function. The question I would like to pose is how would one implement quantum relative entropy in CVX?