Changeset 9b8070 in git for libpolys

Timestamp:

Mar 3, 2022, 9:44:34 AM (2 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

64968f0d7471e7e021e850057f1dd83c900ed890

Parents:

fe4754c768f6e7c9aa8e0f5b9875e689a44095d2

Message:

opt: InpAdd, InpMult for Z/n

File:

• libpolys/coeffs/rmodulon.cc (modified) (3 diffs)

libpolys/coeffs/rmodulon.cc

@@ -207 +207 @@
 }
 
+static void nrnInpMult(number &a, number b, const coeffs r)
+{
+  mpz_mul((mpz_ptr)a, (mpz_ptr)a, (mpz_ptr) b);
+  mpz_mod((mpz_ptr)a, (mpz_ptr)a, r->modNumber);
+}
+
 static void nrnPower(number a, int i, number * result, const coeffs r)
 {
@@ -222 +228 @@
   mpz_mod(erg, erg, r->modNumber);
   return (number) erg;
+}
+
+static void nrnInpAdd(number &a, number b, const coeffs r)
+{
+  mpz_add((mpz_ptr)a, (mpz_ptr)a, (mpz_ptr) b);
+  mpz_mod((mpz_ptr)a, (mpz_ptr)a, r->modNumber);
 }
 
@@ -1006 +1018 @@
   r->cfInt         = nrnInt;
   r->cfAdd         = nrnAdd;
+  r->cfInpAdd      = nrnInpAdd;
   r->cfSub         = nrnSub;
   r->cfMult        = nrnMult;
+  r->cfInpMult     = nrnInpMult;
   r->cfDiv         = nrnDiv;
   r->cfAnn         = nrnAnn;