1axes3 = subplot(4, 4, 3);
2hold(axes3,'on');
3
4axes7 = subplot(4, 4, 7);
5hold(axes7,'on');
6
7axes11 = subplot(4, 4, 11);
8hold(axes11,'on');
9
10axes15 = subplot(4, 4, 15);
11hold(axes15,'on');
12
13addpath(pathname_blis)
14
15if(plot_s)
16    axes(axes3);
17    output_mt_shemm_asm_blis
18    plot(data_mt_shemm_asm_blis(:,1), data_mt_shemm_asm_blis(:,3), 'LineWidth', 1.25,'Color', [0 0 1]);
19end
20
21if(plot_d)
22    axes(axes7);
23    output_mt_dhemm_asm_blis
24    plot(data_mt_dhemm_asm_blis(:,1), data_mt_dhemm_asm_blis(:,3), 'LineWidth', 1.25,'Color', [0 0 1]);
25end
26
27if(plot_c)
28    axes(axes11);
29    output_mt_chemm_1m_blis
30    plot(data_mt_chemm_1m_blis(:,1), data_mt_chemm_1m_blis(:,3), 'LineWidth', 1.25,'Color', [0 0 1]);
31end
32
33if(plot_z)
34    axes(axes15);
35    output_mt_zhemm_1m_blis
36    plot(data_mt_zhemm_1m_blis(:,1), data_mt_zhemm_1m_blis(:,3), 'LineWidth', 1.25,'Color', [0 0 1]);
37end
38
39clear *hemm*
40rmpath(pathname_blis)
41
42addpath(pathname_openblas)
43
44if(plot_s)
45    axes(axes3);
46    output_mt_shemm_openblas
47    plot(data_mt_shemm_openblas(:,1), data_mt_shemm_openblas(:,3), 'LineWidth', 1.25,'Color', [0 1 0]);
48end
49
50if(plot_d)
51    axes(axes7);
52    output_mt_dhemm_openblas
53    plot(data_mt_dhemm_openblas(:,1), data_mt_dhemm_openblas(:,3), 'LineWidth', 1.25,'Color', [0 1 0]);
54end
55
56if(plot_c)
57    axes(axes11);
58    output_mt_chemm_openblas
59    plot(data_mt_chemm_openblas(:,1), data_mt_chemm_openblas(:,3), 'LineWidth', 1.25,'Color', [0 1 0]);
60end
61
62if(plot_z)
63    axes(axes15);
64    output_mt_zhemm_openblas
65    plot(data_mt_zhemm_openblas(:,1), data_mt_zhemm_openblas(:,3), 'LineWidth', 1.25,'Color', [0 1 0]);
66end
67
68clear *hemm*
69rmpath(pathname_openblas)
70
71addpath(pathname_armpl);
72
73if(plot_s)
74    axes(axes3);
75    output_mt_shemm_armpl
76    plot(data_mt_shemm_armpl(:,1), data_mt_shemm_armpl(:,3), '--', 'LineWidth', 1.25,'Color', [1 0 1]);
77end
78
79if(plot_d)
80    axes(axes7);
81    output_mt_dhemm_armpl
82    plot(data_mt_dhemm_armpl(:,1), data_mt_dhemm_armpl(:,3), '--', 'LineWidth', 1.25,'Color', [1 0 1]);
83end
84
85if(plot_c)
86    axes(axes11);
87    output_mt_chemm_armpl
88    plot(data_mt_chemm_armpl(:,1), data_mt_chemm_armpl(:,3), '--', 'LineWidth', 1.25,'Color', [1 0 1]);
89end
90
91if(plot_z)
92    axes(axes15);
93    output_mt_zhemm_armpl
94    plot(data_mt_zhemm_armpl(:,1), data_mt_zhemm_armpl(:,3), '--', 'LineWidth', 1.25,'Color', [1 0 1]);
95end
96
97clear *hemm*
98rmpath(pathname_armpl)
99
100% SSYMM multi threaded
101
102axes(axes3);
103ylabel( 'GFLOPS', 'FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue');
104%xlabel( 'matrix dimension m=n=k', 'FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue' );
105title('SSYMM (multi-threaded)','FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue');
106box(axes3,'on');
107set(axes3,'FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue');
108v = axis;     % extract the current ranges
109axis( [ 0 xmax_mt 0 speak*numcores ] )
110
111
112axes(axes7);
113ylabel( 'GFLOPS', 'FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue');
114%xlabel( 'matrix dimension m=n=k', 'FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue' );
115title('DSYMM (multi-threaded)','FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue');
116box(axes7,'on');
117set(axes7,'FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue');
118%legend({'BLIS', 'OpenBLAS', 'ARMPL'},'FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue', 'Location', 'best');
119v = axis;     % extract the current ranges
120axis( [ 0 xmax_mt 0 dpeak*numcores ] )
121
122% CHEMM multi threaded
123
124axes(axes11);
125ylabel( 'GFLOPS', 'FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue');
126%xlabel( 'matrix dimension m=n=k', 'FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue' );
127title('CHEMM (multi-threaded)','FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue');
128box(axes11,'on');
129set(axes11,'FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue');
130v = axis;     % extract the current ranges
131axis( [ 0 xmax_mt 0 speak*numcores ] )
132
133
134% ZHEMM multi threaded
135axes(axes15);
136ylabel( 'GFLOPS', 'FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue');
137xlabel( 'matrix dimension m=n=k', 'FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue' );
138title('ZHEMM (multi-threaded)','FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue');
139box(axes15,'on');
140set(axes15,'FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue');
141% legend({'BLIS', 'BLIS (AVX2)', 'OpenBLAS', 'MKL'},'FontSize', fontsize, 'FontWeight', 'bold', 'FontName', 'Helvetica Neue', 'Location', 'South');
142v = axis;     % extract the current ranges
143axis( [ 0 xmax_mt 0 dpeak*numcores ] )
144
145
146
147