Here are my results (Fedora Core 19 x64):
first the shell output:
load_function_package_from_file:
'/home/helio/gretl/functions/matrix_perf.gfn' is already loaded
detect blas: confused, found too many libs!
-----
linux-vdso.so.1 => (0x00007fffb9bfe000)
libgtksourceview-3.0.so.1 => /lib64/libgtksourceview-3.0.so.1 (0x00000033d9200000)
libgtk-3.so.0 => /lib64/libgtk-3.so.0 (0x00000033d7800000)
libgdk-3.so.0 => /lib64/libgdk-3.so.0 (0x00000033d7400000)
libatk-1.0.so.0 => /lib64/libatk-1.0.so.0 (0x00000035a7e00000)
libgio-2.0.so.0 => /lib64/libgio-2.0.so.0 (0x00000035a6600000)
libpangocairo-1.0.so.0 => /lib64/libpangocairo-1.0.so.0 (0x00000033d8600000)
libcairo-gobject.so.2 => /lib64/libcairo-gobject.so.2 (0x00000033d7000000)
libpango-1.0.so.0 => /lib64/libpango-1.0.so.0 (0x0000003098c00000)
libcairo.so.2 => /lib64/libcairo.so.2 (0x00000033d8a00000)
libgretl-1.0.so.10 => /usr/local/lib/libgretl-1.0.so.10 (0x00007f87cce3e000)
liblapack.so.3 => /usr/lib64/atlas/liblapack.so.3 (0x00007f87cc5da000)
libblas.so.3 => /lib64/libblas.so.3 (0x000000316e600000)
libgfortran.so.3 => /lib64/libgfortran.so.3 (0x000000327f600000)
libm.so.6 => /lib64/libm.so.6 (0x000000327d600000)
libdl.so.2 => /lib64/libdl.so.2 (0x000000327ca00000)
libz.so.1 => /lib64/libz.so.1 (0x000000327d200000)
libxml2.so.2 => /lib64/libxml2.so.2 (0x0000003283200000)
libgmp.so.10 => /lib64/libgmp.so.10 (0x0000003e4b000000)
libfftw3.so.3 => /lib64/libfftw3.so.3 (0x00000036ca600000)
libcurl.so.4 => /lib64/libcurl.so.4 (0x00000036be600000)
libgdk_pixbuf-2.0.so.0 => /lib64/libgdk_pixbuf-2.0.so.0 (0x00000033d8200000)
libgobject-2.0.so.0 => /lib64/libgobject-2.0.so.0 (0x00000035a5e00000)
libglib-2.0.so.0 => /lib64/libglib-2.0.so.0 (0x00000035a5a00000)
libgomp.so.1 => /lib64/libgomp.so.1 (0x0000003285e00000)
libpthread.so.0 => /lib64/libpthread.so.0 (0x000000327ce00000)
libc.so.6 => /lib64/libc.so.6 (0x000000327c600000)
libgmodule-2.0.so.0 => /lib64/libgmodule-2.0.so.0 (0x00000035a6200000)
libX11.so.6 => /lib64/libX11.so.6 (0x0000003280a00000)
libXi.so.6 => /lib64/libXi.so.6 (0x0000003284600000)
libXfixes.so.3 => /lib64/libXfixes.so.3 (0x0000003284e00000)
libatk-bridge-2.0.so.0 => /lib64/libatk-bridge-2.0.so.0 (0x00000035a9e00000)
libpangoft2-1.0.so.0 => /lib64/libpangoft2-1.0.so.0 (0x0000003099000000)
libfontconfig.so.1 => /lib64/libfontconfig.so.1 (0x0000003284200000)
libXinerama.so.1 => /lib64/libXinerama.so.1 (0x0000003282a00000)
libXrandr.so.2 => /lib64/libXrandr.so.2 (0x0000003282e00000)
libXcursor.so.1 => /lib64/libXcursor.so.1 (0x0000003286200000)
libXcomposite.so.1 => /lib64/libXcomposite.so.1 (0x000000328b200000)
libXdamage.so.1 => /lib64/libXdamage.so.1 (0x0000003289a00000)
libwayland-client.so.0 => /lib64/libwayland-client.so.0 (0x000000328a600000)
libxkbcommon.so.0 => /lib64/libxkbcommon.so.0 (0x000000328ce00000)
libwayland-cursor.so.0 => /lib64/libwayland-cursor.so.0 (0x000000328e200000)
libXext.so.6 => /lib64/libXext.so.6 (0x0000003280e00000)
libffi.so.6 => /lib64/libffi.so.6 (0x000000327f200000)
libselinux.so.1 => /lib64/libselinux.so.1 (0x000000337d200000)
libresolv.so.2 => /lib64/libresolv.so.2 (0x000000327ee00000)
libgthread-2.0.so.0 => /lib64/libgthread-2.0.so.0 (0x00000035a6a00000)
libharfbuzz.so.0 => /lib64/libharfbuzz.so.0 (0x0000003179a00000)
libfreetype.so.6 => /lib64/libfreetype.so.6 (0x0000003281e00000)
libpixman-1.so.0 => /lib64/libpixman-1.so.0 (0x0000003287200000)
libEGL.so.1 => /lib64/libEGL.so.1 (0x000000337e600000)
libpng15.so.15 => /lib64/libpng15.so.15 (0x00000033d6c00000)
libxcb-shm.so.0 => /lib64/libxcb-shm.so.0 (0x0000003288600000)
libxcb-render.so.0 => /lib64/libxcb-render.so.0 (0x0000003288200000)
libxcb.so.1 => /lib64/libxcb.so.1 (0x0000003280600000)
libXrender.so.1 => /lib64/libXrender.so.1 (0x0000003282200000)
libGL.so.1 => /lib64/libGL.so.1 (0x000000337de00000)
librt.so.1 => /lib64/librt.so.1 (0x000000327da00000)
/lib64/ld-linux-x86-64.so.2 (0x000000327c200000)
libf77blas.so.3 => /usr/lib64/atlas/libf77blas.so.3 (0x00007f87cc3ac000)
libcblas.so.3 => /usr/lib64/atlas/libcblas.so.3 (0x000000364d000000)
libquadmath.so.0 => /lib64/libquadmath.so.0 (0x00007f87cc16f000)
libgcc_s.so.1 => /lib64/libgcc_s.so.1 (0x000000327ea00000)
liblzma.so.5 => /lib64/liblzma.so.5 (0x0000003281600000)
libidn.so.11 => /lib64/libidn.so.11 (0x000000329d000000)
libssh2.so.1 => /lib64/libssh2.so.1 (0x00000036be200000)
libssl3.so => /lib64/libssl3.so (0x0000003c22c00000)
libsmime3.so => /lib64/libsmime3.so (0x0000003c22800000)
libnss3.so => /lib64/libnss3.so (0x0000003c22000000)
libnssutil3.so => /lib64/libnssutil3.so (0x0000003c21c00000)
libplds4.so => /lib64/libplds4.so (0x0000003c21400000)
libplc4.so => /lib64/libplc4.so (0x0000003c21000000)
libnspr4.so => /lib64/libnspr4.so (0x0000003c21800000)
libgssapi_krb5.so.2 => /lib64/libgssapi_krb5.so.2 (0x00000036bd200000)
libkrb5.so.3 => /lib64/libkrb5.so.3 (0x00000036bca00000)
libk5crypto.so.3 => /lib64/libk5crypto.so.3 (0x00000036bda00000)
libcom_err.so.2 => /lib64/libcom_err.so.2 (0x0000003290200000)
liblber-2.4.so.2 => /lib64/liblber-2.4.so.2 (0x000000398c800000)
libldap-2.4.so.2 => /lib64/libldap-2.4.so.2 (0x0000003c24000000)
libatspi.so.0 => /lib64/libatspi.so.0 (0x00000035a9600000)
libdbus-1.so.3 => /lib64/libdbus-1.so.3 (0x0000003281200000)
libexpat.so.1 => /lib64/libexpat.so.1 (0x0000003283600000)
libpcre.so.1 => /lib64/libpcre.so.1 (0x000000337ce00000)
libgraphite2.so.3 => /lib64/libgraphite2.so.3 (0x000000328aa00000)
libX11-xcb.so.1 => /lib64/libX11-xcb.so.1 (0x0000003285600000)
libxcb-dri2.so.0 => /lib64/libxcb-dri2.so.0 (0x0000003289600000)
libxcb-xfixes.so.0 => /lib64/libxcb-xfixes.so.0 (0x0000003286600000)
libxcb-shape.so.0 => /lib64/libxcb-shape.so.0 (0x000000328be00000)
libgbm.so.1 => /lib64/libgbm.so.1 (0x000000337e200000)
libwayland-server.so.0 => /lib64/libwayland-server.so.0 (0x0000003287e00000)
libglapi.so.0 => /lib64/libglapi.so.0 (0x000000337ee00000)
libudev.so.1 => /lib64/libudev.so.1 (0x0000003281a00000)
libdrm.so.2 => /lib64/libdrm.so.2 (0x0000003286e00000)
libXau.so.6 => /lib64/libXau.so.6 (0x0000003280200000)
libxcb-glx.so.0 => /lib64/libxcb-glx.so.0 (0x0000003287600000)
libXxf86vm.so.1 => /lib64/libXxf86vm.so.1 (0x0000003283e00000)
libatlas.so.3 => /usr/lib64/atlas/libatlas.so.3 (0x000000364c200000)
libssl.so.10 => /lib64/libssl.so.10 (0x00000036bde00000)
libcrypto.so.10 => /lib64/libcrypto.so.10 (0x00000039a7a00000)
libkrb5support.so.0 => /lib64/libkrb5support.so.0 (0x00000036bd600000)
libkeyutils.so.1 => /lib64/libkeyutils.so.1 (0x00000036bce00000)
libsasl2.so.3 => /lib64/libsasl2.so.3 (0x0000003c23c00000)
libcrypt.so.1 => /lib64/libcrypt.so.1 (0x0000003c22400000)
libfreebl3.so => /lib64/libfreebl3.so (0x0000003c20c00000)
----
gretl versão 1.9.91cvs
Sessão atual: 2014-06-16 22:25
? include matrix_perf.gfn
/home/helio/gretl/functions/matrix_perf.gfn
? matrix_perf(1234)
dgemm experiment 1, variant 1, speed in Gflops
m n k vanilla openmp sysblas
128 128 128 0,76280 3,8224 1,2408
128 128 256 1,4744 4,4971 2,3319
128 128 512 1,5409 4,6394 2,2827
128 128 1024 1,6307 4,7384 2,4370
128 128 2048 1,6366 4,6429 2,3655
result: openmp dominates
dgemm experiment 1, variant 2, speed in Gflops
m n k vanilla openmp sysblas
128 128 128 1,2369 4,2107 2,3447
256 256 128 1,7027 4,2914 2,4852
512 512 128 1,7378 4,3442 2,4599
1024 1024 128 1,7530 4,5483 2,3357
2048 2048 128 1,6727 4,3896 2,4288
result: openmp dominates
dgemm experiment 1, variant 3, speed in Gflops
m n k vanilla openmp sysblas
128 128 128 1,3331 1,6034 1,0660
256 256 256 1,7022 4,6465 2,2175
512 512 512 1,7230 3,6227 2,6052
1024 1024 1024 1,5251 4,3004 2,0975
2048 2048 2048 1,5605 4,1371 2,0563
result: openmp dominates
dgemm experiment 2, variant 1, speed in Gflops
m n k vanilla openmp sysblas
8 8 8 0,54613 0,19646 0,48416
16 8 8 0,81994 0,38813 0,68170
32 8 8 1,0980 0,62199 0,85589
64 8 8 1,3154 0,82638 1,0962
128 8 8 1,4332 1,0507 1,2646
256 8 8 1,4836 1,2024 1,4424
512 8 8 1,5283 1,2608 1,5017
1024 8 8 1,5241 1,3197 1,4720
2048 8 8 1,5281 1,3507 1,5226
4096 8 8 1,5382 1,3684 1,5155
result: vanilla dominates
dgemm experiment 2, variant 2, speed in Gflops
m n k vanilla openmp sysblas
10 2 1000 1,2576 1,1964 1,6339
20 2 1000 1,3761 1,9025 1,8154
40 2 1000 1,2545 1,9819 1,8016
80 2 1000 1,6691 2,2706 2,1917
160 2 1000 1,5171 2,3899 2,4151
320 2 1000 1,6823 2,5837 2,3722
640 2 1000 1,5027 2,8140 2,1026
1280 2 1000 1,4856 2,8152 2,1165
2560 2 1000 1,5764 2,3430 2,0931
5120 2 1000 1,5695 3,4026 2,1107
result: openmp dominates for mnk >= 640000
dgemm experiment 2, variant 3, speed in Gflops
m n k vanilla openmp sysblas
10 10 1000 1,1356 2,6000 1,7212
20 10 1000 1,3239 3,8221 1,5229
40 10 1000 1,2920 3,7631 1,8438
80 10 1000 1,4216 4,1526 2,2323
160 10 1000 1,6549 4,3422 2,4932
320 10 1000 1,6679 4,3360 2,4248
result: openmp dominates
Operating system: Linux (64-bit)
BLAS library: sysblas
Number of processors: 4
OpenMP enabled: yes
Performance summary:
vanilla -
dominates outright in 1 out of 6 tests
openmp -
dominates outright in 4 out of 6 tests
dominates in 1 test(s) for mnk >= 640000
sysblas -
dominates outright in 0 out of 6 tests