这学期当数据结构的助教,使用了一些自动化的方法来改作业。其中一道题需要读入学生的输出(n个浮点数),跟标准答案进行比较,误差在一定范围内都算对。有一个学生直接在代码里输出了n个nan,然后惊奇地发现居然通过了测试。
检查我的代码之后发现问题出在这一句:
if(fabs(rx[i]-x[i])/fabs(rx[i])>0.1)
return 0;
else
return 1;
其中rx是标准答案,x是学生的输出。如果读入的x是一个正常的浮点数,那么当然没有问题;但如果x是nan,那么if语句的条件永远不会满足,因为nan没法跟其他数进行比较(这点跟inf不同)。
定位了问题之后解决方法也很简单,每次读入x之后先检查一下是否是nan即可,C语言提供了一些很好用的函数,包括:
- int fpclassify(x) 用来查看浮点数x的情况,fpclassify可以用任何浮点数表达式作为参数,fpclassify的返回值有以下几种情况。
- FP_NAN:x是一个“not a number”。
- FP_INFINITE: x是正、负无穷。
- FP_ZERO: x是0。
- FP_SUBNORMAL: x太小,以至于不能用浮点数的规格化形