Average Error: 31.3 → 4.5
Time: 2.7m
Precision: 64
Internal Precision: 128
\[e^{\log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.re - \tan^{-1}_* \frac{x.im}{x.re} \cdot y.im} \cdot \cos \left(\log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.im + \tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right)\]
\[\begin{array}{l} \mathbf{if}\;x.re \le -2.8299828101670785 \cdot 10^{+24}:\\ \;\;\;\;e^{\log \left(\cos \left(y.im \cdot \log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) - \sin \left(y.im \cdot \log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) \cdot \sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right)\right) - \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot y.im - y.re \cdot \log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right)}\\ \mathbf{else}:\\ \;\;\;\;\frac{\cos \left((y.im \cdot \left(\log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) + \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right))_*\right)}{\frac{{\left(e^{y.im}\right)}^{\left(\tan^{-1}_* \frac{x.im}{x.re}\right)}}{{\left(\sqrt{x.re^2 + x.im^2}^*\right)}^{y.re}}}\\ \end{array}\]

Error

Bits error versus x.re

Bits error versus x.im

Bits error versus y.re

Bits error versus y.im

Derivation

  1. Split input into 2 regimes
  2. if x.re < -2.8299828101670785e+24

    1. Initial program 41.7

      \[e^{\log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.re - \tan^{-1}_* \frac{x.im}{x.re} \cdot y.im} \cdot \cos \left(\log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.im + \tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right)\]
    2. Simplified5.1

      \[\leadsto \color{blue}{\frac{\cos \left((y.im \cdot \left(\log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) + \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right))_*\right)}{\frac{{\left(e^{y.im}\right)}^{\left(\tan^{-1}_* \frac{x.im}{x.re}\right)}}{{\left(\sqrt{x.re^2 + x.im^2}^*\right)}^{y.re}}}}\]
    3. Using strategy rm
    4. Applied add-exp-log5.1

      \[\leadsto \frac{\cos \left((y.im \cdot \left(\log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) + \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right))_*\right)}{\frac{{\left(e^{y.im}\right)}^{\left(\tan^{-1}_* \frac{x.im}{x.re}\right)}}{{\color{blue}{\left(e^{\log \left(\sqrt{x.re^2 + x.im^2}^*\right)}\right)}}^{y.re}}}\]
    5. Applied pow-exp5.1

      \[\leadsto \frac{\cos \left((y.im \cdot \left(\log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) + \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right))_*\right)}{\frac{{\left(e^{y.im}\right)}^{\left(\tan^{-1}_* \frac{x.im}{x.re}\right)}}{\color{blue}{e^{\log \left(\sqrt{x.re^2 + x.im^2}^*\right) \cdot y.re}}}}\]
    6. Applied pow-exp5.1

      \[\leadsto \frac{\cos \left((y.im \cdot \left(\log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) + \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right))_*\right)}{\frac{\color{blue}{e^{y.im \cdot \tan^{-1}_* \frac{x.im}{x.re}}}}{e^{\log \left(\sqrt{x.re^2 + x.im^2}^*\right) \cdot y.re}}}\]
    7. Applied div-exp0.1

      \[\leadsto \frac{\cos \left((y.im \cdot \left(\log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) + \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right))_*\right)}{\color{blue}{e^{y.im \cdot \tan^{-1}_* \frac{x.im}{x.re} - \log \left(\sqrt{x.re^2 + x.im^2}^*\right) \cdot y.re}}}\]
    8. Applied add-exp-log14.8

      \[\leadsto \frac{\color{blue}{e^{\log \left(\cos \left((y.im \cdot \left(\log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) + \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right))_*\right)\right)}}}{e^{y.im \cdot \tan^{-1}_* \frac{x.im}{x.re} - \log \left(\sqrt{x.re^2 + x.im^2}^*\right) \cdot y.re}}\]
    9. Applied div-exp14.8

      \[\leadsto \color{blue}{e^{\log \left(\cos \left((y.im \cdot \left(\log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) + \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right))_*\right)\right) - \left(y.im \cdot \tan^{-1}_* \frac{x.im}{x.re} - \log \left(\sqrt{x.re^2 + x.im^2}^*\right) \cdot y.re\right)}}\]
    10. Using strategy rm
    11. Applied fma-udef14.8

      \[\leadsto e^{\log \left(\cos \color{blue}{\left(y.im \cdot \log \left(\sqrt{x.re^2 + x.im^2}^*\right) + \tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right)}\right) - \left(y.im \cdot \tan^{-1}_* \frac{x.im}{x.re} - \log \left(\sqrt{x.re^2 + x.im^2}^*\right) \cdot y.re\right)}\]
    12. Applied cos-sum14.6

      \[\leadsto e^{\log \color{blue}{\left(\cos \left(y.im \cdot \log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) \cdot \cos \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right) - \sin \left(y.im \cdot \log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) \cdot \sin \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right)\right)} - \left(y.im \cdot \tan^{-1}_* \frac{x.im}{x.re} - \log \left(\sqrt{x.re^2 + x.im^2}^*\right) \cdot y.re\right)}\]
    13. Taylor expanded around 0 4.9

      \[\leadsto e^{\log \left(\cos \left(y.im \cdot \log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) \cdot \color{blue}{1} - \sin \left(y.im \cdot \log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) \cdot \sin \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right)\right) - \left(y.im \cdot \tan^{-1}_* \frac{x.im}{x.re} - \log \left(\sqrt{x.re^2 + x.im^2}^*\right) \cdot y.re\right)}\]

    if -2.8299828101670785e+24 < x.re

    1. Initial program 21.9

      \[e^{\log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.re - \tan^{-1}_* \frac{x.im}{x.re} \cdot y.im} \cdot \cos \left(\log \left(\sqrt{x.re \cdot x.re + x.im \cdot x.im}\right) \cdot y.im + \tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right)\]
    2. Simplified4.1

      \[\leadsto \color{blue}{\frac{\cos \left((y.im \cdot \left(\log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) + \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot y.re\right))_*\right)}{\frac{{\left(e^{y.im}\right)}^{\left(\tan^{-1}_* \frac{x.im}{x.re}\right)}}{{\left(\sqrt{x.re^2 + x.im^2}^*\right)}^{y.re}}}}\]
  3. Recombined 2 regimes into one program.
  4. Final simplification4.5

    \[\leadsto \begin{array}{l} \mathbf{if}\;x.re \le -2.8299828101670785 \cdot 10^{+24}:\\ \;\;\;\;e^{\log \left(\cos \left(y.im \cdot \log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) - \sin \left(y.im \cdot \log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) \cdot \sin \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right)\right) - \left(\tan^{-1}_* \frac{x.im}{x.re} \cdot y.im - y.re \cdot \log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right)}\\ \mathbf{else}:\\ \;\;\;\;\frac{\cos \left((y.im \cdot \left(\log \left(\sqrt{x.re^2 + x.im^2}^*\right)\right) + \left(y.re \cdot \tan^{-1}_* \frac{x.im}{x.re}\right))_*\right)}{\frac{{\left(e^{y.im}\right)}^{\left(\tan^{-1}_* \frac{x.im}{x.re}\right)}}{{\left(\sqrt{x.re^2 + x.im^2}^*\right)}^{y.re}}}\\ \end{array}\]

Reproduce

herbie shell --seed 2019010 +o rules:numerics
(FPCore (x.re x.im y.re y.im)
  :name "powComplex, real part"
  (* (exp (- (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) (* (atan2 x.im x.re) y.im))) (cos (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re)))))

Details

Time bar (total: 2.7m)Debug log

sample4.1s

Algorithm
intervals

simplify119.0ms

Counts
1 → 1
Calls
1 calls:
Slowest
118.0ms
(* (exp (- (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.re) (* (atan2 x.im x.re) y.im))) (cos (+ (* (log (sqrt (+ (* x.re x.re) (* x.im x.im)))) y.im) (* (atan2 x.im x.re) y.re))))

prune24.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 1.5b

localize36.0ms

Local error

Found 4 expressions with local error:

22.5b
(cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))
0.2b
(* (atan2 x.im x.re) y.re)
0.1b
(fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))
0.0b
(/ (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) (/ (pow (exp y.im) (atan2 x.im x.re)) (pow (hypot x.re x.im) y.re)))

rewrite51.0ms

Algorithm
rewrite-expression-head
Rules
264×times-frac
146×add-cube-cbrt
146×*-un-lft-identity
146×add-sqr-sqrt
144×unpow-prod-down
40×associate-/r*
24×add-exp-log
19×div-exp
add-cbrt-cube
pow-exp
pow-to-exp
div-inv
add-log-exp
log1p-expm1-u
pow1
expm1-log1p-u
associate-/l*
fma-udef
associate-*l*
cbrt-undiv
cos-sum
*-commutative
associate-/r/
div-sub
frac-2neg
clear-num
Counts
4 → 222
Calls
4 calls:
Slowest
20.0ms
(/ (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) (/ (pow (exp y.im) (atan2 x.im x.re)) (pow (hypot x.re x.im) y.re)))
15.0ms
(cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))
5.0ms
(* (atan2 x.im x.re) y.re)
0.0ms
(fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))

series574.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
236.0ms
(/ (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) (/ (pow (exp y.im) (atan2 x.im x.re)) (pow (hypot x.re x.im) y.re)))
162.0ms
(fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))
160.0ms
(cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))
15.0ms
(* (atan2 x.im x.re) y.re)

simplify11.7s

Counts
311 → 234
Calls
311 calls:
Slowest
360.0ms
(/ (* (* (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))) (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))) (/ (* (* (pow (exp y.im) (atan2 x.im x.re)) (pow (exp y.im) (atan2 x.im x.re))) (pow (exp y.im) (atan2 x.im x.re))) (* (* (pow (hypot x.re x.im) y.re) (pow (hypot x.re x.im) y.re)) (pow (hypot x.re x.im) y.re))))
355.0ms
(/ (* (cbrt (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))) (cbrt (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))))) (/ (* (cbrt (pow (exp y.im) (atan2 x.im x.re))) (cbrt (pow (exp y.im) (atan2 x.im x.re)))) (* (cbrt (pow (hypot x.re x.im) y.re)) (cbrt (pow (hypot x.re x.im) y.re)))))
293.0ms
(/ (* (cbrt (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))) (cbrt (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))))) (/ (* (cbrt (pow (exp y.im) (atan2 x.im x.re))) (cbrt (pow (exp y.im) (atan2 x.im x.re)))) (pow 1 y.re)))
220.0ms
(/ (* (cos (- (* (atan2 x.im x.re) y.re) (* y.im (log (/ 1 x.re))))) (exp (* -1 (* y.re (log (/ 1 x.re)))))) (exp (* y.im (atan2 x.im x.re))))
213.0ms
(/ (* (* (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))) (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))) (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))) (* (* (/ (pow (exp y.im) (atan2 x.im x.re)) (pow (hypot x.re x.im) y.re)) (/ (pow (exp y.im) (atan2 x.im x.re)) (pow (hypot x.re x.im) y.re))) (/ (pow (exp y.im) (atan2 x.im x.re)) (pow (hypot x.re x.im) y.re))))

prune4.6s

Pruning

5 alts after pruning (4 fresh and 1 done)

Merged error: 1.0b

localize27.0ms

Local error

Found 4 expressions with local error:

22.6b
(cos (fma (- y.im) (log (/ -1 x.re)) (* y.re (atan2 x.im x.re))))
0.2b
(* y.re (atan2 x.im x.re))
0.1b
(fma (- y.im) (log (/ -1 x.re)) (* y.re (atan2 x.im x.re)))
0.0b
(/ (cos (fma (- y.im) (log (/ -1 x.re)) (* y.re (atan2 x.im x.re)))) (* (pow (/ -1 x.re) y.re) (pow (exp y.im) (atan2 x.im x.re))))

rewrite28.0ms

Algorithm
rewrite-expression-head
Rules
24×add-exp-log
10×div-exp
prod-exp
add-cbrt-cube
add-cube-cbrt
*-un-lft-identity
add-sqr-sqrt
pow-exp
pow-to-exp
add-log-exp
log1p-expm1-u
pow1
expm1-log1p-u
associate-/l*
fma-udef
associate-*r*
times-frac
cos-sum
cbrt-undiv
div-inv
cbrt-unprod
*-commutative
div-sub
frac-2neg
associate-/r*
clear-num
Counts
4 → 65
Calls
4 calls:
Slowest
20.0ms
(/ (cos (fma (- y.im) (log (/ -1 x.re)) (* y.re (atan2 x.im x.re)))) (* (pow (/ -1 x.re) y.re) (pow (exp y.im) (atan2 x.im x.re))))
4.0ms
(cos (fma (- y.im) (log (/ -1 x.re)) (* y.re (atan2 x.im x.re))))
2.0ms
(* y.re (atan2 x.im x.re))
0.0ms
(fma (- y.im) (log (/ -1 x.re)) (* y.re (atan2 x.im x.re)))

series496.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
237.0ms
(/ (cos (fma (- y.im) (log (/ -1 x.re)) (* y.re (atan2 x.im x.re)))) (* (pow (/ -1 x.re) y.re) (pow (exp y.im) (atan2 x.im x.re))))
139.0ms
(cos (fma (- y.im) (log (/ -1 x.re)) (* y.re (atan2 x.im x.re))))
105.0ms
(fma (- y.im) (log (/ -1 x.re)) (* y.re (atan2 x.im x.re)))
14.0ms
(* y.re (atan2 x.im x.re))

simplify6.6s

Counts
38 → 77
Calls
38 calls:
Slowest
670.0ms
(/ (* (* (cos (fma (- y.im) (log (/ -1 x.re)) (* y.re (atan2 x.im x.re)))) (cos (fma (- y.im) (log (/ -1 x.re)) (* y.re (atan2 x.im x.re))))) (cos (fma (- y.im) (log (/ -1 x.re)) (* y.re (atan2 x.im x.re))))) (* (* (* (pow (/ -1 x.re) y.re) (pow (exp y.im) (atan2 x.im x.re))) (* (pow (/ -1 x.re) y.re) (pow (exp y.im) (atan2 x.im x.re)))) (* (pow (/ -1 x.re) y.re) (pow (exp y.im) (atan2 x.im x.re)))))
374.0ms
(- (log (cos (fma (- y.im) (log (/ -1 x.re)) (* y.re (atan2 x.im x.re))))) (+ (log (pow (/ -1 x.re) y.re)) (* (log (exp y.im)) (atan2 x.im x.re))))
373.0ms
(cos (- (* (atan2 x.im x.re) y.re) (+ (* y.im (log -1)) (* y.im (log (/ 1 x.re))))))
361.0ms
(- (log (cos (fma (- y.im) (log (/ -1 x.re)) (* y.re (atan2 x.im x.re))))) (+ (log (pow (/ -1 x.re) y.re)) (* y.im (atan2 x.im x.re))))
355.0ms
(- (log (cos (fma (- y.im) (log (/ -1 x.re)) (* y.re (atan2 x.im x.re))))) (+ (* (log (/ -1 x.re)) y.re) (* (log (exp y.im)) (atan2 x.im x.re))))

prune1.6s

Pruning

5 alts after pruning (3 fresh and 2 done)

Merged error: 1.0b

localize15.0ms

Local error

Found 4 expressions with local error:

22.5b
(cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))
3.8b
(log (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))))
0.3b
(* (log (hypot x.re x.im)) y.re)
0.3b
(* y.im (atan2 x.im x.re))

rewrite8.0ms

Algorithm
rewrite-expression-head
Rules
add-cube-cbrt
*-un-lft-identity
add-sqr-sqrt
pow1
add-exp-log
add-log-exp
log1p-expm1-u
associate-*l*
add-cbrt-cube
expm1-log1p-u
associate-*r*
log-prod
*-commutative
log-pow
fma-udef
rem-log-exp
cos-sum
Counts
4 → 51
Calls
4 calls:
Slowest
3.0ms
(* (log (hypot x.re x.im)) y.re)
2.0ms
(cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))
2.0ms
(log (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))))
1.0ms
(* y.im (atan2 x.im x.re))

series485.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
203.0ms
(log (cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re))))
160.0ms
(cos (fma y.im (log (hypot x.re x.im)) (* (atan2 x.im x.re) y.re)))
113.0ms
(* (log (hypot x.re x.im)) y.re)
8.0ms
(* y.im (atan2 x.im x.re))

simplify940.0ms

Counts
22 → 63
Calls
22 calls:
Slowest
288.0ms
(- (+ (* y.im (* (atan2 x.im x.re) (* y.re (log x.im)))) (+ (* 1/2 (* (pow (atan2 x.im x.re) 2) (pow y.re 2))) (* 1/2 (* (pow y.im 2) (pow (log x.im) 2))))))
170.0ms
(- 1 (+ (* y.im (* (atan2 x.im x.re) (* y.re (log x.im)))) (* 1/2 (* (pow y.im 2) (pow (log x.im) 2)))))
118.0ms
(* -1 (* y.re (log (/ -1 x.re))))
109.0ms
(log (cos (- (* (atan2 x.im x.re) y.re) (* y.im (log (/ -1 x.re))))))
77.0ms
(cos (- (* (atan2 x.im x.re) y.re) (* y.im (log (/ -1 x.re)))))

prune1.1s

Pruning

3 alts after pruning (2 fresh and 1 done)

Merged error: 0.5b

localize33.0ms

Local error

Found 4 expressions with local error:

22.8b
(cos (* (atan2 x.im x.re) y.re))
22.4b
(sin (* (atan2 x.im x.re) y.re))
6.7b
(cos (* y.im (log (hypot x.re x.im))))
6.5b
(sin (* y.im (log (hypot x.re x.im))))

rewrite17.0ms

Algorithm
rewrite-expression-head
Rules
12×log-prod
add-cube-cbrt
*-un-lft-identity
add-sqr-sqrt
sin-sum
distribute-lft-in
cos-sum
distribute-rgt-in
add-log-exp
log1p-expm1-u
add-exp-log
add-cbrt-cube
pow1
expm1-log1p-u
Counts
4 → 48
Calls
4 calls:
Slowest
4.0ms
(cos (* (atan2 x.im x.re) y.re))
4.0ms
(sin (* (atan2 x.im x.re) y.re))
4.0ms
(sin (* y.im (log (hypot x.re x.im))))
4.0ms
(cos (* y.im (log (hypot x.re x.im))))

series300.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
110.0ms
(cos (* y.im (log (hypot x.re x.im))))
105.0ms
(sin (* y.im (log (hypot x.re x.im))))
45.0ms
(sin (* (atan2 x.im x.re) y.re))
40.0ms
(cos (* (atan2 x.im x.re) y.re))

simplify726.0ms

Counts
24 → 60
Calls
24 calls:
Slowest
289.0ms
(- 1 (* 1/2 (* (pow y.im 2) (pow (log x.im) 2))))
105.0ms
(cos (* -1 (* y.im (log (/ -1 x.re)))))
75.0ms
(sin (* -1 (* y.im (log (/ 1 x.re)))))
67.0ms
(sin (* -1 (* y.im (log (/ -1 x.re)))))
64.0ms
(cos (* -1 (* y.im (log (/ 1 x.re)))))

prune1.7s

Pruning

3 alts after pruning (1 fresh and 2 done)

Merged error: 0.2b

regimes262.0ms

Accuracy

7.7% (2.9b remaining)

Error of 4.5b against oracle of 1.6b and baseline of 4.7b

bsearch994.0ms

end0.0ms

sample2.1m

Algorithm
intervals