Average Error: 33.7 → 33.7
Time: 17.2s
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(\sqrt[3]{{\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right)}^{3}}\right) \bmod a\right)} \cdot \sqrt[3]{\left(\left(\sqrt[3]{{\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right)}^{3}}\right) \bmod a\right)}\right) \cdot \sqrt[3]{\left(\left(\sqrt[3]{\sqrt[3]{{\left({\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right)}^{3}\right)}^{3}}}\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(\sqrt[3]{{\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right)}^{3}}\right) \bmod a\right)} \cdot \sqrt[3]{\left(\left(\sqrt[3]{{\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right)}^{3}}\right) \bmod a\right)}\right) \cdot \sqrt[3]{\left(\left(\sqrt[3]{\sqrt[3]{{\left({\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right)}^{3}\right)}^{3}}}\right) \bmod a\right)}\right|
double f(double a) {
        double r7007 = a;
        double r7008 = expm1(r7007);
        double r7009 = sin(r7008);
        double r7010 = expm1(r7009);
        double r7011 = atan(r7007);
        double r7012 = atan2(r7010, r7011);
        double r7013 = fmod(r7012, r7007);
        double r7014 = fabs(r7013);
        return r7014;
}
↓
double f(double a) {
        double r7015 = a;
        double r7016 = expm1(r7015);
        double r7017 = sin(r7016);
        double r7018 = expm1(r7017);
        double r7019 = atan(r7015);
        double r7020 = atan2(r7018, r7019);
        double r7021 = 3.0;
        double r7022 = pow(r7020, r7021);
        double r7023 = cbrt(r7022);
        double r7024 = fmod(r7023, r7015);
        double r7025 = cbrt(r7024);
        double r7026 = r7025 * r7025;
        double r7027 = pow(r7022, r7021);
        double r7028 = cbrt(r7027);
        double r7029 = cbrt(r7028);
        double r7030 = fmod(r7029, r7015);
        double r7031 = cbrt(r7030);
        double r7032 = r7026 * r7031;
        double r7033 = fabs(r7032);
        return r7033;
}

Error

Bits error versus a

Derivation

  1. Initial program 33.7

    \[\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-cbrt-cube33.8

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

  4. Simplified33.8

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

  6. Applied add-cube-cbrt33.8

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

  8. Applied add-cbrt-cube33.7

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

  9. Simplified33.7

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

  10. Final simplification33.7

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

Reproduce

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