\frac{x \cdot x}{y \cdot y} + \frac{z \cdot z}{t \cdot t}\mathsf{hypot}\left(\frac{z}{t}, \frac{x}{y}\right) \cdot \left(\sqrt{\mathsf{hypot}\left(\frac{z}{t}, \frac{x}{y}\right)} \cdot \sqrt{\mathsf{hypot}\left(\frac{z}{t}, \frac{x}{y}\right)}\right)double f(double x, double y, double z, double t) {
double r611996 = x;
double r611997 = r611996 * r611996;
double r611998 = y;
double r611999 = r611998 * r611998;
double r612000 = r611997 / r611999;
double r612001 = z;
double r612002 = r612001 * r612001;
double r612003 = t;
double r612004 = r612003 * r612003;
double r612005 = r612002 / r612004;
double r612006 = r612000 + r612005;
return r612006;
}
double f(double x, double y, double z, double t) {
double r612007 = z;
double r612008 = t;
double r612009 = r612007 / r612008;
double r612010 = x;
double r612011 = y;
double r612012 = r612010 / r612011;
double r612013 = hypot(r612009, r612012);
double r612014 = sqrt(r612013);
double r612015 = r612014 * r612014;
double r612016 = r612013 * r612015;
return r612016;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 33.8 |
|---|---|
| Target | 0.4 |
| Herbie | 0.6 |
Initial program 33.8
Simplified20.0
rmApplied add-sqr-sqrt20.1
Simplified20.0
Simplified0.4
rmApplied add-sqr-sqrt0.6
Final simplification0.6
herbie shell --seed 2020089 +o rules:numerics
(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))))