\frac{\left(\frac{\left(-b_2\right)}{\left(\sqrt{\left(\left(b_2 \cdot b_2\right) - \left(a \cdot c\right)\right)}\right)}\right)}{a}\frac{\sqrt{\frac{b_2 \cdot b_2 + c \cdot a}{\frac{b_2 \cdot b_2 + c \cdot a}{b_2 \cdot b_2 - c \cdot a}}} - b_2}{a}double f(double a, double b_2, double c) {
double r743371 = b_2;
double r743372 = -r743371;
double r743373 = r743371 * r743371;
double r743374 = a;
double r743375 = c;
double r743376 = r743374 * r743375;
double r743377 = r743373 - r743376;
double r743378 = sqrt(r743377);
double r743379 = r743372 + r743378;
double r743380 = r743379 / r743374;
return r743380;
}
double f(double a, double b_2, double c) {
double r743381 = b_2;
double r743382 = r743381 * r743381;
double r743383 = c;
double r743384 = a;
double r743385 = r743383 * r743384;
double r743386 = r743382 + r743385;
double r743387 = r743382 - r743385;
double r743388 = r743386 / r743387;
double r743389 = r743386 / r743388;
double r743390 = sqrt(r743389);
double r743391 = r743390 - r743381;
double r743392 = r743391 / r743384;
return r743392;
}



Bits error versus a



Bits error versus b_2



Bits error versus c
Initial program 1.7
Simplified1.7
rmApplied p16-flip--2.7
rmApplied difference-of-squares2.6
Applied associate-/l*1.7
Final simplification1.7
herbie shell --seed 2019121
(FPCore (a b_2 c)
:name "quad2p (problem 3.2.1, positive)"
(/.p16 (+.p16 (neg.p16 b_2) (sqrt.p16 (-.p16 (*.p16 b_2 b_2) (*.p16 a c)))) a))