3

在进行算术运算或尝试打印（调试）这么大的数字时，我会失去精度：1234567890.123456789

我认为我的问题在于 $d（算术结果）和 $e 的格式化打印。我怎样才能强制长双打？我的 Perl 版本（SUN 上的 5.8.4）说这是可能的。sprintf 有一个长双精度的大小选项（q 或 L 或 ll），但我还没有弄清楚如何使用它，也不知道它是否适用于 printf。

编辑：我添加了 BigFloat，它有效！但我仍然想强制长双打。

尝试添加 1234567890 + 0.123456789 并减去 1234567890 - 0.123456789。

use Config;
use Math::BigFloat;
$a = 1234567890;
$b = 123456789;
$c = $b/1e9;                   # 0.123456789
$d = $a + $c;                  # not enough precision (32-bit or double?)
$e = sprintf("%d.%.9d",$a,$b); # combine as strings
$f = 1234567890.123456789;     # for reference (not enough precision)

# Use BigFloat to bypass lack of longdbl
$aBig = Math::BigFloat->new("$a");
$dSum = $aBig->fadd("$c");         # $dSum = $a + $c
$aBig = Math::BigFloat->new("$a"); # <-- Need a new one for every operation?
$dDif = $aBig->fsub(abs("$c"));    # $dDif = $a - $c

print "a $a\n";             # 1234567890  
print "c $c\n";             # 0.123456789
print "d=a+c $d\n";         # 1234567890.12346  <-- **Problem**
print "dSum=a+c $dSum\n";   # 1234567890.123456789  <-- Solution
print "dDif=a-c $dDif\n";   # 1234567890.876543211  <-- Solution
print "e $e\n";             # 1234567890.123456789
print "f $f\n";             # 1234567890.12346  <-- double, 52-bit, not longdbl? 
printf ("printf    e 20.9f %20.9f\n",$e);    # 1234567890.123456717 <-- **Problem**
printf ("printf dSum 20.9f %20.9f\n",$dSum); # 1234567890.123456717 <-- **Problem**
printf ("printf dSum 20s %20s\n",$dSum);     # 1234567890.123456789 
printf ("printf dDif 20.9f %20.9f\n",$dDif); # 1234567890.876543283 <-- **Problem**
printf ("printf dDif 20s %20s\n",$dDif);     # 1234567890.876543211 

print "uselongdouble $Config{uselongdouble}\n"; # empty. No long doubles by default
print "d_longdbl $Config{d_longdbl}\n";         # "define". Supports long doubles
print "size double longdbl $Config{doublesize} $Config{longdblsize}\n"; # Ans 8 16

我还使用此代码来尝试理解类型，但没有太大帮助。有没有人用它来解释这样的问题？

use Devel::Peek 'Dump';

Dump ($dSum); # Wow, it's complicated
Dump ($f);
4

2 回答 2

3

Perl 有一种浮点数，称为NV. an 的大小NV是在构建 Perl 时决定的。

$ perl -V:nvsize
nvsize='8';

在 Perl 程序中也可以通过 Config 模块访问此信息。

$ perl -MConfig -E'say $Config{nvsize}'
8

NVPerl 构建后无法更改大小。

在构建 Perl 时，您可以强制构建 Perl 以使用long double浮点数，如下所示：

sh Configure -Duselongdouble ...
  -or-
perlbrew install -Duselongdouble ...
  -or-
perlbrew install --ld ...

如果您不想重建perl或制作新的，则需要使用模块。我推荐Math::LongDouble，因为它提供了对原生long double浮点数的访问，并且它尽可能透明。

另一种选择是使用任意精度的库，例如Math::BigFloat，但如果您只需要一个long double.

于 2015-01-15T03:56:49.450 回答
1

bignum将重载当前范围内的所有运算符以使用任意精度整数和浮点运算。

use bignum;

my $f = 123456789.123456789;
print "$f\n";         # 123456789.123456789
print $f + $f, "\n";  # 246913578.246913578

在幕后，bignum 将所有数值常量适当地转换为 Math::BigInt 和 Math::BigNum 对象。

重要的是要注意 bignum 是词法范围的，不会影响整个程序。例如...

{
    use bignum;
    $f = 123456789.123456789;  # $f is a Math::BigNum object
}

$g = 123456789.123456789;      # $g is a regular NV

print "$f\n";     # 123456789.123456789
print "$g\n";     # 123456789.123457

# This will use Math::BigNum's addition method, but $g has already lost precision.
print $f + $g, "\n";  # 246913578.246913789

通过使用Math::BigInt::GMP插件来使用GNU 多精度算术库进行某些操作，您可以获得更好的性能。您必须首先使用正常的 CPAN 安装过程安装该模块（您不需要 GMP）。然后告诉 bignum 使用 GMP。

use bignum lib => "GMP";
于 2015-01-15T00:24:12.943 回答