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 格式打印 float 。这会导致未定义的行为。使用:

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

你会发现你得到了正确的答案。