Average Error: 33.6 → 33.6
Time: 17.4s
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(\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|\]
\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(\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|
double f(double a) {
        double r6027 = a;
        double r6028 = expm1(r6027);
        double r6029 = sin(r6028);
        double r6030 = expm1(r6029);
        double r6031 = atan(r6027);
        double r6032 = atan2(r6030, r6031);
        double r6033 = fmod(r6032, r6027);
        double r6034 = fabs(r6033);
        return r6034;
}
↓
double f(double a) {
        double r6035 = a;
        double r6036 = expm1(r6035);
        double r6037 = sin(r6036);
        double r6038 = expm1(r6037);
        double r6039 = atan(r6035);
        double r6040 = atan2(r6038, r6039);
        double r6041 = fmod(r6040, r6035);
        double r6042 = cbrt(r6041);
        double r6043 = 3.0;
        double r6044 = pow(r6040, r6043);
        double r6045 = cbrt(r6044);
        double r6046 = fmod(r6045, r6035);
        double r6047 = cbrt(r6046);
        double r6048 = r6042 * r6047;
        double r6049 = r6048 * r6047;
        double r6050 = fabs(r6049);
        return r6050;
}

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-cbrt-cube33.7

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

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

    \[\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. Taylor expanded around 0 33.6

    \[\leadsto \left|\left(\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)}} \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|\]

  8. Final simplification33.6

    \[\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(\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|\]

Reproduce

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