Average Error: 36.9 → 0.4
Time: 6.2s
Precision: 64
\[\sin \left(x + \varepsilon\right) - \sin x\]
↓
\[\mathsf{fma}\left(\sin x, \log \left(e^{\frac{{\left(\cos \varepsilon\right)}^{3} - 1}{\mathsf{fma}\left(\cos \varepsilon, \cos \varepsilon + 1, 1\right)}}\right), \cos x \cdot \sin \varepsilon\right)\]
\sin \left(x + \varepsilon\right) - \sin x
↓
\mathsf{fma}\left(\sin x, \log \left(e^{\frac{{\left(\cos \varepsilon\right)}^{3} - 1}{\mathsf{fma}\left(\cos \varepsilon, \cos \varepsilon + 1, 1\right)}}\right), \cos x \cdot \sin \varepsilon\right)
double f(double x, double eps) {
        double r137186 = x;
        double r137187 = eps;
        double r137188 = r137186 + r137187;
        double r137189 = sin(r137188);
        double r137190 = sin(r137186);
        double r137191 = r137189 - r137190;
        return r137191;
}
↓
double f(double x, double eps) {
        double r137192 = x;
        double r137193 = sin(r137192);
        double r137194 = eps;
        double r137195 = cos(r137194);
        double r137196 = 3.0;
        double r137197 = pow(r137195, r137196);
        double r137198 = 1.0;
        double r137199 = r137197 - r137198;
        double r137200 = r137195 + r137198;
        double r137201 = fma(r137195, r137200, r137198);
        double r137202 = r137199 / r137201;
        double r137203 = exp(r137202);
        double r137204 = log(r137203);
        double r137205 = cos(r137192);
        double r137206 = sin(r137194);
        double r137207 = r137205 * r137206;
        double r137208 = fma(r137193, r137204, r137207);
        return r137208;
}

Error

Bits error versus x

Bits error versus eps

Target

Original36.9
Target15.0
Herbie0.4
\[2 \cdot \left(\cos \left(x + \frac{\varepsilon}{2}\right) \cdot \sin \left(\frac{\varepsilon}{2}\right)\right)\]

Derivation

  1. Initial program 36.9

    \[\sin \left(x + \varepsilon\right) - \sin x\]
  2. Using strategy rm

  3. Applied sin-sum21.7

    \[\leadsto \color{blue}{\left(\sin x \cdot \cos \varepsilon + \cos x \cdot \sin \varepsilon\right)} - \sin x\]

  4. Taylor expanded around inf 21.7

    \[\leadsto \color{blue}{\left(\sin \varepsilon \cdot \cos x + \sin x \cdot \cos \varepsilon\right) - \sin x}\]

  5. Simplified0.4

    \[\leadsto \color{blue}{\mathsf{fma}\left(\sin x, \cos \varepsilon - 1, \cos x \cdot \sin \varepsilon\right)}\]
  6. Using strategy rm

  7. Applied add-log-exp0.4

    \[\leadsto \mathsf{fma}\left(\sin x, \cos \varepsilon - \color{blue}{\log \left(e^{1}\right)}, \cos x \cdot \sin \varepsilon\right)\]

  8. Applied add-log-exp0.4

    \[\leadsto \mathsf{fma}\left(\sin x, \color{blue}{\log \left(e^{\cos \varepsilon}\right)} - \log \left(e^{1}\right), \cos x \cdot \sin \varepsilon\right)\]

  9. Applied diff-log0.4

    \[\leadsto \mathsf{fma}\left(\sin x, \color{blue}{\log \left(\frac{e^{\cos \varepsilon}}{e^{1}}\right)}, \cos x \cdot \sin \varepsilon\right)\]

  10. Simplified0.4

    \[\leadsto \mathsf{fma}\left(\sin x, \log \color{blue}{\left(e^{\cos \varepsilon - 1}\right)}, \cos x \cdot \sin \varepsilon\right)\]
  11. Using strategy rm

  12. Applied flip3--0.4

    \[\leadsto \mathsf{fma}\left(\sin x, \log \left(e^{\color{blue}{\frac{{\left(\cos \varepsilon\right)}^{3} - {1}^{3}}{\cos \varepsilon \cdot \cos \varepsilon + \left(1 \cdot 1 + \cos \varepsilon \cdot 1\right)}}}\right), \cos x \cdot \sin \varepsilon\right)\]

  13. Simplified0.4

    \[\leadsto \mathsf{fma}\left(\sin x, \log \left(e^{\frac{\color{blue}{{\left(\cos \varepsilon\right)}^{3} - 1}}{\cos \varepsilon \cdot \cos \varepsilon + \left(1 \cdot 1 + \cos \varepsilon \cdot 1\right)}}\right), \cos x \cdot \sin \varepsilon\right)\]

  14. Simplified0.4

    \[\leadsto \mathsf{fma}\left(\sin x, \log \left(e^{\frac{{\left(\cos \varepsilon\right)}^{3} - 1}{\color{blue}{\mathsf{fma}\left(\cos \varepsilon, \cos \varepsilon + 1, 1\right)}}}\right), \cos x \cdot \sin \varepsilon\right)\]

  15. Final simplification0.4

    \[\leadsto \mathsf{fma}\left(\sin x, \log \left(e^{\frac{{\left(\cos \varepsilon\right)}^{3} - 1}{\mathsf{fma}\left(\cos \varepsilon, \cos \varepsilon + 1, 1\right)}}\right), \cos x \cdot \sin \varepsilon\right)\]

Reproduce

herbie shell --seed 2020001 +o rules:numerics
(FPCore (x eps)
  :name "2sin (example 3.3)"
  :precision binary64

  :herbie-target
  (* 2 (* (cos (+ x (/ eps 2))) (sin (/ eps 2))))

  (- (sin (+ x eps)) (sin x)))