\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{\mathsf{log1p}\left(\mathsf{expm1}\left(\sqrt{\frac{8}{3}} \cdot \left(\sqrt{\frac{8}{3}} \cdot \sin \left(x \cdot 0.5\right)\right)\right)\right)}{\frac{\sin x}{\sin \left(0.5 \cdot x\right)}}double f(double x) {
double r616512 = 8.0;
double r616513 = 3.0;
double r616514 = r616512 / r616513;
double r616515 = x;
double r616516 = 0.5;
double r616517 = r616515 * r616516;
double r616518 = sin(r616517);
double r616519 = r616514 * r616518;
double r616520 = r616519 * r616518;
double r616521 = sin(r616515);
double r616522 = r616520 / r616521;
return r616522;
}
double f(double x) {
double r616523 = 8.0;
double r616524 = 3.0;
double r616525 = r616523 / r616524;
double r616526 = sqrt(r616525);
double r616527 = x;
double r616528 = 0.5;
double r616529 = r616527 * r616528;
double r616530 = sin(r616529);
double r616531 = r616526 * r616530;
double r616532 = r616526 * r616531;
double r616533 = expm1(r616532);
double r616534 = log1p(r616533);
double r616535 = sin(r616527);
double r616536 = r616528 * r616527;
double r616537 = sin(r616536);
double r616538 = r616535 / r616537;
double r616539 = r616534 / r616538;
return r616539;
}




Bits error versus x
Results
| Original | 15.3 |
|---|---|
| Target | 0.3 |
| Herbie | 0.5 |
Initial program 15.3
rmApplied associate-/l*0.5
Simplified0.5
rmApplied add-sqr-sqrt0.5
Applied associate-*l*0.4
rmApplied log1p-expm1-u0.5
Final simplification0.5
herbie shell --seed 2020046 +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)))