#include "libm.h"#if FLT_EVAL_METHOD==0 || FLT_EVAL_METHOD==1#define EPS DBL_EPSILON#elif FLT_EVAL_METHOD==2#define EPS LDBL_EPSILON#endifstatic const double_t toint = 1/EPS;double floor(double x){union {double f; uint64_t i;} u = {x};int e = u.i >> 52 & 0x7ff;double_t y;if (e >= 0x3ff+52 || x == 0)return x;/* y = int(x) - x, where int(x) is an integer neighbor of x */if (u.i >> 63)y = x - toint + toint - x;elsey = x + toint - toint - x;/* special case because of non-nearest rounding modes */if (e <= 0x3ff-1) {FORCE_EVAL(y);return u.i >> 63 ? -1 : 0;}if (y > 0)return x + y - 1;return x + y;}
此处可能存在不合适展示的内容,页面不予展示。您可通过相关编辑功能自查并修改。
如您确认内容无涉及 不当用语 / 纯广告导流 / 暴力 / 低俗色情 / 侵权 / 盗版 / 虚假 / 无价值内容或违法国家有关法律法规的内容,可点击提交进行申诉,我们将尽快为您处理。