Details

simplify48.5m

Calls

16870 calls:

Slowest
2.4s
(/ (* (* (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)) (* (* (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)) (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592))) (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)))) (* (* (* (fma 0.3275911 (fabs x) 1) (fma 0.3275911 (fabs x) 1)) (fma 0.3275911 (fabs x) 1)) (* (* (fma 0.3275911 (fabs x) 1) (fma 0.3275911 (fabs x) 1)) (fma 0.3275911 (fabs x) 1))))
2.4s
(/ (* (* (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)) (* (* (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)) (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592))) (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)))) (* (* (* (fma 0.3275911 (fabs x) 1) (fma 0.3275911 (fabs x) 1)) (fma 0.3275911 (fabs x) 1)) (* (* (fma 0.3275911 (fabs x) 1) (fma 0.3275911 (fabs x) 1)) (fma 0.3275911 (fabs x) 1))))
2.4s
(/ (* (* (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)) (* (* (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)) (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592))) (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)))) (* (* (* (fma 0.3275911 (fabs x) 1) (fma 0.3275911 (fabs x) 1)) (* (fma 0.3275911 (fabs x) 1) (fma 0.3275911 (fabs x) 1))) (* (fma 0.3275911 (fabs x) 1) (fma 0.3275911 (fabs x) 1))))
2.3s
(/ (* (* (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)) (* (* (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)) (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592))) (* (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (fma (/ 1 (fma 0.3275911 (fabs x) 1)) (+ (/ 1.061405429 (fma 0.3275911 (fabs x) 1)) -1.453152027) 1.421413741) -0.284496736) 0.254829592)))) (* (* (* (fma 0.3275911 (fabs x) 1) (fma 0.3275911 (fabs x) 1)) (* (fma 0.3275911 (fabs x) 1) (fma 0.3275911 (fabs x) 1))) (* (fma 0.3275911 (fabs x) 1) (fma 0.3275911 (fabs x) 1))))
2.3s
(/ (* (* (* (fma (+ (+ alpha beta) i) i (* beta alpha)) (* (+ (+ alpha beta) i) i)) (* (fma (+ (+ alpha beta) i) i (* beta alpha)) (* (+ (+ alpha beta) i) i))) (* (fma (+ (+ alpha beta) i) i (* beta alpha)) (* (+ (+ alpha beta) i) i))) (* (* (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta))) (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta)))) (* (fma 2 i (+ alpha beta)) (fma 2 i (+ alpha beta)))))

sample9.4m

Algorithm
88×intervals
Results
6.5m101386×body10240exit
1.0m384320×body80valid
39.4s25800×body1280valid
20.3s19826×body640valid
11.6s70468×body80nan
6.4s10236×body320valid
4.9s1554×body1280nan
4.6s9185×body160valid
3.0s1306×body640nan
1.7s1473×body2560valid
1.6s645×body320nan
1.5s33536×pre80true
546.0ms567×body5120valid
325.0ms354×body160nan

series5.7m

Calls

457 calls:

Slowest
4.4m
(* (/ (/ (+ alpha beta) 1) (sqrt (sqrt (+ (+ (+ alpha beta) (* 2 i)) 2.0)))) (/ (/ (- beta alpha) (+ (+ alpha beta) (* 2 i))) (sqrt (+ (+ (+ alpha beta) (* 2 i)) 2.0))))
2.0s
(/ 1 (* (pow -1 1.0) (* (pow -2 1.0) (pow (log (/ 1 u1)) 1.0))))
1.5s
(* 6 (pow (/ 1 (* (pow -1 1.0) (* (pow -2 1.0) (pow (log (/ 1 u1)) 1.0)))) 0.5))
1.5s
(/ 1 (* 6 (pow (/ 1 (* (pow -1 1.0) (* (pow -2 1.0) (pow (log (/ 1 u1)) 1.0)))) 0.5)))
1.4s
(/ 6 (pow (* -2 (log u1)) 0.5))

rewrite5.4m

Algorithm
149×rewrite-expression-head
Calls

457 calls:

Slowest
1.3m
(* (/ 1 (exp (+ (- (- 1 z) -6) 0.5))) (+ (+ (/ 771.3234287776531 (+ (- 1 z) 2)) (+ (+ (+ (/ -1259.1392167224028 (+ (- 1 z) 1)) 0.9999999999998099) (/ 676.5203681218851 (- 1 z))) (/ -176.6150291621406 (- 5 (+ z 1))))) (+ (+ (/ -0.13857109526572012 (- 7 (+ z 1))) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6))))))
1.3m
(* (/ 1 (exp (+ (- (- 1 z) -6) 0.5))) (+ (+ (/ 771.3234287776531 (+ (- 1 z) 2)) (+ (+ (+ (/ -1259.1392167224028 (+ (- 1 z) 1)) 0.9999999999998099) (/ 676.5203681218851 (- 1 z))) (/ -176.6150291621406 (- 5 (+ z 1))))) (+ (+ (/ -0.13857109526572012 (- 7 (+ z 1))) (/ 12.507343278686905 (- 6 (+ z 1)))) (+ (/ 1.5056327351493116e-07 (- 9 (+ z 1))) (/ 9.984369578019572e-06 (- (- 1 z) -6))))))
36.6s
(* (+ (/ 12.507343278686905 (+ z 4)) (+ (+ (+ (/ 771.3234287776531 (+ z 2)) (+ (/ 676.5203681218851 z) (+ 0.9999999999998099 (/ -1259.1392167224028 (- z -1))))) (/ -176.6150291621406 (+ z 3))) (/ -0.13857109526572012 (- z -5)))) (* (* (pow (+ (+ 7 (- z 1)) 0.5) (+ (- z 1) 0.5)) (sqrt (* PI 2))) (exp (- (+ (+ 7 (- z 1)) 0.5)))))
16.1s
(* (/ (+ (+ (+ (+ (+ 1 (* 0.1049934947 (* x x))) (* 0.0424060604 (* (* x x) (* x x)))) (* 0.0072644182 (* (* (* x x) (* x x)) (* x x)))) (* 0.0005064034 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))) (* 0.0001789971 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (+ (+ (+ (+ (+ (+ 1 (* 0.7715471019 (* x x))) (* 0.2909738639 (* (* x x) (* x x)))) (* 0.0694555761 (* (* (* x x) (* x x)) (* x x)))) (* 0.0140005442 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))) (* 0.0008327945 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (* (* 2 0.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))) x)
15.3s
(* (/ (+ (+ (+ (+ (+ 1 (* 0.1049934947 (* x x))) (* 0.0424060604 (* (* x x) (* x x)))) (* 0.0072644182 (* (* (* x x) (* x x)) (* x x)))) (* 0.0005064034 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))) (* 0.0001789971 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (+ (+ (+ (+ (+ (+ 1 (* (* x x) 0.7715471019)) (* 0.2909738639 (* (* x x) (* x x)))) (* 0.0694555761 (* (* (* x x) (* x x)) (* x x)))) (* 0.0140005442 (* (* (* (* x x) (* x x)) (* x x)) (* x x)))) (* 0.0008327945 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (* (* 2 0.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))) x)
Rules
25406×*-un-lft-identity
20257×times-frac
10167×add-sqr-sqrt
9925×distribute-lft-out
8640×add-cube-cbrt
7104×frac-add
5155×associate-*l/
4268×flip-+ flip3-+
4229×frac-times
3998×add-exp-log
3052×associate-*r/
2640×pow1
2118×prod-exp
1977×add-cbrt-cube
1855×cube-div
1702×frac-sub
1509×prod-diff
1447×sqrt-prod
1274×flip3-- flip--
1207×associate-*l*
1157×associate-*r*
1061×div-exp
901×associate-/l*
837×pow-exp
663×cbrt-unprod
655×associate-/r/
640×pow-prod-down
627×add-log-exp
595×div-inv
569×pow-prod-up
545×associate-/r*
512×distribute-rgt-in distribute-lft-in
460×log1p-expm1-u
458×expm1-log1p-u
457×insert-posit16
453×pow-to-exp
362×unswap-sqr
347×1-exp
337×cbrt-undiv
331×rec-exp
264×pow-plus
237×difference-of-squares
224×sub-div
218×log-pow
211×cbrt-div
200×swap-sqr
193×distribute-lft-out--
184×fma-neg
183×sqrt-div
182×sqrt-pow1
171×unpow-prod-down
165×fma-def
158×*-commutative
145×pow-sqr
121×exp-sum
120×cbrt-prod
118×pow1/2
90×pow2
87×frac-2neg clear-num
82×fma-udef
80×associate-+l+
74×sub-neg
70×sqr-pow
69×associate--r+
63×sum-log
59×associate-/l/
53×log-prod
50×pow1/3
47×associate-+l-
44×pow-flip pow-sub
40×+-commutative
39×exp-neg
32×cube-prod pow-div
31×cube-unmult
30×diff-log
29×exp-diff inv-pow
26×neg-sub0 distribute-rgt1-in distribute-rgt-out
24×associate-+r+ rem-sqrt-square
21×pow-unpow
18×pow-pow
17×un-div-inv
14×hypot-udef
12×div-sub hypot-def exp-prod
11×sin-sum
10×pow3
8×sqrt-unprod
5×unpow-prod-up rem-square-sqrt
4×expm1-log1p expm1-udef rem-cube-cbrt log1p-expm1 unpow3 cube-mult rem-log-exp rem-cbrt-cube log1p-udef
3×associate--l+
2×associate--r- log-div
1×sqrt-pow2 sum-cubes distribute-rgt-out-- cos-sum sin-mult

prune5.1m

regimes8.7s

Accuracy

Total 14.6b remaining (10.7%)

3.1b76.3%Octave 3.8, jcobi/1
2.6b0%_divideComplex, real part
2.1b0%_divideComplex, imaginary part
1.5b39.8%Octave 3.8, jcobi/4
1.0b0%Jmat.Real.lambertw, newton loop step

localize8.2s

bsearch1.4s

end1.0ms