\left(\left(x \cdot \log y - y\right) - z\right) + \log t
\left(\left(\left(\left(2 \cdot \log \left(\sqrt[3]{y}\right)\right) \cdot x + x \cdot \log \left({\left(\frac{1}{y}\right)}^{\frac{-1}{3}}\right)\right) - y\right) - z\right) + \log tdouble f(double x, double y, double z, double t) {
double r83673 = x;
double r83674 = y;
double r83675 = log(r83674);
double r83676 = r83673 * r83675;
double r83677 = r83676 - r83674;
double r83678 = z;
double r83679 = r83677 - r83678;
double r83680 = t;
double r83681 = log(r83680);
double r83682 = r83679 + r83681;
return r83682;
}
double f(double x, double y, double z, double t) {
double r83683 = 2.0;
double r83684 = y;
double r83685 = cbrt(r83684);
double r83686 = log(r83685);
double r83687 = r83683 * r83686;
double r83688 = x;
double r83689 = r83687 * r83688;
double r83690 = 1.0;
double r83691 = r83690 / r83684;
double r83692 = -0.3333333333333333;
double r83693 = pow(r83691, r83692);
double r83694 = log(r83693);
double r83695 = r83688 * r83694;
double r83696 = r83689 + r83695;
double r83697 = r83696 - r83684;
double r83698 = z;
double r83699 = r83697 - r83698;
double r83700 = t;
double r83701 = log(r83700);
double r83702 = r83699 + r83701;
return r83702;
}



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
Taylor expanded around inf 0.1
Final simplification0.1
herbie shell --seed 2020047
(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)))