\left(\left(x \cdot \log y - y\right) - z\right) + \log t
\mathsf{fma}\left(\log y, x, \log \left(\sqrt{t}\right) + \left(\log \left(\sqrt{t}\right) - \left(y + z\right)\right)\right)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 fma(log(y), x, (log(sqrt(t)) + (log(sqrt(t)) - (y + z))));
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t
Results
Initial program 0.1
Simplified0.1
rmApplied add-sqr-sqrt0.1
Applied log-prod0.1
Applied associate--l+0.1
Final simplification0.1
herbie shell --seed 2020058 +o rules:numerics
(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)))