* [Ada] Minor fixes mostly in comments of runtime arithmetic unit
@ 2019-09-17 8:07 Pierre-Marie de Rodat
0 siblings, 0 replies; only message in thread
From: Pierre-Marie de Rodat @ 2019-09-17 8:07 UTC (permalink / raw)
To: gcc-patches; +Cc: Yannick Moy
[-- Attachment #1: Type: text/plain, Size: 804 bytes --]
Multiple comments in functions Double_Divide and Scaled_Divide were
incorrect. Now fixed.
Also change the expression (if Zhi /= 0 then Ylo * Zhi else 0) to the
simpler equivalent (Ylo * Zhi) in Double_Divide.
Also add a comment explaining why the implementation of Algorithm D for
multiple-precision division from the 2nd Edition of The Art of Computer
Programming does not suffer from two bugs discovered on that version.
There is no impact on execution, hence no test.
Tested on x86_64-pc-linux-gnu, committed on trunk
2019-09-17 Yannick Moy <moy@adacore.com>
gcc/ada/
* libgnat/s-arit64.adb (Double_Divide): Simplify needlessly
complex computation. Fix comments.
(Scaled_Divide): Fix comments. Explain why implementation does
not suffer from bugs in Algorithm D from 2nd Edition of TAOCP.
[-- Attachment #2: patch.diff --]
[-- Type: text/x-diff, Size: 2976 bytes --]
--- gcc/ada/libgnat/s-arit64.adb
+++ gcc/ada/libgnat/s-arit64.adb
@@ -161,7 +161,7 @@ package body System.Arith_64 is
end if;
else
- T2 := (if Zhi /= 0 then Ylo * Zhi else 0);
+ T2 := Ylo * Zhi;
end if;
T1 := Ylo * Zlo;
@@ -179,7 +179,7 @@ package body System.Arith_64 is
Den_Pos := (Y < 0) = (Z < 0);
- -- Check overflow case of largest negative number divided by 1
+ -- Check overflow case of largest negative number divided by -1
if X = Int64'First and then Du = 1 and then not Den_Pos then
Raise_Error;
@@ -404,15 +404,16 @@ package body System.Arith_64 is
Ru := T2 rem Zlo;
end if;
- -- If divisor is double digit and too large, raise error
+ -- If divisor is double digit and dividend is too large, raise error
elsif (D (1) & D (2)) >= Zu then
Raise_Error;
-- This is the complex case where we definitely have a double digit
-- divisor and a dividend of at least three digits. We use the classical
- -- multiple division algorithm (see section (4.3.1) of Knuth's "The Art
- -- of Computer Programming", Vol. 2 for a description (algorithm D).
+ -- multiple-precision division algorithm (see section (4.3.1) of Knuth's
+ -- "The Art of Computer Programming", Vol. 2 for a description
+ -- (algorithm D).
else
-- First normalize the divisor so that it has the leading bit on.
@@ -450,7 +451,7 @@ package body System.Arith_64 is
-- Note that when we scale up the dividend, it still fits in four
-- digits, since we already tested for overflow, and scaling does
- -- not change the invariant that (D (1) & D (2)) >= Zu.
+ -- not change the invariant that (D (1) & D (2)) < Zu.
T1 := Shift_Left (D (1) & D (2), Scale);
D (1) := Hi (T1);
@@ -485,6 +486,19 @@ package body System.Arith_64 is
-- Adjust quotient digit if it was too high
+ -- We use the version of the algorithm in the 2nd Edition of
+ -- "The Art of Computer Programming". This had a bug not
+ -- discovered till 1995, see Vol 2 errata:
+ -- http://www-cs-faculty.stanford.edu/~uno/err2-2e.ps.gz.
+ -- Under rare circumstances the expression in the test could
+ -- overflow. This version was further corrected in 2005, see
+ -- Vol 2 errata:
+ -- http://www-cs-faculty.stanford.edu/~uno/all2-pre.ps.gz.
+ -- This implementation is not impacted by these bugs, due to the
+ -- use of a word-size comparison done in function Le3 instead of
+ -- a comparison on two-word integer quantities in the original
+ -- algorithm.
+
loop
exit when Le3 (S1, S2, S3, D (J + 1), D (J + 2), D (J + 3));
Qd (J + 1) := Qd (J + 1) - 1;
^ permalink raw reply [flat|nested] only message in thread
only message in thread, other threads:[~2019-09-17 8:06 UTC | newest]
Thread overview: (only message) (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2019-09-17 8:07 [Ada] Minor fixes mostly in comments of runtime arithmetic unit Pierre-Marie de Rodat
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).