\[\sin \left(x + \varepsilon\right) - \sin x\]
Test:
NMSE example 3.3
Bits:
128 bits
Bits error versus x
Bits error versus eps
Time: 21.3 s
Input Error: 35.5
Output Error: 1.8
Log:
Profile: 🕒
\(\begin{cases} \left(\sin x \cdot \cos \varepsilon + \cos x \cdot \sin \varepsilon\right) - \sin x & \text{when } \varepsilon \le -5.952085980018888 \cdot 10^{-41} \\ \varepsilon - \frac{1}{2} \cdot \left(\left(\varepsilon + x\right) \cdot \left(x \cdot \varepsilon\right)\right) & \text{when } \varepsilon \le 7.687581306962853 \cdot 10^{-33} \\ \left(\sin x \cdot \cos \varepsilon + \cos x \cdot \sin \varepsilon\right) - \sin x & \text{otherwise} \end{cases}\)

    if eps < -5.952085980018888e-41 or 7.687581306962853e-33 < eps

    1. Started with
      \[\sin \left(x + \varepsilon\right) - \sin x\]
      29.3
    2. Using strategy rm
      29.3
    3. Applied sin-sum to get
      \[\color{red}{\sin \left(x + \varepsilon\right)} - \sin x \leadsto \color{blue}{\left(\sin x \cdot \cos \varepsilon + \cos x \cdot \sin \varepsilon\right)} - \sin x\]
      2.9

    if -5.952085980018888e-41 < eps < 7.687581306962853e-33

    1. Started with
      \[\sin \left(x + \varepsilon\right) - \sin x\]
      44.4
    2. Applied taylor to get
      \[\sin \left(x + \varepsilon\right) - \sin x \leadsto \varepsilon - \left(\frac{1}{2} \cdot \left(\varepsilon \cdot {x}^2\right) + \frac{1}{2} \cdot \left({\varepsilon}^2 \cdot x\right)\right)\]
      10.0
    3. Taylor expanded around 0 to get
      \[\color{red}{\varepsilon - \left(\frac{1}{2} \cdot \left(\varepsilon \cdot {x}^2\right) + \frac{1}{2} \cdot \left({\varepsilon}^2 \cdot x\right)\right)} \leadsto \color{blue}{\varepsilon - \left(\frac{1}{2} \cdot \left(\varepsilon \cdot {x}^2\right) + \frac{1}{2} \cdot \left({\varepsilon}^2 \cdot x\right)\right)}\]
      10.0
    4. Applied simplify to get
      \[\color{red}{\varepsilon - \left(\frac{1}{2} \cdot \left(\varepsilon \cdot {x}^2\right) + \frac{1}{2} \cdot \left({\varepsilon}^2 \cdot x\right)\right)} \leadsto \color{blue}{\varepsilon - \frac{1}{2} \cdot \left(\left(\varepsilon + x\right) \cdot \left(x \cdot \varepsilon\right)\right)}\]
      0.1

  1. Removed slow pow expressions

Original test:


(lambda ((x default) (eps default))
  #:name "NMSE example 3.3"
  (- (sin (+ x eps)) (sin x))
  #:target
  (* 2 (* (cos (+ x (/ eps 2))) (sin (/ eps 2)))))