Changeset 4c2cc1 in git for libpolys

Timestamp:

Dec 1, 2016, 4:30:01 PM (7 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

72f5181b0615d9c0f91fee1c8cc55d57b70d9b58

Parents:

83f824ccb8af0c09308f49e50e6acf5ea6c667c7

git-author:

Hans Schoenemann <hannes@mathematik.uni-kl.de>2016-12-01 16:30:01+01:00

git-committer:

Hans Schoenemann <hannes@mathematik.uni-kl.de>2016-12-01 16:33:16+01:00

Message:

add: nrzXExtGcd, nrzQuotRem (rintegers.cc)

File:

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

libpolys/coeffs/rintegers.cc

@@ -129 +129 @@
 }
 
+number  nrzXExtGcd (number a, number b, number *s, number *t, number *u, number *v, const coeffs )
+{
+  mpz_ptr erg = (mpz_ptr) omAllocBin(gmp_nrz_bin);
+  mpz_ptr bs = (mpz_ptr) omAllocBin(gmp_nrz_bin);
+  mpz_ptr bt = (mpz_ptr) omAllocBin(gmp_nrz_bin);
+  mpz_init(erg);
+  mpz_init(bs);
+  mpz_init(bt);
+
+  mpz_gcdext(erg, bs, bt, (mpz_ptr)a, (mpz_ptr)b);
+
+  mpz_ptr bu = (mpz_ptr) omAllocBin(gmp_nrz_bin);
+  mpz_ptr bv = (mpz_ptr) omAllocBin(gmp_nrz_bin);
+
+  mpz_init_set(bu, (mpz_ptr) b);
+  mpz_init_set(bv, (mpz_ptr) a);
+
+  assume(mpz_cmp_si(erg, 0));
+
+  mpz_div(bu, bu, erg);
+  mpz_div(bv, bv, erg);
+
+  mpz_mul_si(bu, bu, -1);
+  *u = (number) bu;
+  *v = (number) bv;
+
+  *s = (number) bs;
+  *t = (number) bt;
+  return (number) erg;
+}
+
 void nrzPower (number a, int i, number * result, const coeffs)
 {
@@ -279 +310 @@
   mpz_tdiv_q(erg, (mpz_ptr) a, (mpz_ptr) b);
   return (number) erg;
+}
+
+number nrzQuotRem (number a, number b, number * r, const coeffs )
+{
+  mpz_ptr qq = (mpz_ptr) omAllocBin(gmp_nrz_bin);
+  mpz_init(qq);
+  mpz_init((mpz_ptr)(*r));
+  mpz_fdiv_qr(qq, (mpz_ptr)(*r), (mpz_ptr) a, (mpz_ptr) b);
+  return (number) qq;
 }
 
@@ -509 +549 @@
   r->cfSize  = nrzSize;
   r->cfInt  = nrzInt;
-  //#ifdef HAVE_RINGS
-  r->cfDivComp = nrzDivComp; // only for ring stuff
-  r->cfIsUnit = nrzIsUnit; // only for ring stuff
-  r->cfGetUnit = nrzGetUnit; // only for ring stuff
-  r->cfExtGcd = nrzExtGcd; // only for ring stuff
-  r->cfDivBy = nrzDivBy; // only for ring stuff
-  //#endif
+  r->cfDivComp = nrzDivComp;
+  r->cfIsUnit = nrzIsUnit;
+  r->cfGetUnit = nrzGetUnit;
+  r->cfExtGcd = nrzExtGcd;
+  r->cfXExtGcd = nrzXExtGcd;
+  r->cfDivBy = nrzDivBy;
   r->cfInpNeg   = nrzNeg;
   r->cfInvers= nrzInvers;