Average Error: 33.2 → 33.3
Time: 48.1s
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 r1019704 = a;
        double r1019705 = expm1(r1019704);
        double r1019706 = sin(r1019705);
        double r1019707 = expm1(r1019706);
        double r1019708 = atan(r1019704);
        double r1019709 = atan2(r1019707, r1019708);
        double r1019710 = fmod(r1019709, r1019704);
        double r1019711 = fabs(r1019710);
        return r1019711;
}
↓
double f(double a) {
        double r1019712 = a;
        double r1019713 = expm1(r1019712);
        double r1019714 = sin(r1019713);
        double r1019715 = expm1(r1019714);
        double r1019716 = atan(r1019712);
        double r1019717 = atan2(r1019715, r1019716);
        double r1019718 = cbrt(r1019717);
        double r1019719 = r1019718 * r1019718;
        double r1019720 = r1019719 * r1019718;
        double r1019721 = fmod(r1019720, r1019712);
        double r1019722 = fabs(r1019721);
        return r1019722;
}

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 2019163 
(FPCore (a)
  :name "Random Jason Timeout Test 006"
  (fabs (fmod (atan2 (expm1 (sin (expm1 a))) (atan a)) a)))