\left(\left(x \cdot \log y - y\right) - z\right) + \log t
\left(\left(x \cdot \log y - y\right) - z\right) + \log t
double code(double x, double y, double z, double t) {
return ((((x * log(y)) - y) - z) + log(t));
}
double code(double x, double y, double z, double t) {
return ((((x * log(y)) - y) - z) + log(t));
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t
Results
Initial program 0.1
Final simplification0.1
herbie shell --seed 2020091
(FPCore (x y z t)
:name "Numeric.SpecFunctions:incompleteGamma from math-functions-0.1.5.2, A"
:precision binary64
(+ (- (- (* x (log y)) y) z) (log t)))