diff options
| author | Kaz Kylheku <kaz@kylheku.com> | 2011-12-10 20:41:06 -0800 |
|---|---|---|
| committer | Kaz Kylheku <kaz@kylheku.com> | 2011-12-10 20:41:06 -0800 |
| commit | 858859cf6fd6becef038cf4d9c1e650c3f2a9cc6 (patch) | |
| tree | 081517aa9628b45d86e7fb1e139e3da996bdc89b /mpi-patches/fix-mult-bug | |
| parent | c940739b8f863f3331f00b36c23ad1526ab1147e (diff) | |
| download | txr-858859cf6fd6becef038cf4d9c1e650c3f2a9cc6.tar.gz txr-858859cf6fd6becef038cf4d9c1e650c3f2a9cc6.tar.bz2 txr-858859cf6fd6becef038cf4d9c1e650c3f2a9cc6.zip | |
Bignum support in mult function.
* arith.c: Regenerated.
* arith.txr (CNUM_BIT): New constant.
(bignum, bignum_dbl_ipt): New static functions.
(@{add-fname}): Use bignum function.
(mul): New functions, rewrite of mul from lib.c.
* lib.c (mul): Function removed.
* mpi-patches/add-mp-set-intptr (mp_set_intptr): Revised patch.
Local variable v should be int_ptr_t not unsigned long.
Also, the mp_set interface doesn't set the sign; it's an unsigned
interface. We must do that ourselves.
* mpi-patches/fix-mult-bug: The main multiplication function is
also broken in the same way, requiring the cast.
* mpi-patches/mpi-set-double-intptr: Fixed use of wrong type for
local variable v.
Diffstat (limited to 'mpi-patches/fix-mult-bug')
| -rw-r--r-- | mpi-patches/fix-mult-bug | 13 |
1 files changed, 11 insertions, 2 deletions
diff --git a/mpi-patches/fix-mult-bug b/mpi-patches/fix-mult-bug index e86d0363..bb8b0f0d 100644 --- a/mpi-patches/fix-mult-bug +++ b/mpi-patches/fix-mult-bug @@ -1,7 +1,7 @@ Index: mpi-1.8.6/mpi.c =================================================================== ---- mpi-1.8.6.orig/mpi.c 2011-12-10 12:05:39.000000000 -0800 -+++ mpi-1.8.6/mpi.c 2011-12-10 12:05:43.000000000 -0800 +--- mpi-1.8.6.orig/mpi.c 2011-12-10 19:41:00.000000000 -0800 ++++ mpi-1.8.6/mpi.c 2011-12-10 19:43:09.000000000 -0800 @@ -3263,7 +3263,7 @@ } @@ -11,3 +11,12 @@ Index: mpi-1.8.6/mpi.c dp[ix] = ACCUM(w); k = CARRYOUT(w); } +@@ -3480,7 +3480,7 @@ + pa = DIGITS(a); + for(jx = 0; jx < ua; ++jx, ++pa) { + pt = pbt + ix + jx; +- w = *pb * *pa + k + *pt; ++ w = *pb * (mp_word) *pa + k + *pt; + *pt = ACCUM(w); + k = CARRYOUT(w); + } |