Average Error: 0.0 → 0.2
Time: 13.8s
Precision: 64
\[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b\]
↓
\[\left(\left(x - \left(\sqrt[3]{y - 1} \cdot \sqrt[3]{y - 1}\right) \cdot \left(\sqrt[3]{y - 1} \cdot z\right)\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b\]
\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b
↓
\left(\left(x - \left(\sqrt[3]{y - 1} \cdot \sqrt[3]{y - 1}\right) \cdot \left(\sqrt[3]{y - 1} \cdot z\right)\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b
double f(double x, double y, double z, double t, double a, double b) {
        double r32998 = x;
        double r32999 = y;
        double r33000 = 1.0;
        double r33001 = r32999 - r33000;
        double r33002 = z;
        double r33003 = r33001 * r33002;
        double r33004 = r32998 - r33003;
        double r33005 = t;
        double r33006 = r33005 - r33000;
        double r33007 = a;
        double r33008 = r33006 * r33007;
        double r33009 = r33004 - r33008;
        double r33010 = r32999 + r33005;
        double r33011 = 2.0;
        double r33012 = r33010 - r33011;
        double r33013 = b;
        double r33014 = r33012 * r33013;
        double r33015 = r33009 + r33014;
        return r33015;
}
↓
double f(double x, double y, double z, double t, double a, double b) {
        double r33016 = x;
        double r33017 = y;
        double r33018 = 1.0;
        double r33019 = r33017 - r33018;
        double r33020 = cbrt(r33019);
        double r33021 = r33020 * r33020;
        double r33022 = z;
        double r33023 = r33020 * r33022;
        double r33024 = r33021 * r33023;
        double r33025 = r33016 - r33024;
        double r33026 = t;
        double r33027 = r33026 - r33018;
        double r33028 = a;
        double r33029 = r33027 * r33028;
        double r33030 = r33025 - r33029;
        double r33031 = r33017 + r33026;
        double r33032 = 2.0;
        double r33033 = r33031 - r33032;
        double r33034 = b;
        double r33035 = r33033 * r33034;
        double r33036 = r33030 + r33035;
        return r33036;
}

Error

Bits error versus x

Bits error versus y

Bits error versus z

Bits error versus t

Bits error versus a

Bits error versus b

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Initial program 0.0

    \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b\]
  2. Using strategy rm

  3. Applied add-cube-cbrt0.2

    \[\leadsto \left(\left(x - \color{blue}{\left(\left(\sqrt[3]{y - 1} \cdot \sqrt[3]{y - 1}\right) \cdot \sqrt[3]{y - 1}\right)} \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b\]

  4. Applied associate-*l*0.2

    \[\leadsto \left(\left(x - \color{blue}{\left(\sqrt[3]{y - 1} \cdot \sqrt[3]{y - 1}\right) \cdot \left(\sqrt[3]{y - 1} \cdot z\right)}\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b\]

  5. Final simplification0.2

    \[\leadsto \left(\left(x - \left(\sqrt[3]{y - 1} \cdot \sqrt[3]{y - 1}\right) \cdot \left(\sqrt[3]{y - 1} \cdot z\right)\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b\]

Reproduce

herbie shell --seed 2019305 +o rules:numerics
(FPCore (x y z t a b)
  :name "Statistics.Distribution.Beta:$centropy from math-functions-0.1.5.2"
  :precision binary64
  (+ (- (- x (* (- y 1) z)) (* (- t 1) a)) (* (- (+ y t) 2) b)))