\left(\left(\left(\left(x \cdot \log y + z\right) + t\right) + a\right) + \left(b - 0.5\right) \cdot \log c\right) + y \cdot i
\left(\left(\left(\left(\left(x \cdot \left(2 \cdot \log \left(\sqrt[3]{y}\right)\right) + \log \left(\mathsf{expm1}\left(\mathsf{log1p}\left(\sqrt[3]{y}\right)\right)\right) \cdot x\right) + z\right) + t\right) + a\right) + \left(b - 0.5\right) \cdot \log c\right) + y \cdot idouble f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r82800 = x;
double r82801 = y;
double r82802 = log(r82801);
double r82803 = r82800 * r82802;
double r82804 = z;
double r82805 = r82803 + r82804;
double r82806 = t;
double r82807 = r82805 + r82806;
double r82808 = a;
double r82809 = r82807 + r82808;
double r82810 = b;
double r82811 = 0.5;
double r82812 = r82810 - r82811;
double r82813 = c;
double r82814 = log(r82813);
double r82815 = r82812 * r82814;
double r82816 = r82809 + r82815;
double r82817 = i;
double r82818 = r82801 * r82817;
double r82819 = r82816 + r82818;
return r82819;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r82820 = x;
double r82821 = 2.0;
double r82822 = y;
double r82823 = cbrt(r82822);
double r82824 = log(r82823);
double r82825 = r82821 * r82824;
double r82826 = r82820 * r82825;
double r82827 = log1p(r82823);
double r82828 = expm1(r82827);
double r82829 = log(r82828);
double r82830 = r82829 * r82820;
double r82831 = r82826 + r82830;
double r82832 = z;
double r82833 = r82831 + r82832;
double r82834 = t;
double r82835 = r82833 + r82834;
double r82836 = a;
double r82837 = r82835 + r82836;
double r82838 = b;
double r82839 = 0.5;
double r82840 = r82838 - r82839;
double r82841 = c;
double r82842 = log(r82841);
double r82843 = r82840 * r82842;
double r82844 = r82837 + r82843;
double r82845 = i;
double r82846 = r82822 * r82845;
double r82847 = r82844 + r82846;
return r82847;
}



Bits error versus x



Bits error versus y



Bits error versus z



Bits error versus t



Bits error versus a



Bits error versus b



Bits error versus c



Bits error versus i
Results
Initial program 0.1
rmApplied add-cube-cbrt0.1
Applied log-prod0.1
Applied distribute-lft-in0.1
Simplified0.1
Simplified0.1
rmApplied expm1-log1p-u0.1
Final simplification0.1
herbie shell --seed 2019323 +o rules:numerics
(FPCore (x y z t a b c i)
:name "Numeric.SpecFunctions:logBeta from math-functions-0.1.5.2, B"
:precision binary64
(+ (+ (+ (+ (+ (* x (log y)) z) t) a) (* (- b 0.5) (log c))) (* y i)))