Details

simplify46.2m

Calls

15599 calls:

Slowest
2.1s
(/ (* (* (* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0)))) (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0)))) (* (* (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i)))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))))) (* (* (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0)) (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))))
1.8s
(* (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (- (* 0.254829592 0.254829592) (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))) (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429)))))))))) (exp (- (* (fabs x) (fabs x))))) (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (- (* 0.254829592 0.254829592) (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))) (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429)))))))))) (exp (- (* (fabs x) (fabs x))))))
1.7s
(* (* (* (sqrt (/ 1 PI)) (sqrt (/ 1 PI))) (sqrt (/ 1 PI))) (* (* (+ (* (/ 1 (pow (fabs x) (/ 7 2))) (/ 15/8 (pow (fabs x) (/ 7 2)))) (+ (+ (/ 1/2 (* (* (fabs x) (fabs x)) (fabs x))) (/ 3/4 (pow (fabs x) 5))) (/ 1 (fabs x)))) (+ (* (/ 1 (pow (fabs x) (/ 7 2))) (/ 15/8 (pow (fabs x) (/ 7 2)))) (+ (+ (/ 1/2 (* (* (fabs x) (fabs x)) (fabs x))) (/ 3/4 (pow (fabs x) 5))) (/ 1 (fabs x))))) (+ (* (/ 1 (pow (fabs x) (/ 7 2))) (/ 15/8 (pow (fabs x) (/ 7 2)))) (+ (+ (/ 1/2 (* (* (fabs x) (fabs x)) (fabs x))) (/ 3/4 (pow (fabs x) 5))) (/ 1 (fabs x))))))
1.7s
(* (* (* (* (/ 1 (sqrt PI)) (sqrt (+ (+ (+ (* 2 (fabs x)) (* (/ 2 3) (* (* (fabs x) (fabs x)) (fabs x)))) (* (/ 1 5) (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)))) (* (/ 1 21) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))))) (* (/ 1 (sqrt PI)) (sqrt (+ (+ (+ (* 2 (fabs x)) (* (/ 2 3) (* (* (fabs x) (fabs x)) (fabs x)))) (* (/ 1 5) (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)))) (* (/ 1 21) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x))))))) (* (/ 1 (sqrt PI)) (sqrt (+ (+ (+ (* 2 (fabs x)) (* (/ 2 3) (* (* (fabs x) (fabs x)) (fabs x)))) (* (/ 1 5) (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)))) (* (/ 1 21) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x))))))) (* (* (sqrt (+ (+ (+ (* 2 (fabs x)) (* (/ 2 3) (* (* (fabs x) (fabs x)) (fabs x)))) (* (/ 1 5) (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)))) (* (/ 1 21) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x))))) (sqrt (+ (+ (+ (* 2 (fabs x)) (* (/ 2 3) (* (* (fabs x) (fabs x)) (fabs x)))) (* (/ 1 5) (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)))) (* (/ 1 21) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))))) (sqrt (+ (+ (+ (* 2 (fabs x)) (* (/ 2 3) (* (* (fabs x) (fabs x)) (fabs x)))) (* (/ 1 5) (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)))) (* (/ 1 21) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)))))))
1.6s
(/ (* (* (* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i)))) (* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))))) (* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))))) (* (* (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0)) (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))))

sample9.2m

Algorithm
87×intervals
Results
6.6m101844×body10240exit
1.1m388126×body80valid
26.3s22603×body1280valid
13.4s17309×body640valid
10.7s68966×body80nan
6.2s1598×body1280nan
4.5s9003×body320valid
4.3s1360×body640nan
1.8s4767×body160valid
1.6s1463×body2560valid
1.5s33536×pre80true
1.2s656×body320nan
618.0ms607×body5120valid
452.0ms338×body160nan

prune4.0m

rewrite3.4m

Algorithm
154×rewrite-expression-head
Calls

485 calls:

Slowest
25.5s
(* (+ (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (+ z 7))) (+ (+ (/ -0.13857109526572012 (- z -5)) (/ 12.507343278686905 (+ z 4))) (+ (+ (/ 771.3234287776531 (+ z 2)) (+ (/ 676.5203681218851 z) (+ 0.9999999999998099 (/ -1259.1392167224028 (- z -1))))) (/ -176.6150291621406 (+ z 3))))) (exp (- (+ (- z -6) 0.5))))
19.1s
(* (* (pow (+ (- z -6) 0.5) (- z (- 1 0.5))) (sqrt (* PI 2))) (* (+ (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (+ z 7))) (+ (+ (/ -0.13857109526572012 (- z -5)) (/ 12.507343278686905 (+ z 4))) (+ (+ (/ 771.3234287776531 (+ z 2)) (+ (/ 676.5203681218851 z) (+ 0.9999999999998099 (/ -1259.1392167224028 (- z -1))))) (/ -176.6150291621406 (+ z 3))))) (exp (- (+ (- z -6) 0.5)))))
18.5s
(* (* (/ (sqrt (+ (+ (+ (+ (+ 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) (/ (sqrt (+ (+ (+ (+ (+ 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)
18.2s
(* (/ (sqrt (+ (+ (+ (+ (+ 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) (/ (sqrt (+ (+ (+ (+ (+ 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))))))
16.4s
(* (/ (+ (+ (+ (+ (+ 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)
Rules
21088×*-un-lft-identity
17898×times-frac
9209×add-sqr-sqrt
7766×distribute-lft-out
7155×add-cube-cbrt
3660×add-exp-log
3239×pow1
3131×frac-add
3028×sqrt-prod
2910×add-cbrt-cube
2444×frac-times
2135×flip-+ flip3-+
1953×prod-exp
1832×associate-*r*
1580×associate-*l/
1468×associate-*l*
1402×cbrt-unprod
1320×associate-*r/
1281×pow-prod-down
876×associate-/r/
807×associate-/r*
704×add-log-exp
617×div-exp
613×associate-/l*
554×exp-neg
528×difference-of-squares
509×exp-diff
507×neg-sub0
501×pow-prod-up
485×insert-posit16
470×div-inv
454×distribute-lft-out--
451×sqrt-div
418×cbrt-prod
388×cbrt-undiv
355×unswap-sqr
284×cbrt-div
245×pow-sub
236×pow-plus
190×flip3-- flip--
175×sqrt-pow1
171×*-commutative
159×unpow-prod-down
156×log-pow
147×pow-sqr
135×swap-sqr pow1/3 1-exp
106×associate-/l/
105×pow1/2
99×rec-exp
95×log-prod
86×pow2
85×distribute-rgt-in sqr-pow distribute-lft-in
84×frac-2neg clear-num
63×pow-exp
62×sum-log sub-neg
57×diff-log
55×un-div-inv
54×frac-sub
51×+-commutative
50×pow-to-exp
47×exp-prod
36×rem-sqrt-square
33×cube-unmult
29×associate--l+
18×pow-flip div-sub
16×log-div
15×associate-+r+
12×pow-unpow pow-pow inv-pow
10×pow-div associate--r+
associate-+l+ pow3 associate--l- difference-of-sqr-1 rem-log-exp
sin-sum exp-sum sqrt-unprod rem-exp-log
distribute-rgt-out distribute-rgt1-in
sin-mult sub-div rem-square-sqrt
unpow-prod-up distribute-lft-neg-out distribute-lft-neg-in difference-cubes distribute-rgt-neg-in unpow1/3 associate--r- distribute-rgt-out--
sum-cubes associate-+l-
sqr-sin neg-mul-1 remove-posit16 distribute-neg-in rem-cbrt-cube cos-sum

series1.7m

Calls

485 calls:

Slowest
8.0s
(- (* 1 1) (/ (* (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (- (* 0.254829592 0.254829592) (* (/ (- (+ (/ 1.061405429 (* (+ 1 (* (fabs x) 0.3275911)) (* (+ 1 (* (fabs x) 0.3275911)) (+ 1 (* (fabs x) 0.3275911))))) (/ 1.421413741 (+ 1 (* (fabs x) 0.3275911)))) (+ 0.284496736 (/ 1.453152027 (* (+ 1 (* (fabs x) 0.3275911)) (+ 1 (* (fabs x) 0.3275911)))))) (* (+ 1 (* (fabs x) 0.3275911)) (+ 1 (* (fabs x) 0.3275911)))) (- (+ (/ 1.061405429 (* (+ 1 (* (fabs x) 0.3275911)) (* (+ 1 (* (fabs x) 0.3275911)) (+ 1 (* (fabs x) 0.3275911))))) (/ 1.421413741 (+ 1 (* (fabs x) 0.3275911)))) (+ 0.284496736 (/ 1.453152027 (* (+ 1 (* (fabs x) 0.3275911)) (+ 1 (* (fabs x) 0.3275911))))))))) 1) (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))) (* (- 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429)))))))) (exp (* (fabs x) (fabs x))))))
8.0s
(- (* 1 1) (/ (* (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (- (* 0.254829592 0.254829592) (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))) (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429)))))))))) 1) (* (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x)))))) (* (- 0.254829592 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1 (+ 1 (* 0.3275911 (fabs x)))) 1.061405429)))))))) (exp (* (fabs x) (fabs x))))))
3.3s
(/ (* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i)))) (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0)))
3.1s
(* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (/ (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i))) (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))))
2.9s
(/ (/ (* (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2 i))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2 i)))) (+ (+ (+ alpha beta) (* 2 i)) (sqrt 1.0))) (- (+ (+ alpha beta) (* 2 i)) (sqrt 1.0)))

bsearch16.1s

regimes11.6s

Accuracy

Total 73.5b remaining (31.6%)

9.9b58.3%math.log10 on complex, real part
9.6b59.2%math.log/1 on complex, real part
9.6b56.6%math.abs on complex
9.4b59.3%math.log/2 on complex, real part
8.6b48%math.sqrt on complex, real part

localize8.1s

end1.0ms