\sqrt{16 \cdot {x}^{2} + 1} - 4 \cdot x\sqrt{16 \cdot {x}^{2} + 1} - 4 \cdot xdouble code(double x) {
return ((double) (((double) sqrt(((double) (((double) (16.0 * ((double) pow(x, 2.0)))) + 1.0)))) - ((double) (4.0 * x))));
}
double code(double x) {
return ((double) (((double) sqrt(((double) (((double) (16.0 * ((double) pow(x, 2.0)))) + 1.0)))) - ((double) (4.0 * x))));
}



Bits error versus x
Results
Initial program 23.7
Final simplification23.7
herbie shell --seed 2020153
(FPCore (x)
:name "(- (sqrt (+ (* 16 (pow x 2)) 1)) (* 4 x))"
:precision binary64
(- (sqrt (+ (* 16.0 (pow x 2.0)) 1.0)) (* 4.0 x)))