Random Jason Timeout Test 006

Average Error: 33.8 → 33.8

Time: 18.4s

Precision: 64

• could not determine a ground truth for program body

  1. a = 7.41237110397388e+37

Math

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

C

double f(double a) {
        double r9301 = a;
        double r9302 = expm1(r9301);
        double r9303 = sin(r9302);
        double r9304 = expm1(r9303);
        double r9305 = atan(r9301);
        double r9306 = atan2(r9304, r9305);
        double r9307 = fmod(r9306, r9301);
        double r9308 = fabs(r9307);
        return r9308;
}

↓

double f(double a) {
        double r9309 = a;
        double r9310 = expm1(r9309);
        double r9311 = sin(r9310);
        double r9312 = cbrt(r9311);
        double r9313 = r9312 * r9312;
        double r9314 = r9313 * r9312;
        double r9315 = expm1(r9314);
        double r9316 = atan(r9309);
        double r9317 = atan2(r9315, r9316);
        double r9318 = 3.0;
        double r9319 = pow(r9317, r9318);
        double r9320 = cbrt(r9319);
        double r9321 = fmod(r9320, r9309);
        double r9322 = cbrt(r9321);
        double r9323 = r9322 * r9322;
        double r9324 = expm1(r9311);
        double r9325 = atan2(r9324, r9316);
        double r9326 = fmod(r9325, r9309);
        double r9327 = cbrt(r9326);
        double r9328 = r9323 * r9327;
        double r9329 = fabs(r9328);
        return r9329;
}

Error

Bits error versus a

Derivation

1. Initial program 33.8

   \[\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-cube 33.9

   \[\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. Simplified 33.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-cbrt 33.9

   \[\leadsto \left|\left(\left(\sqrt[3]{{\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)}^{3}}\right) \bmod a\right)\right|\]
7. Using strategy rm

8. Applied add-cube-cbrt 33.9

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

9. Taylor expanded around 0 33.8

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

10. Final simplification 33.8

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

Reproduce

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