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

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 30.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 simplification0.4

    \[\leadsto \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. Final simplification0.4

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

Reproduce

herbie shell --seed 2018362 +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: 10.3s)Debug log

start1.4s

Algorithm
intervals

setup131.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 0.3b

localize18.0ms

Local error

Found 2 expressions with local error:

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

rewrite5.0ms

Algorithm
rewrite-expression-head
Counts
2 → 43
Calls

2 calls. Slowest were:

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

series132.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

104.0ms
(/ (log (hypot re im)) (log base))
28.0ms
(hypot re im)

simplify577.0ms

Counts
37 → 49
Calls

37 calls. Slowest were:

265.0ms
(* -1 (/ (log (/ -1 re)) (- (log -1) (log (/ -1 base)))))
146.0ms
(/ (* (* (log (hypot re im)) (log (hypot re im))) (log (hypot re im))) (* (* (log base) (log base)) (log base)))
34.0ms
(/ (* (cbrt (log (hypot re im))) (cbrt (log (hypot re im)))) (* (cbrt (log base)) (cbrt (log base))))

prune632.0ms

Pruning

13 alts after pruning (13 fresh and 0 done)

Merged error: 0.0b

localize9.0ms

Local error

Found 2 expressions with local error:

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

rewrite13.0ms

Algorithm
rewrite-expression-head
Counts
2 → 43
Calls

2 calls. Slowest were:

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

series113.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

90.0ms
(/ (log (* 1 (hypot re im))) (log base))
23.0ms
(hypot re im)

simplify1.0s

Counts
37 → 49
Calls

37 calls. Slowest were:

366.0ms
(/ (* (* (log (* 1 (hypot re im))) (log (* 1 (hypot re im)))) (log (* 1 (hypot re im)))) (* (* (log base) (log base)) (log base)))
263.0ms
(* -1 (/ (log (/ -1 re)) (- (log -1) (log (/ -1 base)))))
28.0ms
(/ (* (cbrt (log (* 1 (hypot re im)))) (cbrt (log (* 1 (hypot re im))))) 1)

prune633.0ms

Pruning

13 alts after pruning (12 fresh and 1 done)

Merged error: 0.0b

localize7.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)

rewrite8.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)

series373.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

223.0ms
(/ 1 (log base))
111.0ms
(* (log (hypot re im)) (/ 1 (log base)))
38.0ms
(hypot re im)

simplify757.0ms

Counts
25 → 61
Calls

25 calls. Slowest were:

366.0ms
(* (* (* (log (hypot re im)) (log (hypot re im))) (log (hypot re im))) (* (* (/ 1 (log base)) (/ 1 (log base))) (/ 1 (log base))))
229.0ms
(* -1 (/ (log (/ -1 re)) (- (log -1) (log (/ -1 base)))))
45.0ms
(/ 1 (- (log -1) (log (/ -1 base))))

prune852.0ms

Pruning

14 alts after pruning (13 fresh and 1 done)

Merged error: 0.0b

localize6.0ms

Local error

Found 3 expressions with local error:

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

rewrite6.0ms

Algorithm
rewrite-expression-head
Counts
3 → 76
Calls

3 calls. Slowest were:

3.0ms
(/ (log base) (log (hypot re im)))
2.0ms
(/ 1 (/ (log base) (log (hypot re im))))
0.0ms
(hypot re im)

series326.0ms

Counts
3 → 9
Calls

3 calls. Slowest were:

152.0ms
(/ (log base) (log (hypot re im)))
135.0ms
(/ 1 (/ (log base) (log (hypot re im))))
38.0ms
(hypot re im)

simplify1.7s

Counts
59 → 85
Calls

59 calls. Slowest were:

1.1s
(* -1 (/ (- (log -1) (log (/ -1 base))) (log (/ -1 re))))
228.0ms
(* -1 (/ (log (/ -1 re)) (- (log -1) (log (/ -1 base)))))
152.0ms
(/ (* (* (log base) (log base)) (log base)) (* (* (log (hypot re im)) (log (hypot re im))) (log (hypot re im))))

prune966.0ms

Pruning

13 alts after pruning (12 fresh and 1 done)

Merged error: 0b

regimes660.0ms

Accuracy

0% (0.4b remaining)

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

bsearch1.0ms