Average Error: 33.3 → 33.4
Time: 38.7s
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(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right) \bmod a\right)}\right) \cdot \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) \bmod a\right)}} \cdot \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) \bmod a\right)}}\right) \cdot \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) \bmod a\right)}}\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(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right) \bmod a\right)}\right) \cdot \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) \bmod a\right)}} \cdot \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) \bmod a\right)}}\right) \cdot \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) \bmod a\right)}}\right)\right|
double f(double a) {
        double r13479 = a;
        double r13480 = expm1(r13479);
        double r13481 = sin(r13480);
        double r13482 = expm1(r13481);
        double r13483 = atan(r13479);
        double r13484 = atan2(r13482, r13483);
        double r13485 = fmod(r13484, r13479);
        double r13486 = fabs(r13485);
        return r13486;
}
↓
double f(double a) {
        double r13487 = a;
        double r13488 = expm1(r13487);
        double r13489 = sin(r13488);
        double r13490 = expm1(r13489);
        double r13491 = atan(r13487);
        double r13492 = atan2(r13490, r13491);
        double r13493 = fmod(r13492, r13487);
        double r13494 = cbrt(r13493);
        double r13495 = r13494 * r13494;
        double r13496 = cbrt(r13494);
        double r13497 = r13496 * r13496;
        double r13498 = r13497 * r13496;
        double r13499 = r13495 * r13498;
        double r13500 = fabs(r13499);
        return r13500;
}

Error

Bits error versus a

Derivation

  1. Initial program 33.3

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

  6. Final simplification33.4

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

Reproduce

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