Changeset d1d27f1 in git for libpolys/coeffs/longrat.cc

Timestamp:

Aug 3, 2012, 7:33:11 PM (12 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

6993c8371c2a089de1e05f7e0a7e9ab0f8c76145

Parents:

69c542f9869e7a3e82b1cc481139d89110b7afd3

Message:

add: nlClearContent

File:

• libpolys/coeffs/longrat.cc (modified) (2 diffs)

libpolys/coeffs/longrat.cc

@@ -2660 +2660 @@
   } while (numberCollectionEnumerator.MoveNext() ); 
  
+  cand=nlCopy(cand,cf);
   // part 2: compute gcd(cand,all coeffs)
   numberCollectionEnumerator.Reset();
@@ -2683 +2684 @@
 }
 
-static void nlClearDenominators(ICoeffsEnumerator& /*numberCollectionEnumerator*/, number& c, const coeffs cf)
+static void nlClearDenominators(ICoeffsEnumerator& numberCollectionEnumerator, number& c, const coeffs cf)
 {
   assume(cf != NULL);
   assume(getCoeffType(cf) == ID);
-  // all coeffs are given by integers!!!
-
-  c = n_Init(1, cf);
-  assume(FALSE); // TODO: NOT YET IMPLEMENTED!!!
+  // all coeffs are given by integers after returning from this routine
+
+  // part 1, collect product of all denominators /gcds
+  number cand,cand1;
+  cand=ALLOC_RNUMBER();
+  #if defined(LDEBUG)
+  cand->debug=123456;
+  #endif
+  cand->s=3;
+
+  int s=0;
+  mpz_t tmp;
+  mpz_init(tmp);
+  numberCollectionEnumerator.Reset();
+  do
+  {
+    if (!(SR_HDL(cand1)&SR_INT))
+    {
+      nlNormalize(numberCollectionEnumerator.Current(),cf);
+      cand1= numberCollectionEnumerator.Current();
+      if ((!(SR_HDL(cand1)&SR_INT)) // not a short int
+      && (cand1->s==1))             // and is rational
+      {
+        if (s==0) // first denom, we meet
+        {
+          mpz_init_set(cand->z,cand1->n); 
+          s=1; 
+        }
+        else // we have already something
+        {
+          mpz_gcd(tmp,cand->z,cand1->n);
+          if (mpz_cmp_si(tmp,1)!=0)
+          {
+            mpz_div(cand->z,cand->z,tmp);
+          }
+          mpz_mul(cand->z,cand->z,cand1->n);
+        }
+      }
+    }
+  } while (numberCollectionEnumerator.MoveNext() ); 
+
+  if (s==0) // nothing to do, all coeffs are already integers
+  {
+    mpz_clear(tmp);
+    nlDelete(&cand,cf); 
+    c=nlInit(1,cf); 
+    return;
+  }
+  cand=nlShort3(cand);
+ 
+  // part2: all coeffs = all coeffs * cand
+  c=cand;
+  numberCollectionEnumerator.Reset();
+  do
+  {
+    number t=nlMult(numberCollectionEnumerator.Current(),cand,cf);
+    nlDelete(&numberCollectionEnumerator.Current(),cf);
+    nlNormalize(t,cf);
+    numberCollectionEnumerator.Current()=t;
+  } while (numberCollectionEnumerator.MoveNext() ); 
+ 
 }