Date:Tuesday, February 25th, 2020
Commit:7c13946c on findroot
Hostname:warfa with Racket 7.5
Seed:2020056
Parameters:256 points for 4 iterations
Flags:
setup:simplifyrules:arithmeticrules:polynomialsrules:fractionsrules:exponentsrules:trigonometryrules:hyperbolicrules:complexrules:specialrules:boolsrules:branchesgenerate:rrgenerate:taylorgenerate:simplifyreduce:regimesreduce:avg-errorreduce:binary-searchreduce:branch-expressionsprecision:doubleprecision:fallback
default

Details

rewrite12.1m (46.5%)

Algorithm
162×rewrite-expression-head
Calls

470 calls:

9.7m
(/ (* (* (sqrt (* PI 2.0)) (pow (+ (+ (- z 1.0) 7.0) 0.5) (+ (- z 1.0) 0.5))) (+ (+ (+ (+ (+ (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ -1259.1392167224028 (+ (- z 1.0) 2.0))) (/ 771.3234287776531 (+ (- z 1.0) 3.0))) (/ -176.6150291621406 (+ (- z 1.0) 4.0))) (/ 12.507343278686905 (+ (- z 1.0) 5.0))) (/ -0.13857109526572012 (+ (- z 1.0) 6.0))) (/ 9.984369578019572e-06 (+ (- z 1.0) 7.0))) (/ 1.5056327351493116e-07 (+ (- z 1.0) 8.0)))) (exp (+ (+ (- z 1.0) 7.0) 0.5)))
1.2m
(* (* (* (sqrt (+ (+ (+ (+ (+ 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))))) (sqrt (/ 1 (* (+ (+ (+ (* (* (pow x 2) (* x (pow x 3))) 0.0694555761) 1.0) (* (pow x 2) (+ 0.7715471019 (* 0.2909738639 (pow x 2))))) (+ (* 0.0140005442 (pow x 8)) (* (pow x 2) (+ (* 0.0008327945 (* (* (* (* x x) (* x x)) (* x x)) (* x x))) (* (* 2.0 0.0001789971) (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x))))))) 1)))) (* (sqrt (+ (+ (+ (+ (+ 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))))) (sqrt (/ 1 (* (+ (+ (+ (* (* (pow x 2) (* x (pow x 3))) 0.0694555761) 1.0) (* (pow x 2) (+ 0.7715471019 (* 0.2909738639 (pow x 2))))) (+ (* 0.0140005442 (pow x 8)) (* (pow x 2) (+ (* 0.0008327945 (* (* (* (* x x) (* x x)) (* x x)) (* x x))) (* (* 2.0 0.0001789971) (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x))))))) 1))))) x)
10.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 (* (+ (+ (+ (* (* (pow x 2) (* x (pow x 3))) 0.0694555761) 1.0) (* (pow x 2) (+ 0.7715471019 (* 0.2909738639 (pow x 2))))) (+ (* 0.0140005442 (pow x 8)) (* (pow x 2) (+ (* 0.0008327945 (* (* (* (* x x) (* x x)) (* x x)) (* x x))) (* (* 2.0 0.0001789971) (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x))))))) 1))) x)
6.6s
(* (/ (+ (+ (+ (+ (+ 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)
4.1s
(* (+ (+ (+ (+ (+ 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 (* (+ (+ (+ (* (* (pow x 2) (* x (pow x 3))) 0.0694555761) 1.0) (* (pow x 2) (+ 0.7715471019 (* 0.2909738639 (pow x 2))))) (+ (* 0.0140005442 (pow x 8)) (* (pow x 2) (+ (* 0.0008327945 (* (* (* (* x x) (* x x)) (* x x)) (* x x))) (* (* 2.0 0.0001789971) (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x))))))) 1)))
Rules
16958×times-frac
12689×*-un-lft-identity
10127×add-sqr-sqrt
7683×add-cube-cbrt
3944×add-exp-log
2625×add-cbrt-cube
2286×pow1
2202×prod-exp
1957×distribute-lft-out
1898×sqrt-prod
1747×associate-*l*
1200×associate-*r*
1130×cbrt-unprod
1104×flip-+ flip3-+
934×associate-/r/
915×pow-prod-down
812×associate-*r/
736×add-log-exp
725×frac-times
722×difference-of-squares
691×associate-/l*
669×sqrt-div
606×div-inv
571×div-exp
554×distribute-lft-out--
550×associate-/r*
503×associate-*l/
482×cbrt-div
470×pow-exp
429×cbrt-undiv
366×unswap-sqr
318×unpow-prod-down
282×cbrt-prod
251×log-pow
249×pow-to-exp
236×frac-add
233×sqr-pow
202×flip3-- flip--
198×cube-prod
168×pow-unpow
161×*-commutative
146×sqrt-pow1
137×distribute-rgt-in associate-/l/ distribute-lft-in
115×pow1/2 swap-sqr
103×sub-neg
98×pow-prod-up log-prod
91×diff-log
86×clear-num frac-2neg
74×pow1/3
70×pow-sqr
60×1-exp
58×sum-log
55×rem-sqrt-square
54×pow-plus pow-pow
51×rec-exp
39×associate--l+
37×pow2 +-commutative
29×exp-neg exp-prod
28×exp-diff
27×neg-sub0
26×pow-flip
25×inv-pow
24×cube-div distribute-neg-frac
20×unpow3 associate-+r+ cube-mult
17×associate-+l+
16×associate--l- distribute-rgt-neg-in
15×rem-cube-cbrt div-sub
12×un-div-inv distribute-lft-neg-in sin-sum associate--r+
9×unpow2
7×associate-+l-
6×rem-exp-log distribute-neg-in
5×exp-sum rem-log-exp log-div
4×sqrt-unprod neg-mul-1 frac-sub rem-cbrt-cube
3×sqrt-pow2 difference-of-sqr-1 distribute-rgt-neg-out
2×sqrt-undiv difference-cubes pow3 rem-square-sqrt neg-log
1×cos-sum unpow-prod-up pow-sub distribute-frac-neg associate-+r- log-rec cos-mult associate--r- exp-to-pow

series10.8m (41.6%)

Calls

463 calls:

20.0s
(- (pow 1.0 3) (/ (pow (* (* 1.0 (- (* 0.254829592 0.254829592) (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))) (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429)))))))))) (exp (- (* (fabs x) (fabs x))))) 3) (pow (* (+ 1.0 (* 0.3275911 (fabs x))) (- 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) 3)))
3.2s
(- (pow 1.0 3) (pow (* (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 0.254829592 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -0.284496736 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ 1.421413741 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) (+ -1.453152027 (* (/ 1.0 (+ 1.0 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))) 3))
2.2s
(/ (* (/ (/ i (/ (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) i))) (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2.0 i)))) (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0)))
2.2s
(/ (* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2.0 i))) (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2.0 i)))) (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0)))
2.2s
(* (/ (/ (* i (+ (+ alpha beta) i)) (+ (+ alpha beta) (* 2.0 i))) (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (/ (+ (* beta alpha) (* i (+ (+ alpha beta) i))) (+ (+ alpha beta) (* 2.0 i))) (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))))

sample1.6m (6.3%)

Algorithm
97×intervals
Results
44.4s11865×body10240exit
16.9s434049×body80valid
11.0s27353×body1280valid
5.1s21229×body640valid
2.5s69933×body80nan
1.7s1629×body1280nan
1.6s10546×body320valid
975.0ms1344×body640nan
733.0ms5579×body160valid
656.0ms49536×pre80true
342.0ms659×body320nan
341.0ms1394×body2560valid
145.0ms593×body5120valid
112.0ms326×body160nan

prune1.0m (3.9%)

simplify14.5s (0.9%)

Algorithm
256×egg-herbie

regimes6.4s (0.4%)

Accuracy

Total 116.1b remaining (48.3%)

Threshold costs 3.6b (1.5%)

23.9b35.3%math.sqrt on complex, imaginary part, im greater than 0 branch
11.0b56.1%math.log10 on complex, real part
10.9b56.2%math.log/2 on complex, real part
10.8b56.5%math.log/1 on complex, real part
10.3b57.3%math.abs on complex

bsearch2.7s (0.2%)

localize2.6s (0.2%)

end1.0ms (0.0%)