\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 r75075 = x;
double r75076 = y;
double r75077 = r75075 * r75076;
double r75078 = z;
double r75079 = r75077 + r75078;
double r75080 = r75079 * r75076;
double r75081 = 27464.7644705;
double r75082 = r75080 + r75081;
double r75083 = r75082 * r75076;
double r75084 = 230661.510616;
double r75085 = r75083 + r75084;
double r75086 = r75085 * r75076;
double r75087 = t;
double r75088 = r75086 + r75087;
double r75089 = a;
double r75090 = r75076 + r75089;
double r75091 = r75090 * r75076;
double r75092 = b;
double r75093 = r75091 + r75092;
double r75094 = r75093 * r75076;
double r75095 = c;
double r75096 = r75094 + r75095;
double r75097 = r75096 * r75076;
double r75098 = i;
double r75099 = r75097 + r75098;
double r75100 = r75088 / r75099;
return r75100;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r75101 = x;
double r75102 = y;
double r75103 = r75101 * r75102;
double r75104 = z;
double r75105 = r75103 + r75104;
double r75106 = r75105 * r75102;
double r75107 = 27464.7644705;
double r75108 = r75106 + r75107;
double r75109 = r75108 * r75102;
double r75110 = 230661.510616;
double r75111 = r75109 + r75110;
double r75112 = r75111 * r75102;
double r75113 = t;
double r75114 = r75112 + r75113;
double r75115 = a;
double r75116 = r75102 + r75115;
double r75117 = r75116 * r75102;
double r75118 = b;
double r75119 = r75117 + r75118;
double r75120 = r75119 * r75102;
double r75121 = c;
double r75122 = r75120 + r75121;
double r75123 = r75122 * r75102;
double r75124 = i;
double r75125 = r75123 + r75124;
double r75126 = r75114 / r75125;
return r75126;
}



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 29.5
Final simplification29.5
herbie shell --seed 2020046 +o rules:numerics
(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)))