| 149× | rewrite-expression-head |
430 calls:
| 15.9s | (* (/ (+ (+ (+ (+ (+ 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) |
| 9.2s | (* (/ 1.0 21.0) (* (* (* (* (* (* (fabs x) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x)) (fabs x))) |
| 8.9s | (/ (+ (+ (+ (+ (+ 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))))) |
| 4.6s | (* (/ (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))))) 1) (/ (* x (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)))))) (fma (* 2.0 0.0001789971) (* (pow x 3) (* (pow x 3) (* (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)))))))) |
| 1.7s | (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))))) |
| 17444× | times-frac |
| 15001× | *-un-lft-identity |
| 10507× | add-sqr-sqrt |
| 8097× | add-cube-cbrt |
| 3023× | add-exp-log |
| 2595× | distribute-lft-out |
| 1960× | add-cbrt-cube |
| 1537× | fma-def |
| 1526× | pow1 |
| 1346× | prod-diff |
| 1065× | prod-exp |
| 1027× | sqrt-prod |
| 957× | div-exp |
| 899× | difference-of-squares |
| 898× | distribute-lft-out-- |
| 800× | associate-*l* |
| 779× | associate-/l* |
| 708× | associate-*r* |
| 649× | associate-/r* |
| 636× | associate-/r/ |
| 620× | add-log-exp |
| 598× | cbrt-unprod |
| 485× | div-inv |
| 451× | cbrt-undiv |
| 434× | log1p-expm1-u |
| 430× | expm1-log1p-u |
| 425× | pow-prod-down |
| 397× | cbrt-prod |
| 381× | flip-+ flip3-+ |
| 345× | unpow-prod-down |
| 273× | 1-exp |
| 250× | associate-*r/ |
| 242× | log-pow |
| 229× | unswap-sqr |
| 225× | sqrt-pow1 |
| 218× | distribute-rgt-in distribute-lft-in |
| 212× | sqrt-div |
| 193× | fma-neg |
| 175× | flip3-- flip-- |
| 156× | associate-*l/ |
| 155× | associate--l+ |
| 136× | cube-prod |
| 129× | rec-exp |
| 125× | exp-sum |
| 124× | sqr-pow |
| 123× | pow1/2 |
| 116× | *-commutative |
| 115× | cbrt-div |
| 112× | log-prod |
| 109× | frac-times |
| 96× | pow-prod-up |
| 88× | swap-sqr |
| 84× | pow-exp |
| 83× | frac-2neg clear-num |
| 81× | pow1/3 sub-neg |
| 75× | associate-/l/ |
| 74× | diff-log pow-unpow |
| 71× | exp-prod associate--r+ |
| 64× | pow-sqr |
| 55× | pow-to-exp |
| 52× | fma-udef |
| 50× | frac-add |
| 38× | rem-sqrt-square |
| 36× | distribute-lft-neg-in distribute-rgt-neg-in |
| 35× | pow-plus |
| 33× | pow-flip |
| 32× | inv-pow sum-log |
| 31× | associate-+l+ |
| 29× | neg-mul-1 |
| 28× | pow2 |
| 19× | pow-pow unpow3 cube-mult |
| 18× | +-commutative |
| 17× | hypot-def |
| 15× | associate--l- |
| 14× | cube-div |
| 13× | rem-cube-cbrt |
| 12× | rem-exp-log hypot-udef div-sub rem-log-exp |
| 10× | un-div-inv sqrt-unprod |
| 6× | neg-log |
| 5× | rem-cbrt-cube log-div unpow2 |
| 4× | sin-sum expm1-log1p rem-square-sqrt frac-sub expm1-udef associate-+l- |
| 3× | log1p-udef log1p-expm1 neg-sub0 log-rec |
| 2× | exp-diff associate-+r+ distribute-rgt-neg-out |
| 1× | sqrt-undiv sqrt-pow2 exp-neg distribute-neg-in unpow-prod-up pow-sub distribute-frac-neg cos-sum associate--r- |
| 99× | intervals |
| 42.9s | 12064× | body | 10240 | exit |
| 21.8s | 495844× | body | 80 | valid |
| 11.2s | 27503× | body | 1280 | valid |
| 5.1s | 21253× | body | 640 | valid |
| 2.5s | 68070× | body | 80 | nan |
| 1.8s | 1658× | body | 1280 | nan |
| 1.6s | 10646× | body | 320 | valid |
| 1.0s | 1391× | body | 640 | nan |
| 983.0ms | 9277× | body | 160 | valid |
| 616.0ms | 49536× | pre | 80 | true |
| 373.0ms | 1552× | body | 2560 | valid |
| 337.0ms | 646× | body | 320 | nan |
| 136.0ms | 587× | body | 5120 | valid |
| 105.0ms | 301× | body | 160 | nan |
430 calls:
| 4.9s | (fma (* (cbrt (fma (sqrt 1.0) (sqrt 1.0) (- (* (cbrt z) (* (cbrt z) (cbrt z)))))) (cbrt (fma (sqrt 1.0) (sqrt 1.0) (- (* (cbrt z) (* (cbrt z) (cbrt z))))))) (cbrt (fma (sqrt 1.0) (sqrt 1.0) (- (* (cbrt z) (* (cbrt z) (cbrt z)))))) (fma -1 z (- z 1.0))) |
| 3.3s | (pow (* (/ (+ alpha beta) (+ (+ alpha beta) (fma 2.0 i 2.0))) (/ (- beta alpha) (+ (fma 2.0 i beta) alpha))) 3) |
| 2.4s | (/ (* (* (/ i (+ (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0))) (/ (+ (+ alpha beta) i) (- (+ (+ alpha beta) (* 2.0 i)) (sqrt 1.0)))) (sqrt (fma beta alpha (* i (+ (+ alpha beta) i))))) (/ (fma i 2.0 (+ alpha beta)) (/ (sqrt (fma beta alpha (* i (+ (+ alpha beta) i)))) (fma i 2.0 (+ alpha beta))))) |
| 2.1s | (/ (* (- (* (+ (+ 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)))) |
| 2.1s | (/ (* (* (- (* (+ (+ 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)))) |
Total 43.3b remaining (34.6%)
Threshold costs 5.6b (4.5%)
| 9.0b | 26.4% | _divideComplex, imaginary part |
| 8.2b | 24.7% | _divideComplex, real part |
| 6.3b | 0% | Octave 3.8, jcobi/2 |
| 4.7b | 19.5% | math.sqrt on complex, real part |
| 4.4b | 26.5% | math.sqrt on complex, imaginary part, im greater than 0 branch |