| 158× | rewrite-expression-head |
477 calls:
| 2.3m | (* (+ (+ (+ (+ (+ (/ 676.5203681218851 (- 1.0 z)) (/ -1259.1392167224028 (+ 2.0 (- z)))) (/ 771.3234287776531 (+ 3.0 (- z)))) 0.9999999999998099) (/ 12.507343278686905 (+ (- z) 5.0))) (+ (/ -176.6150291621406 (+ (- z) 4.0)) (+ (+ (/ -0.13857109526572012 (+ 6.0 (- z))) (/ 1.5056327351493116e-07 (+ 8.0 (- z)))) (/ 9.984369578019572e-06 (+ 7.0 (- z)))))) (sqrt (* 2.0 PI))) |
| 1.7m | (* (+ (+ (+ (+ (+ (/ 676.5203681218851 (- 1.0 z)) (/ -1259.1392167224028 (+ 2.0 (- z)))) (/ 771.3234287776531 (+ 3.0 (- z)))) 0.9999999999998099) (/ 12.507343278686905 (+ (- z) 5.0))) (+ (/ -176.6150291621406 (+ (- z) 4.0)) (+ (+ (/ -0.13857109526572012 (+ 6.0 (- z))) (/ 1.5056327351493116e-07 (+ 8.0 (- z)))) (/ 9.984369578019572e-06 (+ 7.0 (- z)))))) (sqrt (* 2.0 PI))) |
| 38.9s | (* (/ (sqrt (+ (* (* (* x x) (* x x)) (+ (+ (* (* (* x x) (* (* x x) (* x x))) 0.0001789971) (* 0.0005064034 (* (* x x) (* x x)))) 0.0424060604)) (+ (* (* x x) 0.1049934947) (+ (* 0.0072644182 (* (* x x) (* (* x x) (* x x)))) 1.0)))) (/ (+ 1.0 (+ (* (* x x) (+ (+ (* (* x x) 0.2909738639) (* 0.0694555761 (* (* x x) (* x x)))) 0.7715471019)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* x x) (+ 0.0008327945 (* 0.0003579942 (* x x)))) 0.0140005442)))) (sqrt (+ (* (* (* x x) (* x x)) (+ (+ (* (* (* x x) (* (* x x) (* x x))) 0.0001789971) (* 0.0005064034 (* (* x x) (* x x)))) 0.0424060604)) (+ (* (* x x) 0.1049934947) (+ (* 0.0072644182 (* (* x x) (* (* x x) (* x x)))) 1.0)))))) x) |
| 23.7s | (/ (sqrt (+ (* (* (* x x) (* x x)) (+ (+ (* (* (* x x) (* (* x x) (* x x))) 0.0001789971) (* 0.0005064034 (* (* x x) (* x x)))) 0.0424060604)) (+ (* (* x x) 0.1049934947) (+ (* 0.0072644182 (* (* x x) (* (* x x) (* x x)))) 1.0)))) (/ (+ 1.0 (+ (* (* x x) (+ (+ (* (* x x) 0.2909738639) (* 0.0694555761 (* (* x x) (* x x)))) 0.7715471019)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* x x) (+ 0.0008327945 (* 0.0003579942 (* x x)))) 0.0140005442)))) (sqrt (+ (* (* (* x x) (* x x)) (+ (+ (* (* (* x x) (* (* x x) (* x x))) 0.0001789971) (* 0.0005064034 (* (* x x) (* x x)))) 0.0424060604)) (+ (* (* x x) 0.1049934947) (+ (* 0.0072644182 (* (* x x) (* (* x x) (* x x)))) 1.0)))))) |
| 21.7s | (* (/ (+ (* (* (* x x) (* x x)) (+ (+ (* (* (* x x) (* (* x x) (* x x))) 0.0001789971) (* 0.0005064034 (* (* x x) (* x x)))) 0.0424060604)) (+ (* (* x x) 0.1049934947) (+ (* 0.0072644182 (* (* x x) (* (* x x) (* x x)))) 1.0))) (+ 1.0 (+ (* (* x x) (+ (+ (* (* x x) 0.2909738639) (* 0.0694555761 (* (* x x) (* x x)))) 0.7715471019)) (* (* (* (* x x) (* x x)) (* (* x x) (* x x))) (+ (* (* x x) (+ 0.0008327945 (* (* (* x x) 2.0) 0.0001789971))) 0.0140005442))))) x) |
| 24755× | frac-times |
| 16185× | associate-*r/ |
| 13149× | times-frac |
| 11405× | *-un-lft-identity |
| 9706× | associate-*l/ |
| 8672× | add-sqr-sqrt |
| 7847× | exp-neg |
| 7645× | flip-+ flip3-+ |
| 6751× | frac-add |
| 6735× | exp-diff |
| 6730× | neg-sub0 |
| 6489× | add-cube-cbrt |
| 4017× | sqrt-prod |
| 3143× | add-exp-log |
| 2682× | pow1 |
| 2535× | associate-*r* |
| 2352× | add-cbrt-cube |
| 1994× | distribute-lft-out |
| 1419× | prod-exp |
| 1354× | associate-*l* |
| 1128× | un-div-inv |
| 847× | cbrt-unprod associate-/r/ |
| 762× | associate-/l* |
| 740× | difference-of-squares |
| 738× | add-log-exp |
| 644× | pow-prod-up |
| 601× | distribute-lft-out-- |
| 595× | pow-prod-down |
| 586× | div-exp |
| 453× | associate-/r* |
| 441× | cbrt-undiv |
| 439× | div-inv |
| 336× | sqrt-div |
| 321× | pow-plus |
| 320× | unswap-sqr |
| 294× | cbrt-prod |
| 249× | log-pow |
| 206× | swap-sqr |
| 177× | flip3-- flip-- |
| 176× | cbrt-div |
| 174× | *-commutative |
| 162× | unpow-prod-down |
| 161× | log-prod sqrt-pow1 |
| 144× | pow-sqr |
| 139× | associate-/l/ |
| 136× | pow1/2 |
| 124× | pow-exp |
| 109× | pow2 |
| 103× | sum-log |
| 102× | distribute-rgt-in distribute-lft-in |
| 94× | distribute-rgt-neg-in |
| 89× | frac-2neg clear-num |
| 85× | sqr-pow |
| 79× | pow1/3 |
| 66× | cube-prod pow-to-exp |
| 58× | sub-neg |
| 57× | +-commutative |
| 52× | cube-unmult |
| 49× | rem-sqrt-square |
| 47× | diff-log |
| 36× | 1-exp |
| 32× | exp-sum div-sub |
| 31× | rec-exp |
| 25× | exp-prod |
| 24× | cube-div |
| 23× | pow3 |
| 22× | frac-sub |
| 21× | pow-unpow |
| 20× | associate--l+ |
| 19× | log-div |
| 16× | associate--r+ |
| 13× | difference-cubes pow-pow |
| 12× | associate-+r+ rem-log-exp |
| 11× | associate-+r- |
| 10× | unpow-prod-up |
| 9× | rem-cube-cbrt |
| 8× | distribute-rgt-out |
| 7× | associate-+l+ pow-div |
| 6× | sin-sum pow-sub |
| 5× | distribute-rgt1-in |
| 4× | pow-flip unpow3 cube-mult sum-cubes inv-pow unsub-neg |
| 3× | associate-+l- associate--r- |
| 2× | sqrt-unprod rem-square-sqrt rem-cbrt-cube |
| 1× | associate--l- distribute-rgt-out-- rem-exp-log cos-sum count-2 |
| 98× | intervals |
| 52.3s | 12154× | body | 10240 | exit |
| 46.0s | 437147× | body | 80 | valid |
| 25.9s | 27121× | body | 1280 | valid |
| 11.8s | 21329× | body | 640 | valid |
| 6.8s | 67962× | body | 80 | nan |
| 5.3s | 10658× | body | 320 | valid |
| 3.0s | 1648× | body | 1280 | nan |
| 2.0s | 9643× | body | 160 | valid |
| 1.8s | 1340× | body | 640 | nan |
| 1.5s | 49536× | pre | 80 | true |
| 865.0ms | 1496× | body | 2560 | valid |
| 833.0ms | 625× | body | 5120 | valid |
| 524.0ms | 669× | body | 320 | nan |
| 175.0ms | 339× | body | 160 | nan |
474 calls:
| 21.2s | (- (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))) |
| 4.8s | (- (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)) |
| 4.0s | (+ (* (* (* (/ 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))))) (* (* (/ 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)))))) (* 1.0 (* (* (/ 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.3s | (* (/ (* 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)))) |
| 3.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)))) |
Total 99.9b remaining (39.4%)
Threshold costs 3.0b (1.2%)
| 10.4b | 58.8% | math.log/2 on complex, real part |
| 10.3b | 58.7% | math.log10 on complex, real part |
| 10.1b | 59.3% | math.log/1 on complex, real part |
| 10.0b | 58.8% | math.abs on complex |
| 9.4b | 38.4% | Octave 3.8, jcobi/4 |