\frac{x \cdot y - \left(z \cdot 9\right) \cdot t}{a \cdot 2}\frac{x \cdot y - 9 \cdot \left(t \cdot z\right)}{a \cdot 2}double f(double x, double y, double z, double t, double a) {
double r673188 = x;
double r673189 = y;
double r673190 = r673188 * r673189;
double r673191 = z;
double r673192 = 9.0;
double r673193 = r673191 * r673192;
double r673194 = t;
double r673195 = r673193 * r673194;
double r673196 = r673190 - r673195;
double r673197 = a;
double r673198 = 2.0;
double r673199 = r673197 * r673198;
double r673200 = r673196 / r673199;
return r673200;
}
double f(double x, double y, double z, double t, double a) {
double r673201 = x;
double r673202 = y;
double r673203 = r673201 * r673202;
double r673204 = 9.0;
double r673205 = t;
double r673206 = z;
double r673207 = r673205 * r673206;
double r673208 = r673204 * r673207;
double r673209 = r673203 - r673208;
double r673210 = a;
double r673211 = 2.0;
double r673212 = r673210 * r673211;
double r673213 = r673209 / r673212;
return r673213;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 7.7 |
|---|---|
| Target | 5.6 |
| Herbie | 7.6 |
Initial program 7.7
Taylor expanded around inf 7.6
Final simplification7.6
herbie shell --seed 2020060 +o rules:numerics
(FPCore (x y z t a)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, I"
:precision binary64
:herbie-target
(if (< a -2.090464557976709e+86) (- (* 0.5 (/ (* y x) a)) (* 4.5 (/ t (/ a z)))) (if (< a 2.144030707833976e+99) (/ (- (* x y) (* z (* 9 t))) (* a 2)) (- (* (/ y a) (* x 0.5)) (* (/ t a) (* z 4.5)))))
(/ (- (* x y) (* (* z 9) t)) (* a 2)))