Average Error: 33.6 → 33.6
Time: 47.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(\tan^{-1}_* \frac{\left(\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)} \cdot \sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}\right) \cdot \left(\left(\sqrt[3]{\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}} \cdot \sqrt[3]{\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}}\right) \cdot \sqrt[3]{\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\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(\tan^{-1}_* \frac{\left(\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)} \cdot \sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}\right) \cdot \left(\left(\sqrt[3]{\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}} \cdot \sqrt[3]{\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}}\right) \cdot \sqrt[3]{\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}}\right)}{\tan^{-1} a}\right) \bmod a\right)\right|
double f(double a) {
        double r30849 = a;
        double r30850 = expm1(r30849);
        double r30851 = sin(r30850);
        double r30852 = expm1(r30851);
        double r30853 = atan(r30849);
        double r30854 = atan2(r30852, r30853);
        double r30855 = fmod(r30854, r30849);
        double r30856 = fabs(r30855);
        return r30856;
}
↓
double f(double a) {
        double r30857 = a;
        double r30858 = expm1(r30857);
        double r30859 = sin(r30858);
        double r30860 = expm1(r30859);
        double r30861 = cbrt(r30860);
        double r30862 = r30861 * r30861;
        double r30863 = cbrt(r30861);
        double r30864 = r30863 * r30863;
        double r30865 = r30864 * r30863;
        double r30866 = r30862 * r30865;
        double r30867 = atan(r30857);
        double r30868 = atan2(r30866, r30867);
        double r30869 = fmod(r30868, r30857);
        double r30870 = fabs(r30869);
        return r30870;
}

Error

Bits error versus a

Derivation

  1. Initial program 33.6

    \[\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.6

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

  5. Applied add-cube-cbrt33.6

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

  6. Final simplification33.6

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

Reproduce

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