Average Error: 33.4 → 33.5
Time: 20.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(\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\left(\left(\left(\sqrt[3]{\sin \left(\mathsf{expm1}\left(a\right)\right)} \cdot \sqrt[3]{\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]{\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]{\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|\]
\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(\left(\sqrt[3]{\sin \left(\mathsf{expm1}\left(a\right)\right)} \cdot \sqrt[3]{\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]{\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]{\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|
double f(double a) {
        double r26111 = a;
        double r26112 = expm1(r26111);
        double r26113 = sin(r26112);
        double r26114 = expm1(r26113);
        double r26115 = atan(r26111);
        double r26116 = atan2(r26114, r26115);
        double r26117 = fmod(r26116, r26111);
        double r26118 = fabs(r26117);
        return r26118;
}

double f(double a) {
        double r26119 = a;
        double r26120 = expm1(r26119);
        double r26121 = sin(r26120);
        double r26122 = cbrt(r26121);
        double r26123 = r26122 * r26122;
        double r26124 = cbrt(r26123);
        double r26125 = r26122 * r26124;
        double r26126 = cbrt(r26122);
        double r26127 = r26126 * r26126;
        double r26128 = cbrt(r26127);
        double r26129 = r26125 * r26128;
        double r26130 = cbrt(r26126);
        double r26131 = r26129 * r26130;
        double r26132 = r26131 * r26122;
        double r26133 = expm1(r26132);
        double r26134 = atan(r26119);
        double r26135 = atan2(r26133, r26134);
        double r26136 = fmod(r26135, r26119);
        double r26137 = fabs(r26136);
        return r26137;
}

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

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

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

    \[\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(\sqrt[3]{\sqrt[3]{\sin \left(\mathsf{expm1}\left(a\right)\right)} \cdot \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]{\sin \left(\mathsf{expm1}\left(a\right)\right)}\right)}{\tan^{-1} a}\right) \bmod a\right)\right|\]
  7. Applied associate-*r*33.5

    \[\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 \sqrt[3]{\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]{\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-cbrt33.5

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

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

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

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

Reproduce

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