Average Error: 39.5 → 0.9
Time: 26.4s
Precision: 64
Internal precision: 2176
\[\cos \left(x + \varepsilon\right) - \cos x\]
⬇
\[\begin{array}{l}
\mathbf{if}\;\varepsilon \le -5.5164922447739564 \cdot 10^{-14}:\\
\;\;\;\;\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right) - \cos x\\
\mathbf{if}\;\varepsilon \le 1.5989598071946224 \cdot 10^{-07}:\\
\;\;\;\;-2 \cdot \left(\sin \left(\frac{\varepsilon}{2}\right) \cdot \sin \left(\frac{\varepsilon + \left(x + x\right)}{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 < -5.5164922447739564e-14
Initial program 31.9
\[\cos \left(x + \varepsilon\right) - \cos x\]
- Using strategy
rm
Applied cos-sum 2.0
\[\leadsto \color{blue}{\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right)} - \cos x\]
if -5.5164922447739564e-14 < eps < 1.5989598071946224e-07
Initial program 48.8
\[\cos \left(x + \varepsilon\right) - \cos x\]
- Using strategy
rm
Applied diff-cos 37.2
\[\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 0.3
\[\leadsto -2 \cdot \color{blue}{\left(\sin \left(\frac{\varepsilon}{2}\right) \cdot \sin \left(\frac{\varepsilon + \left(x + x\right)}{2}\right)\right)}\]
if 1.5989598071946224e-07 < eps
Initial program 30.0
\[\cos \left(x + \varepsilon\right) - \cos x\]
- Using strategy
rm
Applied cos-sum 1.1
\[\leadsto \color{blue}{\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right)} - \cos x\]
Applied associate--l- 1.1
\[\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 '#(1064651971 495577305 2200811460 13024471 864198081 231948279)'
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
(- (cos (+ x eps)) (cos x)))