logq (problem 3.4.3)

Average Error: 58.4 → 0.0

Time: 13.6s

Precision: 64

Math

\[\log \left(\frac{1 - \varepsilon}{1 + \varepsilon}\right)\]

↓

\[\mathsf{log1p}\left(-\varepsilon\right) - \mathsf{log1p}\left(\varepsilon\right)\]

C

double f(double eps) {
        double r1657824 = 1.0;
        double r1657825 = eps;
        double r1657826 = r1657824 - r1657825;
        double r1657827 = r1657824 + r1657825;
        double r1657828 = r1657826 / r1657827;
        double r1657829 = log(r1657828);
        return r1657829;
}

↓

double f(double eps) {
        double r1657830 = eps;
        double r1657831 = -r1657830;
        double r1657832 = log1p(r1657831);
        double r1657833 = log1p(r1657830);
        double r1657834 = r1657832 - r1657833;
        return r1657834;
}

Error

Bits error versus eps

Target

Original	58.4
Target	0.3
Herbie	0.0

\[-2 \cdot \left(\left(\varepsilon + \frac{{\varepsilon}^{3}}{3}\right) + \frac{{\varepsilon}^{5}}{5}\right)\]

Derivation

1. Initial program 58.4

   \[\log \left(\frac{1 - \varepsilon}{1 + \varepsilon}\right)\]
2. Using strategy rm

3. Applied log-div 58.4

   \[\leadsto \color{blue}{\log \left(1 - \varepsilon\right) - \log \left(1 + \varepsilon\right)}\]

4. Simplified 50.4

   \[\leadsto \log \left(1 - \varepsilon\right) - \color{blue}{\mathsf{log1p}\left(\varepsilon\right)}\]
5. Using strategy rm

6. Applied sub-neg 50.4

   \[\leadsto \log \color{blue}{\left(1 + \left(-\varepsilon\right)\right)} - \mathsf{log1p}\left(\varepsilon\right)\]

7. Applied log1p-def 0.0

   \[\leadsto \color{blue}{\mathsf{log1p}\left(-\varepsilon\right)} - \mathsf{log1p}\left(\varepsilon\right)\]

8. Final simplification 0.0

   \[\leadsto \mathsf{log1p}\left(-\varepsilon\right) - \mathsf{log1p}\left(\varepsilon\right)\]

Reproduce

herbie shell --seed 2019151 +o rules:numerics
(FPCore (eps)
  :name "logq (problem 3.4.3)"

  :herbie-target
  (* -2 (+ (+ eps (/ (pow eps 3) 3)) (/ (pow eps 5) 5)))

  (log (/ (- 1 eps) (+ 1 eps))))