\frac{\left(\frac{8}{3} \cdot \sin \left(x \cdot 0.5\right)\right) \cdot \sin \left(x \cdot 0.5\right)}{\sin x}\left(8 \cdot \frac{\sin \left(0.5 \cdot x\right)}{3}\right) \cdot \mathsf{log1p}\left(\mathsf{expm1}\left(\frac{\sin \left(0.5 \cdot x\right)}{\sin x}\right)\right)double f(double x) {
double r592977 = 8.0;
double r592978 = 3.0;
double r592979 = r592977 / r592978;
double r592980 = x;
double r592981 = 0.5;
double r592982 = r592980 * r592981;
double r592983 = sin(r592982);
double r592984 = r592979 * r592983;
double r592985 = r592984 * r592983;
double r592986 = sin(r592980);
double r592987 = r592985 / r592986;
return r592987;
}
double f(double x) {
double r592988 = 8.0;
double r592989 = 0.5;
double r592990 = x;
double r592991 = r592989 * r592990;
double r592992 = sin(r592991);
double r592993 = 3.0;
double r592994 = r592992 / r592993;
double r592995 = r592988 * r592994;
double r592996 = sin(r592990);
double r592997 = r592992 / r592996;
double r592998 = expm1(r592997);
double r592999 = log1p(r592998);
double r593000 = r592995 * r592999;
return r593000;
}




Bits error versus x
Results
| Original | 14.9 |
|---|---|
| Target | 0.3 |
| Herbie | 0.6 |
Initial program 14.9
rmApplied *-un-lft-identity14.9
Applied times-frac0.5
Simplified0.5
Simplified0.5
rmApplied div-inv0.5
Applied associate-*l*0.5
Simplified0.3
rmApplied log1p-expm1-u0.6
Final simplification0.6
herbie shell --seed 2020036 +o rules:numerics
(FPCore (x)
:name "Graphics.Rasterific.Svg.PathConverter:segmentToBezier from rasterific-svg-0.2.3.1, A"
:precision binary64
:herbie-target
(/ (/ (* 8 (sin (* x 0.5))) 3) (/ (sin x) (sin (* x 0.5))))
(/ (* (* (/ 8 3) (sin (* x 0.5))) (sin (* x 0.5))) (sin x)))