x=0.1:0.1:16;

loglog(x, x*12.8e9)
hold on

xx = ones(1, length(x));
belp = xx.*16e9;

loglog(x, belp, '--')
grid on


%% Plot points
gemm_threaded = 2*0.184;
gemm_threaded_tiled = 2*0.153;

plot(gemm_threaded, gemm_threaded*12.8e9, 'ok')
plot(gemm_threaded_tiled, gemm_threaded_tiled*12.8e9, 'or')
legend('DRAM GB/s', 'Peak Flop/s', 'gemm threaded', 'gemm threaded tiled')


title('Roofline Model')
xlabel('Flop/byte')
ylabel('Flop/s')
xlim([0.1 16])
ylim([0 10e10])

hold off