Average Error: 58.6 → 0.0
Time: 11.3s
Precision: 64
Internal Precision: 128
\[\frac{1}{2} \cdot \log \left(\frac{1 + x}{1 - x}\right)\]
↓
\[(\left(\log_* (1 + \left(-x\right))\right) \cdot \frac{-1}{2} + \left(\frac{\log_* (1 + x)}{2}\right))_*\]

Error

Bits error versus x

Derivation

  1. Initial program 58.6

    \[\frac{1}{2} \cdot \log \left(\frac{1 + x}{1 - x}\right)\]

  2. Simplified50.5

    \[\leadsto \color{blue}{(\left(\log \left(1 - x\right)\right) \cdot \frac{-1}{2} + \left(\frac{\log_* (1 + x)}{2}\right))_*}\]
  3. Using strategy rm

  4. Applied sub-neg50.5

    \[\leadsto (\left(\log \color{blue}{\left(1 + \left(-x\right)\right)}\right) \cdot \frac{-1}{2} + \left(\frac{\log_* (1 + x)}{2}\right))_*\]

  5. Applied log1p-def0.0

    \[\leadsto (\color{blue}{\left(\log_* (1 + \left(-x\right))\right)} \cdot \frac{-1}{2} + \left(\frac{\log_* (1 + x)}{2}\right))_*\]

  6. Final simplification0.0

    \[\leadsto (\left(\log_* (1 + \left(-x\right))\right) \cdot \frac{-1}{2} + \left(\frac{\log_* (1 + x)}{2}\right))_*\]

Reproduce

herbie shell --seed 2019022 +o rules:numerics
(FPCore (x)
  :name "Hyperbolic arc-(co)tangent"
  (* (/ 1 2) (log (/ (+ 1 x) (- 1 x)))))

Details

Time bar (total: 10.6s)Debug log

sample208.0ms

Algorithm
intervals
Results
92.0ms115×body1280valid
65.0ms83×body640valid
22.0ms246×body80nan
13.0ms37×body320valid
4.0ms16×body160valid
1.0ms5×body80valid

simplify62.0ms

Counts
1 → 1
Calls
1 calls:
Slowest
62.0ms
(* (/ 1 2) (log (/ (+ 1 x) (- 1 x))))

prune23.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 51.0b

localize22.0ms

Local error

Found 2 expressions with local error:

6.2b
(log (- 1 x))
0.0b
(fma (log (- 1 x)) -1/2 (/ (log1p x) 2))

rewrite1.0ms

Algorithm
rewrite-expression-head
Rules
3×add-cube-cbrt
3×add-exp-log
3×*-un-lft-identity
3×log-prod
3×pow1
3×add-sqr-sqrt
2×add-log-exp
2×log1p-expm1-u
2×add-cbrt-cube
2×log1p-def
2×log-div
2×expm1-log1p-u
1×flip--
1×fma-udef
1×log-pow
1×rem-log-exp
1×flip3--
1×sub-neg
Counts
2 → 28
Calls
2 calls:
Slowest
1.0ms
(log (- 1 x))
0.0ms
(fma (log (- 1 x)) -1/2 (/ (log1p x) 2))

series185.0ms

Counts
2 → 6
Calls
2 calls:
Slowest
162.0ms
(fma (log (- 1 x)) -1/2 (/ (log1p x) 2))
23.0ms
(log (- 1 x))

simplify992.0ms

Counts
14 → 34
Calls
14 calls:
Slowest
277.0ms
(- (log -1) (+ (log (/ 1 x)) (+ (* 1/2 (/ 1 (pow x 2))) (/ 1 x))))
250.0ms
(- (+ x (+ (* 1/3 (pow x 3)) (* 1/2 (pow x 2)))))
230.0ms
(- (+ (log (/ -1 x)) (+ (* 1/2 (/ 1 (pow x 2))) (/ 1 x))))
75.0ms
(- (+ (* 1/3 (/ 1 (pow x 3))) (/ 1 x)) (* 1/2 (log -1)))
64.0ms
(+ (/ 1 x) (+ (* 1/3 (/ 1 (pow x 3))) (* 1/2 (log -1))))

prune186.0ms

Pruning

2 alts after pruning (2 fresh and 0 done)

Merged error: 0b

localize11.0ms

Local error

Found 1 expressions with local error:

0.0b
(fma (log1p (- x)) -1/2 (/ (log1p x) 2))

rewrite0.0ms

Algorithm
rewrite-expression-head
Rules
1×add-log-exp
1×log1p-expm1-u
1×fma-udef
1×add-cube-cbrt
1×add-exp-log
1×add-cbrt-cube
1×*-un-lft-identity
1×pow1
1×expm1-log1p-u
1×add-sqr-sqrt
Counts
1 → 10
Calls
1 calls:
Slowest
0.0ms
(fma (log1p (- x)) -1/2 (/ (log1p x) 2))

series130.0ms

Counts
1 → 3
Calls
1 calls:
Slowest
130.0ms
(fma (log1p (- x)) -1/2 (/ (log1p x) 2))

simplify220.0ms

Counts
3 → 13
Calls
3 calls:
Slowest
94.0ms
(- (+ (* 1/3 (/ 1 (pow x 3))) (/ 1 x)) (* 1/2 (log -1)))
65.0ms
(+ (/ 1 x) (+ (* 1/3 (/ 1 (pow x 3))) (* 1/2 (log -1))))
60.0ms
(+ x (+ (* 1/3 (pow x 3)) (* 1/5 (pow x 5))))

prune98.0ms

Pruning

2 alts after pruning (1 fresh and 1 done)

Merged error: 0b

localize35.0ms

Local error

Found 3 expressions with local error:

0.1b
(* 1/3 (pow x 3))
0.1b
(* 1/5 (pow x 5))
0.0b
(+ (* 1/3 (pow x 3)) (* 1/5 (pow x 5)))

rewrite12.0ms

Algorithm
rewrite-expression-head
Rules
8×associate-*r*
5×add-log-exp
5×add-cube-cbrt
5×*-un-lft-identity
5×add-sqr-sqrt
3×log1p-expm1-u
3×add-exp-log
3×add-cbrt-cube
3×pow1
3×expm1-log1p-u
2×*-commutative
1×flip-+
1×unpow3
1×flip3-+
1×fma-def
1×sum-log
1×+-commutative
1×cube-mult
Counts
3 → 42
Calls
3 calls:
Slowest
6.0ms
(+ (* 1/3 (pow x 3)) (* 1/5 (pow x 5)))
2.0ms
(* 1/5 (pow x 5))
2.0ms
(* 1/3 (pow x 3))

series126.0ms

Counts
3 → 9
Calls
3 calls:
Slowest
73.0ms
(+ (* 1/3 (pow x 3)) (* 1/5 (pow x 5)))
27.0ms
(* 1/3 (pow x 3))
25.0ms
(* 1/5 (pow x 5))

simplify140.0ms

Counts
17 → 51
Calls
17 calls:
Slowest
35.0ms
(* (exp (* 1/3 (pow x 3))) (exp (* 1/5 (pow x 5))))
22.0ms
(+ (* 1/3 (pow x 3)) (* 1/5 (pow x 5)))
13.0ms
(+ (* 1/3 (pow x 3)) (* 1/5 (pow x 5)))
12.0ms
(+ (* 1/3 (pow x 3)) (* 1/5 (pow x 5)))
10.0ms
(* 1/3 (pow x 3))

prune562.0ms

Pruning

2 alts after pruning (0 fresh and 2 done)

Merged error: 0b

regimes33.0ms

Accuracy

0% (0.0b remaining)

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

bsearch1.0ms

end0.0ms

sample7.6s

Algorithm
intervals
Results
4.0s3409×body1280valid
1.6s2484×body640valid
922.0ms7987×body80nan
597.0ms1300×body320valid
202.0ms604×body160valid
24.0ms203×body80valid