1 USING: math math.combinatorics math.matrices.simd
2 prettyprint sequences typed ;
3 IN: benchmark.matrix-exponential-simd
5 TYPED:: e^m4 ( m: matrix4 iterations: fixnum -- e^m: matrix4 )
7 iterations <iota> [| i |
8 m i m4^n i factorial >float m4/n m4+
11 :: matrix-exponential-simd-benchmark ( -- )
14 identity-matrix4 20 e^m4 result!
18 MAIN: matrix-exponential-simd-benchmark