Average Error: 38.9 → 1.3
Time: 1.7m
Precision: 64
Internal precision: 2432
\[\cos \left(x + \varepsilon\right) - \cos x\]
⬇
\[\begin{array}{l}
\mathbf{if}\;\varepsilon \le -3.645937152382937 \cdot 10^{+19}:\\
\;\;\;\;\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right) - \cos x\\
\mathbf{if}\;\varepsilon \le 4.729663737457019 \cdot 10^{-05}:\\
\;\;\;\;-2 \cdot \left(\sin \left(\frac{\varepsilon}{2}\right) \cdot \sin \left(\frac{\left(x + \varepsilon\right) + x}{2}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\cos x \cdot \cos \varepsilon - \left(\sin x \cdot \sin \varepsilon + \cos x\right)\\
\end{array}\]
Derivation
- Split input into 3 regimes.
-
if eps < -3.645937152382937e+19
Initial program 29.2
\[\cos \left(x + \varepsilon\right) - \cos x\]
- Using strategy
rm
Applied cos-sum 0.8
\[\leadsto \color{blue}{\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right)} - \cos x\]
if -3.645937152382937e+19 < eps < 4.729663737457019e-05
Initial program 48.4
\[\cos \left(x + \varepsilon\right) - \cos x\]
- Using strategy
rm
Applied diff-cos 37.7
\[\leadsto \color{blue}{-2 \cdot \left(\sin \left(\frac{\left(x + \varepsilon\right) - x}{2}\right) \cdot \sin \left(\frac{\left(x + \varepsilon\right) + x}{2}\right)\right)}\]
Applied simplify 1.7
\[\leadsto -2 \cdot \color{blue}{\left(\sin \left(\frac{\varepsilon}{2}\right) \cdot \sin \left(\frac{\left(x + \varepsilon\right) + x}{2}\right)\right)}\]
if 4.729663737457019e-05 < eps
Initial program 29.3
\[\cos \left(x + \varepsilon\right) - \cos x\]
- Using strategy
rm
Applied cos-sum 0.9
\[\leadsto \color{blue}{\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right)} - \cos x\]
Applied associate--l- 0.9
\[\leadsto \color{blue}{\cos x \cdot \cos \varepsilon - \left(\sin x \cdot \sin \varepsilon + \cos x\right)}\]
- Recombined 3 regimes into one program.
- Removed slow pow expressions
Runtime
Please include this information when filing a bug report:
herbie shell --seed '#(772101555 1905824529 294602591 2478279198 2123125427 4197813737)'
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
(- (cos (+ x eps)) (cos x)))