\frac{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.7644705\right) \cdot y + 230661.510616\right) \cdot y + t}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}\left(t + \left(y \cdot \left(y \cdot \left(z + x \cdot y\right) + 27464.7644705\right) + 230661.510616\right) \cdot y\right) \cdot \frac{1}{i + y \cdot \left(\left(\left(a + y\right) \cdot y + b\right) \cdot y + c\right)}double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r3054097 = x;
double r3054098 = y;
double r3054099 = r3054097 * r3054098;
double r3054100 = z;
double r3054101 = r3054099 + r3054100;
double r3054102 = r3054101 * r3054098;
double r3054103 = 27464.7644705;
double r3054104 = r3054102 + r3054103;
double r3054105 = r3054104 * r3054098;
double r3054106 = 230661.510616;
double r3054107 = r3054105 + r3054106;
double r3054108 = r3054107 * r3054098;
double r3054109 = t;
double r3054110 = r3054108 + r3054109;
double r3054111 = a;
double r3054112 = r3054098 + r3054111;
double r3054113 = r3054112 * r3054098;
double r3054114 = b;
double r3054115 = r3054113 + r3054114;
double r3054116 = r3054115 * r3054098;
double r3054117 = c;
double r3054118 = r3054116 + r3054117;
double r3054119 = r3054118 * r3054098;
double r3054120 = i;
double r3054121 = r3054119 + r3054120;
double r3054122 = r3054110 / r3054121;
return r3054122;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i) {
double r3054123 = t;
double r3054124 = y;
double r3054125 = z;
double r3054126 = x;
double r3054127 = r3054126 * r3054124;
double r3054128 = r3054125 + r3054127;
double r3054129 = r3054124 * r3054128;
double r3054130 = 27464.7644705;
double r3054131 = r3054129 + r3054130;
double r3054132 = r3054124 * r3054131;
double r3054133 = 230661.510616;
double r3054134 = r3054132 + r3054133;
double r3054135 = r3054134 * r3054124;
double r3054136 = r3054123 + r3054135;
double r3054137 = 1.0;
double r3054138 = i;
double r3054139 = a;
double r3054140 = r3054139 + r3054124;
double r3054141 = r3054140 * r3054124;
double r3054142 = b;
double r3054143 = r3054141 + r3054142;
double r3054144 = r3054143 * r3054124;
double r3054145 = c;
double r3054146 = r3054144 + r3054145;
double r3054147 = r3054124 * r3054146;
double r3054148 = r3054138 + r3054147;
double r3054149 = r3054137 / r3054148;
double r3054150 = r3054136 * r3054149;
return r3054150;
}



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