#include "libm.h"#if FLT_EVAL_METHOD==2#undef sqrt#define sqrt sqrtl#endif/* acosh(x) = log(x + sqrt(x*x-1)) */double acosh(double x){union {double f; uint64_t i;} u = {.f = x};unsigned e = u.i >> 52 & 0x7ff;/* x < 1 domain error is handled in the called functions */if (e < 0x3ff + 1)/* |x| < 2, up to 2ulp error in [1,1.125] */return log1p(x-1 + sqrt((x-1)*(x-1)+2*(x-1)));if (e < 0x3ff + 26)/* |x| < 0x1p26 */return log(2*x - 1/(x+sqrt(x*x-1)));/* |x| >= 0x1p26 or nan */return log(x) + 0.693147180559945309417232121458176568;}
此处可能存在不合适展示的内容,页面不予展示。您可通过相关编辑功能自查并修改。
如您确认内容无涉及 不当用语 / 纯广告导流 / 暴力 / 低俗色情 / 侵权 / 盗版 / 虚假 / 无价值内容或违法国家有关法律法规的内容,可点击提交进行申诉,我们将尽快为您处理。