Average Error: 33.4 → 33.3
Time: 46.0s
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(\sqrt[3]{\left(\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right) \bmod a\right)} \cdot \sqrt[3]{\left(\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) \cdot \sqrt[3]{\left(\left(\log \left(\log \left(e^{e^{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}}\right)\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(\sqrt[3]{\left(\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right) \bmod a\right)} \cdot \sqrt[3]{\left(\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) \cdot \sqrt[3]{\left(\left(\log \left(\log \left(e^{e^{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}}\right)\right)\right) \bmod a\right)}\right|
double f(double a) {
        double r938487 = a;
        double r938488 = expm1(r938487);
        double r938489 = sin(r938488);
        double r938490 = expm1(r938489);
        double r938491 = atan(r938487);
        double r938492 = atan2(r938490, r938491);
        double r938493 = fmod(r938492, r938487);
        double r938494 = fabs(r938493);
        return r938494;
}
↓
double f(double a) {
        double r938495 = a;
        double r938496 = expm1(r938495);
        double r938497 = sin(r938496);
        double r938498 = expm1(r938497);
        double r938499 = atan(r938495);
        double r938500 = atan2(r938498, r938499);
        double r938501 = fmod(r938500, r938495);
        double r938502 = cbrt(r938501);
        double r938503 = exp(r938500);
        double r938504 = log(r938503);
        double r938505 = fmod(r938504, r938495);
        double r938506 = cbrt(r938505);
        double r938507 = r938502 * r938506;
        double r938508 = exp(r938503);
        double r938509 = log(r938508);
        double r938510 = log(r938509);
        double r938511 = fmod(r938510, r938495);
        double r938512 = cbrt(r938511);
        double r938513 = r938507 * r938512;
        double r938514 = fabs(r938513);
        return r938514;
}

Error

Bits error versus a

Derivation

  1. Initial program 33.4

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

    \[\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-cube-cbrt33.4

    \[\leadsto \left|\color{blue}{\left(\sqrt[3]{\left(\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)} \cdot \sqrt[3]{\left(\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) \cdot \sqrt[3]{\left(\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|\]
  6. Using strategy rm

  7. Applied add-log-exp33.4

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

  8. Taylor expanded around 0 33.3

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

  9. Final simplification33.3

    \[\leadsto \left|\left(\sqrt[3]{\left(\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right) \bmod a\right)} \cdot \sqrt[3]{\left(\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) \cdot \sqrt[3]{\left(\left(\log \left(\log \left(e^{e^{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}}\right)\right)\right) \bmod a\right)}\right|\]

Reproduce

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