\frac{x \cdot x}{y \cdot y} + \frac{z \cdot z}{t \cdot t}\left(\frac{x}{y} \cdot \frac{\sqrt[3]{x} \cdot \left(\sqrt[3]{\sqrt[3]{x} \cdot \sqrt[3]{x}} \cdot \sqrt[3]{\sqrt[3]{x}}\right)}{\sqrt[3]{y} \cdot \sqrt[3]{y}}\right) \cdot \frac{\sqrt[3]{x}}{\sqrt[3]{y}} + \frac{z}{t} \cdot \frac{z}{t}double f(double x, double y, double z, double t) {
double r648468 = x;
double r648469 = r648468 * r648468;
double r648470 = y;
double r648471 = r648470 * r648470;
double r648472 = r648469 / r648471;
double r648473 = z;
double r648474 = r648473 * r648473;
double r648475 = t;
double r648476 = r648475 * r648475;
double r648477 = r648474 / r648476;
double r648478 = r648472 + r648477;
return r648478;
}
double f(double x, double y, double z, double t) {
double r648479 = x;
double r648480 = y;
double r648481 = r648479 / r648480;
double r648482 = cbrt(r648479);
double r648483 = r648482 * r648482;
double r648484 = cbrt(r648483);
double r648485 = cbrt(r648482);
double r648486 = r648484 * r648485;
double r648487 = r648482 * r648486;
double r648488 = cbrt(r648480);
double r648489 = r648488 * r648488;
double r648490 = r648487 / r648489;
double r648491 = r648481 * r648490;
double r648492 = r648482 / r648488;
double r648493 = r648491 * r648492;
double r648494 = z;
double r648495 = t;
double r648496 = r648494 / r648495;
double r648497 = r648496 * r648496;
double r648498 = r648493 + r648497;
return r648498;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 33.6 |
|---|---|
| Target | 0.4 |
| Herbie | 0.9 |
Initial program 33.6
rmApplied times-frac18.7
rmApplied times-frac0.4
rmApplied add-cube-cbrt0.8
Applied add-cube-cbrt0.9
Applied times-frac0.9
Applied associate-*r*0.9
rmApplied add-cube-cbrt0.9
Applied cbrt-prod0.9
Final simplification0.9
herbie shell --seed 2019322
(FPCore (x y z t)
:name "Graphics.Rasterific.Svg.PathConverter:arcToSegments from rasterific-svg-0.2.3.1"
:precision binary64
:herbie-target
(+ (pow (/ x y) 2) (pow (/ z t) 2))
(+ (/ (* x x) (* y y)) (/ (* z z) (* t t))))