\frac{x \cdot y - \left(z \cdot 9\right) \cdot t}{a \cdot 2}0.5 \cdot \frac{x}{\frac{a}{y}} - \frac{4.5 \cdot \left(t \cdot z\right)}{a}double f(double x, double y, double z, double t, double a) {
double r437652 = x;
double r437653 = y;
double r437654 = r437652 * r437653;
double r437655 = z;
double r437656 = 9.0;
double r437657 = r437655 * r437656;
double r437658 = t;
double r437659 = r437657 * r437658;
double r437660 = r437654 - r437659;
double r437661 = a;
double r437662 = 2.0;
double r437663 = r437661 * r437662;
double r437664 = r437660 / r437663;
return r437664;
}
double f(double x, double y, double z, double t, double a) {
double r437665 = 0.5;
double r437666 = x;
double r437667 = a;
double r437668 = y;
double r437669 = r437667 / r437668;
double r437670 = r437666 / r437669;
double r437671 = r437665 * r437670;
double r437672 = 4.5;
double r437673 = t;
double r437674 = z;
double r437675 = r437673 * r437674;
double r437676 = r437672 * r437675;
double r437677 = r437676 / r437667;
double r437678 = r437671 - r437677;
return r437678;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 7.3 |
|---|---|
| Target | 5.5 |
| Herbie | 8.0 |
if (* x y) < -inf.0Initial program 64.0
Taylor expanded around 0 64.0
rmApplied associate-*r/64.0
rmApplied *-un-lft-identity64.0
Applied times-frac5.8
Simplified5.8
if -inf.0 < (* x y) < -1.8798820468816307e+17Initial program 4.3
Taylor expanded around 0 4.3
rmApplied associate-/l*1.9
if -1.8798820468816307e+17 < (* x y) < 2.7107991116943106e+188Initial program 3.8
Taylor expanded around 0 3.7
rmApplied div-inv3.7
if 2.7107991116943106e+188 < (* x y) Initial program 27.9
Taylor expanded around 0 27.7
rmApplied associate-*r/27.7
rmApplied associate-/l*6.6
Final simplification8.0
herbie shell --seed 2019294
(FPCore (x y z t a)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, I"
:precision binary64
:herbie-target
(if (< a -2.090464557976709e86) (- (* 0.5 (/ (* y x) a)) (* 4.5 (/ t (/ a z)))) (if (< a 2.14403070783397609e99) (/ (- (* x y) (* z (* 9 t))) (* a 2)) (- (* (/ y a) (* x 0.5)) (* (/ t a) (* z 4.5)))))
(/ (- (* x y) (* (* z 9) t)) (* a 2)))