% Plot confusion matrix (true vs inferred errors) as heatmap.
% Saves PNG in images/.

thisdir = fileparts(mfilename('fullpath'));
outdir  = fullfile(thisdir,'..','images');
if ~exist(outdir,'dir'), mkdir(outdir);
end

rng(1);
thisdir = fileparts(mfilename('fullpath'));
addpath(genpath(fullfile(thisdir,'..','src')));
warning('off','all');
p = 0.12;
N = 2000;
C = confusion_bitflip(N, p, 0);

figure('visible','off');
imagesc(C);
axis image;
colorbar;
set(gca,'XTick',1:4,'XTickLabel',{'I','X1','X2','X3'});
set(gca,'YTick',1:5,'YTickLabel',{'I','X1','X2','X3','multi'});
xlabel('inferred');
ylabel('true');
title(sprintf('Bit-flip syndrome decoder confusion matrix, p=%.2f, N=%d', p, N));

if ~exist(outdir,'dir'), mkdir(outdir);
end
print(fullfile(outdir,'bitflip_decoder_confusion_matrix.png'),'-dpng');