Changes in directory llvm/lib/Support:
APInt.cpp updated: 1.66 -> 1.67
---
Log message:
Use a better algorithm for rounding sqrt results. Change the FIXME about
this to a NOTE: because pari/gp results start to get rounded incorrectly
after 192 bits of precision. APInt and pari/gp never differ by more than
1, but APInt is more accurate because it does not lose precision after 192
bits as does pari/gp.
---
Diffs of the changes: (+12 -8)
APInt.cpp | 20 ++++++++++++--------
1 files changed, 12 insertions(+), 8 deletions(-)
Index: llvm/lib/Support/APInt.cpp
diff -u llvm/lib/Support/APInt.cpp:1.66 llvm/lib/Support/APInt.cpp:1.67
--- llvm/lib/Support/APInt.cpp:1.66 Thu Mar 1 20:01:34 2007
+++ llvm/lib/Support/APInt.cpp Thu Mar 1 22:21:55 2007
@@ -1239,19 +1239,23 @@
}
// Make sure we return the closest approximation
- // FIXME: This still has an off-by-one error in it. Test case:
- // 190 bits: sqrt(694114394047834196220892040454508646882614255319893124270)
=
- // 26346050824513229049493703285 (not 26346050824513229049493703284)
+ // NOTE: The rounding calculation below is correct. It will produce an
+ // off-by-one discrepancy with results from pari/gp. That discrepancy has
been
+ // determined to be a rounding issue with pari/gp as it begins to use a
+ // floating point representation after 192 bits. There are no discrepancies
+ // between this algorithm and pari/gp for bit widths < 192 bits.
APInt square(x_old * x_old);
APInt nextSquare((x_old + 1) * (x_old +1));
if (this->ult(square))
return x_old;
- else if (this->ule(nextSquare))
- if ((nextSquare - *this).ult(*this - square))
- return x_old + 1;
- else
+ else if (this->ule(nextSquare)) {
+ APInt midpoint((nextSquare - square).udiv(two));
+ APInt offset(*this - square);
+ if (offset.ult(midpoint))
return x_old;
- else
+ else
+ return x_old + 1;
+ } else
assert(0 && "Error in APInt::sqrt computation");
return x_old + 1;
}
_______________________________________________
llvm-commits mailing list
llvm-commits@cs.uiuc.edu
http://lists.cs.uiuc.edu/mailman/listinfo/llvm-commits
