Statistical benchmarking in J.
See the documentation or browse the files here.
Using pacman
install 'github:jitwit/bonsai@a'
Inverting a random 50 by 50 matrix.
load 'stats/bonsai' bonsai '%. ? 50 50 $ 0' ββββββββββ¬ββββββββββ¬βββββββββββ¬βββββββββββ βN = 2000βlower βestimate βupper β ββββββββββΌββββββββββΌβββββββββββΌβββββββββββ€ βΞΌ β9.3306e_5β9.3794e_5 β9.42814e_5β ββββββββββΌββββββββββΌβββββββββββΌβββββββββββ€ βΟ β9.0145e_6β1.09376e_5β1.26696e_5β ββββββββββ΄ββββββββββ΄βββββββββββ΄βββββββββββ
Comparing two ways to get the 1000 by 1000 identity matrix.
'e. i. 1000' bonsai '= i. 1000' ββββββββββββ¬ββββββββββββ¬ββββββββββββ¬ββββββββββββ βcomparisonβlower βestimate βupper β ββββββββββββΌββββββββββββΌββββββββββββΌββββββββββββ€ β- & ΞΌ β0.00147287 β0.00148659 β0.00150041 β ββββββββββββΌββββββββββββΌββββββββββββΌββββββββββββ€ βΞΌ(x) β0.00165997 β0.00167314 β0.00168693 β ββββββββββββΌββββββββββββΌββββββββββββΌββββββββββββ€ βΞΌ(y) β0.000185358β0.000186555β0.000187817β ββββββββββββΌββββββββββββΌββββββββββββΌββββββββββββ€ βx (-~%[) yβ_0.888336 β_0.8885 β_0.888663 β ββββββββββββ΄ββββββββββββ΄ββββββββββββ΄ββββββββββββ