Average Error: 33.6 → 33.7
Time: 39.8s
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(\log \left(\sqrt{e^{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}}\right) + \log \left(\sqrt{e^{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}}\right)\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(\log \left(\sqrt{e^{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}}\right) + \log \left(\sqrt{e^{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}}\right)\right) \bmod a\right)\right|
double f(double a) {
        double r19062 = a;
        double r19063 = expm1(r19062);
        double r19064 = sin(r19063);
        double r19065 = expm1(r19064);
        double r19066 = atan(r19062);
        double r19067 = atan2(r19065, r19066);
        double r19068 = fmod(r19067, r19062);
        double r19069 = fabs(r19068);
        return r19069;
}
↓
double f(double a) {
        double r19070 = a;
        double r19071 = expm1(r19070);
        double r19072 = sin(r19071);
        double r19073 = expm1(r19072);
        double r19074 = atan(r19070);
        double r19075 = atan2(r19073, r19074);
        double r19076 = exp(r19075);
        double r19077 = sqrt(r19076);
        double r19078 = log(r19077);
        double r19079 = r19078 + r19078;
        double r19080 = fmod(r19079, r19070);
        double r19081 = fabs(r19080);
        return r19081;
}

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-log-exp33.7

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

  5. Applied add-sqr-sqrt33.7

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

  6. Applied log-prod33.7

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

  7. Final simplification33.7

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

Reproduce

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