double code(double x) {
return ((double) (((double) (1.0 - ((double) cos(x)))) / ((double) sin(x))));
}
double code(double x) {
double VAR;
if ((((double) (((double) (1.0 - ((double) cos(x)))) / ((double) sin(x)))) <= -0.003356608539896645)) {
VAR = ((double) log(((double) exp(((double) (((double) (1.0 - ((double) cos(x)))) / ((double) sin(x))))))));
} else {
double VAR_1;
if ((((double) (((double) (1.0 - ((double) cos(x)))) / ((double) sin(x)))) <= 4.6661405104372516e-07)) {
VAR_1 = ((double) (((double) (0.041666666666666664 * ((double) pow(x, 3.0)))) + ((double) (((double) (0.004166666666666667 * ((double) pow(x, 5.0)))) + ((double) (x * 0.5))))));
} else {
VAR_1 = ((double) (((double) (((double) log(((double) exp(((double) (((double) pow(1.0, 3.0)) - ((double) pow(((double) cos(x)), 3.0)))))))) / ((double) (((double) (1.0 * 1.0)) + ((double) (((double) cos(x)) * ((double) (1.0 + ((double) cos(x)))))))))) / ((double) sin(x))));
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x
Results
| Original | 30.3 |
|---|---|
| Target | 0.0 |
| Herbie | 0.6 |
if (/ (- 1.0 (cos x)) (sin x)) < -0.00335660853989664483Initial program 0.9
rmApplied add-log-exp1.1
if -0.00335660853989664483 < (/ (- 1.0 (cos x)) (sin x)) < 4.6661405104372516e-7Initial program 60.1
Taylor expanded around 0 0.1
Simplified0.1
if 4.6661405104372516e-7 < (/ (- 1.0 (cos x)) (sin x)) Initial program 1.2
rmApplied flip3--1.3
Simplified1.3
rmApplied add-log-exp1.3
Applied add-log-exp1.3
Applied diff-log1.4
Simplified1.3
Final simplification0.6
herbie shell --seed 2020191
(FPCore (x)
:name "tanhf (example 3.4)"
:precision binary64
:herbie-expected 2
:herbie-target
(tan (/ x 2.0))
(/ (- 1.0 (cos x)) (sin x)))