Calculating 10 * 1.000001 ** 6e7 should minimize the propagation of rounding errors, if exponentiation is sufficiently well implemented.
I have found a case where that's not so for the calculation in question - namely, Windows 7, perl-5.20.0, nvtype is double.
But the same perl version and configuration on Ubuntu supports your assertion. So it might just be that the Windows exponentiation (or something else there) is buggy.
I haven't experimented much at all, and don't really intend to.
What's the reasoning behind your assertion ?
Cheers, Rob | [reply] |
> What's the reasoning behind your assertion ?
Exponentiation laws.
E.g. x^10 = x^8 * x^2 and x^8 = ((x^2)^2)^2
Just because 10 is 1010 binary, hence < 4*2 multiplications.
6e7 should have a 26 bit representation so <52 multiplications needed.
Of course there are even faster implementations which are less exact.
I think that's what you are observing.
edit
Of course a binary calculation should be done with integer x=1000001 without fraction, dividing later thru 10^(6*60) should be easy enough.
update
proof of concept -> here
| [reply] [d/l] [select] |
proof of concept -> here
I don't think there's any guarantee that you'll get a more accurate result using that exponentiation approach.
I ran the below script (on perl-5.16 with NV of "double") and found that sometimes your approach gave the best approximation, other times the "left-to-right" variant of your approach gave the best approximation, and other times just doing the series of multiplications gave the best approximation. (See the comments in the script.)
I used the mpfr library as my reference for the correct 15 digit value - calculated with 1000 bits of precision (which is way overkill).
It's true that your approach uses fewer calculations, but some of those squarings stand to really amplify the errors produced by the roundings.
#!perl -l
use strict;
use warnings;
use Math::MPFR qw(:mpfr);
my $e=107;
my $v = "5.34";
Rmpfr_set_default_prec(1000);
mpfr($v, $e);
print lanx($v, $e); # Wins for 5.34 ** 107 & 5.34 ** 100.
print menezes($v, $e); # Wins for 5.231 ** 107 & 5.231 ** 100.
print by_mul($v, $e); # Wins for 5.35 ** 107 & 5.35 ** 100.
sub mpfr {
my $val = shift;
my $exp = shift;
my $ret = Math::MPFR->new($val);
$ret **= $e;
Rmpfr_out_str($ret, 10, 15, MPFR_RNDN);
printf "\n";
}
sub lanx {
my $val = shift;
$val += 0;
my $exp = shift;
my $ret = 1;
for my $bit (reverse split //,sprintf '%b',$exp) {
$ret *= $val if $bit;
$val **= 2;
}
return $ret;
}
sub menezes {
my $val = shift;
$val += 0;
my $exp = shift;
my $ret = 1;
for my $bit (split //,sprintf '%b',$exp) {
$ret **= 2;
$ret *= $val if $bit;
}
return $ret;
}
sub by_mul {
my $val = shift;
$val += 0;
my $exp = shift;
my $ret = 1;
$ret *= $val for 1..$exp;
return $ret;
}
I could be wrong, of course. (It's happened before ;-)
Cheers, Rob | [reply] [d/l] |