Average Error: 30.6 → 17.6
Time: 57.6s
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 -1.339088896496117 \cdot 10^{+154}:\\ \;\;\;\;\frac{\left(\sqrt[3]{\log \left(-re\right)} \cdot \sqrt[3]{\log \left(-re\right)}\right) \cdot \left(\log base \cdot \sqrt[3]{\log \left(-re\right)}\right)}{\log base \cdot \log base}\\ \mathbf{elif}\;im \le -5.717705561230724 \cdot 10^{-79}:\\ \;\;\;\;\frac{1}{\frac{\log base \cdot \log base}{\log \left(\sqrt{re \cdot re + im \cdot im}\right) \cdot \log base}}\\ \mathbf{elif}\;im \le 2.9780377269530606 \cdot 10^{-08}:\\ \;\;\;\;\frac{-1}{\log base} \cdot \log \left(\frac{-1}{re}\right)\\ \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 < -1.339088896496117e+154

    1. Initial program 62.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. Simplified62.0

      \[\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 -inf 51.6

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

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

      \[\leadsto \frac{\color{blue}{\left(\left(\sqrt[3]{\log \left(-re\right)} \cdot \sqrt[3]{\log \left(-re\right)}\right) \cdot \sqrt[3]{\log \left(-re\right)}\right)} \cdot \log base}{\log base \cdot \log base}\]
    7. Applied associate-*l*51.7

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

    if -1.339088896496117e+154 < im < -5.717705561230724e-79

    1. Initial program 13.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. Simplified13.9

      \[\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 clear-num13.9

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

    if -5.717705561230724e-79 < im < 2.9780377269530606e-08

    1. Initial program 24.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. Simplified24.0

      \[\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 -inf 62.8

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

      \[\leadsto \color{blue}{\frac{-1}{\log base} \cdot \log \left(\frac{-1}{re}\right)}\]
    5. Taylor expanded around -inf 62.8

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

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

    if 2.9780377269530606e-08 < im

    1. Initial program 38.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. Simplified38.9

      \[\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 14.1

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

    \[\leadsto \begin{array}{l} \mathbf{if}\;im \le -1.339088896496117 \cdot 10^{+154}:\\ \;\;\;\;\frac{\left(\sqrt[3]{\log \left(-re\right)} \cdot \sqrt[3]{\log \left(-re\right)}\right) \cdot \left(\log base \cdot \sqrt[3]{\log \left(-re\right)}\right)}{\log base \cdot \log base}\\ \mathbf{elif}\;im \le -5.717705561230724 \cdot 10^{-79}:\\ \;\;\;\;\frac{1}{\frac{\log base \cdot \log base}{\log \left(\sqrt{re \cdot re + im \cdot im}\right) \cdot \log base}}\\ \mathbf{elif}\;im \le 2.9780377269530606 \cdot 10^{-08}:\\ \;\;\;\;\frac{-1}{\log base} \cdot \log \left(\frac{-1}{re}\right)\\ \mathbf{else}:\\ \;\;\;\;\frac{\log im}{\log base}\\ \end{array}\]

Reproduce

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

sample1.1s

Algorithm
intervals
Results
781.0ms252×body10240exit
146.0ms448×body80nan
60.0ms257×body80valid

simplify12.0ms

Counts
1 → 1
Calls
1 calls:
Slowest
12.0ms
(/ (+ (* (log (sqrt (+ (* re re) (* im im)))) (log base)) (* (atan2 im re) 0)) (+ (* (log base) (log base)) (* 0 0)))

prune22.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 30.8b

localize66.0ms

Local error

Found 4 expressions with local error:

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

rewrite34.0ms

Algorithm
rewrite-expression-head
Rules
20×add-exp-log
20×add-cbrt-cube
12×pow1
10×add-sqr-sqrt
add-cube-cbrt
*-un-lft-identity
associate-*l*
cbrt-unprod
prod-exp
associate-*r*
add-log-exp
div-exp
cbrt-undiv
sqrt-prod
*-commutative
log-pow
pow-prod-down
sqrt-div
pow1/2
associate-/l*
div-inv
flip-+
pow-prod-up
pow-plus
times-frac
flip3-+
frac-2neg
associate-/r*
clear-num
rem-sqrt-square
pow2
Counts
4 → 74
Calls
4 calls:
Slowest
13.0ms
(/ (* (log (sqrt (+ (* re re) (* im im)))) (log base)) (* (log base) (log base)))
9.0ms
(* (log (sqrt (+ (* re re) (* im im)))) (log base))
8.0ms
(* (log base) (log base))
4.0ms
(sqrt (+ (* re re) (* im im)))

series385.0ms

Counts
4 → 12
Calls
4 calls:
Slowest
181.0ms
(* (log base) (log base))
109.0ms
(/ (* (log (sqrt (+ (* re re) (* im im)))) (log base)) (* (log base) (log base)))
76.0ms
(* (log (sqrt (+ (* re re) (* im im)))) (log base))
19.0ms
(sqrt (+ (* re re) (* im im)))

simplify3.3s

Counts
49 → 86
Calls
49 calls:
Slowest
426.0ms
(* -1 (* (log (/ -1 re)) (- (log -1) (log (/ -1 base)))))
378.0ms
(- (+ (log (log (sqrt (+ (* re re) (* im im))))) (log (log base))) (+ (log (log base)) (log (log base))))
341.0ms
(- (+ (log (log (sqrt (+ (* re re) (* im im))))) (log (log base))) (log (* (log base) (log base))))
305.0ms
(- (log (* (log (sqrt (+ (* re re) (* im im)))) (log base))) (log (* (log base) (log base))))
257.0ms
(/ (* (* (* (log (sqrt (+ (* re re) (* im im)))) (log (sqrt (+ (* re re) (* im im))))) (log (sqrt (+ (* re re) (* im im))))) (* (* (log base) (log base)) (log base))) (* (* (* (log base) (log base)) (* (log base) (log base))) (* (log base) (log base))))

prune1.1s

Pruning

14 alts after pruning (14 fresh and 0 done)

Merged error: 6.3b

localize24.0ms

Local error

Found 3 expressions with local error:

0.3b
(* (/ -1 (log base)) (log (/ -1 re)))
0.2b
(/ -1 (log base))
0.0b
(log (/ -1 re))

rewrite11.0ms

Algorithm
rewrite-expression-head
Rules
12×log-prod
add-cube-cbrt
*-un-lft-identity
add-sqr-sqrt
pow1
add-exp-log
div-inv
add-cbrt-cube
associate-*l*
associate-*r*
distribute-lft-in
distribute-rgt-in
add-log-exp
associate-/r*
log-pow
cbrt-unprod
*-commutative
prod-exp
associate-*l/
pow-prod-down
rem-log-exp
frac-2neg
clear-num
log-div
Counts
3 → 55
Calls
3 calls:
Slowest
8.0ms
(* (/ -1 (log base)) (log (/ -1 re)))
1.0ms
(log (/ -1 re))
1.0ms
(/ -1 (log base))

series535.0ms

Counts
3 → 9
Calls
3 calls:
Slowest
250.0ms
(/ -1 (log base))
167.0ms
(log (/ -1 re))
118.0ms
(* (/ -1 (log base)) (log (/ -1 re)))

simplify1.1s

Counts
36 → 64
Calls
36 calls:
Slowest
441.0ms
(* (* (* (/ -1 (log base)) (/ -1 (log base))) (/ -1 (log base))) (* (* (log (/ -1 re)) (log (/ -1 re))) (log (/ -1 re))))
125.0ms
(* -1 (/ (log (/ -1 re)) (- (log -1) (log (/ -1 base)))))
114.0ms
(* -1 (/ (- (log -1) (log re)) (log base)))
83.0ms
(* (log (/ 1 re)) (/ -1 (log base)))
50.0ms
(* (log (/ -1 re)) (/ -1 (log base)))

prune693.0ms

Pruning

21 alts after pruning (21 fresh and 0 done)

Merged error: 6.2b

localize1.0ms

Local error

Found 3 expressions with local error:

0.3b
(* (/ -1 (log base)) (log (/ -1 re)))
0.2b
(/ -1 (log base))
0.0b
(log (/ -1 re))

rewrite6.0ms

Algorithm
rewrite-expression-head
Rules
12×log-prod
add-cube-cbrt
*-un-lft-identity
add-sqr-sqrt
pow1
add-exp-log
div-inv
add-cbrt-cube
associate-*l*
associate-*r*
distribute-lft-in
distribute-rgt-in
add-log-exp
associate-/r*
log-pow
cbrt-unprod
*-commutative
prod-exp
associate-*l/
pow-prod-down
rem-log-exp
frac-2neg
clear-num
log-div
Counts
3 → 55
Calls
3 calls:
Slowest
4.0ms
(* (/ -1 (log base)) (log (/ -1 re)))
0.0ms
(log (/ -1 re))
0.0ms
(/ -1 (log base))

series506.0ms

Counts
3 → 9
Calls
3 calls:
Slowest
224.0ms
(/ -1 (log base))
189.0ms
(log (/ -1 re))
92.0ms
(* (/ -1 (log base)) (log (/ -1 re)))

simplify1.3s

Counts
36 → 64
Calls
36 calls:
Slowest
492.0ms
(* (* (* (/ -1 (log base)) (/ -1 (log base))) (/ -1 (log base))) (* (* (log (/ -1 re)) (log (/ -1 re))) (log (/ -1 re))))
141.0ms
(* -1 (/ (log (/ -1 re)) (- (log -1) (log (/ -1 base)))))
125.0ms
(* -1 (/ (- (log -1) (log re)) (log base)))
56.0ms
(* (log (/ 1 re)) (/ -1 (log base)))
48.0ms
(* (/ -1 (log base)) (log (/ -1 re)))

prune713.0ms

Pruning

21 alts after pruning (20 fresh and 1 done)

Merged error: 6.2b

localize11.0ms

Local error

Found 3 expressions with local error:

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

rewrite17.0ms

Algorithm
rewrite-expression-head
Rules
19×add-exp-log
19×add-cbrt-cube
11×pow1
add-cube-cbrt
associate-*l*
*-un-lft-identity
add-sqr-sqrt
cbrt-unprod
prod-exp
associate-*r*
div-exp
cbrt-undiv
add-log-exp
*-commutative
pow-prod-down
associate-/l*
div-inv
pow-prod-up
log-pow
pow-plus
times-frac
frac-2neg
associate-/r*
clear-num
pow2
Counts
3 → 59
Calls
3 calls:
Slowest
8.0ms
(/ (* (log (- re)) (log base)) (* (log base) (log base)))
4.0ms
(* (log (- re)) (log base))
4.0ms
(* (log base) (log base))

series433.0ms

Counts
3 → 9
Calls
3 calls:
Slowest
227.0ms
(* (log base) (log base))
120.0ms
(/ (* (log (- re)) (log base)) (* (log base) (log base)))
86.0ms
(* (log (- re)) (log base))

simplify3.5s

Counts
40 → 68
Calls
40 calls:
Slowest
652.0ms
(pow (- (log -1) (log (/ -1 base))) 2)
482.0ms
(- (+ (log (log (- re))) (log (log base))) (+ (log (log base)) (log (log base))))
350.0ms
(* -1 (* (log (/ -1 re)) (- (log -1) (log (/ -1 base)))))
293.0ms
(- (log (* (log (- re)) (log base))) (log (* (log base) (log base))))
269.0ms
(- (log (* (log (- re)) (log base))) (+ (log (log base)) (log (log base))))

prune741.0ms

Pruning

25 alts after pruning (24 fresh and 1 done)

Merged error: 6.2b

regimes885.0ms

Accuracy

52.4% (10.3b remaining)

Error of 17.6b against oracle of 7.3b and baseline of 29.0b

bsearch2.5s

end0.0ms

sample35.9s

Algorithm
intervals
Results
27.3s8052×body10240exit
5.1s16054×body80nan
2.6s8000×body80valid