Dear Allin,

Running on my Mac Mini (2.5 GHz Intel Core i5, 8 GB 1333 MHz DDR3 / OS X 10.9.3):

gretl versão 1.9.91cvs
Sessão atual: 2014-06-16 19:42

? include matrix_perf.gfn
 /Users/henrique/Library/Application Support/gretl/functions/matrix_perf.gfn
? matrix_perf(1234)
Aviso: foram gerados valores não-finitos
dgemm experiment 1, variant 1, speed in Gflops

         m         n         k   vanilla    veclib
       128       128       128    1,6428    19,713
       128       128       256    1,9294    19,294
       128       128       512    2,0505    18,455
       128       128      1024    2,0972       nan
       128       128      2048    1,9174       nan

result: vanilla dominates for mnk >= 16777216
  veclib dominates for mnk < 16777216

dgemm experiment 1, variant 2, speed in Gflops

         m         n         k   vanilla    veclib
       128       128       128    1,9713    9,8566
       256       256       128    2,0505    18,455
       512       512       128    1,9174    13,422
      1024      1024       128    2,0649    13,422
      2048      2048       128    2,0649    17,896

result: veclib dominates

dgemm experiment 1, variant 3, speed in Gflops

         m         n         k   vanilla    veclib
       128       128       128    1,9713    3,9426
       256       256       256    2,0972    16,777
       512       512       512    1,7896    26,844
      1024      1024      1024    1,7602    30,678
      2048      2048      2048    1,8276    35,791

result: veclib dominates

dgemm experiment 2, variant 1, speed in Gflops

         m         n         k   vanilla    veclib
         8         8         8   0,29257   0,30118
        16         8         8   0,48762   0,56889
        32         8         8   0,68267   0,93091
        64         8         8    1,0240    1,4629
       128         8         8    1,1378    1,7067
       256         8         8    1,2800    2,5600
       512         8         8    1,4629    2,0480
      1024         8         8    1,4629    1,7067
      2048         8         8    1,4629    2,5600
      4096         8         8    1,4629    2,5600

result: veclib dominates

dgemm experiment 2, variant 2, speed in Gflops

         m         n         k   vanilla    veclib
        10         2      1000    1,6260    3,6364
        20         2      1000    1,7857    4,6512
        40         2      1000    1,6667    3,6364
        80         2      1000    1,8519    4,0816
       160         2      1000    1,9802    5,8824
       320         2      1000    1,8182    4,0816
       640         2      1000    1,7391    3,3333
      1280         2      1000    1,6529    3,2787
      2560         2      1000    1,7544    3,7037
      5120         2      1000    1,7544    3,7736

result: veclib dominates

dgemm experiment 2, variant 3, speed in Gflops

         m         n         k   vanilla    veclib
        10        10      1000    1,7391    6,6667
        20        10      1000    1,9048    6,6667
        40        10      1000    1,6000    8,0000
        80        10      1000    1,4815    10,000
       160        10      1000    1,6000    13,333
       320        10      1000    1,6000    13,333

result: veclib dominates

Operating system: Mac OS X (64-bit)
BLAS library: Apple VecLib
Number of processors: 4
OpenMP enabled: no

Performance summary:

vanilla -
  dominates outright in 0 out of 6 tests
  dominates in 1 test(s) for mnk >= 16777216

veclib -
  dominates outright in 5 out of 6 tests
  dominates in 1 test(s) for mnk < 16777216