\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{8 \cdot \frac{\sin \left(0.5 \cdot x\right)}{3}}{\mathsf{log1p}\left(\mathsf{expm1}\left(\frac{\sin x}{\sin \left(0.5 \cdot x\right)}\right)\right)}double f(double x) {
double r335027 = 8.0;
double r335028 = 3.0;
double r335029 = r335027 / r335028;
double r335030 = x;
double r335031 = 0.5;
double r335032 = r335030 * r335031;
double r335033 = sin(r335032);
double r335034 = r335029 * r335033;
double r335035 = r335034 * r335033;
double r335036 = sin(r335030);
double r335037 = r335035 / r335036;
return r335037;
}
double f(double x) {
double r335038 = 8.0;
double r335039 = 0.5;
double r335040 = x;
double r335041 = r335039 * r335040;
double r335042 = sin(r335041);
double r335043 = 3.0;
double r335044 = r335042 / r335043;
double r335045 = r335038 * r335044;
double r335046 = sin(r335040);
double r335047 = r335046 / r335042;
double r335048 = expm1(r335047);
double r335049 = log1p(r335048);
double r335050 = r335045 / r335049;
return r335050;
}




Bits error versus x
Results
| Original | 14.8 |
|---|---|
| Target | 0.3 |
| Herbie | 0.4 |
Initial program 14.8
rmApplied associate-/l*0.5
Simplified0.5
rmApplied div-inv0.5
Applied associate-*l*0.5
Simplified0.3
rmApplied log1p-expm1-u0.4
Final simplification0.4
herbie shell --seed 2019322 +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)))