\left(\left(x \cdot \log y - y\right) - z\right) + \log t
\left(\left(-1 \cdot \left(\log \left(\frac{1}{y}\right) \cdot x\right) - y\right) - z\right) + \log tdouble f(double x, double y, double z, double t) {
double r119884 = x;
double r119885 = y;
double r119886 = log(r119885);
double r119887 = r119884 * r119886;
double r119888 = r119887 - r119885;
double r119889 = z;
double r119890 = r119888 - r119889;
double r119891 = t;
double r119892 = log(r119891);
double r119893 = r119890 + r119892;
return r119893;
}
double f(double x, double y, double z, double t) {
double r119894 = -1.0;
double r119895 = 1.0;
double r119896 = y;
double r119897 = r119895 / r119896;
double r119898 = log(r119897);
double r119899 = x;
double r119900 = r119898 * r119899;
double r119901 = r119894 * r119900;
double r119902 = r119901 - r119896;
double r119903 = z;
double r119904 = r119902 - r119903;
double r119905 = t;
double r119906 = log(r119905);
double r119907 = r119904 + r119906;
return r119907;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t
Results
Initial program 0.1
rmApplied add-cube-cbrt0.1
Applied log-prod0.1
Applied distribute-lft-in0.1
Simplified0.1
rmApplied fma-def0.1
Taylor expanded around inf 0.2
Simplified0.1
Final simplification0.1
herbie shell --seed 2020046 +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)))