Average Error: 31.0 → 17.3
Time: 1.2m
Precision: 64
Internal Precision: 128
\[\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right) \cdot \log base + \tan^{-1}_* \frac{im}{re} \cdot 0}{\log base \cdot \log base + 0 \cdot 0}\]
\[\begin{array}{l} \mathbf{if}\;im \le -4.135402362840663 \cdot 10^{+167}:\\ \;\;\;\;\frac{\log \left(-re\right)}{\log base}\\ \mathbf{elif}\;im \le -1.8502929772009478 \cdot 10^{-155}:\\ \;\;\;\;\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log base}\\ \mathbf{elif}\;im \le 1.66338613518387 \cdot 10^{-124}:\\ \;\;\;\;\frac{\log \left(-re\right)}{\log base}\\ \mathbf{elif}\;im \le 6.266361841679746 \cdot 10^{+28}:\\ \;\;\;\;\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log base}\\ \mathbf{else}:\\ \;\;\;\;\frac{\log im}{\log base}\\ \end{array}\]

Error

Bits error versus re

Bits error versus im

Bits error versus base

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 3 regimes
  2. if im < -4.135402362840663e+167 or -1.8502929772009478e-155 < im < 1.66338613518387e-124

    1. Initial program 38.6

      \[\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right) \cdot \log base + \tan^{-1}_* \frac{im}{re} \cdot 0}{\log base \cdot \log base + 0 \cdot 0}\]
    2. Initial simplification38.6

      \[\leadsto \frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right) \cdot \log base}{\log base \cdot \log base}\]
    3. Taylor expanded around -inf 20.5

      \[\leadsto \frac{\log \color{blue}{\left(-1 \cdot re\right)} \cdot \log base}{\log base \cdot \log base}\]
    4. Simplified20.5

      \[\leadsto \frac{\log \color{blue}{\left(-re\right)} \cdot \log base}{\log base \cdot \log base}\]
    5. Using strategy rm
    6. Applied times-frac20.4

      \[\leadsto \color{blue}{\frac{\log \left(-re\right)}{\log base} \cdot \frac{\log base}{\log base}}\]
    7. Simplified20.4

      \[\leadsto \frac{\log \left(-re\right)}{\log base} \cdot \color{blue}{1}\]

    if -4.135402362840663e+167 < im < -1.8502929772009478e-155 or 1.66338613518387e-124 < im < 6.266361841679746e+28

    1. Initial program 17.3

      \[\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right) \cdot \log base + \tan^{-1}_* \frac{im}{re} \cdot 0}{\log base \cdot \log base + 0 \cdot 0}\]
    2. Initial simplification17.3

      \[\leadsto \frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right) \cdot \log base}{\log base \cdot \log base}\]
    3. Using strategy rm
    4. Applied times-frac17.2

      \[\leadsto \color{blue}{\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log base} \cdot \frac{\log base}{\log base}}\]
    5. Simplified17.2

      \[\leadsto \frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log base} \cdot \color{blue}{1}\]

    if 6.266361841679746e+28 < im

    1. Initial program 40.9

      \[\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right) \cdot \log base + \tan^{-1}_* \frac{im}{re} \cdot 0}{\log base \cdot \log base + 0 \cdot 0}\]
    2. Initial simplification40.9

      \[\leadsto \frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right) \cdot \log base}{\log base \cdot \log base}\]
    3. Taylor expanded around 0 12.1

      \[\leadsto \color{blue}{\frac{\log im}{\log base}}\]
  3. Recombined 3 regimes into one program.
  4. Final simplification17.3

    \[\leadsto \begin{array}{l} \mathbf{if}\;im \le -4.135402362840663 \cdot 10^{+167}:\\ \;\;\;\;\frac{\log \left(-re\right)}{\log base}\\ \mathbf{elif}\;im \le -1.8502929772009478 \cdot 10^{-155}:\\ \;\;\;\;\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log base}\\ \mathbf{elif}\;im \le 1.66338613518387 \cdot 10^{-124}:\\ \;\;\;\;\frac{\log \left(-re\right)}{\log base}\\ \mathbf{elif}\;im \le 6.266361841679746 \cdot 10^{+28}:\\ \;\;\;\;\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log base}\\ \mathbf{else}:\\ \;\;\;\;\frac{\log im}{\log base}\\ \end{array}\]

Reproduce

herbie shell --seed 2018360 
(FPCore (re im base)
  :name "math.log/2 on complex, real part"
  (/ (+ (* (log (sqrt (+ (* re re) (* im im)))) (log base)) (* (atan2 im re) 0)) (+ (* (log base) (log base)) (* 0 0))))

Details

Time bar (total: 18.3s)Debug log

start1.7s

Algorithm
intervals

setup25.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 30.5b

localize43.0ms

Local error

Found 4 expressions with local error:

29.1b
(sqrt (+ (* re re) (* im im)))
0.5b
(* (log base) (log base))
0.4b
(* (log (sqrt (+ (* re re) (* im im)))) (log base))
0.4b
(/ (* (log (sqrt (+ (* re re) (* im im)))) (log base)) (* (log base) (log base)))

rewrite20.0ms

Algorithm
rewrite-expression-head
Counts
4 → 74
Calls

4 calls. Slowest were:

9.0ms
(/ (* (log (sqrt (+ (* re re) (* im im)))) (log base)) (* (log base) (log base)))
4.0ms
(* (log (sqrt (+ (* re re) (* im im)))) (log base))
3.0ms
(* (log base) (log base))

series451.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

210.0ms
(* (log base) (log base))
124.0ms
(/ (* (log (sqrt (+ (* re re) (* im im)))) (log base)) (* (log base) (log base)))
96.0ms
(* (log (sqrt (+ (* re re) (* im im)))) (log base))
21.0ms
(sqrt (+ (* re re) (* im im)))

simplify2.6s

Counts
49 → 86
Calls

49 calls. Slowest were:

317.0ms
(- (+ (log (log (sqrt (+ (* re re) (* im im))))) (log (log base))) (log (* (log base) (log base))))
296.0ms
(pow (- (log -1) (log (/ -1 base))) 2)
257.0ms
(- (+ (log (log (sqrt (+ (* re re) (* im im))))) (log (log base))) (+ (log (log base)) (log (log base))))

prune1.1s

Pruning

13 alts after pruning (13 fresh and 0 done)

Merged error: 7.4b

localize11.0ms

Local error

Found 2 expressions with local error:

0.4b
(* (/ -1 (log base)) (log (/ -1 re)))
0.3b
(/ -1 (log base))

rewrite5.0ms

Algorithm
rewrite-expression-head
Counts
2 → 41
Calls

2 calls. Slowest were:

4.0ms
(* (/ -1 (log base)) (log (/ -1 re)))
0.0ms
(/ -1 (log base))

series308.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

195.0ms
(/ -1 (log base))
113.0ms
(* (/ -1 (log base)) (log (/ -1 re)))

simplify1.2s

Counts
27 → 47
Calls

27 calls. Slowest were:

365.0ms
(* (* (* (/ -1 (log base)) (/ -1 (log base))) (/ -1 (log base))) (* (* (log (/ -1 re)) (log (/ -1 re))) (log (/ -1 re))))
128.0ms
(* -1 (/ (- (log -1) (log re)) (log base)))
111.0ms
(* -1 (/ (log (/ -1 re)) (- (log -1) (log (/ -1 base)))))

prune642.0ms

Pruning

21 alts after pruning (20 fresh and 1 done)

Merged error: 7.3b

localize13.0ms

Local error

Found 3 expressions with local error:

0.5b
(* (log base) (log base))
0.4b
(* (log (- re)) (log base))
0.4b
(/ (* (log (- re)) (log base)) (* (log base) (log base)))

rewrite19.0ms

Algorithm
rewrite-expression-head
Counts
3 → 59
Calls

3 calls. Slowest were:

8.0ms
(/ (* (log (- re)) (log base)) (* (log base) (log base)))
4.0ms
(* (log (- re)) (log base))
3.0ms
(* (log base) (log base))

series453.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

204.0ms
(* (log base) (log base))
146.0ms
(/ (* (log (- re)) (log base)) (* (log base) (log base)))
103.0ms
(* (log (- re)) (log base))

simplify2.6s

Counts
40 → 68
Calls

40 calls. Slowest were:

316.0ms
(- (+ (log (log (- re))) (log (log base))) (log (* (log base) (log base))))
312.0ms
(pow (- (log -1) (log (/ -1 base))) 2)
252.0ms
(* -1 (* (log (/ 1 base)) (- (log -1) (log (/ 1 re)))))

prune1.2s

Pruning

21 alts after pruning (20 fresh and 1 done)

Merged error: 7.3b

localize11.0ms

Local error

Found 1 expressions with local error:

0.3b
(/ (log (- re)) (log base))

rewrite12.0ms

Algorithm
rewrite-expression-head
Counts
1 → 31
Calls

1 calls. Slowest were:

11.0ms
(/ (log (- re)) (log base))

series90.0ms

Counts
1 → 3
Calls

1 calls. Slowest were:

89.0ms
(/ (log (- re)) (log base))

simplify633.0ms

Counts
34 → 34
Calls

34 calls. Slowest were:

182.0ms
(* -1 (/ (- (log -1) (log (/ 1 re))) (log (/ 1 base))))
162.0ms
(/ (* (* (log (- re)) (log (- re))) (log (- re))) (* (* (log base) (log base)) (log base)))
139.0ms
(* -1 (/ (log (/ -1 re)) (- (log -1) (log (/ -1 base)))))

prune364.0ms

Pruning

21 alts after pruning (19 fresh and 2 done)

Merged error: 7.3b

regimes661.0ms

Accuracy

54.3% (10.1b remaining)

Error of 17.3b against oracle of 7.2b and baseline of 29.3b

bsearch4.1s