| 159× | rewrite-expression-head |
456 calls:
| 9.6m | (/ (* (* (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))) |
| 3.4m | (/ (/ (* (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))))) (hypot (pow (* 0.0001789971 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x))) 3/2) (pow (+ (+ (+ (+ 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)))) 3/2))) (fma (* 2.0 0.0001789971) (* (pow x 2) (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3)))))) (fma 0.0008327945 (* (pow x 2) (* (pow x 2) (* (pow x 2) (* x (pow x 3))))) (fma 0.0140005442 (* (pow x 2) (* (pow x 2) (* x (pow x 3)))) (fma 0.0694555761 (pow x 6) (fma 0.2909738639 (pow x 4) (fma (* 0.7715471019 x) x 1.0))))))) (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)))) (+ (+ (+ (+ 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))) (* 0.0001789971 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (* (+ (+ (+ (+ 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)))))))) |
| 34.6s | (- (* (* 0.0001789971 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x))) (* 0.0001789971 (* (* (* (* (* x x) (* x x)) (* x x)) (* x x)) (* x x)))) (* (+ (+ (+ (+ 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))))) |
| 7.3s | (* (/ (+ (+ (+ (+ (+ 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.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))))) |
| 13211× | times-frac |
| 11447× | *-un-lft-identity |
| 10060× | add-sqr-sqrt |
| 7890× | add-cube-cbrt |
| 2617× | sqrt-prod |
| 2521× | add-exp-log |
| 2266× | prod-diff |
| 1940× | add-cbrt-cube |
| 1471× | associate-*r* |
| 1430× | pow1 |
| 1366× | distribute-lft-out |
| 1098× | associate-*l* |
| 827× | prod-exp |
| 734× | associate-/l* |
| 711× | add-log-exp |
| 709× | unpow-prod-down |
| 663× | div-exp |
| 651× | associate-/r/ |
| 647× | cbrt-prod |
| 616× | difference-of-squares |
| 571× | cbrt-unprod |
| 563× | associate-/r* |
| 534× | fma-def |
| 522× | cube-prod |
| 515× | div-inv |
| 476× | distribute-lft-out-- |
| 451× | log1p-expm1-u |
| 450× | expm1-log1p-u |
| 439× | cbrt-undiv |
| 357× | pow-prod-down |
| 356× | flip-+ flip3-+ |
| 350× | log-pow |
| 337× | distribute-rgt-in distribute-lft-in |
| 315× | fma-neg |
| 288× | unswap-sqr |
| 215× | sqr-pow |
| 207× | distribute-rgt-neg-in |
| 198× | associate--l+ |
| 197× | exp-prod |
| 195× | distribute-lft-neg-in |
| 184× | flip3-- flip-- |
| 183× | associate-*r/ |
| 180× | sqrt-div |
| 174× | sqrt-pow1 |
| 151× | frac-times |
| 131× | *-commutative |
| 129× | associate-*l/ |
| 122× | exp-sum |
| 116× | unpow3 cube-mult |
| 113× | 1-exp log-prod |
| 108× | pow-exp |
| 100× | cube-div |
| 94× | sub-neg |
| 93× | associate-/l/ clear-num frac-2neg |
| 91× | pow-unpow |
| 88× | diff-log |
| 87× | pow1/2 |
| 76× | swap-sqr |
| 75× | rec-exp neg-mul-1 |
| 73× | associate--r+ |
| 62× | pow1/3 |
| 59× | pow-to-exp cbrt-div |
| 57× | exp-neg |
| 55× | pow-prod-up |
| 54× | neg-sub0 pow-pow |
| 51× | exp-diff |
| 45× | rem-sqrt-square |
| 43× | fma-udef |
| 42× | exp-to-pow |
| 41× | pow-sqr |
| 36× | sum-log |
| 34× | rem-cube-cbrt |
| 25× | un-div-inv |
| 24× | pow-plus |
| 22× | frac-add +-commutative |
| 21× | associate--l- |
| 20× | hypot-def |
| 18× | distribute-neg-in |
| 17× | associate-+l+ |
| 16× | pow2 |
| 12× | associate-+r+ div-sub |
| 10× | inv-pow pow-flip |
| 9× | unpow2 |
| 7× | rem-log-exp |
| 5× | rem-exp-log associate-+l- neg-log |
| 4× | pow-div sin-sum expm1-log1p frac-sub expm1-udef |
| 3× | sqrt-unprod distribute-lft-neg-out rem-cbrt-cube |
| 2× | sqrt-undiv difference-cubes sqrt-pow2 log1p-udef pow3 hypot-udef rem-square-sqrt distribute-rgt-neg-out log-div |
| 1× | log1p-expm1 unpow-prod-up pow-sub unsub-neg distribute-rgt-out cos-sum associate--r- distribute-rgt-out-- |
450 calls:
| 5.7m | (- (pow (pow 1.0 3) 3) (pow (log1p (expm1 (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))) 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.0s | (/ (* (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0) (fma i 2.0 (+ alpha beta))) (/ (fma beta alpha (* i (+ (+ alpha beta) i))) (fma i 2.0 (+ alpha beta)))) |
| 1.9s | (/ (* (* (- (* (+ (+ alpha beta) (* 2.0 i)) (+ (+ alpha beta) (* 2.0 i))) 1.0) (fma i 2.0 (+ alpha beta))) (fma i 2.0 (+ alpha beta))) (fma beta alpha (* i (+ (+ alpha beta) i)))) |
| 987.0ms | (- (pow 1.0 3) (log1p (expm1 (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)))) |
| 98× | intervals |
| 44.8s | 7687× | body | 10240 | exit |
| 33.5s | 454135× | body | 80 | valid |
| 14.1s | 27267× | body | 1280 | valid |
| 6.9s | 21054× | body | 640 | valid |
| 3.8s | 68309× | body | 80 | nan |
| 2.7s | 1644× | body | 1280 | nan |
| 2.3s | 10373× | body | 320 | valid |
| 1.7s | 9489× | body | 160 | valid |
| 1.6s | 1308× | body | 640 | nan |
| 969.0ms | 49536× | pre | 80 | true |
| 545.0ms | 1427× | body | 2560 | valid |
| 521.0ms | 601× | body | 320 | nan |
| 259.0ms | 4306× | body | 80 | overflowed |
| 247.0ms | 598× | body | 5120 | valid |
| 197.0ms | 344× | body | 160 | nan |
| 255× | egg-herbie |
Total 44.3b remaining (40.4%)
Threshold costs 5.5b (5.1%)
| 8.2b | 34.5% | _divideComplex, imaginary part |
| 7.0b | 36.3% | _divideComplex, real part |
| 5.5b | 14.9% | Octave 3.8, jcobi/2 |
| 5.1b | 56.7% | Octave 3.8, jcobi/4 |
| 4.4b | 23.1% | math.sqrt on complex, imaginary part, im greater than 0 branch |