Does anyone know why? Thanks!

```
using LinearAlgebra, Convex, MosekTools
rho = HermitianSemidefinite(4)
add_constraint!(rho, tr(rho) == 1)
p = satisfy(tr(rho) == 1)
solve!(p, Mosek.Optimizer)
p.status
rho_opt = round.(evaluate(rho))
tr(rho_opt) # should be 1
eigen(rho_opt)
```