Details

rewrite12.4m (52.2%)

Algorithm
157×rewrite-expression-head
Calls

480 calls:

8.9m
(/ (+ (/ (/ 5688119651184367.0 3.777893186295716e+22) (+ (- z 1.0) 8.0)) (+ (/ (/ 5893731530977871.0 5.902958103587057e+20) (+ (- z 1.0) 7.0)) (+ (/ (/ -4992549864024149.0 3.602879701896397e+16) (+ (- z 1.0) 6.0)) (+ (/ (/ 7041008316161153.0 562949953421312.0) (+ (- z 1.0) 5.0)) (+ (/ (/ -6214088902520669.0 35184372088832.0) (+ (- z 1.0) 4.0)) (+ (/ (/ 848079078717083.0 1099511627776.0) (+ (- z 1.0) 3.0)) (+ (/ (/ -5537752839100187.0 4398046511104.0) (+ (- z 1.0) 2.0)) (+ (/ 5950736089418501.0 (* (- z 0) 8796093022208.0)) (/ 562949953421205.0 562949953421312.0))))))))) (/ (exp (+ (+ (- z 1.0) 7.0) (/ 1.0 2.0))) (* (sqrt (* PI 2.0)) (pow (+ (+ (- z 1.0) 7.0) (/ 1.0 2.0)) (+ (- z 1.0) (/ 1.0 2.0))))))
30.1s
(/ (* x (+ (+ (* (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3))))) (/ 825478423409049.0 4.611686018427388e+18)) (+ (+ (/ (* 8375303961237363.0 (* (pow x 2) (* x (pow x 3)))) 1.152921504606847e+18) 1.0) (* (pow x 2) (+ (/ 3782789308857969.0 3.602879701896397e+16) (* (/ 763919671262763.0 18014398509481984.0) (pow x 2)))))) (/ (* 583843369866023.0 (* (pow x 2) (* (pow x 2) (* x (pow x 3))))) 1.152921504606847e+18))) (+ (+ (* (* (pow x 2) (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3)))))) (* 2.0 (/ 825478423409049.0 4.611686018427388e+18))) (+ (+ (/ (* 2502400853142105.0 (* (pow x 2) (* x (pow x 3)))) 3.602879701896397e+16) 1.0) (* (pow x 2) (+ (/ 1737369620307813.0 2251799813685248.0) (* (/ 2620859570069187.0 9007199254740992.0) (pow x 2)))))) (* (pow x 2) (+ (* (/ 2017691060547333.0 1.4411518807585587e+17) (* (* (* x x) x) (* (* x x) x))) (* (/ 7681173503746455.0 9.223372036854776e+18) (* (* (* (* x x) (* x x)) (* x x)) (* x x)))))))
23.7s
(* (/ (+ (+ (+ (+ (+ 1.0 (* 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 (* 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 0.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x))))) x)
13.8s
(/ (+ (+ (+ (+ (+ 1.0 (* 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 (* 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 0.0001789971) (* (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)) (* x x)))))
3.8s
(- (sqrt (+ 1.0 (+ (* (/ 5125037576950889.0 18014398509481984.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 2))) (* (/ 6544414927309881.0 4503599627370496.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 4)))))) (sqrt (+ (* (/ 4780145094533421.0 4503599627370496.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 5))) (+ (* (/ 6401478394306903.0 4503599627370496.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 3))) (* (/ 2295300911148351.0 9007199254740992.0) (/ (exp (- (pow (fabs x) 2))) (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0)))))))
Rules
7814×times-frac
7042×*-un-lft-identity
6523×add-sqr-sqrt
4774×add-cube-cbrt
3982×sqrt-prod
2497×add-exp-log
1812×add-cbrt-cube
1567×pow1
1378×frac-add
1267×associate-*l*
1190×associate-*r/
1134×frac-times
1072×associate-*l/
987×distribute-lft-out
985×add-log-exp
979×prod-exp
965×associate-*r*
952×flip-+ flip3-+
922×cbrt-div
878×frac-sub
637×cbrt-unprod
632×cbrt-prod
609×associate-/l*
507×div-exp
460×associate-/r*
454×sqrt-div
435×pow-prod-down
427×associate-/r/
399×unpow-prod-down
398×difference-of-squares flip3-- flip--
315×div-inv unswap-sqr
304×cbrt-undiv
303×distribute-lft-out--
250×log-pow
185×sum-log
170×swap-sqr
167×log-prod sqrt-pow1
160×sqr-pow
152×associate-/l/
142×pow-exp
140×cube-prod
138×diff-log
135×pow1/3
128×*-commutative
125×pow-unpow
124×pow-to-exp
107×pow1/2
97×distribute-rgt-neg-in distribute-lft-neg-in sub-neg
95×distribute-rgt-in distribute-lft-in
84×frac-2neg clear-num
56×rem-sqrt-square
54×exp-prod
53×exp-sum
50×pow-prod-up 1-exp
49×rec-exp associate--l+
48×+-commutative
46×pow-pow
42×pow-sqr
39×unpow3 cube-mult
28×log-div
26×cube-div
21×inv-pow pow-flip pow-plus associate-+l- pow2
20×div-sub
18×associate-+l+
17×un-div-inv
16×pow-sub
15×rem-exp-log associate--l-
14×exp-diff
10×associate--r+ rem-log-exp
rem-cube-cbrt sin-sum associate-+r+ neg-mul-1
rem-cbrt-cube
sqrt-undiv unpow2 unpow-prod-up neg-sub0
sqrt-pow2 sqrt-unprod exp-neg rem-square-sqrt distribute-frac-neg
cos-sum exp-sqrt associate-+r- distribute-rgt-out-- unpow1/2 neg-log sum-cubes

simplify4.4m (18.4%)

sample3.2m (13.5%)

Algorithm
95×intervals
Results
1.2m11786×body10240exit
43.7s415411×body80valid
19.6s22424×body1280valid
9.1s17815×body640valid
6.7s69676×body80nan
3.8s1551×body1280nan
3.7s9332×body320valid
2.7s1286×body640nan
2.1s8759×body160valid
1.9s49536×pre80true
831.0ms669×body320nan
387.0ms747×body2560valid
358.0ms300×body160nan
200.0ms272×body5120valid

prune2.0m (8.2%)

series1.5m (6.2%)

Calls

478 calls:

1.7s
(cbrt (- (+ 1.0 (+ (* (/ 5125037576950889.0 18014398509481984.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 2))) (* (/ 6544414927309881.0 4503599627370496.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 4))))) (+ (* (/ 4780145094533421.0 4503599627370496.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 5))) (+ (* (/ 6401478394306903.0 4503599627370496.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 3))) (* (/ 2295300911148351.0 9007199254740992.0) (/ (exp (- (pow (fabs x) 2))) (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0)))))))
1.7s
(cbrt (- (+ 1.0 (+ (* (/ 5125037576950889.0 18014398509481984.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 2))) (* (/ 6544414927309881.0 4503599627370496.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 4))))) (+ (* (/ 4780145094533421.0 4503599627370496.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 5))) (+ (* (/ 6401478394306903.0 4503599627370496.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 3))) (* (/ 2295300911148351.0 9007199254740992.0) (/ (exp (- (pow (fabs x) 2))) (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0)))))))
1.6s
(/ 6.0 (pow (* -2.0 (log u1)) (/ 1.0 2.0)))
1.5s
(- (sqrt (+ 1.0 (+ (* (/ 5125037576950889.0 18014398509481984.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 2))) (* (/ 6544414927309881.0 4503599627370496.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 4)))))) (sqrt (+ (* (/ 4780145094533421.0 4503599627370496.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 5))) (+ (* (/ 6401478394306903.0 4503599627370496.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 3))) (* (/ 2295300911148351.0 9007199254740992.0) (/ (exp (- (pow (fabs x) 2))) (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0)))))))
1.3s
(- (+ 1.0 (+ (* (/ 5125037576950889.0 18014398509481984.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 2))) (* (/ 6544414927309881.0 4503599627370496.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 4))))) (+ (* (/ 4780145094533421.0 4503599627370496.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 5))) (+ (* (/ 6401478394306903.0 4503599627370496.0) (/ (exp (- (pow (fabs x) 2))) (pow (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0) 3))) (* (/ 2295300911148351.0 9007199254740992.0) (/ (exp (- (pow (fabs x) 2))) (+ (* (/ 1475339155889891.0 4503599627370496.0) (fabs x)) 1.0))))))

regimes10.8s (0.8%)

Accuracy

Total 74.5b remaining (33.2%)

Threshold costs 5.0b (2.2%)

10.4b56.6%math.log/2 on complex, real part
10.2b57.7%math.log10 on complex, real part
10.0b58.5%math.log/1 on complex, real part
9.6b58.6%math.abs on complex
8.7b59%math.sqrt on complex, real part

localize6.2s (0.4%)

bsearch3.2s (0.2%)

end1.0ms (0.0%)