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);

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