\left(x \cdot y - z \cdot y\right) \cdot t
\left(t \cdot y\right) \cdot \left(x - z\right)
double f(double x, double y, double z, double t) {
double r398314 = x;
double r398315 = y;
double r398316 = r398314 * r398315;
double r398317 = z;
double r398318 = r398317 * r398315;
double r398319 = r398316 - r398318;
double r398320 = t;
double r398321 = r398319 * r398320;
return r398321;
}
double f(double x, double y, double z, double t) {
double r398322 = t;
double r398323 = y;
double r398324 = r398322 * r398323;
double r398325 = x;
double r398326 = z;
double r398327 = r398325 - r398326;
double r398328 = r398324 * r398327;
return r398328;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.2 |
|---|---|
| Target | 3.1 |
| Herbie | 7.1 |
if (- (* x y) (* z y)) < -1.7752227888120415e+293 or 3.1698261024189057e+198 < (- (* x y) (* z y)) Initial program 38.2
rmApplied add-cube-cbrt38.6
Applied associate-*l*38.6
Taylor expanded around inf 38.2
Simplified0.7
if -1.7752227888120415e+293 < (- (* x y) (* z y)) < -5.7409949950736745e-146 or 6.1708799165572e-321 < (- (* x y) (* z y)) < 3.1698261024189057e+198Initial program 0.3
if -5.7409949950736745e-146 < (- (* x y) (* z y)) < 6.1708799165572e-321Initial program 9.4
rmApplied distribute-rgt-out--9.4
Applied associate-*l*1.8
Final simplification7.1
herbie shell --seed 2019308
(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))