主页 / 问题 / 1149924
Accepted
Olga
Asked:2020-07-07 16:59:32 +0000 UTC2020-07-07 16:59:32 +0000 UTC

Fortran Visual Studio 2013 上数字乘法结果的准确性

  • 772

在 Visual Studio 2013 中创建了两个项目，第一个项目用于一个内核，第二个项目使用 OpenMP 指令用于四个内核。正在解决相同的问题（所有输入数据都相同）。结果是不同的。检查时发现，乘法结果的准确性因项目而异。将两个变量传递给要相乘的代码的非并行部分。

print*,CT1,TA(I,J,K),(CT1 * TA(I,J,K))
результат на одном ядре:  -3,5е-6    12,89648   -4,5137684е-5
результат на 4-х ядрах:   -3,5е-6    12,89648   -4,5137665е-5

Python中的乘法

a=-3.5e-6
b=12.89648
c=a*b
print("{:.10e}".format(c))
-4.5137680000e-05

变量声明相同，工作室设置相同。浮点模型设置为精确。更改这些设置无济于事。请告诉我可能是什么问题。

visual-studio-2013

1 个回答

  1. Best Answer

    正如 Captain Obvious 所指出的，在某种意义上，差异构成了最后一位数字的单位（例如，在一种情况下是结果，在另一种情况下是输入数据）：

      program t
    implicit none
    real*4 r4, a, b
    integer*4 i4
    equivalence (r4, i4)

    a = -3.5e-6
    b = 12.89648
    r4 = a*b
    print '(z8,es15.7)', i4, r4
    r4 = r4*(1.e0 + epsilon(1.e0))
    print '(z8,es15.7)', i4, r4
    a = a*(1.e0 - epsilon(1.e0))
    b = b*(1.e0 - epsilon(1.e0))
    r4 = a*b
    print '(z8,es15.7)', i4, r4
  end program t

    获取您的所有号码：

    B83D5237 -4.5137680E-05
B83D5238 -4.5137684E-05
B83D5233 -4.5137665E-05
    • -1

相关问题