From e2cb1eff0a954a83e0661ea1e7f70a47ed54e893 Mon Sep 17 00:00:00 2001
From: Alberto Ruiz <aruiz@um.es>
Date: Mon, 8 Jun 2015 10:09:39 +0200
Subject: modular C matrix product

---
 packages/base/src/Internal/C/lapack-aux.c | 38 ++++++++++++++-----------------
 1 file changed, 17 insertions(+), 21 deletions(-)

(limited to 'packages/base/src/Internal/C/lapack-aux.c')

diff --git a/packages/base/src/Internal/C/lapack-aux.c b/packages/base/src/Internal/C/lapack-aux.c
index 7da6f6a..1601bef 100644
--- a/packages/base/src/Internal/C/lapack-aux.c
+++ b/packages/base/src/Internal/C/lapack-aux.c
@@ -1290,29 +1290,25 @@ int multiplyQ(int ta, int tb, KQMAT(a),KQMAT(b),QMAT(r)) {
 }
 
 
-int multiplyI(KOIMAT(a), KOIMAT(b), OIMAT(r)) {
-    { TRAV(r,i,j) {
-        int k;
-        AT(r,i,j) = 0;
-        for (k=0;k<ac;k++) {
-            AT(r,i,j) += AT(a,i,k) * AT(b,k,j);
-        }
-      }
+#define MULT_IMP_VER(OP)        \
+    { TRAV(r,i,j) {             \
+        int k;                  \
+        AT(r,i,j) = 0;          \
+        for (k=0;k<ac;k++) {    \
+            OP                  \
+        }                       \
+      }                         \
     }
-    OK
-}
 
-int multiplyL(KOLMAT(a), KOLMAT(b), OLMAT(r)) {
-    { TRAV(r,i,j) {
-        int k;
-        AT(r,i,j) = 0;
-        for (k=0;k<ac;k++) {
-            AT(r,i,j) += AT(a,i,k) * AT(b,k,j);
-        }
-      }
-    }
-    OK
-}
+#define MULT_IMP {                                                                   \
+    if (m==1) {                                                                      \
+        MULT_IMP_VER( AT(r,i,j) += AT(a,i,k) * AT(b,k,j); )                          \
+    } else {                                                                         \
+        MULT_IMP_VER( AT(r,i,j) = (AT(r,i,j) + (AT(a,i,k) * AT(b,k,j)) % m) % m ; )  \
+    } OK }
+
+int multiplyI(int     m, KOIMAT(a), KOIMAT(b), OIMAT(r)) MULT_IMP
+int multiplyL(int32_t m, KOLMAT(a), KOLMAT(b), OLMAT(r)) MULT_IMP
 
 
 ////////////////// sparse matrix-product ///////////////////////////////////////
-- 
cgit v1.2.3