Random Jason Timeout Test 006

Average Error: 33.5 → 33.6

Time: 40.9s

Precision: 64

• could not determine a ground truth for program body

  1. a = 7.953652505345186e+78

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

C

double f(double a) {
        double r24496 = a;
        double r24497 = expm1(r24496);
        double r24498 = sin(r24497);
        double r24499 = expm1(r24498);
        double r24500 = atan(r24496);
        double r24501 = atan2(r24499, r24500);
        double r24502 = fmod(r24501, r24496);
        double r24503 = fabs(r24502);
        return r24503;
}

↓

double f(double a) {
        double r24504 = a;
        double r24505 = expm1(r24504);
        double r24506 = sin(r24505);
        double r24507 = expm1(r24506);
        double r24508 = cbrt(r24507);
        double r24509 = r24508 * r24508;
        double r24510 = cbrt(r24508);
        double r24511 = r24510 * r24510;
        double r24512 = r24511 * r24510;
        double r24513 = r24509 * r24512;
        double r24514 = atan(r24504);
        double r24515 = atan2(r24513, r24514);
        double r24516 = fmod(r24515, r24504);
        double r24517 = fabs(r24516);
        return r24517;
}

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-cube-cbrt 33.5

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

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

6. Final simplification 33.6

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

Reproduce

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