Average Error: 33.5 → 33.5
Time: 17.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(\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(\log \left(e^{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}\right)\right) \bmod a\right)}\right) \cdot \sqrt[3]{\left(\left(\log \left(\log \left(e^{e^{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}}\right)\right)\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(\log \left(e^{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}\right)\right) \bmod a\right)}\right) \cdot \sqrt[3]{\left(\left(\log \left(\log \left(e^{e^{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}}\right)\right)\right) \bmod a\right)}\right|
double f(double a) {
        double r5612 = a;
        double r5613 = expm1(r5612);
        double r5614 = sin(r5613);
        double r5615 = expm1(r5614);
        double r5616 = atan(r5612);
        double r5617 = atan2(r5615, r5616);
        double r5618 = fmod(r5617, r5612);
        double r5619 = fabs(r5618);
        return r5619;
}
↓
double f(double a) {
        double r5620 = a;
        double r5621 = expm1(r5620);
        double r5622 = sin(r5621);
        double r5623 = expm1(r5622);
        double r5624 = atan(r5620);
        double r5625 = atan2(r5623, r5624);
        double r5626 = fmod(r5625, r5620);
        double r5627 = cbrt(r5626);
        double r5628 = exp(r5625);
        double r5629 = log(r5628);
        double r5630 = fmod(r5629, r5620);
        double r5631 = cbrt(r5630);
        double r5632 = r5627 * r5631;
        double r5633 = exp(r5628);
        double r5634 = log(r5633);
        double r5635 = log(r5634);
        double r5636 = fmod(r5635, r5620);
        double r5637 = cbrt(r5636);
        double r5638 = r5632 * r5637;
        double r5639 = fabs(r5638);
        return r5639;
}

Error

Bits error versus a

Derivation

  1. Initial program 33.5

    \[\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-log-exp33.5

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

  5. Applied add-cube-cbrt33.5

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

  7. Applied add-log-exp33.5

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

  8. Taylor expanded around 0 33.5

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

  9. Final simplification33.5

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

Reproduce

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