Average Error: 30.9 → 17.3
Time: 4.6s
Precision: 64
Internal Precision: 128
\[\log \left(\sqrt{re \cdot re + im \cdot im}\right)\]
\[\begin{array}{l} \mathbf{if}\;re \le -1.747299818609954 \cdot 10^{+153}:\\ \;\;\;\;\log \left(-re\right)\\ \mathbf{elif}\;re \le 1.0209258191409414 \cdot 10^{+147}:\\ \;\;\;\;\log \left(\sqrt{e^{\log \left(im \cdot im + re \cdot re\right)}}\right)\\ \mathbf{else}:\\ \;\;\;\;\log re\\ \end{array}\]

Error

Bits error versus re

Bits error versus im

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Split input into 3 regimes
  2. if re < -1.747299818609954e+153

    1. Initial program 61.6

      \[\log \left(\sqrt{re \cdot re + im \cdot im}\right)\]
    2. Using strategy rm
    3. Applied add-exp-log61.6

      \[\leadsto \log \left(\sqrt{\color{blue}{e^{\log \left(re \cdot re + im \cdot im\right)}}}\right)\]
    4. Taylor expanded around -inf 6.8

      \[\leadsto \log \color{blue}{\left(-1 \cdot re\right)}\]
    5. Simplified6.8

      \[\leadsto \log \color{blue}{\left(-re\right)}\]

    if -1.747299818609954e+153 < re < 1.0209258191409414e+147

    1. Initial program 20.8

      \[\log \left(\sqrt{re \cdot re + im \cdot im}\right)\]
    2. Using strategy rm
    3. Applied add-exp-log20.8

      \[\leadsto \log \left(\sqrt{\color{blue}{e^{\log \left(re \cdot re + im \cdot im\right)}}}\right)\]

    if 1.0209258191409414e+147 < re

    1. Initial program 60.3

      \[\log \left(\sqrt{re \cdot re + im \cdot im}\right)\]
    2. Using strategy rm
    3. Applied add-exp-log60.3

      \[\leadsto \log \left(\sqrt{\color{blue}{e^{\log \left(re \cdot re + im \cdot im\right)}}}\right)\]
    4. Taylor expanded around inf 6.6

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

    \[\leadsto \begin{array}{l} \mathbf{if}\;re \le -1.747299818609954 \cdot 10^{+153}:\\ \;\;\;\;\log \left(-re\right)\\ \mathbf{elif}\;re \le 1.0209258191409414 \cdot 10^{+147}:\\ \;\;\;\;\log \left(\sqrt{e^{\log \left(im \cdot im + re \cdot re\right)}}\right)\\ \mathbf{else}:\\ \;\;\;\;\log re\\ \end{array}\]

Reproduce

herbie shell --seed 2019005 
(FPCore (re im)
  :name "math.log/1 on complex, real part"
  (log (sqrt (+ (* re re) (* im im)))))

Details

Time bar (total: 4.0s)Debug log

sample27.0ms

Algorithm
intervals

simplify3.0ms

Counts
1 → 1
Calls

1 calls. Slowest were:

3.0ms
(log (sqrt (+ (* re re) (* im im))))

prune5.0ms

Pruning

1 alts after pruning (1 fresh and 0 done)

Merged error: 30.3b

localize18.0ms

Local error

Found 2 expressions with local error:

29.0b
(sqrt (+ (* re re) (* im im)))
0.0b
(+ (* re re) (* im im))

rewrite8.0ms

Algorithm
rewrite-expression-head
Counts
2 → 25
Calls

2 calls. Slowest were:

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

series56.0ms

Counts
2 → 6
Calls

2 calls. Slowest were:

32.0ms
(sqrt (+ (* re re) (* im im)))
24.0ms
(+ (* re re) (* im im))

simplify277.0ms

Counts
13 → 31
Calls

13 calls. Slowest were:

157.0ms
(sqrt (- (* (* re re) (* re re)) (* (* im im) (* im im))))
71.0ms
(sqrt (+ (pow (* re re) 3) (pow (* im im) 3)))
11.0ms
(* (exp (* re re)) (exp (* im im)))

prune311.0ms

Pruning

4 alts after pruning (4 fresh and 0 done)

Merged error: 7.9b

localize21.0ms

Local error

Found 4 expressions with local error:

30.3b
(log (+ (* re re) (* im im)))
29.0b
(sqrt (exp (log (+ (* re re) (* im im)))))
3.0b
(exp (log (+ (* re re) (* im im))))
0.0b
(+ (* re re) (* im im))

rewrite10.0ms

Algorithm
rewrite-expression-head
Counts
4 → 60
Calls

4 calls. Slowest were:

3.0ms
(sqrt (exp (log (+ (* re re) (* im im)))))
3.0ms
(exp (log (+ (* re re) (* im im))))
2.0ms
(+ (* re re) (* im im))

series86.0ms

Counts
4 → 12
Calls

4 calls. Slowest were:

39.0ms
(log (+ (* re re) (* im im)))
18.0ms
(sqrt (exp (log (+ (* re re) (* im im)))))
15.0ms
(exp (log (+ (* re re) (* im im))))
14.0ms
(+ (* re re) (* im im))

simplify752.0ms

Counts
39 → 72
Calls

39 calls. Slowest were:

163.0ms
(log (- (* (* re re) (* re re)) (* (* im im) (* im im))))
151.0ms
(sqrt (exp (log (- (* (* re re) (* re re)) (* (* im im) (* im im))))))
103.0ms
(exp (log (- (* (* re re) (* re re)) (* (* im im) (* im im)))))

prune771.0ms

Pruning

4 alts after pruning (3 fresh and 1 done)

Merged error: 7.9b

localize8.0ms

Local error

Found 0 expressions with local error:

rewrite0.0ms

Algorithm
rewrite-expression-head
Counts
0 → 0
Calls

0 calls. Slowest were:

series0.0ms

Counts
0 → 0
Calls

0 calls. Slowest were:

simplify0.0ms

Counts
0 → 0
Calls

0 calls. Slowest were:

prune6.0ms

Pruning

4 alts after pruning (2 fresh and 2 done)

Merged error: 7.9b

localize5.0ms

Local error

Found 0 expressions with local error:

rewrite0.0ms

Algorithm
rewrite-expression-head
Counts
0 → 0
Calls

0 calls. Slowest were:

series0.0ms

Counts
0 → 0
Calls

0 calls. Slowest were:

simplify0.0ms

Counts
0 → 0
Calls

0 calls. Slowest were:

prune0.0ms

Pruning

4 alts after pruning (1 fresh and 3 done)

Merged error: 7.9b

regimes349.0ms

Accuracy

57.5% (10.0b remaining)

Error of 17.3b against oracle of 7.3b and baseline of 30.9b

bsearch99.0ms

end0.0ms

sample1.2s

Algorithm
intervals