Changeset f55f39 in git for factory/facFactorize.cc

Timestamp:

Jul 31, 2012, 3:31:51 PM (11 years ago)

Author:

Martin Lee <martinlee84@…>

Branches:

(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', '2234726c50d679d6664181a5c72f75a6fd64a787')

Children:

97686ecde8614d72ee9d26e96dc9344df8a7da32

Parents:

3ace5b6cacf35176c78dd144d681a504076ac119

git-author:

Martin Lee <martinlee84@web.de>2012-07-31 15:31:51+02:00

git-committer:

Martin Lee <martinlee84@web.de>2012-09-04 17:25:37+02:00

Message:

chg: multiply by common denominator

File:

• factory/facFactorize.cc (modified) (2 diffs)

factory/facFactorize.cc

@@ -658 +658 @@
   }
 
+  A *= bCommonDen (A);
   // check main variable
   CFList Aeval, list, evaluation, bufEvaluation, bufAeval;
@@ -946 +947 @@
   bool noOneToOne= false;
 
+  CFList commonDenominators;
+  for (CFListIterator iter=biFactors; iter.hasItem(); iter++)
+    commonDenominators.append (bCommonDen (iter.getItem()));
+  CanonicalForm tmp1, tmp2, tmp3=1;
+  CFListIterator iter1, iter2;
+  for (int i= 0; i < A.level() - 2; i++)
+  {
+    iter2= commonDenominators;
+    for (iter1= leadingCoeffs2[i]; iter1.hasItem(); iter1++, iter2++)
+    {
+      tmp1= bCommonDen (iter1.getItem());
+      Off (SW_RATIONAL);
+      iter2.getItem()= lcm (iter2.getItem(), tmp1);
+      On (SW_RATIONAL);
+    }
+  }
+  tmp1= prod (commonDenominators);
+  for (iter1= Aeval; iter1.hasItem(); iter1++)
+  {
+    tmp2= bCommonDen (iter1.getItem());
+    Off (SW_RATIONAL);
+    tmp3= lcm (tmp2,tmp3);
+    On (SW_RATIONAL);
+  }
+  CanonicalForm multiplier;
+  multiplier= tmp3/tmp1;
+  iter2= commonDenominators;
+  for (iter1=biFactors; iter1.hasItem(); iter1++, iter2++)
+    iter1.getItem() *= iter2.getItem()*multiplier;
+
+  for (iter1= Aeval; iter1.hasItem(); iter1++)
+    iter1.getItem() *= tmp3*power (multiplier, biFactors.length() - 1);
+
+  for (int i= 0; i < A.level() - 2; i++)
+  {
+    iter2= commonDenominators;
+    for (iter1= leadingCoeffs2[i]; iter1.hasItem(); iter1++, iter2++)
+      iter1.getItem() *= iter2.getItem()*multiplier;
+  }
+
+
   factors= nonMonicHenselLift (Aeval, biFactors, leadingCoeffs2, diophant,
                                Pi, liftBounds, liftBoundsLength, noOneToOne);