\frac{\left(\frac{8}{3} \cdot \sin \left(x \cdot 0.5\right)\right) \cdot \sin \left(x \cdot 0.5\right)}{\sin x}\frac{\sin \left(x \cdot 0.5\right)}{\frac{\mathsf{expm1}\left(\mathsf{log1p}\left(3 \cdot \frac{\sin x}{8}\right)\right)}{\sin \left(x \cdot 0.5\right)}}double f(double x) {
double r375737 = 8.0;
double r375738 = 3.0;
double r375739 = r375737 / r375738;
double r375740 = x;
double r375741 = 0.5;
double r375742 = r375740 * r375741;
double r375743 = sin(r375742);
double r375744 = r375739 * r375743;
double r375745 = r375744 * r375743;
double r375746 = sin(r375740);
double r375747 = r375745 / r375746;
return r375747;
}
double f(double x) {
double r375748 = x;
double r375749 = 0.5;
double r375750 = r375748 * r375749;
double r375751 = sin(r375750);
double r375752 = 3.0;
double r375753 = sin(r375748);
double r375754 = 8.0;
double r375755 = r375753 / r375754;
double r375756 = r375752 * r375755;
double r375757 = log1p(r375756);
double r375758 = expm1(r375757);
double r375759 = r375758 / r375751;
double r375760 = r375751 / r375759;
return r375760;
}




Bits error versus x
Results
| Original | 15.2 |
|---|---|
| Target | 0.3 |
| Herbie | 0.4 |
Initial program 15.2
Simplified15.1
rmApplied associate-/l*0.3
Simplified0.4
rmApplied expm1-log1p-u0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2019174 +o rules:numerics
(FPCore (x)
:name "Graphics.Rasterific.Svg.PathConverter:segmentToBezier from rasterific-svg-0.2.3.1, A"
:herbie-target
(/ (/ (* 8.0 (sin (* x 0.5))) 3.0) (/ (sin x) (sin (* x 0.5))))
(/ (* (* (/ 8.0 3.0) (sin (* x 0.5))) (sin (* x 0.5))) (sin x)))