summaryrefslogtreecommitdiff
path: root/sys/src/libmp/port/mpvecadd.c
diff options
context:
space:
mode:
authorTaru Karttunen <taruti@taruti.net>2011-03-30 15:46:40 +0300
committerTaru Karttunen <taruti@taruti.net>2011-03-30 15:46:40 +0300
commite5888a1ffdae813d7575f5fb02275c6bb07e5199 (patch)
treed8d51eac403f07814b9e936eed0c9a79195e2450 /sys/src/libmp/port/mpvecadd.c
Import sources from 2011-03-30 iso image
Diffstat (limited to 'sys/src/libmp/port/mpvecadd.c')
-rwxr-xr-xsys/src/libmp/port/mpvecadd.c35
1 files changed, 35 insertions, 0 deletions
diff --git a/sys/src/libmp/port/mpvecadd.c b/sys/src/libmp/port/mpvecadd.c
new file mode 100755
index 000000000..98fdcc901
--- /dev/null
+++ b/sys/src/libmp/port/mpvecadd.c
@@ -0,0 +1,35 @@
+#include "os.h"
+#include <mp.h>
+#include "dat.h"
+
+// prereq: alen >= blen, sum has at least blen+1 digits
+void
+mpvecadd(mpdigit *a, int alen, mpdigit *b, int blen, mpdigit *sum)
+{
+ int i, carry;
+ mpdigit x, y;
+
+ carry = 0;
+ for(i = 0; i < blen; i++){
+ x = *a++;
+ y = *b++;
+ x += carry;
+ if(x < carry)
+ carry = 1;
+ else
+ carry = 0;
+ x += y;
+ if(x < y)
+ carry++;
+ *sum++ = x;
+ }
+ for(; i < alen; i++){
+ x = *a++ + carry;
+ if(x < carry)
+ carry = 1;
+ else
+ carry = 0;
+ *sum++ = x;
+ }
+ *sum = carry;
+}