奇怪的 C 程序行为

Question

我有以下 C 程序：

#include <stdio.h>
#include <stdlib.h>
#include <assert.h>
#include <math.h>

int main() {
    const int opt_count = 2;

    int oc = 30;

    int c = 900;

    printf("%d %f\n", c, pow(oc, opt_count));
    assert(c == (int)(pow(oc, opt_count)));
}

我在 Windows 8.1 上运行 MinGW。 Gcc 版本 4.9.3。我编译我的程序：

gcc program.c -o program.exe

当我运行它时，我得到这个输出：

$ program
900 900.000000
Assertion failed: c == (int)(pow(oc, opt_count)), file program.c, line 16

This application has requested the Runtime to terminate it in an unusual way.
Please contact the application's support team for more information.

发生了什么事？我希望断言能够通过，因为 900 == 30^2。

谢谢！

编辑

我没有使用任何分数或小数。我只使用整数。

Answer 1

当pow 的实现是通过时会发生这种情况

pow(x,y) = exp(log(x)*y)

其他库实现首先通过整数幂来减少指数，从而避免这种小的浮点错误。

---

更多涉及的实现包含如下步骤

pow(x,y) {
    if(y<0) return 1/pow(x, -y);

    n = (int)round(y);
    y = y-n;
    px = x; powxn = 1;
    while(n>0) {
        if(n%2==1) powxn *= px;
        n /=2; px *= px;
    }
    return powxn * exp(log(x)*y);
}

与通常的分而治之分别。整数幂 powxn 的减半平方方法。

Answer 2

@LutzL 给你一个很好的答案（和解决方案），另一个解决方案是将差异与 epsilon 进行比较，例如：0.00001，这样你就可以使用 math.h 中包含的标准函数 pow

#define EPSILON 0.0001
#define EQ(a, b) (fabs(a - b) < EPSILON)

assert(EQ((double)c, pow(oc, opt_count)));