skyscraper_plot.m

d = 4;

% choi_ground     = randomCPTP_quasi_pure(d,0.9);
choi_ground     = randomCPTP(d,d*d);
choi_ground_vec = reshape(choi_ground,[],1);

A = PM_minimal(d);

p               = real(A*choi_ground_vec);
n               = p; % noiseless scenario

[choi_ml_vec, solution, costs] = gdapB(A,n);

choi_ml = reshape(choi_ml_vec,[],d*d);

trace_dist(choi_ground/trace(choi_ground),choi_ml/trace(choi_ml))

figure('Position',[1 0 300 300]);
% subplot(2,1,1)

b1 = bar3(real(choi_ml));

for k = 1:length(b1)
    zdata = b1(k).ZData;
    b1(k).CData = zdata;
    b1(k).FaceColor = 'interp';
end


xticks([1,d*d/2,d*d])
yticks([1,d*d/2,d*d])

xticklabels([])
yticklabels([])

ylim([0,d*d+0.5])
xlim([0,d*d+0.5])
zlim([-1,1])
zlabel('Re$C_{\mathcal{E}}$','Interpreter','latex')
set(gca,'fontsize',18)
saveas(gcf,'./plots/FULLRANKssREAL.png')
figure('Position',[1 0 300 300]);
b2 = bar3(imag(choi_ml))

for k = 1:length(b2)
    zdata = b2(k).ZData;
    b2(k).CData = zdata;
    b2(k).FaceColor = 'interp';
end

% xlabel('i')
% ylabel('j')
xticks([1,d*d/2,d*d])
yticks([1,d*d/2,d*d])
xticklabels([])
yticklabels([])
ylim([0,d*d+0.5])
xlim([0,d*d+0.5])
zlim([-1,1])
zlabel('Im$C_{\mathcal{E}}$','Interpreter','latex')

set(gca,'fontsize',18)
saveas(gcf,'./plots/FULLRANKssIMAG.png')

% saveas(gcf,'./plots/ss.eps','epsc')