Average Error: 33.9 → 34.0
Time: 26.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{\mathsf{expm1}\left(\left({\left(\sqrt[3]{\sqrt[3]{\sin \left(\mathsf{expm1}\left(a\right)\right)}}\right)}^{5} \cdot \sqrt[3]{\sqrt[3]{\sin \left(\mathsf{expm1}\left(a\right)\right)}}\right) \cdot \left({\left(\sqrt[3]{\sqrt[3]{\sin \left(\mathsf{expm1}\left(a\right)\right)}}\right)}^{2} \cdot \sqrt[3]{\sqrt[3]{\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{\mathsf{expm1}\left(\left({\left(\sqrt[3]{\sqrt[3]{\sin \left(\mathsf{expm1}\left(a\right)\right)}}\right)}^{5} \cdot \sqrt[3]{\sqrt[3]{\sin \left(\mathsf{expm1}\left(a\right)\right)}}\right) \cdot \left({\left(\sqrt[3]{\sqrt[3]{\sin \left(\mathsf{expm1}\left(a\right)\right)}}\right)}^{2} \cdot \sqrt[3]{\sqrt[3]{\sin \left(\mathsf{expm1}\left(a\right)\right)}}\right)\right)}{\tan^{-1} a}\right) \bmod a\right)\right|
double f(double a) {
        double r21599 = a;
        double r21600 = expm1(r21599);
        double r21601 = sin(r21600);
        double r21602 = expm1(r21601);
        double r21603 = atan(r21599);
        double r21604 = atan2(r21602, r21603);
        double r21605 = fmod(r21604, r21599);
        double r21606 = fabs(r21605);
        return r21606;
}

double f(double a) {
        double r21607 = a;
        double r21608 = expm1(r21607);
        double r21609 = sin(r21608);
        double r21610 = cbrt(r21609);
        double r21611 = cbrt(r21610);
        double r21612 = 5.0;
        double r21613 = pow(r21611, r21612);
        double r21614 = r21613 * r21611;
        double r21615 = 2.0;
        double r21616 = pow(r21611, r21615);
        double r21617 = r21616 * r21611;
        double r21618 = r21614 * r21617;
        double r21619 = expm1(r21618);
        double r21620 = atan(r21607);
        double r21621 = atan2(r21619, r21620);
        double r21622 = fmod(r21621, r21607);
        double r21623 = fabs(r21622);
        return r21623;
}

Error

Bits error versus a

Derivation

  1. Initial program 33.9

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

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

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

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

    \[\leadsto \left|\left(\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\left(\color{blue}{{\left(\sqrt[3]{\sqrt[3]{\sin \left(\mathsf{expm1}\left(a\right)\right)}}\right)}^{5}} \cdot \sqrt[3]{\sqrt[3]{\sin \left(\mathsf{expm1}\left(a\right)\right)}}\right) \cdot \sqrt[3]{\sin \left(\mathsf{expm1}\left(a\right)\right)}\right)}{\tan^{-1} a}\right) \bmod a\right)\right|\]
  8. Using strategy rm
  9. Applied add-cube-cbrt34.0

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

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

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

Reproduce

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