The important code is listed here:

cvx_begin sdp

cvx_solver mosek

cvx_precision best

variable G(48,48) hermitian ;

dual variables y{262}

E_ph= 0.5+0.5*(trace(G’*H_1))/4.0/P_pass;

maximize( E_ph)

subject to

G >= 0;

for k=1:262,

trace(G’*matcell{k}) == vecb(k) : y{k};

end

cvx_end

Now,i wang to compute this problem more accurately,But nomattaer how I adjust the options of “cvx_precison”, it always give me a tolrence of “1e-8”