\left(x \cdot 2 - \left(\left(y \cdot 9\right) \cdot z\right) \cdot t\right) + \left(a \cdot 27\right) \cdot b
\begin{array}{l}
\mathbf{if}\;t \le -2.171737053366690492496454059018873464008 \cdot 10^{128}:\\
\;\;\;\;\left(2 \cdot x - 9 \cdot \left(t \cdot \left(z \cdot y\right)\right)\right) + \left(a \cdot 27\right) \cdot b\\
\mathbf{elif}\;t \le 4.404072813174863545657687298340743366284 \cdot 10^{58}:\\
\;\;\;\;a \cdot \left(27 \cdot b\right) + \left(2 \cdot x - y \cdot \left(\left(9 \cdot z\right) \cdot t\right)\right)\\
\mathbf{else}:\\
\;\;\;\;27 \cdot \left(a \cdot b\right) + \left(x \cdot 2 - \left(y \cdot \left(z \cdot 9\right)\right) \cdot t\right)\\
\end{array}double f(double x, double y, double z, double t, double a, double b) {
double r558133 = x;
double r558134 = 2.0;
double r558135 = r558133 * r558134;
double r558136 = y;
double r558137 = 9.0;
double r558138 = r558136 * r558137;
double r558139 = z;
double r558140 = r558138 * r558139;
double r558141 = t;
double r558142 = r558140 * r558141;
double r558143 = r558135 - r558142;
double r558144 = a;
double r558145 = 27.0;
double r558146 = r558144 * r558145;
double r558147 = b;
double r558148 = r558146 * r558147;
double r558149 = r558143 + r558148;
return r558149;
}
double f(double x, double y, double z, double t, double a, double b) {
double r558150 = t;
double r558151 = -2.1717370533666905e+128;
bool r558152 = r558150 <= r558151;
double r558153 = 2.0;
double r558154 = x;
double r558155 = r558153 * r558154;
double r558156 = 9.0;
double r558157 = z;
double r558158 = y;
double r558159 = r558157 * r558158;
double r558160 = r558150 * r558159;
double r558161 = r558156 * r558160;
double r558162 = r558155 - r558161;
double r558163 = a;
double r558164 = 27.0;
double r558165 = r558163 * r558164;
double r558166 = b;
double r558167 = r558165 * r558166;
double r558168 = r558162 + r558167;
double r558169 = 4.4040728131748635e+58;
bool r558170 = r558150 <= r558169;
double r558171 = r558164 * r558166;
double r558172 = r558163 * r558171;
double r558173 = r558156 * r558157;
double r558174 = r558173 * r558150;
double r558175 = r558158 * r558174;
double r558176 = r558155 - r558175;
double r558177 = r558172 + r558176;
double r558178 = r558163 * r558166;
double r558179 = r558164 * r558178;
double r558180 = r558154 * r558153;
double r558181 = r558157 * r558156;
double r558182 = r558158 * r558181;
double r558183 = r558182 * r558150;
double r558184 = r558180 - r558183;
double r558185 = r558179 + r558184;
double r558186 = r558170 ? r558177 : r558185;
double r558187 = r558152 ? r558168 : r558186;
return r558187;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a




Bits error versus b
Results
| Original | 3.5 |
|---|---|
| Target | 2.4 |
| Herbie | 1.0 |
if t < -2.1717370533666905e+128Initial program 1.4
Taylor expanded around inf 1.4
if -2.1717370533666905e+128 < t < 4.4040728131748635e+58Initial program 4.5
rmApplied associate-*l*4.5
rmApplied pow14.5
Applied pow14.5
Applied pow14.5
Applied pow14.5
Applied pow-prod-down4.5
Applied pow-prod-down4.5
Applied pow-prod-down4.5
Simplified1.1
if 4.4040728131748635e+58 < t Initial program 0.6
rmApplied associate-*l*0.5
rmApplied associate-*l*0.6
Simplified0.6
rmApplied pow10.6
Applied pow10.6
Applied pow-prod-down0.6
Applied pow10.6
Applied pow-prod-down0.6
Simplified0.5
Final simplification1.0
herbie shell --seed 2019209
(FPCore (x y z t a b)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, A"
:precision binary64
:herbie-target
(if (< y 7.590524218811189e-161) (+ (- (* x 2) (* (* (* y 9) z) t)) (* a (* 27 b))) (+ (- (* x 2) (* 9 (* y (* t z)))) (* (* a 27) b)))
(+ (- (* x 2) (* (* (* y 9) z) t)) (* (* a 27) b)))