\frac{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.764470499998\right) \cdot y + 230661.510616000014\right) \cdot y + t}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}\frac{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.764470499998\right) \cdot y + 230661.510616000014\right) \cdot y + t}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r59780 = x;
double r59781 = y;
double r59782 = r59780 * r59781;
double r59783 = z;
double r59784 = r59782 + r59783;
double r59785 = r59784 * r59781;
double r59786 = 27464.7644705;
double r59787 = r59785 + r59786;
double r59788 = r59787 * r59781;
double r59789 = 230661.510616;
double r59790 = r59788 + r59789;
double r59791 = r59790 * r59781;
double r59792 = t;
double r59793 = r59791 + r59792;
double r59794 = a;
double r59795 = r59781 + r59794;
double r59796 = r59795 * r59781;
double r59797 = b;
double r59798 = r59796 + r59797;
double r59799 = r59798 * r59781;
double r59800 = c;
double r59801 = r59799 + r59800;
double r59802 = r59801 * r59781;
double r59803 = i;
double r59804 = r59802 + r59803;
double r59805 = r59793 / r59804;
return r59805;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r59806 = x;
double r59807 = y;
double r59808 = r59806 * r59807;
double r59809 = z;
double r59810 = r59808 + r59809;
double r59811 = r59810 * r59807;
double r59812 = 27464.7644705;
double r59813 = r59811 + r59812;
double r59814 = r59813 * r59807;
double r59815 = 230661.510616;
double r59816 = r59814 + r59815;
double r59817 = r59816 * r59807;
double r59818 = t;
double r59819 = r59817 + r59818;
double r59820 = a;
double r59821 = r59807 + r59820;
double r59822 = r59821 * r59807;
double r59823 = b;
double r59824 = r59822 + r59823;
double r59825 = r59824 * r59807;
double r59826 = c;
double r59827 = r59825 + r59826;
double r59828 = r59827 * r59807;
double r59829 = i;
double r59830 = r59828 + r59829;
double r59831 = r59819 / r59830;
return r59831;
}



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 28.7
Final simplification28.7
herbie shell --seed 2020027
(FPCore (x y z t a b c i)
:name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2"
:precision binary64
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))