2-ancestry mixing, zero discriminant

Average Error: 15.2 → 0.8

Time: 18.1s

Precision: 64

Math

\[\sqrt[3]{\frac{g}{2 \cdot a}}\]

↓

\[\frac{\sqrt[3]{g}}{\sqrt[3]{2 \cdot a}}\]

C

double f(double g, double a) {
        double r23376501 = g;
        double r23376502 = 2.0;
        double r23376503 = a;
        double r23376504 = r23376502 * r23376503;
        double r23376505 = r23376501 / r23376504;
        double r23376506 = cbrt(r23376505);
        return r23376506;
}

↓

double f(double g, double a) {
        double r23376507 = g;
        double r23376508 = cbrt(r23376507);
        double r23376509 = 2.0;
        double r23376510 = a;
        double r23376511 = r23376509 * r23376510;
        double r23376512 = cbrt(r23376511);
        double r23376513 = r23376508 / r23376512;
        return r23376513;
}

Error

Bits error versus g

Bits error versus a

Derivation

1. Initial program 15.2

   \[\sqrt[3]{\frac{g}{2 \cdot a}}\]
2. Using strategy rm

3. Applied cbrt-div 0.8

   \[\leadsto \color{blue}{\frac{\sqrt[3]{g}}{\sqrt[3]{2 \cdot a}}}\]

4. Final simplification 0.8

   \[\leadsto \frac{\sqrt[3]{g}}{\sqrt[3]{2 \cdot a}}\]

Reproduce

herbie shell --seed 2019128 
(FPCore (g a)
  :name "2-ancestry mixing, zero discriminant"
  (cbrt (/ g (* 2 a))))