Average Error: 30.9 → 18.3
Time: 55.5s
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.83216196722596 \cdot 10^{+96}:\\ \;\;\;\;\frac{\log \left(-re\right)}{\log base}\\ \mathbf{elif}\;im \le -4.822064075697262 \cdot 10^{-20}:\\ \;\;\;\;\frac{\log \left({\left(\sqrt{re \cdot re + im \cdot im}\right)}^{\frac{1}{3}} \cdot \left(\sqrt[3]{\sqrt{\sqrt{re \cdot re + im \cdot im}} \cdot \sqrt{\sqrt{re \cdot re + im \cdot im}}} \cdot \sqrt[3]{\sqrt{re \cdot re + im \cdot im}}\right)\right)}{\log base}\\ \mathbf{elif}\;im \le 2.0164282822791598 \cdot 10^{-135}:\\ \;\;\;\;\frac{\log \left(-re\right)}{\log base}\\ \mathbf{elif}\;im \le 3.401373965950837 \cdot 10^{+105}:\\ \;\;\;\;\frac{\log \left(\left(\sqrt[3]{\sqrt{re \cdot re + im \cdot im}} \cdot \sqrt[3]{\sqrt{re \cdot re + im \cdot im}}\right) \cdot \sqrt[3]{\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 4 regimes
  2. if im < -4.83216196722596e+96 or -4.822064075697262e-20 < im < 2.0164282822791598e-135

    1. Initial program 33.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. Simplified33.3

      \[\leadsto \color{blue}{\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-frac33.2

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

      \[\leadsto \frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log base} \cdot \color{blue}{1}\]
    6. Taylor expanded around -inf 22.7

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

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

    if -4.83216196722596e+96 < im < -4.822064075697262e-20

    1. Initial program 15.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. Simplified15.6

      \[\leadsto \color{blue}{\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-frac15.5

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

      \[\leadsto \frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log base} \cdot \color{blue}{1}\]
    6. Using strategy rm
    7. Applied add-cube-cbrt15.5

      \[\leadsto \frac{\log \color{blue}{\left(\left(\sqrt[3]{\sqrt{re \cdot re + im \cdot im}} \cdot \sqrt[3]{\sqrt{re \cdot re + im \cdot im}}\right) \cdot \sqrt[3]{\sqrt{re \cdot re + im \cdot im}}\right)}}{\log base} \cdot 1\]
    8. Using strategy rm
    9. Applied pow1/315.5

      \[\leadsto \frac{\log \left(\left(\sqrt[3]{\sqrt{re \cdot re + im \cdot im}} \cdot \sqrt[3]{\sqrt{re \cdot re + im \cdot im}}\right) \cdot \color{blue}{{\left(\sqrt{re \cdot re + im \cdot im}\right)}^{\frac{1}{3}}}\right)}{\log base} \cdot 1\]
    10. Using strategy rm
    11. Applied add-sqr-sqrt15.5

      \[\leadsto \frac{\log \left(\left(\sqrt[3]{\sqrt{re \cdot re + im \cdot im}} \cdot \sqrt[3]{\sqrt{\color{blue}{\sqrt{re \cdot re + im \cdot im} \cdot \sqrt{re \cdot re + im \cdot im}}}}\right) \cdot {\left(\sqrt{re \cdot re + im \cdot im}\right)}^{\frac{1}{3}}\right)}{\log base} \cdot 1\]
    12. Applied sqrt-prod15.5

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

    if 2.0164282822791598e-135 < im < 3.401373965950837e+105

    1. Initial program 15.0

      \[\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. Simplified15.0

      \[\leadsto \color{blue}{\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-frac14.9

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

      \[\leadsto \frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log base} \cdot \color{blue}{1}\]
    6. Using strategy rm
    7. Applied add-cube-cbrt14.9

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

    if 3.401373965950837e+105 < im

    1. Initial program 51.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. Simplified51.3

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

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

    \[\leadsto \begin{array}{l} \mathbf{if}\;im \le -4.83216196722596 \cdot 10^{+96}:\\ \;\;\;\;\frac{\log \left(-re\right)}{\log base}\\ \mathbf{elif}\;im \le -4.822064075697262 \cdot 10^{-20}:\\ \;\;\;\;\frac{\log \left({\left(\sqrt{re \cdot re + im \cdot im}\right)}^{\frac{1}{3}} \cdot \left(\sqrt[3]{\sqrt{\sqrt{re \cdot re + im \cdot im}} \cdot \sqrt{\sqrt{re \cdot re + im \cdot im}}} \cdot \sqrt[3]{\sqrt{re \cdot re + im \cdot im}}\right)\right)}{\log base}\\ \mathbf{elif}\;im \le 2.0164282822791598 \cdot 10^{-135}:\\ \;\;\;\;\frac{\log \left(-re\right)}{\log base}\\ \mathbf{elif}\;im \le 3.401373965950837 \cdot 10^{+105}:\\ \;\;\;\;\frac{\log \left(\left(\sqrt[3]{\sqrt{re \cdot re + im \cdot im}} \cdot \sqrt[3]{\sqrt{re \cdot re + im \cdot im}}\right) \cdot \sqrt[3]{\sqrt{re \cdot re + im \cdot im}}\right)}{\log base}\\ \mathbf{else}:\\ \;\;\;\;\frac{\log im}{\log base}\\ \end{array}\]

Reproduce

herbie shell --seed 2019005 
(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: 53.8s)Debug log

sample1.1s

Algorithm
intervals

simplify12.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

12.0ms
(/ (+ (* (log (sqrt (+ (* re re) (* im im)))) (log base)) (* (atan2 im re) 0)) (+ (* (log base) (log base)) (* 0 0)))

prune11.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 29.8b

localize38.0ms

Local error

Found 4 expressions with local error:

28.6b
(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)))

rewrite21.0ms

Algorithm
rewrite-expression-head
Counts
4 → 74
Calls

4 calls. Slowest were:

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

series321.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

161.0ms
(* (log base) (log base))
72.0ms
(/ (* (log (sqrt (+ (* re re) (* im im)))) (log base)) (* (log base) (log base)))
68.0ms
(* (log (sqrt (+ (* re re) (* im im)))) (log base))
19.0ms
(sqrt (+ (* re re) (* im im)))

simplify3.1s

Counts
49 → 86
Calls

49 calls. Slowest were:

387.0ms
(- (+ (log (log (sqrt (+ (* re re) (* im im))))) (log (log base))) (log (* (log base) (log base))))
361.0ms
(* -1 (* (log (/ -1 re)) (- (log -1) (log (/ -1 base)))))
310.0ms
(- (+ (log (log (sqrt (+ (* re re) (* im im))))) (log (log base))) (+ (log (log base)) (log (log base))))

prune1.0s

Pruning

12 alts after pruning (12 fresh and 0 done)

Merged error: 8.2b

localize6.0ms

Local error

Found 2 expressions with local error:

28.6b
(sqrt (+ (* re re) (* im im)))
0.3b
(/ (log (sqrt (+ (* re re) (* im im)))) (log base))

rewrite11.0ms

Algorithm
rewrite-expression-head
Counts
2 → 51
Calls

2 calls. Slowest were:

8.0ms
(/ (log (sqrt (+ (* re re) (* im im)))) (log base))
2.0ms
(sqrt (+ (* re re) (* im im)))

series125.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

102.0ms
(/ (log (sqrt (+ (* re re) (* im im)))) (log base))
23.0ms
(sqrt (+ (* re re) (* im im)))

simplify957.0ms

Counts
51 → 57
Calls

51 calls. Slowest were:

165.0ms
(/ (log (sqrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))) (log base))
158.0ms
(sqrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))
117.0ms
(* -1 (/ (log (/ -1 re)) (- (log -1) (log (/ -1 base)))))

prune660.0ms

Pruning

13 alts after pruning (13 fresh and 0 done)

Merged error: 8.2b

localize38.0ms

Local error

Found 4 expressions with local error:

28.6b
(sqrt (+ (* re re) (* im im)))
28.6b
(sqrt (+ (* re re) (* im im)))
28.6b
(sqrt (+ (* re re) (* im im)))
0.5b
(cbrt (sqrt (+ (* re re) (* im im))))

rewrite15.0ms

Algorithm
rewrite-expression-head
Counts
4 → 58
Calls

4 calls. Slowest were:

4.0ms
(sqrt (+ (* re re) (* im im)))
4.0ms
(sqrt (+ (* re re) (* im im)))
4.0ms
(sqrt (+ (* re re) (* im im)))

series216.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

117.0ms
(cbrt (sqrt (+ (* re re) (* im im))))
37.0ms
(sqrt (+ (* re re) (* im im)))
31.0ms
(sqrt (+ (* re re) (* im im)))
31.0ms
(sqrt (+ (* re re) (* im im)))

simplify866.0ms

Counts
38 → 70
Calls

38 calls. Slowest were:

153.0ms
(sqrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))
142.0ms
(sqrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))
129.0ms
(cbrt (sqrt (- (* (* re re) (* re re)) (* (* im im) (* im im)))))

prune1.4s

Pruning

14 alts after pruning (13 fresh and 1 done)

Merged error: 8.2b

localize20.0ms

Local error

Found 4 expressions with local error:

28.6b
(sqrt (+ (* re re) (* im im)))
28.6b
(sqrt (+ (* re re) (* im im)))
28.6b
(sqrt (+ (* re re) (* im im)))
5.0b
(pow (sqrt (+ (* re re) (* im im))) 1/3)

rewrite8.0ms

Algorithm
rewrite-expression-head
Counts
4 → 60
Calls

4 calls. Slowest were:

3.0ms
(pow (sqrt (+ (* re re) (* im im))) 1/3)
2.0ms
(sqrt (+ (* re re) (* im im)))
2.0ms
(sqrt (+ (* re re) (* im im)))

series164.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

88.0ms
(pow (sqrt (+ (* re re) (* im im))) 1/3)
28.0ms
(sqrt (+ (* re re) (* im im)))
28.0ms
(sqrt (+ (* re re) (* im im)))
19.0ms
(sqrt (+ (* re re) (* im im)))

simplify730.0ms

Counts
39 → 72
Calls

39 calls. Slowest were:

160.0ms
(sqrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))
138.0ms
(sqrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))
124.0ms
(sqrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))

prune1.1s

Pruning

13 alts after pruning (12 fresh and 1 done)

Merged error: 8.2b

regimes595.0ms

Accuracy

53.5% (10.8b remaining)

Error of 18.3b against oracle of 7.5b and baseline of 30.8b

bsearch3.0s

end0.0ms

sample38.2s

Algorithm
intervals