set.seed(4987) library(microbenchmark) library(Matrix) A <- matrix(rnorm(1e6, 1e3, 1e2), ncol = 1e3) B <- matrix(rnorm(1e6, 1e3, 1e2), ncol = 1e3) A <- crossprod(A, A) A <- A * 1000 / mean(A) colnames(A) <- colnames(B) <- NULL options(scipen=4) BLAS <- microbenchmark( sort(c(as.vector(A), as.vector(B))), det(A), A %*% B, t(A) %*% B, crossprod(A, B), solve(A), solve(A, t(B)), solve(B), chol(A), chol(B, pivot = TRUE), qr(A, LAPACK = TRUE), svd(A), eigen(A, symmetric = TRUE), eigen(A, symmetric = FALSE), eigen(B, symmetric = FALSE), lu(A), fft(A), Hilbert(3000), toeplitz(A[1:500, 1]), princomp(A), times=25L, unit='ms', control = list(order = 'block') ) NotBLAS <- microbenchmark( A + 2, A - 2, A * 2, A / 2, A * 0.5, A ^ 2, sqrt(A[1:10000]), sin(A[1:10000]), A + B, A - B, A * B, A / B, A[1:1e5] %% B[1:1e5], A[1:1e5] %/% B[1:1e5], times = 5000L, unit='ms', control = list(order = 'block') ) library(dplyr) summarize(group_by(BLAS, expr), Min = min(time/1e6), LQ = quantile(time/1e6, .25), Median=median(time/1e6), UQ = quantile(time/1e6, .75), Max = max(time/1e6), Mean=mean(time/1e6), SD = sd(time/1e6), CV = SD/Mean, n = n()) summarize(group_by(NotBLAS, expr), Min = min(time/1e6), LQ = quantile(time/1e6, .25), Median=median(time/1e6), UQ = quantile(time/1e6, .75), Max = max(time/1e6), Mean=mean(time/1e6), SD = sd(time/1e6), CV = SD/Mean, n = n()) Test Platform: i7-3740QM @ 2.7Ghz with 8MB RAM; Win7 64 R-devel_2015-09-10 Unit: milliseconds ----------------------- 463-SJLJ-Ref-TG BLAS: expr Min LQ Median UQ Max Mean SD CV n (fctr) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (int) 1 sort(c(as.vector(A), as.vector(B))) 346.670472 411.416125 413.877411 417.024694 491.331660 410.553484 28.352892 0.069060168 25 2 det(A) 219.487502 219.701329 220.149146 224.544383 233.087167 222.413028 3.516384 0.015810154 25 3 A %*% B 667.596926 672.010045 680.406347 686.980939 694.779515 680.316311 8.845379 0.013001862 25 4 t(A) %*% B 677.693186 683.155272 690.410152 696.095576 714.086292 692.379455 10.734744 0.015504135 25 5 crossprod(A, B) 1174.316860 1179.612297 1188.170300 1197.340106 1224.832402 1191.142853 13.777636 0.011566737 25 6 solve(A) 1057.255641 1065.229236 1074.788647 1088.213294 1149.432714 1080.882305 21.247271 0.019657340 25 7 solve(A, t(B)) 1466.743904 1476.934521 1488.949511 1508.065670 1647.901686 1501.146277 40.739752 0.027139096 25 8 solve(B) 1055.135259 1067.052468 1071.265457 1079.953582 1095.306092 1074.366153 11.042234 0.010277906 25 9 chol(A) 198.156984 198.435111 200.645284 202.481072 274.923715 203.694160 14.975968 0.073521832 25 10 chol(B, pivot = TRUE) 3.565423 3.608036 3.617548 3.625537 9.268729 4.435963 1.936535 0.436553454 25 11 qr(A, LAPACK = TRUE) 686.699769 688.311456 694.458776 699.672032 704.654339 694.455397 6.273625 0.009033878 25 12 svd(A) 3523.933341 3552.889321 3570.097782 3590.092836 3656.964578 3577.719134 36.856300 0.010301619 25 13 eigen(A, symmetric = TRUE) 1537.113724 1544.136514 1549.912111 1563.911654 1632.498193 1557.346331 19.538132 0.012545785 25 14 eigen(A, symmetric = FALSE) 5879.346346 5923.193319 5934.092380 5951.791272 6039.831070 5938.558658 30.899473 0.005203194 25 15 eigen(B, symmetric = FALSE) 6599.224214 6688.869231 6701.224747 6724.945022 6780.983907 6703.843732 36.325394 0.005418592 25 16 lu(A) 224.395617 227.229007 231.011681 233.769357 301.897044 240.613051 26.705430 0.110989116 25 17 fft(A) 158.840316 159.487121 162.078530 163.045695 168.644751 161.924910 2.556407 0.015787610 25 18 Hilbert(3000) 175.193091 250.913139 251.308831 252.707073 477.203516 258.187361 52.921503 0.204973251 25 19 toeplitz(A[1:500, 1]) 5.452194 5.653084 5.685044 5.720809 12.367691 6.423225 2.083354 0.324347012 25 20 princomp(A) 2912.846707 2929.681924 2957.535290 2991.234258 3038.586159 2961.546004 35.907469 0.012124569 25 NotBLAS: expr Min LQ Median UQ Max Mean SD CV n (fctr) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (int) 1 A + 2 2.468897 2.601682 2.622228 2.666363 85.564434 3.3402080 2.90901905 0.87090955 5000 2 A - 2 2.474984 2.603585 2.627935 2.696040 82.606628 3.4400976 2.97548597 0.86494232 5000 3 A * 2 2.486399 2.604345 2.625652 2.668361 80.452384 3.3729374 2.95901869 0.87728243 5000 4 A/2 4.319523 4.427197 4.450025 4.545905 83.112279 5.2318349 2.96751643 0.56720376 5000 5 A * 0.5 2.475745 2.600541 2.620135 2.673972 81.964007 3.4033950 2.95926873 0.86950493 5000 6 A^2 2.488301 2.599400 2.617662 2.652761 80.574516 3.3406635 2.94784791 0.88241390 5000 7 sqrt(A[1:10000]) 0.183389 0.184530 0.184911 0.186052 6.568124 0.1923734 0.10229228 0.53173813 5000 8 sin(A[1:10000]) 0.633870 0.635012 0.635772 0.637294 1.828179 0.6412265 0.04838875 0.07546281 5000 9 A + B 1.577827 1.615874 1.643838 1.676369 75.259674 1.8073003 2.09311187 1.15814281 5000 10 A - B 1.575924 1.597992 1.620820 1.654682 73.819199 1.7873970 2.04293905 1.14296880 5000 11 A * B 1.575924 1.601796 1.628429 1.667999 75.471978 1.7970439 2.07679488 1.15567288 5000 12 A/B 5.232660 5.244835 5.268045 5.369250 78.904235 5.4518620 2.05531369 0.37699297 5000 13 A[1:100000]%%B[1:100000] 3.852300 3.886543 3.893392 3.965301 77.128562 3.9986451 1.47445036 0.36873749 5000 14 A[1:100000]%/%B[1:100000] 3.517484 3.532321 3.596622 3.614123 76.134384 3.6596916 1.44511068 0.39487226 5000 ---------------------- 493-SEH-POSIX-Ref-TG BLAS: expr Min LQ Median UQ Max Mean SD CV n (fctr) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (int) 1 sort(c(as.vector(A), as.vector(B))) 338.381844 410.412054 411.837310 415.040141 478.423316 406.999795 28.648678 0.070389907 25 2 det(A) 211.491079 211.808774 213.916981 216.401477 230.526576 214.665343 4.055187 0.018890741 25 3 A %*% B 635.262718 639.341400 644.261310 646.274398 661.466347 644.272359 7.081330 0.010991205 25 4 t(A) %*% B 647.980444 650.184150 653.953887 663.895294 668.152038 656.219762 6.835029 0.010415761 25 5 crossprod(A, B) 1170.194042 1172.998897 1175.121562 1181.986836 1198.825858 1178.317285 8.579228 0.007280915 25 6 solve(A) 1050.296389 1055.599818 1061.713656 1075.142488 1139.137466 1068.768847 20.938019 0.019590783 25 7 solve(A, t(B)) 1464.983449 1475.457521 1483.305179 1492.519120 1566.337937 1487.362661 20.597782 0.013848527 25 8 solve(B) 1061.938516 1067.720581 1076.397292 1081.991782 1091.797360 1076.639790 9.391143 0.008722642 25 9 chol(A) 198.071377 198.686604 202.318990 202.661037 277.394134 204.139787 15.387098 0.075375300 25 10 chol(B, pivot = TRUE) 3.519767 3.635050 3.646464 3.762889 9.120725 4.501496 1.970007 0.437633786 25 11 qr(A, LAPACK = TRUE) 680.090553 682.158050 684.750981 687.508276 697.853366 685.297798 4.329356 0.006317481 25 12 svd(A) 3402.099456 3429.376404 3445.659552 3476.605791 3527.316517 3456.487936 36.443381 0.010543471 25 13 eigen(A, symmetric = TRUE) 1497.170796 1508.003653 1515.501275 1524.413880 1593.282731 1518.044851 18.657349 0.012290380 25 14 eigen(A, symmetric = FALSE) 5577.349618 5601.437812 5622.938782 5650.597726 5710.294868 5627.537161 35.746762 0.006352115 25 15 eigen(B, symmetric = FALSE) 6269.312103 6370.354603 6392.572385 6401.074078 6460.657460 6385.649158 36.090058 0.005651745 25 16 lu(A) 216.243199 218.889777 219.535822 224.236960 293.879314 231.447522 27.371226 0.118261044 25 17 fft(A) 158.338090 160.060496 161.838452 163.557053 167.644866 161.902706 2.656454 0.016407715 25 18 Hilbert(3000) 172.142450 246.993495 247.904730 249.638550 462.391280 253.838923 50.756220 0.199954442 25 19 toeplitz(A[1:500, 1]) 5.668683 5.860442 5.888598 5.925884 12.543850 6.617860 2.049864 0.309747276 25 20 princomp(A) 2811.828559 2848.070617 2858.973481 2878.443481 2977.561162 2871.749801 40.471198 0.014092870 25 NotBLAS: expr Min LQ Median UQ Max Mean SD CV n (fctr) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (int) 1 A + 2 2.418674 2.573812 2.590268 2.613096 87.401743 3.3176997 3.04139885 0.91671914 5000 2 A - 2 2.443785 2.589127 2.607769 2.646578 81.559183 3.4098016 3.06170737 0.89791363 5000 3 A * 2 2.439981 2.590553 2.608530 2.643534 87.271621 3.4043742 3.10091151 0.91086095 5000 4 A/2 2.858122 2.972264 3.014497 3.047979 81.612069 3.7626603 2.98407959 0.79307707 5000 5 A * 0.5 2.434654 2.586843 2.603965 2.621847 81.566412 3.3707782 3.07975952 0.91366423 5000 6 A^2 2.443405 2.589507 2.607009 2.638303 81.988739 3.3955817 3.07268185 0.90490587 5000 7 sqrt(A[1:10000]) 0.175399 0.180345 0.181868 0.183389 7.468706 0.1892085 0.11420929 0.60361592 5000 8 sin(A[1:10000]) 0.617510 0.627022 0.630826 0.645665 1.770346 0.6375912 0.04872665 0.07642303 5000 9 A + B 1.626908 1.661150 1.673325 1.707188 75.651562 1.8341317 2.06798396 1.12750028 5000 10 A - B 1.625386 1.664194 1.679413 1.715177 75.467413 1.8418374 2.06893026 1.12329690 5000 11 A * B 1.625766 1.657726 1.667238 1.681315 74.191303 1.8239461 2.05043314 1.12417421 5000 12 A/B 2.628315 2.635924 2.645056 2.705932 76.369136 2.8169433 2.07280333 0.73583424 5000 13 A[1:100000]%%B[1:100000] 4.322947 4.334742 4.338927 4.407888 76.947076 4.4387060 1.45612013 0.32805059 5000 14 A[1:100000]%/%B[1:100000] 3.939048 3.949702 3.954268 3.979474 74.364799 4.0555842 1.42448375 0.35124008 5000 ---------------------- 493-SEH-POSIX-OPB-TG BLAS: expr Min LQ Median UQ Max Mean SD CV n (fctr) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (int) 1 sort(c(as.vector(A), as.vector(B))) 346.307880 409.971084 412.729522 415.560629 481.249097 407.873714 27.013964 0.06623120 25 2 det(A) 25.179774 25.313320 25.499754 26.092912 37.814938 27.056227 3.272499 0.12095177 25 3 A %*% B 32.317469 33.685653 39.154206 46.913214 57.920329 40.802831 7.227388 0.17712957 25 4 t(A) %*% B 43.173915 43.719135 48.276072 58.702584 73.592056 52.603995 9.716491 0.18471013 25 5 crossprod(A, B) 29.180460 29.459347 36.169388 45.822015 57.051327 39.234015 9.886814 0.25199598 25 6 solve(A) 68.554579 71.990641 75.369250 84.929803 166.340220 82.552723 19.985428 0.24209289 25 7 solve(A, t(B)) 77.967508 81.424876 85.531713 94.843054 154.065746 90.099335 15.516925 0.17222019 25 8 solve(B) 88.574745 92.768330 99.096755 106.515619 119.114257 99.423613 7.972067 0.08018284 25 9 chol(A) 10.977438 11.635278 11.876118 12.908724 88.026862 15.454507 15.194240 0.98315916 25 10 chol(B, pivot = TRUE) 3.923449 3.957692 3.975574 4.114827 9.395807 4.856174 2.000903 0.41203279 25 11 qr(A, LAPACK = TRUE) 129.071263 130.399118 130.798235 133.491992 151.146369 132.853541 4.864756 0.03661743 25 12 svd(A) 596.404520 604.072975 608.158126 618.875319 689.199483 623.758567 31.733271 0.05087428 25 13 eigen(A, symmetric = TRUE) 279.824221 284.914964 287.559640 290.222578 359.950919 290.351132 14.891298 0.05128721 25 14 eigen(A, symmetric = FALSE) 1345.076666 1350.822585 1355.789294 1362.679679 1434.472092 1361.770864 22.142251 0.01625989 25 15 eigen(B, symmetric = FALSE) 4792.117719 4809.680402 4820.062017 4829.306016 4855.549976 4820.368391 14.864522 0.00308369 25 16 lu(A) 30.719477 32.852414 34.078302 35.978009 110.573755 45.977766 28.040390 0.60986847 25 17 fft(A) 158.454514 160.357266 161.992924 162.896169 166.749230 161.838192 2.153864 0.01330875 25 18 Hilbert(3000) 169.336453 246.676940 247.443595 250.194803 475.146293 254.182628 53.268882 0.20956933 25 19 toeplitz(A[1:500, 1]) 5.634061 5.825439 5.881369 5.913709 12.558688 6.594544 2.040480 0.30941943 25 20 princomp(A) 434.713313 446.791843 453.297188 465.176350 562.749686 471.074398 39.135293 0.08307667 25 NotBLAS: expr Min LQ Median UQ Max Mean SD CV n (fctr) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (int) 1 A + 2 2.416011 2.580376 2.597687 2.631359 86.794126 3.3246622 3.04057878 0.91455270 5000 2 A - 2 2.436557 2.593312 2.615379 2.669406 83.664727 3.4119219 3.05840683 0.89638829 5000 3 A * 2 2.452156 2.595214 2.615189 2.662178 81.755507 3.4092902 3.05264857 0.89539124 5000 4 A/2 2.840620 2.967699 2.999278 3.054066 82.002816 3.7354622 2.94315732 0.78789642 5000 5 A * 0.5 2.442264 2.589888 2.607770 2.649241 82.696420 3.3905433 3.06878468 0.90510117 5000 6 A^2 2.452917 2.587985 2.603965 2.624891 82.177834 3.3809857 3.06701904 0.90713755 5000 7 sqrt(A[1:10000]) 0.175399 0.180345 0.181867 0.182628 7.417723 0.1883772 0.11272926 0.59842300 5000 8 sin(A[1:10000]) 0.617890 0.626261 0.629685 0.635012 1.783663 0.6353244 0.04901841 0.07715493 5000 9 A + B 1.625006 1.660294 1.675988 1.709851 75.279839 1.8392480 2.07033888 1.12564424 5000 10 A - B 1.626908 1.658487 1.671994 1.704144 74.425295 1.8352869 2.06002358 1.12245314 5000 11 A * B 1.624624 1.670186 1.699958 1.724690 75.242552 1.8541044 2.08118712 1.12247568 5000 12 A/B 2.629076 2.636685 2.641632 2.693757 76.012252 2.8111608 2.05562305 0.73123637 5000 13 A[1:100000]%%B[1:100000] 4.324850 4.336263 4.340829 4.417685 75.107104 4.4458514 1.42743953 0.32107225 5000 14 A[1:100000]%/%B[1:100000] 3.939429 3.952365 4.022182 4.045201 74.483126 4.0785691 1.42384261 0.34910347 5000 ---------------------- 493-SEH-POSIX-OPB-TN BLAS: expr Min LQ Median UQ Max Mean SD CV n (fctr) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (int) 1 sort(c(as.vector(A), as.vector(B))) 337.213028 409.064795 410.194043 412.993951 484.426436 405.382111 29.481412 0.07272499 25 2 det(A) 25.287829 25.372294 25.448769 29.880532 37.590839 28.022327 4.245258 0.15149554 25 3 A %*% B 32.075867 32.750828 41.853671 49.000495 63.442530 42.848184 9.647849 0.22516354 25 4 t(A) %*% B 44.049005 44.257886 48.174866 58.273789 74.823650 52.422129 10.207673 0.19472068 25 5 crossprod(A, B) 28.694214 29.117300 33.848876 39.668607 59.277480 35.716943 8.192299 0.22936732 25 6 solve(A) 70.251494 73.326105 77.655900 87.781076 156.925009 83.621231 17.866184 0.21365607 25 7 solve(A, t(B)) 78.658830 83.161740 87.114105 92.341437 169.351292 92.029023 17.823982 0.19367783 25 8 solve(B) 89.515657 92.129514 96.643839 101.953735 112.393182 98.131097 6.888087 0.07019270 25 9 chol(A) 12.309858 13.065100 13.289579 13.591295 92.378724 16.833847 15.826692 0.94017080 25 10 chol(B, pivot = TRUE) 3.893772 4.079063 4.133090 4.301640 9.699806 5.015212 2.018333 0.40244232 25 11 qr(A, LAPACK = TRUE) 128.787810 130.288781 130.915421 133.212343 158.127307 133.122840 5.926681 0.04452039 25 12 svd(A) 603.183427 609.069361 611.918351 642.481072 688.660352 630.878287 31.756437 0.05033687 25 13 eigen(A, symmetric = TRUE) 284.978504 288.197314 291.064187 295.776358 390.315033 296.038383 20.733087 0.07003513 25 14 eigen(A, symmetric = FALSE) 1367.448541 1377.090515 1382.783929 1418.322871 1557.108778 1409.269505 50.194480 0.03561737 25 15 eigen(B, symmetric = FALSE) 4650.988092 4698.182095 4706.879733 4725.462466 4951.773771 4726.326539 68.648561 0.01452472 25 16 lu(A) 30.423468 34.784081 36.049919 38.482289 112.387856 47.798973 28.253927 0.59109904 25 17 fft(A) 158.335426 160.232089 161.947266 164.481605 166.713465 162.145098 2.656089 0.01638094 25 18 Hilbert(3000) 172.198760 246.243580 246.902561 249.914014 462.000153 255.145882 50.443935 0.19770625 25 19 toeplitz(A[1:500, 1]) 5.626071 5.818210 5.849028 5.884032 12.472701 6.569007 2.030083 0.30903950 25 20 princomp(A) 448.094967 456.912834 466.493171 486.485942 550.962600 479.347289 31.503360 0.06572137 25 NotBLAS: expr Min LQ Median UQ Max Mean SD CV n (fctr) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (int) 1 A + 2 2.452917 2.592932 2.610433 2.652286 87.372447 3.3546970 3.02092771 0.90050686 5000 2 A - 2 2.446449 2.588746 2.607009 2.631359 84.036070 3.3791704 3.09526644 0.91598414 5000 3 A * 2 2.450254 2.596355 2.620706 2.667885 82.293498 3.4108719 3.06978745 0.90000080 5000 4 A/2 2.855077 2.963513 2.996234 3.055303 81.036792 3.7466073 2.93954950 0.78458970 5000 5 A * 0.5 2.450634 2.592931 2.611194 2.641252 84.453069 3.3938652 3.09108315 0.91078549 5000 6 A^2 2.443405 2.588366 2.606248 2.630598 81.693490 3.3737386 3.07150401 0.91041554 5000 7 sqrt(A[1:10000]) 0.166267 0.167409 0.167790 0.184911 7.397558 0.1779341 0.11297189 0.63490852 5000 8 sin(A[1:10000]) 0.603812 0.604954 0.605335 0.606856 1.703383 0.6108681 0.04838561 0.07920796 5000 9 A + B 1.631854 1.661150 1.685501 1.714417 75.956322 1.8429324 2.08253557 1.13001191 5000 10 A - B 1.631093 1.667999 1.693491 1.717461 75.607808 1.8500438 2.09123392 1.13036992 5000 11 A * B 1.629952 1.680174 1.701861 1.724785 75.548834 1.8583676 2.09261679 1.12605103 5000 12 A/B 2.630217 2.637827 2.650193 2.709737 76.482518 2.8199917 2.07967625 0.73747602 5000 13 A[1:100000]%%B[1:100000] 4.308109 4.321045 4.330936 4.416163 78.209109 4.4450468 1.46250262 0.32901849 5000 14 A[1:100000]%/%B[1:100000] 3.932961 3.943614 3.948941 4.022372 76.773200 4.0555702 1.45150857 0.35790493 500 ---------------------- 493-SEH-POSIX-OPB-AN BLAS: expr Min LQ Median UQ Max Mean SD CV n (fctr) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (int) 1 sort(c(as.vector(A), as.vector(B))) 338.070236 409.391242 410.171214 415.280600 478.151657 406.139300 27.670445 0.06813043 25 2 det(A) 25.070579 25.203744 25.528289 26.098619 40.153332 26.703725 3.340470 0.12509380 25 3 A %*% B 31.829320 32.168703 36.989309 45.802230 54.174181 38.917978 7.713159 0.19819012 25 4 t(A) %*% B 43.356923 45.510787 55.702926 60.049462 69.068980 53.990610 8.193159 0.15175156 25 5 crossprod(A, B) 28.787049 29.109691 35.639006 48.920215 62.759198 39.510315 11.315157 0.28638488 25 6 solve(A) 68.606705 73.490089 76.856143 86.750372 144.050147 81.875510 15.143255 0.18495464 25 7 solve(A, t(B)) 78.407717 81.909981 87.044478 91.318343 155.031009 90.026832 15.196156 0.16879586 25 8 solve(B) 88.235361 94.640643 100.222198 102.224632 116.761405 99.472161 7.287927 0.07326600 25 9 chol(A) 11.537877 13.067763 13.100864 13.469924 89.405700 16.305719 15.278163 0.93698188 25 10 chol(B, pivot = TRUE) 3.882738 3.926873 4.064985 4.327132 9.266827 4.886840 1.911920 0.39123852 25 11 qr(A, LAPACK = TRUE) 128.206826 129.343302 129.879771 131.059240 158.334666 131.401028 5.829694 0.04436566 25 12 svd(A) 594.256744 598.659971 601.371229 627.951528 683.380512 617.790861 30.879183 0.04998323 25 13 eigen(A, symmetric = TRUE) 280.515923 283.413233 285.617319 289.141651 372.716205 290.123624 17.842270 0.06149885 25 14 eigen(A, symmetric = FALSE) 1418.878743 1427.698893 1431.613972 1443.283872 1512.860024 1440.704060 23.652272 0.01641716 25 15 eigen(B, symmetric = FALSE) 4734.774570 4789.189591 4795.702926 4804.545144 4876.409466 4799.317689 25.957351 0.00540855 25 16 lu(A) 30.834000 33.003843 33.569228 35.529431 107.357228 45.017114 27.507827 0.61105265 25 17 fft(A) 154.031884 155.367348 157.674543 158.455656 163.550965 157.565804 2.447198 0.01553128 25 18 Hilbert(3000) 168.989461 247.013279 248.481528 250.298292 475.175969 254.993570 53.423863 0.20951063 25 19 toeplitz(A[1:500, 1]) 5.809459 5.967356 6.022524 6.045353 12.230339 6.701458 1.937309 0.28908775 25 20 princomp(A) 439.784652 450.565005 456.988168 469.762965 538.406575 469.695819 33.055649 0.07037671 25 NotBLAS: expr Min LQ Median UQ Max Mean SD CV n (fctr) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (int) 1 A + 2 2.504661 2.664080 2.678918 2.698322 87.418483 3.3983430 3.05405466 0.89868934 5000 2 A - 2 2.512652 2.675114 2.694327 2.741031 81.950691 3.5013618 3.06928965 0.87659883 5000 3 A * 2 2.516456 2.670548 2.689571 2.718107 85.239889 3.4743264 3.09382426 0.89048176 5000 4 A/2 2.984058 3.102291 3.155652 3.184949 82.517597 3.8927537 3.02935221 0.77820290 5000 5 A * 0.5 2.514173 2.673211 2.692235 2.723814 83.803219 3.4764812 3.09662027 0.89073407 5000 6 A^2 2.515695 2.672070 2.689952 2.709737 85.336150 3.4585685 3.10400250 0.89748186 5000 7 sqrt(A[1:10000]) 0.874330 0.875852 0.876994 0.894876 7.878096 0.8899799 0.11092121 0.12463338 5000 8 sin(A[1:10000]) 1.233497 1.235019 1.238443 1.255565 2.415250 1.2506302 0.05097291 0.04075778 5000 9 A + B 1.616635 1.653921 1.672945 1.706046 74.484268 1.8330335 2.06861743 1.12852136 5000 10 A - B 1.617396 1.649736 1.669140 1.701861 75.806415 1.8293446 2.08478826 1.13963674 5000 11 A * B 1.616635 1.648975 1.668760 1.698817 76.186129 1.8284716 2.09280167 1.14456337 5000 12 A/B 2.638969 2.648100 2.652666 2.673306 76.799452 2.8162757 2.07349971 0.73625594 5000 13 A[1:100000]%%B[1:100000] 4.318381 4.330937 4.341019 4.425294 76.245102 4.4518243 1.44204810 0.32392296 5000 14 A[1:100000]%/%B[1:100000] 3.961877 3.975574 3.978998 3.995454 75.204885 4.0754915 1.43546816 0.35221964 5000 ---------------------- 493-SEH-POSIX-OPB-TN-LTO BLAS: expr Min LQ Median UQ Max Mean SD CV n (fctr) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (int) 1 sort(c(as.vector(A), as.vector(B))) 340.558536 409.867215 411.248716 413.309744 477.90435 406.254309 27.336146 0.067288262 25 2 det(A) 27.324887 27.515124 27.770803 28.928586 35.90077 28.899959 2.305326 0.079769176 25 3 A %*% B 32.033635 34.582810 40.307424 45.174068 66.20553 40.601926 7.830137 0.192851369 25 4 t(A) %*% B 43.521669 44.060801 48.206446 58.209870 70.86216 52.072306 8.863299 0.170211372 25 5 crossprod(A, B) 28.730739 28.926683 29.297645 35.906099 59.40722 35.318891 10.478419 0.296680305 25 6 solve(A) 77.447020 81.739528 84.994484 95.989423 161.12202 89.614367 16.382909 0.182815648 25 7 solve(A, t(B)) 87.636115 92.239852 100.343569 104.028460 172.51303 101.826520 16.870794 0.165681728 25 8 solve(B) 98.305749 102.707454 105.001712 108.554960 132.09223 106.990451 7.987065 0.074652128 25 9 chol(A) 15.000571 16.397672 16.631283 16.774341 92.56934 19.879755 15.195733 0.764382313 25 10 chol(B, pivot = TRUE) 4.971274 5.067154 5.272610 5.369251 10.38504 5.990002 1.863047 0.311026207 25 11 qr(A, LAPACK = TRUE) 130.530762 130.915421 131.340792 134.512423 156.40718 133.232021 5.234107 0.039285652 25 12 svd(A) 599.521745 604.568353 610.394552 641.227030 696.29532 626.818385 31.515757 0.050278928 25 13 eigen(A, symmetric = TRUE) 283.823004 286.884679 288.655785 290.506031 366.98893 293.611080 16.837928 0.057347728 25 14 eigen(A, symmetric = FALSE) 1402.538143 1414.911160 1420.025111 1428.449949 1487.35190 1425.261105 19.448975 0.013645903 25 15 eigen(B, symmetric = FALSE) 4663.388123 4669.545334 4674.298596 4683.889207 4793.97519 4686.093049 30.756637 0.006563386 25 16 lu(A) 33.084124 33.933722 35.683522 37.658943 109.66404 47.168862 27.664198 0.586492799 25 17 fft(A) 157.868966 158.696876 161.171480 162.357418 169.64426 161.197291 2.877905 0.017853307 25 18 Hilbert(3000) 169.154967 246.404901 247.133128 248.278355 463.81425 252.812130 51.380214 0.203234766 25 19 toeplitz(A[1:500, 1]) 5.669064 5.869193 5.900393 5.930450 12.54728 6.617845 2.042603 0.308650813 25 20 princomp(A) 432.820836 441.650116 450.239319 476.593616 544.67983 467.748964 36.576092 0.078195986 25 NotBLAS: expr Min LQ Median UQ Max Mean SD CV n (fctr) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (int) 1 A + 2 2.440361 2.596736 2.609292 2.626413 86.483660 3.3156484 2.97676420 0.89779248 5000 2 A - 2 2.478408 2.591028 2.610053 2.656851 82.008903 3.3752912 2.96901986 0.87963369 5000 3 A * 2 2.466614 2.585702 2.600161 2.626032 81.941940 3.3251111 2.96781542 0.89254625 5000 4 A/2 2.885516 3.003843 3.037325 3.071567 83.144999 3.7964458 2.95906446 0.77943019 5000 5 A * 0.5 2.481071 2.588365 2.604346 2.643915 83.808166 3.3599627 2.98011193 0.88694791 5000 6 A^2 2.481072 2.586463 2.602443 2.627935 80.749535 3.3274835 2.95750492 0.88881130 5000 7 sqrt(A[1:10000]) 0.165887 0.167028 0.167409 0.167790 6.842446 0.1740123 0.10554791 0.60655443 5000 8 sin(A[1:10000]) 0.603812 0.604954 0.605715 0.608854 1.737245 0.6123042 0.04881917 0.07973026 5000 9 A + B 1.592665 1.636800 1.648214 1.673326 75.685425 1.8065513 2.07547387 1.14885966 5000 10 A - B 1.586958 1.629951 1.644029 1.664574 75.693415 1.8008704 2.07725925 1.15347515 5000 11 A * B 1.593806 1.648595 1.675228 1.700719 75.657649 1.8280364 2.09904485 1.14825113 5000 12 A/B 2.638968 2.651524 2.662178 2.723815 76.903703 2.8353980 2.09631478 0.73933704 5000 13 A[1:100000]%%B[1:100000] 4.383062 4.392954 4.405129 4.488073 75.980672 4.5165001 1.44404570 0.31972670 5000 14 A[1:100000]%/%B[1:100000] 3.930298 3.940190 3.945897 4.023609 78.317164 4.0556036 1.46563175 0.36138437 5000 ----------------------- 493-SEH-POSIX-OPB-AN-LTO BLAS: expr Min LQ Median UQ Max Mean SD CV n (fctr) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (int) 1 sort(c(as.vector(A), as.vector(B))) 337.827494 410.093216 411.384165 414.789028 485.03672 407.029335 28.906795 0.071018949 25 2 det(A) 27.510939 27.720200 27.978922 28.437774 35.59373 28.956238 2.276051 0.078603147 25 3 A %*% B 31.944604 32.177073 32.623749 40.002283 61.34574 37.693293 8.333761 0.221093995 25 4 t(A) %*% B 43.391927 43.841267 44.607541 47.460336 76.44485 49.743409 10.531100 0.211708441 25 5 crossprod(A, B) 29.035499 29.272915 32.178975 39.732147 60.89526 36.703604 9.550389 0.260203028 25 6 solve(A) 78.760416 80.972493 83.411711 93.028955 163.80893 89.787787 17.554056 0.195506057 25 7 solve(A, t(B)) 88.751285 93.062056 105.039760 122.716205 212.39508 110.557593 25.864917 0.233949709 25 8 solve(B) 100.237797 104.213751 105.335008 111.516570 131.88144 108.083264 7.321555 0.067739951 25 9 chol(A) 14.351102 15.072100 15.501275 16.496215 92.45596 18.918419 15.386295 0.813297054 25 10 chol(B, pivot = TRUE) 4.866263 5.019976 5.091124 5.171404 10.52658 5.915048 1.992732 0.336891925 25 11 qr(A, LAPACK = TRUE) 130.162843 130.890690 132.013469 134.329035 162.14397 133.785794 6.294900 0.047052083 25 12 svd(A) 603.576076 611.122779 614.020089 638.721607 723.26827 632.760294 36.667892 0.057949103 25 13 eigen(A, symmetric = TRUE) 285.298862 287.357989 288.944946 290.951566 361.17947 293.482267 16.107772 0.054884992 25 14 eigen(A, symmetric = FALSE) 1471.813340 1474.687821 1479.895370 1488.222045 1550.27432 1488.134658 21.525279 0.014464604 25 15 eigen(B, symmetric = FALSE) 4842.096032 4921.975042 4927.130846 4933.083362 4950.56120 4924.717270 19.570796 0.003973994 25 16 lu(A) 32.917476 33.312408 35.818590 37.641061 112.61995 47.200366 28.019708 0.593633280 25 17 fft(A) 154.482746 155.970019 158.706008 159.758399 163.39117 158.356307 2.583029 0.016311503 25 18 Hilbert(3000) 169.070502 247.656281 249.680782 250.975536 474.84267 257.501519 53.550114 0.207960382 25 19 toeplitz(A[1:500, 1]) 6.073889 6.193737 6.622152 6.826847 12.62451 7.179166 1.954066 0.272185730 25 20 princomp(A) 440.818780 450.165126 456.150744 488.413804 548.81292 474.494738 37.590410 0.079221974 25 NotBLAS: Source: local data frame [14 x 10] expr Min LQ Median UQ Max Mean SD CV n (fctr) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (dbl) (int) 1 A + 2 2.685767 2.839098 2.857361 2.876765 87.453107 3.5864157 3.06287775 0.85402197 5000 2 A - 2 2.662558 2.840620 2.861926 2.899593 82.553742 3.6577826 3.09258632 0.84548117 5000 3 A * 2 2.672070 2.840239 2.862687 2.889701 82.458243 3.6544753 3.09686931 0.84741832 5000 4 A/2 3.127116 3.229369 3.273599 3.321920 81.775292 4.0329116 2.96593802 0.73543342 5000 5 A * 0.5 2.675494 2.836435 2.858121 2.877145 82.497051 3.6318914 3.10344439 0.85449812 5000 6 A^2 2.663700 2.842903 2.864590 2.913290 82.562494 3.6729462 3.08607449 0.84021773 5000 7 sqrt(A[1:10000]) 0.874329 0.875471 0.875852 0.877373 8.139102 0.8838267 0.11404226 0.12903238 5000 8 sin(A[1:10000]) 1.233497 1.235019 1.236922 1.254804 2.411065 1.2490264 0.05153307 0.04125859 5000 9 A + B 1.622342 1.660009 1.679032 1.713275 75.909524 1.8425114 2.09472155 1.13688392 5000 10 A - B 1.625006 1.663813 1.685500 1.717840 75.927787 1.8481592 2.09249212 1.13220339 5000 11 A * B 1.621581 1.656965 1.679413 1.716319 75.992467 1.8420509 2.09863780 1.13929412 5000 12 A/B 2.640491 2.655710 2.787734 2.849371 77.599590 2.9064907 2.09512050 0.72084196 5000 13 A[1:100000]%%B[1:100000] 4.314956 4.326752 4.332839 4.418827 77.936309 4.4452670 1.47134520 0.33099141 5000 14 A[1:100000]%/%B[1:100000] 3.919644 3.929536 4.001827 4.017426 77.226725 4.0562177 1.46632776 0.36150125 5000