Random Jason Timeout Test 006

Average Error: 33.3 → 33.3

Time: 21.3s

Precision: 64

• could not determine a ground truth for program body

  1. a = 3.1717798053764645e+40

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

C

double f(double a) {
        double r24192 = a;
        double r24193 = expm1(r24192);
        double r24194 = sin(r24193);
        double r24195 = expm1(r24194);
        double r24196 = atan(r24192);
        double r24197 = atan2(r24195, r24196);
        double r24198 = fmod(r24197, r24192);
        double r24199 = fabs(r24198);
        return r24199;
}

↓

double f(double a) {
        double r24200 = a;
        double r24201 = expm1(r24200);
        double r24202 = sin(r24201);
        double r24203 = cbrt(r24202);
        double r24204 = cbrt(r24203);
        double r24205 = r24204 * r24204;
        double r24206 = r24205 * r24204;
        double r24207 = r24203 * r24206;
        double r24208 = cbrt(r24207);
        double r24209 = r24208 * r24208;
        double r24210 = r24203 * r24203;
        double r24211 = cbrt(r24210);
        double r24212 = log(r24211);
        double r24213 = exp(r24212);
        double r24214 = r24209 * r24213;
        double r24215 = r24214 * r24203;
        double r24216 = expm1(r24215);
        double r24217 = atan(r24200);
        double r24218 = atan2(r24216, r24217);
        double r24219 = fmod(r24218, r24200);
        double r24220 = fabs(r24219);
        return r24220;
}

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-cbrt 33.4

   \[\leadsto \left|\left(\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) \bmod a\right)\right|\]
4. Using strategy rm

5. Applied add-cube-cbrt 33.4

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

7. Applied add-cube-cbrt 33.4

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

8. Simplified 33.4

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

10. Applied add-exp-log 33.3

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

11. Final simplification 33.3

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

Reproduce

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