Average Error: 33.2 → 33.3
Time: 36.3s
Precision: 64
\[\left|\left(\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right) \bmod a\right)\right|\]
↓
\[\left|\left(\left(\left(\sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}} \cdot \sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}\right) \cdot \sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}\right) \bmod a\right)\right|\]
\left|\left(\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right) \bmod a\right)\right|
↓
\left|\left(\left(\left(\sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}} \cdot \sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}\right) \cdot \sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}\right) \bmod a\right)\right|
double f(double a) {
        double r16558 = a;
        double r16559 = expm1(r16558);
        double r16560 = sin(r16559);
        double r16561 = expm1(r16560);
        double r16562 = atan(r16558);
        double r16563 = atan2(r16561, r16562);
        double r16564 = fmod(r16563, r16558);
        double r16565 = fabs(r16564);
        return r16565;
}
↓
double f(double a) {
        double r16566 = a;
        double r16567 = expm1(r16566);
        double r16568 = sin(r16567);
        double r16569 = expm1(r16568);
        double r16570 = atan(r16566);
        double r16571 = atan2(r16569, r16570);
        double r16572 = cbrt(r16571);
        double r16573 = r16572 * r16572;
        double r16574 = r16573 * r16572;
        double r16575 = fmod(r16574, r16566);
        double r16576 = fabs(r16575);
        return r16576;
}

Error

Bits error versus a

Derivation

  1. Initial program 33.2

    \[\left|\left(\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right) \bmod a\right)\right|\]
  2. Using strategy rm

  3. Applied add-cube-cbrt33.3

    \[\leadsto \left|\left(\color{blue}{\left(\left(\sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}} \cdot \sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}\right) \cdot \sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}\right)} \bmod a\right)\right|\]

  4. Final simplification33.3

    \[\leadsto \left|\left(\left(\left(\sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}} \cdot \sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}\right) \cdot \sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}\right) \bmod a\right)\right|\]

Reproduce

herbie shell --seed 2019325 
(FPCore (a)
  :name "Random Jason Timeout Test 006"
  :precision binary64
  (fabs (fmod (atan2 (expm1 (sin (expm1 a))) (atan a)) a)))