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 │ └──────────┴───────────┴───────────┴───────────┘