Average Error: 31.1 → 0.4
Time: 45.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}\]
\[\frac{1}{\log base} \cdot \log \left(\sqrt[3]{\sqrt{re^2 + im^2}^*} \cdot \left(\sqrt[3]{\sqrt{re^2 + im^2}^*} \cdot \sqrt[3]{\sqrt{re^2 + im^2}^*}\right)\right)\]

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. Initial program 31.1

    \[\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. Simplified0.4

    \[\leadsto \color{blue}{\frac{\log \left(\sqrt{re^2 + im^2}^*\right)}{\log base}}\]
  3. Using strategy rm
  4. Applied div-inv0.4

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

    \[\leadsto \log \color{blue}{\left(\left(\sqrt[3]{\sqrt{re^2 + im^2}^*} \cdot \sqrt[3]{\sqrt{re^2 + im^2}^*}\right) \cdot \sqrt[3]{\sqrt{re^2 + im^2}^*}\right)} \cdot \frac{1}{\log base}\]
  7. Final simplification0.4

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

Reproduce

herbie shell --seed 2019008 +o rules:numerics
(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: 43.7s)Debug log

sample1.1s

Algorithm
intervals

simplify43.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

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

prune15.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 0.3b

localize31.0ms

Local error

Found 2 expressions with local error:

0.4b
(/ (log (hypot re im)) (log base))
0.0b
(hypot re im)

rewrite11.0ms

Algorithm
rewrite-expression-head
Counts
2 → 43
Calls

2 calls. Slowest were:

10.0ms
(/ (log (hypot re im)) (log base))
0.0ms
(hypot re im)

series135.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

103.0ms
(/ (log (hypot re im)) (log base))
32.0ms
(hypot re im)

simplify459.0ms

Counts
37 → 49
Calls

37 calls. Slowest were:

236.0ms
(* -1 (/ (log (/ -1 re)) (- (log -1) (log (/ -1 base)))))
119.0ms
(/ (* (* (log (hypot re im)) (log (hypot re im))) (log (hypot re im))) (* (* (log base) (log base)) (log base)))
15.0ms
(/ (* (cbrt (log (hypot re im))) (cbrt (log (hypot re im)))) (* (cbrt (log base)) (cbrt (log base))))

prune568.0ms

Pruning

15 alts after pruning (15 fresh and 0 done)

Merged error: 0.0b

localize5.0ms

Local error

Found 3 expressions with local error:

0.4b
(* (log (hypot re im)) (/ 1 (log base)))
0.3b
(/ 1 (log base))
0.0b
(hypot re im)

rewrite6.0ms

Algorithm
rewrite-expression-head
Counts
3 → 52
Calls

3 calls. Slowest were:

5.0ms
(* (log (hypot re im)) (/ 1 (log base)))
1.0ms
(/ 1 (log base))
0.0ms
(hypot re im)

series309.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

197.0ms
(/ 1 (log base))
79.0ms
(* (log (hypot re im)) (/ 1 (log base)))
33.0ms
(hypot re im)

simplify605.0ms

Counts
25 → 61
Calls

25 calls. Slowest were:

267.0ms
(* (* (* (log (hypot re im)) (log (hypot re im))) (log (hypot re im))) (* (* (/ 1 (log base)) (/ 1 (log base))) (/ 1 (log base))))
204.0ms
(* -1 (/ (log (/ -1 re)) (- (log -1) (log (/ -1 base)))))
47.0ms
(/ 1 (- (log -1) (log (/ -1 base))))

prune692.0ms

Pruning

13 alts after pruning (13 fresh and 0 done)

Merged error: 0.0b

localize19.0ms

Local error

Found 4 expressions with local error:

0.5b
(cbrt (hypot re im))
0.5b
(cbrt (hypot re im))
0.5b
(cbrt (hypot re im))
0.5b
(* (cbrt (hypot re im)) (cbrt (hypot re im)))

rewrite17.0ms

Algorithm
rewrite-expression-head
Counts
4 → 69
Calls

4 calls. Slowest were:

10.0ms
(* (cbrt (hypot re im)) (cbrt (hypot re im)))
4.0ms
(cbrt (hypot re im))
1.0ms
(cbrt (hypot re im))

series418.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

123.0ms
(cbrt (hypot re im))
110.0ms
(cbrt (hypot re im))
97.0ms
(* (cbrt (hypot re im)) (cbrt (hypot re im)))
87.0ms
(cbrt (hypot re im))

simplify151.0ms

Counts
41 → 81
Calls

41 calls. Slowest were:

29.0ms
(+ (* 1/6 (* (pow re 2) (pow (/ 1 (pow im 5)) 1/3))) (pow im 1/3))
28.0ms
(+ (pow im 2/3) (* 1/3 (* (pow re 2) (pow (/ 1 (pow im 4)) 1/3))))
25.0ms
(+ (* 1/6 (* (pow re 2) (pow (/ 1 (pow im 5)) 1/3))) (pow im 1/3))

prune1.2s

Pruning

13 alts after pruning (12 fresh and 1 done)

Merged error: 0.0b

localize21.0ms

Local error

Found 4 expressions with local error:

5.1b
(pow (hypot re im) 1/3)
0.5b
(cbrt (hypot re im))
0.5b
(cbrt (hypot re im))
0.5b
(* (cbrt (hypot re im)) (pow (hypot re im) 1/3))

rewrite16.0ms

Algorithm
rewrite-expression-head
Counts
4 → 71
Calls

4 calls. Slowest were:

10.0ms
(* (cbrt (hypot re im)) (pow (hypot re im) 1/3))
2.0ms
(pow (hypot re im) 1/3)
1.0ms
(cbrt (hypot re im))

series507.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

139.0ms
(cbrt (hypot re im))
124.0ms
(cbrt (hypot re im))
122.0ms
(pow (hypot re im) 1/3)
121.0ms
(* (cbrt (hypot re im)) (pow (hypot re im) 1/3))

simplify301.0ms

Counts
42 → 83
Calls

42 calls. Slowest were:

49.0ms
(+ (* 1/6 (* (pow re 2) (pow (/ 1 (pow im 5)) 1/3))) (pow im 1/3))
48.0ms
(+ (* 1/6 (* (pow re 2) (pow (/ 1 (pow im 5)) 1/3))) (pow im 1/3))
41.0ms
(* (hypot re im) (* (* (pow (hypot re im) 1/3) (pow (hypot re im) 1/3)) (pow (hypot re im) 1/3)))

prune1.0s

Pruning

13 alts after pruning (11 fresh and 2 done)

Merged error: 0.0b

regimes406.0ms

Accuracy

0% (0.4b remaining)

Error of 0.4b against oracle of 0.0b and baseline of 0.4b

bsearch1.0ms

end0.0ms

sample35.7s

Algorithm
intervals