\log \left(x + \sqrt{x \cdot x + 1}\right)\begin{array}{l}
\mathbf{if}\;x \le -1.021894228612228561203778554045129567385:\\
\;\;\;\;\log \left(\frac{0.125}{{x}^{3}} - \left(\frac{0.0625}{{x}^{5}} + \frac{0.5}{x}\right)\right)\\
\mathbf{elif}\;x \le 0.8867151979513619997064211020187940448523:\\
\;\;\;\;\left(\log \left(\sqrt{1}\right) + \frac{x}{\sqrt{1}}\right) - \frac{\frac{1}{6}}{1} \cdot \frac{{x}^{3}}{\sqrt{1}}\\
\mathbf{else}:\\
\;\;\;\;\log \left(\left(\left(x + \frac{0.5}{x}\right) - \frac{0.125}{{x}^{3}}\right) + x\right)\\
\end{array}double f(double x) {
double r111382 = x;
double r111383 = r111382 * r111382;
double r111384 = 1.0;
double r111385 = r111383 + r111384;
double r111386 = sqrt(r111385);
double r111387 = r111382 + r111386;
double r111388 = log(r111387);
return r111388;
}
double f(double x) {
double r111389 = x;
double r111390 = -1.0218942286122286;
bool r111391 = r111389 <= r111390;
double r111392 = 0.125;
double r111393 = 3.0;
double r111394 = pow(r111389, r111393);
double r111395 = r111392 / r111394;
double r111396 = 0.0625;
double r111397 = 5.0;
double r111398 = pow(r111389, r111397);
double r111399 = r111396 / r111398;
double r111400 = 0.5;
double r111401 = r111400 / r111389;
double r111402 = r111399 + r111401;
double r111403 = r111395 - r111402;
double r111404 = log(r111403);
double r111405 = 0.886715197951362;
bool r111406 = r111389 <= r111405;
double r111407 = 1.0;
double r111408 = sqrt(r111407);
double r111409 = log(r111408);
double r111410 = r111389 / r111408;
double r111411 = r111409 + r111410;
double r111412 = 0.16666666666666666;
double r111413 = r111412 / r111407;
double r111414 = r111394 / r111408;
double r111415 = r111413 * r111414;
double r111416 = r111411 - r111415;
double r111417 = r111389 + r111401;
double r111418 = r111417 - r111395;
double r111419 = r111418 + r111389;
double r111420 = log(r111419);
double r111421 = r111406 ? r111416 : r111420;
double r111422 = r111391 ? r111404 : r111421;
return r111422;
}




Bits error versus x
Results
| Original | 53.1 |
|---|---|
| Target | 45.4 |
| Herbie | 0.2 |
if x < -1.0218942286122286Initial program 63.3
Taylor expanded around -inf 0.1
Simplified0.1
if -1.0218942286122286 < x < 0.886715197951362Initial program 58.7
Taylor expanded around 0 0.3
Simplified0.3
if 0.886715197951362 < x Initial program 31.9
Taylor expanded around inf 0.2
Simplified0.2
Final simplification0.2
herbie shell --seed 2019196
(FPCore (x)
:name "Hyperbolic arcsine"
:herbie-target
(if (< x 0.0) (log (/ -1.0 (- x (sqrt (+ (* x x) 1.0))))) (log (+ x (sqrt (+ (* x x) 1.0)))))
(log (+ x (sqrt (+ (* x x) 1.0)))))