summaryrefslogtreecommitdiff
path: root/src/libbiguint/bu_mul.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/libbiguint/bu_mul.c')
-rw-r--r--src/libbiguint/bu_mul.c32
1 files changed, 0 insertions, 32 deletions
diff --git a/src/libbiguint/bu_mul.c b/src/libbiguint/bu_mul.c
deleted file mode 100644
index 16043fd..0000000
--- a/src/libbiguint/bu_mul.c
+++ /dev/null
@@ -1,32 +0,0 @@
-/* ISC license. */
-
-#include <stdint.h>
-#include <skalibs/uint64.h>
-#include <skalibs/biguint.h>
-
- /* No Karatsuba. Keep it simple, stupid. */
-
-int bu_mul (uint32_t *x, unsigned int xn, uint32_t const *a, unsigned int an, uint32_t const *b, unsigned int bn)
-{
- unsigned int alen = bu_len(a, an) ;
- unsigned int blen = bu_len(b, bn) ;
- uint32_t c[alen + blen] ;
- unsigned int i = 0 ;
- bu_zero(c, alen + blen) ;
- for (; i < alen ; i++)
- {
- uint32_t carry = 0 ;
- unsigned int j = 0 ;
- for (; j < blen ; j++)
- {
- uint64_t t = a[i] ;
- t *= b[j] ;
- t += c[i+j] ;
- t += carry ;
- c[i+j] = (uint32_t)t ;
- carry = (uint32_t)(t >> 32) ;
- }
- c[i+j] += carry ;
- }
- return bu_copy(x, xn, c, alen+blen) ;
-}