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