求一判断浮点数为非零数的C++解决代码
有一次面试的时候,遇到这么一个问题 编程实现判断一个浮点数为非零数
比如:float a = 0.00001
编程实现判断a为非零数
[解决办法]
if(0xFFFFFFFF&a) a是非零数
[解决办法]
1)不管他是否准确,当计算有误差时,这个判断就不准确了。
if(a) //这对于输入数据,和常数等比较有用。
{
。。。。
}
2)判断在某个精度范围,是否可以认定为0
指定范围的数,认定为0;
这是防止除数为零错误,采取的措施,一些数值计算,也经常这样。
const double eps =1e-10; //误差精度
if(fabs(a) < eps )// 或者 if(fabs(a) <= eps )//
{
.....
}
[解决办法]
如果是用C++的话,参考这个 http://msdn.microsoft.com/en-us/library/6x7575x3.aspx
就是楼上自定义的epsilon, 在limits中有标准定义,更好的写法是
#include <limits>
if(fabs(a) < std::numeric_limits<double>::epsilon())
now a can be regarded as 0;
或者
#include <cmath>
#include <limits>
template<class Float>
inline bool IsZero(Float v)
{
return std::abs(v)<=std::numeric_limits<Float>::epsilon();
}
[解决办法]
浮点数的判断只能近似比较,一般都是取差的绝对值,然后判断是否在某个精度范围内(比如if(fabs(x-0.0) < 0.00000000001)就可以认定为是0了)
[解决办法]
不管什么数据类型, 内存中都是以二进制表示的. 只需要判断每个bit位为0, 那么这个数就是0.
假设float是32bit的, 那么你完全可以把它指向的地址以DWORD来判断.
float val = 0.0f;
DWORD &dwVal = *(DWORD*)&val;
if(dwVal == 0)
{
//do something
}
[解决办法]
浮点数0分+0和-0, 都可以直接与0比较.
double f;
//...
bool zero = f == 0;
bool nz = f != 0;
[解决办法]
你如果这么想就进入歧路了,浮点数的位排列不能反映数的大小。
一般正规的用法if(fabs(x)<1e-8) 其中的1e-8根据你的精度要求来定。
[解决办法]
请问:epsilon定义为多少最好呢?分double型和float型讨论。
定义成多少要看你需要表达多大的数值,限制精度的不是数值大小,而是有效数字个数
float的精度大概是10^7分之一,double的精度大概是10^15分之一,如果你用float表示1亿左右的数值,那么比较是否为0的阈值就应该是10左右,如果你用float表示1左右的数值,那么这个阈值就可以是0.0000001左右
[解决办法]
include\float.h
...
#define DBL_DIG 15 /* # of decimal digits of precision */
#define DBL_EPSILON 2.2204460492503131e-016 /* smallest such that 1.0+DBL_EPSILON != 1.0 */
#define DBL_MANT_DIG 53 /* # of bits in mantissa */
#define DBL_MAX 1.7976931348623158e+308 /* max value */
#define DBL_MAX_10_EXP 308 /* max decimal exponent */
#define DBL_MAX_EXP 1024 /* max binary exponent */
#define DBL_MIN 2.2250738585072014e-308 /* min positive value */
#define DBL_MIN_10_EXP (-307) /* min decimal exponent */
#define DBL_MIN_EXP (-1021) /* min binary exponent */
#define _DBL_RADIX 2 /* exponent radix */
#define _DBL_ROUNDS 1 /* addition rounding: near */
#define FLT_DIG 6 /* # of decimal digits of precision */
#define FLT_EPSILON 1.192092896e-07F /* smallest such that 1.0+FLT_EPSILON != 1.0 */
#define FLT_GUARD 0
#define FLT_MANT_DIG 24 /* # of bits in mantissa */
#define FLT_MAX 3.402823466e+38F /* max value */
#define FLT_MAX_10_EXP 38 /* max decimal exponent */
#define FLT_MAX_EXP 128 /* max binary exponent */
#define FLT_MIN 1.175494351e-38F /* min positive value */
#define FLT_MIN_10_EXP (-37) /* min decimal exponent */
#define FLT_MIN_EXP (-125) /* min binary exponent */
#define FLT_NORMALIZE 0
#define FLT_RADIX 2 /* exponent radix */
#define FLT_ROUNDS 1 /* addition rounding: near */
...