Details

rewrite10.6m (49.7%)

Algorithm
144×rewrite-expression-head
Calls

426 calls:

2.2m
(* (/ (* (pow (+ (- z -6) 0.5) (+ (- z 1) 0.5)) (sqrt (* PI 2))) (exp (+ (- z -6) 0.5))) (+ (+ (+ (+ (/ -0.13857109526572012 (+ 6 (- z 1))) (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ 771.3234287776531 (+ z 2))) (/ -1259.1392167224028 (+ z 1)))) (/ -176.6150291621406 (- z -3))) (/ 12.507343278686905 (- z -4))) (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (+ z 7)))))
1.5m
(* (/ (pow (+ (- z -6) 0.5) (+ (- z 1) 0.5)) (/ (exp (+ (- z -6) 0.5)) (sqrt (* PI 2)))) (+ (+ (+ (+ (/ -0.13857109526572012 (+ 6 (- z 1))) (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ 771.3234287776531 (+ z 2))) (/ -1259.1392167224028 (+ z 1)))) (/ -176.6150291621406 (- z -3))) (/ 12.507343278686905 (- z -4))) (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (+ z 7)))))
1.5m
(* (/ (pow (sqrt (+ (- z -6) 0.5)) (+ (- z 1) 0.5)) (/ (exp (+ (- z -6) 0.5)) (sqrt (* PI 2)))) (+ (+ (+ (+ (/ -0.13857109526572012 (+ 6 (- z 1))) (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ 771.3234287776531 (+ z 2))) (/ -1259.1392167224028 (+ z 1)))) (/ -176.6150291621406 (- z -3))) (/ 12.507343278686905 (- z -4))) (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (+ z 7)))))
1.5m
(* (/ (sqrt (* PI 2)) (exp 0.5)) (+ (+ (+ (+ (/ -0.13857109526572012 (+ 6 (- z 1))) (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ 771.3234287776531 (+ z 2))) (/ -1259.1392167224028 (+ z 1)))) (/ -176.6150291621406 (- z -3))) (/ 12.507343278686905 (- z -4))) (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (+ z 7)))))
1.2m
(* (/ (pow (sqrt (+ (- z -6) 0.5)) (+ (- z 1) 0.5)) 1) (* (/ (pow (sqrt (+ (- z -6) 0.5)) (+ (- z 1) 0.5)) (/ (exp (+ (- z -6) 0.5)) (sqrt (* PI 2)))) (+ (+ (+ (+ (/ -0.13857109526572012 (+ 6 (- z 1))) (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 z)) (/ 771.3234287776531 (+ z 2))) (/ -1259.1392167224028 (+ z 1)))) (/ -176.6150291621406 (- z -3))) (/ 12.507343278686905 (- z -4))) (+ (/ 9.984369578019572e-06 (- z -6)) (/ 1.5056327351493116e-07 (+ z 7))))))
Rules
15232×times-frac
12148×*-un-lft-identity
10497×add-sqr-sqrt
7651×add-cube-cbrt
5547×sqrt-prod
2949×add-exp-log
2260×frac-add
2126×associate-*l*
2023×distribute-lft-out
1966×pow1
1786×add-cbrt-cube
1572×associate-*r*
1281×prod-exp
1107×prod-diff
798×flip-+ flip3-+
796×cbrt-prod
787×associate-/l*
755×unpow-prod-down
694×div-exp
619×associate-*r/
570×difference-of-squares
568×frac-times
566×add-log-exp
556×exp-sum
536×cbrt-unprod
466×unswap-sqr
460×pow-prod-up
446×associate-/r*
432×log1p-expm1-u
429×div-inv
426×expm1-log1p-u insert-posit16
420×associate-/r/
415×distribute-lft-out--
377×cbrt-undiv
348×swap-sqr
341×pow-prod-down
245×pow-exp
239×log-pow
237×pow-plus
211×associate-*l/
188×distribute-rgt-in distribute-lft-in
155×flip3-- flip--
142×sqr-pow
140×pow-to-exp
136×sqrt-div
131×*-commutative
127×fma-neg
115×unpow-prod-up
114×1-exp
106×pow-sqr
101×sqrt-pow1
92×sub-neg
86×frac-2neg rec-exp clear-num
84×cbrt-div pow2
80×fma-def pow1/2
77×associate-/l/
75×associate--r+
72×fma-udef
66×log-prod
53×pow1/3
48×associate--l+
45×diff-log
32×exp-prod
30×sum-log
29×rem-sqrt-square
28×associate-+l+
24×cube-unmult
23×+-commutative div-sub
22×frac-sub
20×distribute-rgt1-in distribute-rgt-out
17×un-div-inv
15×pow-flip hypot-udef
12×inv-pow hypot-def
10×exp-diff associate-+l-
pow-unpow cos-sum
expm1-udef pow-div
expm1-log1p pow3 sqrt-unprod rem-log-exp log-div
pow-pow rem-exp-log rem-square-sqrt
sqrt-undiv log1p-expm1 pow-sub log1p-udef
associate--l- sin-sum associate-+r+ difference-of-sqr-1 expm1-def
associate--r-
distribute-rgt-out--

sample4.4m (20.5%)

Algorithm
99×intervals
Results
1.4m12886×body10240exit
1.1m434996×body80valid
37.1s28295×body1280valid
19.8s21277×body640valid
9.3s69653×body80nan
7.4s2109×body1280nan
6.5s11164×body320valid
4.9s1592×body640nan
4.8s9533×body160valid
2.1s49536×pre80true
1.5s869×body320nan
1.4s1416×body2560valid
605.0ms642×body5120valid
572.0ms423×body160nan

simplify2.5m (11.6%)

prune2.2m (10.2%)

series1.4m (6.7%)

Calls

426 calls:

1.7s
(- 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 1) (* (* 0.3275911 (fabs x)) (* 0.3275911 (fabs x))))) (- 1 (* 0.3275911 (fabs x)))) 1.061405429))))))))) (exp (- (* (fabs x) (fabs x))))))
1.5s
(* (pow (* -2 (log u1)) 0.5) (* (pow (* -2 (log u1)) 0.5) 1/216))
1.2s
(pow (* (pow (log u1) 1.0) (pow -2 1.0)) 0.5)
1.2s
(- 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))))))
1.2s
(* (* (pow (* -2 (log u1)) 0.5) (sqrt 1/6)) (sqrt 1/6))

regimes8.6s (0.7%)

Accuracy

Total 30.0b remaining (21.8%)

Threshold costs 5.1b (3.7%)

7.6b4.9%Octave 3.8, jcobi/2
6.8b37.1%_divideComplex, imaginary part
3.2b60.9%Octave 3.8, jcobi/4
3.2b74.6%Octave 3.8, jcobi/1
1.7b0%Jmat.Real.lambertw, newton loop step

localize6.1s (0.5%)

bsearch1.4s (0.1%)

end1.0ms (0.0%)