\frac{1}{2} \cdot \left(x + y \cdot \sqrt{z}\right)x \cdot \frac{1}{2} + \left(\frac{1}{2} \cdot \sqrt{z}\right) \cdot ydouble f(double x, double y, double z) {
double r219604 = 1.0;
double r219605 = 2.0;
double r219606 = r219604 / r219605;
double r219607 = x;
double r219608 = y;
double r219609 = z;
double r219610 = sqrt(r219609);
double r219611 = r219608 * r219610;
double r219612 = r219607 + r219611;
double r219613 = r219606 * r219612;
return r219613;
}
double f(double x, double y, double z) {
double r219614 = x;
double r219615 = 1.0;
double r219616 = 2.0;
double r219617 = r219615 / r219616;
double r219618 = r219614 * r219617;
double r219619 = z;
double r219620 = sqrt(r219619);
double r219621 = r219617 * r219620;
double r219622 = y;
double r219623 = r219621 * r219622;
double r219624 = r219618 + r219623;
return r219624;
}



Bits error versus x



Bits error versus y



Bits error versus z
Results
Initial program 0.2
rmApplied add-sqr-sqrt1.3
Applied associate-*l*1.0
rmApplied distribute-lft-in0.9
Applied distribute-lft-in1.0
Simplified0.5
Simplified0.1
Final simplification0.1
herbie shell --seed 2019235
(FPCore (x y z)
:name "Diagrams.Solve.Polynomial:quadForm from diagrams-solve-0.1, B"
:precision binary64
(* (/ 1 2) (+ x (* y (sqrt z)))))