在 C 中使用 openmp 时的数值差异 [重复]

Question

这个问题在这里已经有了答案：

浮点数学被破坏了吗？ (30 个回答)

4天前关闭。

我正在做一个代码来求解一个积分-微分方程，这个方程有一个参数，我想看看如何影响我的解决方案，当在串行版本中运行 for 循环时，我有一个结果列表，但是当运行代码使用 omp 部分将循环一分为二，前半部分相同，但第二部分在小数点后第三位后更改，但我需要相同的结果。

我能做些什么来解决这个问题？

系列版本

for (int i = 0; i < vetores; i++)
{
   eta = 1.0;
   D   = 1.0;
   K = (double) i*1.0;
    
   K_vetor[i] = calculo(theta,theta_zero,omega,x,t,eta,D,K,g,f);
    
   printf("K %d\n",i);

}

Openmp 版本

#pragma omp parallel default(none) shared(x,omega,K_vetor,vetores) firstprivate(theta_zero,theta,theta2,theta3,t) private(eta,D,K)  //private(i,aux,temp1,temp2) 
{
    #pragma omp sections nowait
    {
        #pragma omp section
        {
            for (int i = 0; i < vetores/2; i++)
            {
                eta = 1.0;
                D   = 1.0;
                K = (double) i*1.0;
                
                K_vetor[i] = calculo(theta,theta_zero,omega,x,t,eta,D,K,g,f);
                
                printf("K %d\n",i);

            }
        }
        #pragma omp section
        {
            for (int i = vetores/2; i < vetores; i++)
            {
                eta = 1.0;
                D   = 1.0;
                K = (double) i*1.0;
                
                K_vetor[i] = calculo(theta2,theta_zero,omega,x,t,eta,D,K,g,f);
                
                printf("K %d\n",i);

            }
        }
    }
}

差速器 output

OMP                                 Serial
0 0.406062                          0 0.406062
1 0.910302                          1 0.910302
2 0.95234                           2 0.95234
3 0.935899                          3 0.935899
4 0.845607                          4 0.845607
5 0.800523                          5 0.800523
6 0.784356                         |6 0.785621
7 0.771353                         |7 0.771651
8 0.742191                         |8 0.742313
9 0.764885                         |9 0.76473
10 0.774907                        |10 0.775604
11 0.776134                        |11 0.77482
12 0.817413                        |12 0.813391

Obs：我所有的变量都是双倍的。

theta 和 theta2 的一些图像

Answer 1

这个问题在这里已经有了答案：

浮点数学被破坏了吗？ (30 个回答)

4天前关闭。

我正在做一个代码来求解一个积分-微分方程，这个方程有一个参数，我想看看如何影响我的解决方案，当在串行版本中运行 for 循环时，我有一个结果列表，但是当运行代码使用 omp 部分将循环一分为二，前半部分相同，但第二部分在小数点后第三位后更改，但我需要相同的结果。

我能做些什么来解决这个问题？

系列版本

for (int i = 0; i < vetores; i++)
{
   eta = 1.0;
   D   = 1.0;
   K = (double) i*1.0;
    
   K_vetor[i] = calculo(theta,theta_zero,omega,x,t,eta,D,K,g,f);
    
   printf("K %d\n",i);

}

Openmp 版本

#pragma omp parallel default(none) shared(x,omega,K_vetor,vetores) firstprivate(theta_zero,theta,theta2,theta3,t) private(eta,D,K)  //private(i,aux,temp1,temp2) 
{
    #pragma omp sections nowait
    {
        #pragma omp section
        {
            for (int i = 0; i < vetores/2; i++)
            {
                eta = 1.0;
                D   = 1.0;
                K = (double) i*1.0;
                
                K_vetor[i] = calculo(theta,theta_zero,omega,x,t,eta,D,K,g,f);
                
                printf("K %d\n",i);

            }
        }
        #pragma omp section
        {
            for (int i = vetores/2; i < vetores; i++)
            {
                eta = 1.0;
                D   = 1.0;
                K = (double) i*1.0;
                
                K_vetor[i] = calculo(theta2,theta_zero,omega,x,t,eta,D,K,g,f);
                
                printf("K %d\n",i);

            }
        }
    }
}

差速器 output

OMP                                 Serial
0 0.406062                          0 0.406062
1 0.910302                          1 0.910302
2 0.95234                           2 0.95234
3 0.935899                          3 0.935899
4 0.845607                          4 0.845607
5 0.800523                          5 0.800523
6 0.784356                         |6 0.785621
7 0.771353                         |7 0.771651
8 0.742191                         |8 0.742313
9 0.764885                         |9 0.76473
10 0.774907                        |10 0.775604
11 0.776134                        |11 0.77482
12 0.817413                        |12 0.813391

Obs：我所有的变量都是双倍的。

theta 和 theta2 的一些图像