Random Jason Timeout Test 006

Average Error: 33.5 → 33.6

Time: 47.3s

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 r32113 = a;
        double r32114 = expm1(r32113);
        double r32115 = sin(r32114);
        double r32116 = expm1(r32115);
        double r32117 = atan(r32113);
        double r32118 = atan2(r32116, r32117);
        double r32119 = fmod(r32118, r32113);
        double r32120 = fabs(r32119);
        return r32120;
}

↓

double f(double a) {
        double r32121 = a;
        double r32122 = expm1(r32121);
        double r32123 = sin(r32122);
        double r32124 = expm1(r32123);
        double r32125 = cbrt(r32124);
        double r32126 = r32125 * r32125;
        double r32127 = cbrt(r32125);
        double r32128 = r32127 * r32127;
        double r32129 = r32128 * r32127;
        double r32130 = r32126 * r32129;
        double r32131 = atan(r32121);
        double r32132 = atan2(r32130, r32131);
        double r32133 = fmod(r32132, r32121);
        double r32134 = fabs(r32133);
        return r32134;
}

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 +o rules:numerics
(FPCore (a)
  :name "Random Jason Timeout Test 006"
  :precision binary64
  (fabs (fmod (atan2 (expm1 (sin (expm1 a))) (atan a)) a)))