Random Jason Timeout Test 015

Average Error: 0.1 → 0.1

Time: 14.3s

Precision: 64

• Falling back on racket egraph because egg-herbie package not installed

• could not determine a ground truth for program body

  1. a = 2.2108085494553437e+120
  2. b = 7.734281002290599e-114

Math

\[\sin \left({\left(\sqrt{\tan^{-1}_* \frac{b}{b}}\right)}^{\left(b - a\right)}\right)\]

↓

\[\sin \left({\left(\sqrt{\sqrt{\tan^{-1}_* \frac{b}{b}}}\right)}^{\left(b - a\right)} \cdot \mathsf{expm1}\left(\mathsf{log1p}\left({\left(\sqrt{\sqrt{\tan^{-1}_* \frac{b}{b}}}\right)}^{\left(b - a\right)}\right)\right)\right)\]

C

double f(double a, double b) {
        double r11033 = b;
        double r11034 = atan2(r11033, r11033);
        double r11035 = sqrt(r11034);
        double r11036 = a;
        double r11037 = r11033 - r11036;
        double r11038 = pow(r11035, r11037);
        double r11039 = sin(r11038);
        return r11039;
}

↓

double f(double a, double b) {
        double r11040 = b;
        double r11041 = atan2(r11040, r11040);
        double r11042 = sqrt(r11041);
        double r11043 = sqrt(r11042);
        double r11044 = a;
        double r11045 = r11040 - r11044;
        double r11046 = pow(r11043, r11045);
        double r11047 = log1p(r11046);
        double r11048 = expm1(r11047);
        double r11049 = r11046 * r11048;
        double r11050 = sin(r11049);
        return r11050;
}

Error

Bits error versus a

Bits error versus b

Derivation

1. Initial program 0.1

   \[\sin \left({\left(\sqrt{\tan^{-1}_* \frac{b}{b}}\right)}^{\left(b - a\right)}\right)\]
2. Using strategy rm

3. Applied add-sqr-sqrt 0.1

   \[\leadsto \sin \left({\left(\sqrt{\color{blue}{\sqrt{\tan^{-1}_* \frac{b}{b}} \cdot \sqrt{\tan^{-1}_* \frac{b}{b}}}}\right)}^{\left(b - a\right)}\right)\]

4. Applied sqrt-prod 0.1

   \[\leadsto \sin \left({\color{blue}{\left(\sqrt{\sqrt{\tan^{-1}_* \frac{b}{b}}} \cdot \sqrt{\sqrt{\tan^{-1}_* \frac{b}{b}}}\right)}}^{\left(b - a\right)}\right)\]

5. Applied unpow-prod-down 0.1

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

7. Applied expm1-log1p-u 0.1

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

8. Final simplification 0.1

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

Reproduce

herbie shell --seed 2020047 +o rules:numerics
(FPCore (a b)
  :name "Random Jason Timeout Test 015"
  :precision binary64
  (sin (pow (sqrt (atan2 b b)) (- b a))))