Random Jason Timeout Test 006

Average Error: 33.7 → 33.7

Time: 17.3s

Precision: 64

• could not determine a ground truth for program body

  1. a = 9.924781450718653e+287

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(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right)}^{3}}\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(\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|\]

C

double f(double a) {
        double r18379 = a;
        double r18380 = expm1(r18379);
        double r18381 = sin(r18380);
        double r18382 = expm1(r18381);
        double r18383 = atan(r18379);
        double r18384 = atan2(r18382, r18383);
        double r18385 = fmod(r18384, r18379);
        double r18386 = fabs(r18385);
        return r18386;
}

↓

double f(double a) {
        double r18387 = a;
        double r18388 = expm1(r18387);
        double r18389 = sin(r18388);
        double r18390 = expm1(r18389);
        double r18391 = atan(r18387);
        double r18392 = atan2(r18390, r18391);
        double r18393 = 3.0;
        double r18394 = pow(r18392, r18393);
        double r18395 = cbrt(r18394);
        double r18396 = fmod(r18395, r18387);
        double r18397 = cbrt(r18396);
        double r18398 = fmod(r18392, r18387);
        double r18399 = cbrt(r18398);
        double r18400 = r18397 * r18399;
        double r18401 = r18392 * r18392;
        double r18402 = r18401 * r18392;
        double r18403 = cbrt(r18402);
        double r18404 = fmod(r18403, r18387);
        double r18405 = cbrt(r18404);
        double r18406 = r18400 * r18405;
        double r18407 = fabs(r18406);
        return r18407;
}

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-cube 33.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. Simplified 33.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-cbrt 33.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. Using strategy rm

8. Applied add-cbrt-cube 33.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}{\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)}}^{3}}\right) \bmod a\right)}\right|\]

9. Applied rem-cube-cbrt 33.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}{\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|\]

10. Taylor expanded around 0 33.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]{\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) \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} \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|\]

11. Final simplification 33.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(\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(\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|\]

Reproduce

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