Changeset bac8611 in git for kernel/rmodulon.cc

Timestamp:

Feb 1, 2008, 4:11:34 PM (16 years ago)

Author:

Oliver Wienand <wienand@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

befecbc5f2fd0929fe1c00f9ed52cc5c128a3612

Parents:

3901ebf059dc5ef7abe8c2726771e6d8edee6d96

Message:

ring.h:
no simple alloc for modn

kutil.cc:
delete generated coeffs in pair creation

number.cc:
copy, delete, size for rmodulon and rintegers

kspoly.cc:
pSetCoeff durch pSetCoeff0 ersetzt

rint*, rmod*: delete, etc.


git-svn-id: file:///usr/local/Singular/svn/trunk@10550 2c84dea3-7e68-4137-9b89-c4e89433aadc

File:

• kernel/rmodulon.cc (modified) (17 diffs)

kernel/rmodulon.cc

@@ -2 +2 @@
 *  Computer Algebra System SINGULAR     *
 ****************************************/
-/* $Id: rmodulon.cc,v 1.14 2008-02-01 10:31:44 wienand Exp $ */
+/* $Id: rmodulon.cc,v 1.15 2008-02-01 15:11:34 wienand Exp $ */
 /*
 * ABSTRACT: numbers modulo n
@@ -23 +23 @@
 
 typedef MP_INT *int_number;
-omBin gmp_bin = omGetSpecBin(sizeof(MP_INT));
+omBin gmp_nrn_bin = omGetSpecBin(sizeof(MP_INT));
 
 int_number nrnModul = NULL;
@@ -35 +35 @@
 number nrnInit (int i)
 {
-  int_number erg = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
+  int_number erg = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
   mpz_init_set_si(erg, i);
   mpz_mod(erg, erg, nrnModul);
@@ -44 +44 @@
 {
   mpz_clear((int_number) *a);
-  omFreeBin((ADDRESS) *a, gmp_bin);
+  omFreeBin((ADDRESS) *a, gmp_nrn_bin);
+  *a = NULL;
+}
+
+number nrnCopy(number a)
+{
+  int_number erg = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
+  mpz_init_set(erg, (int_number) a);
+  return (number) erg;
+}
+
+number cfrnCopy(number a, const ring r)
+{
+  return nrnCopy(a);
+}
+
+int nrnSize(number a)
+{
+  if (a == NULL) return 0;
+  return sizeof(MP_INT);
 }
 
@@ -60 +79 @@
 number nrnMult (number a, number b)
 {
-  int_number erg = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
+  int_number erg = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
   mpz_init(erg);
   mpz_mul(erg, (int_number) a, (int_number) b);
@@ -69 +88 @@
 void nrnPower (number a, int i, number * result)
 {
-  int_number erg = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
+  int_number erg = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
   mpz_init(erg);
   mpz_powm_ui(erg, (int_number) a, i, nrnModul);
@@ -77 +96 @@
 number nrnAdd (number a, number b)
 {
-  int_number erg = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
+  int_number erg = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
   mpz_init(erg);
   mpz_add(erg, (int_number) a, (int_number) b);
@@ -86 +105 @@
 number nrnSub (number a, number b)
 {
-  int_number erg = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
+  int_number erg = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
   mpz_init(erg);
   mpz_sub(erg, (int_number) a, (int_number) b);
@@ -95 +114 @@
 number nrnNeg (number c)
 {
-  int_number erg = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
+  int_number erg = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
   mpz_init(erg);
   mpz_sub(erg, nrnModul, (int_number) c);
@@ -103 +122 @@
 number  nrnInvers (number c)
 {
-  int_number erg = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
+  int_number erg = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
   mpz_init(erg);
   mpz_invert(erg, (int_number) c, nrnModul);
@@ -129 +148 @@
 {
   if ((a == NULL) && (b == NULL)) return nrnInit(0);
-  int_number erg = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
+  int_number erg = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
   mpz_init_set(erg, nrnModul);
   if (a != NULL) mpz_gcd(erg, erg, (int_number) a);
@@ -139 +158 @@
 number nrnGcd3 (number a,number b, number c,ring r)
 {
-  int_number erg = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
+  int_number erg = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
   mpz_init(erg);
   if (a == NULL) a = (number) nrnModul;
@@ -157 +176 @@
 number  nrnExtGcd (number a, number b, number *s, number *t)
 {
-  int_number erg = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
-  int_number bs = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
-  int_number bt = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
+  int_number erg = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
+  int_number bs = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
+  int_number bt = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
   mpz_init(erg);
   mpz_init(bs);
@@ -264 +283 @@
 {
   if (a == NULL) a = (number) nrnModul;
-  int_number erg = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
+  int_number erg = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
   mpz_init(erg);
   if (mpz_divisible_p((int_number) a, (int_number) b))
@@ -296 +315 @@
 number nrnIntDiv (number a,number b)
 {
-  int_number erg = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
+  int_number erg = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
   mpz_init(erg);
   if (a == NULL) a = (number) nrnModul;
@@ -330 +349 @@
     nrnDelete((number*) &nrnMinusOne, NULL);
   }
-  nrnModul = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
+  nrnModul = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
   mpz_init(nrnModul);
   mpz_set_ull(nrnModul, nrnBase);
   mpz_pow_ui(nrnModul, nrnModul, nrnExponent);
 
-  nrnMinusOne = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
+  nrnMinusOne = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
   mpz_init(nrnMinusOne);
   mpz_sub_ui(nrnMinusOne, nrnModul, 1);
@@ -408 +427 @@
 char * nrnRead (char *s, number *a)
 {
-  int_number z = (int_number) omAllocBin(gmp_bin); // evtl. spaeter mit bin
+  int_number z = (int_number) omAllocBin(gmp_nrn_bin); // evtl. spaeter mit bin
   {
     s = nlCPEatLongC(s, z);