\left(x \cdot y - z \cdot y\right) \cdot t
t \cdot \left(y \cdot \left(x - z\right)\right)
double f(double x, double y, double z, double t) {
double r452478 = x;
double r452479 = y;
double r452480 = r452478 * r452479;
double r452481 = z;
double r452482 = r452481 * r452479;
double r452483 = r452480 - r452482;
double r452484 = t;
double r452485 = r452483 * r452484;
return r452485;
}
double f(double x, double y, double z, double t) {
double r452486 = t;
double r452487 = y;
double r452488 = x;
double r452489 = z;
double r452490 = r452488 - r452489;
double r452491 = r452487 * r452490;
double r452492 = r452486 * r452491;
return r452492;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.0 |
|---|---|
| Target | 3.2 |
| Herbie | 7.0 |
if y < -368485883737218.6 or 2.1087438558749025e-26 < y Initial program 14.7
Simplified14.7
rmApplied associate-*r*3.8
if -368485883737218.6 < y < 2.1087438558749025e-26Initial program 2.0
Simplified2.0
rmApplied associate-*r*8.8
rmApplied associate-*l*2.0
Final simplification7.0
herbie shell --seed 2019297
(FPCore (x y z t)
:name "Linear.Projection:inverseInfinitePerspective from linear-1.19.1.3"
:precision binary64
:herbie-target
(if (< t -9.2318795828867769e-80) (* (* y t) (- x z)) (if (< t 2.5430670515648771e83) (* y (* t (- x z))) (* (* y (- x z)) t)))
(* (- (* x y) (* z y)) t))