c - Openmp减少循环错误？

Question

这是测试代码的一部分：

    float a = 0;
    float b = 0;
    int c = 0;
    int d = 0;
    #pragma omp parallel for schedule (dynamic, 1)  reduction(+ : a, b, c, d)
    for(i=0; i<100; i++) {
        a +=1;
        b +=1;
        c +=1;
        d +=1;
    }

    printf("a: %d, b: %d, c: %d, d: %d\n", a, b, c, d);

由于某些原因，我的结果总是：

a: 100, b: 100, c: 0, d: 202
a: 100, b: 100, c: 0, d: 202
a: 100, b: 100, c: 0, d: 202
a: 100, b: 100, c: 0, d: 202
a: 100, b: 100, c: 0, d: 202
a: 100, b: 100, c: 0, d: 202
a: 100, b: 100, c: 0, d: 202

为什么 a、b、c、d 不都等于 100？

Answer 1

您正在使用%d格式来打印浮点数。这会导致未定义的行为。利用：

   printf("a: %f, b: %f, c: %d, d: %d\n", a, b, c, d);

你会看到你得到了正确的答案。