Changeset 50a2aa9 in git


Ignore:
Timestamp:
Jun 21, 2011, 12:42:12 PM (13 years ago)
Author:
Martin Lee <martinlee84@…>
Branches:
(u'spielwiese', '5b153614cbc72bfa198d75b1e9e33dab2645d9fe')
Children:
64e7cb73f028e3866617e41d55fd0265168a395d
Parents:
8c1a84580fcaf92907fa2c80e52eec3de2d8a5c2
Message:
resultantZ needs polys over Z


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

Legend:

Unmodified
Added
Removed
  • factory/facAlgExt.cc

    r8c1a84 r50a2aa9  
    5050  CanonicalForm mipo= getMipo (alpha);
    5151  mipo= mipo (x, alpha);
     52  mipo *= bCommonDen (mipo);
    5253
    5354  int degg= degree (g);
     
    7273      {
    7374        g= F (y - i*alpha, y);
     75        g *= bCommonDen (g);
    7476        if (degg >= 8 || degmipo >= 8)
    7577          norm= resultantZ (g (x, alpha), mipo, x);
     
    8082      {
    8183        g= F (y + i*alpha, y);
     84        g *= bCommonDen (g);
    8285        if (degg >= 8 || degmipo >= 8)
    8386          norm= resultantZ (g (x, alpha), mipo, x);
     
    109112  CanonicalForm norm= sqrfNorm (f, alpha, shift);
    110113  ASSERT (degree (norm, alpha) <= 0, "wrong norm computed");
    111   CFFList normFactors= factorize (norm); //maybe compute norm (not necessarily squarefree), split f and procede recursively
     114  CFFList normFactors= factorize (norm);
    112115  CFList factors;
    113116  if (normFactors.length() <= 2)
     
    151154  CFList sqrfFactors= AlgExtSqrfFactorize (sqrf, alpha);
    152155
    153   if (!isOn (SW_RATIONAL))
    154     On (SW_RATIONAL);
     156  bool save_rat=!isOn (SW_RATIONAL);
     157  On (SW_RATIONAL);
    155158  CanonicalForm buf= F/Lc (F);
    156159  CFFList factors;
     
    169172  factors.insert (CFFactor (Lc(F), 1));
    170173  ASSERT (degree (buf) <= 0, "bug in AlgExtFactorize");
     174  if (save_rat) Off(SW_RATIONAL);
    171175  return factors;
    172176}
Note: See TracChangeset for help on using the changeset viewer.