Changeset dd8668 in git for kernel/longalg.cc


Ignore:
Timestamp:
Oct 11, 2009, 4:19:08 PM (14 years ago)
Author:
Hans Schönemann <hannes@…>
Branches:
(u'spielwiese', 'd0474371d8c5d8068ab70bfb42719c97936b18a6')
Children:
1b83ffaae1ee686f47a1448856f99c2cff0b66dd
Parents:
05998dae5a4237480bafa680febec202b19a3bf6
Message:
*hannes: less Normalize


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

Legend:

Unmodified
Added
Removed

  • kernel/longalg.cc

    r05998da rdd8668  
    22*  Computer Algebra System SINGULAR     *
    33****************************************/
    4 /* $Id: longalg.cc,v 1.62 2009-10-10 17:15:40 Singular Exp $ */
     4/* $Id: longalg.cc,v 1.63 2009-10-11 14:19:08 Singular Exp $ */
    55/*
    66* ABSTRACT:   algebraic numbers
    … …  
    6767static int napExpi(int i, napoly a, napoly b);
    6868       ring nacRing;
     69
     70void naCoefNormalize(number pp);
    6971
    7072#define napCopy(p)       p_Copy(p,nacRing)
    … …  
    947949{
    948950  lnumber l=(lnumber)n;
    949   naNormalize(n);
    950951  if ((l!=NULL)&&(l->n==NULL)&&(p_IsConstant(l->z,r->algring)))
    951952  {
    … …  
    978979  if (p==NULL) return NULL;
    979980  naTest(p);
    980   //naNormalize(p);
    981981  lnumber erg;
    982982  lnumber src = (lnumber)p;
    … …  
    10561056  {
    10571057     number luu=(number)lu;
    1058      naNormalize(luu);
     1058     //if (p_IsConstant(lu->n,nacRing)) naCoefNormalize(luu);
     1059     //else           
     1060                naNormalize(luu);
    10591061     lu=(lnumber)luu;
    10601062  }
    … …  
    11061108    else            x = NULL;
    11071109  }
    1108   //if (x!=NULL)
    1109   //{
    1110   //  if (p_LmIsConstant(x,nacRing))
    1111   //  {
    1112   //    number inv=nacInvers(pGetCoeff(x));
    1113   //    napMultN(lu->z,inv);
    1114   //    n_Delete(&inv,nacRing);
    1115   //    napDelete(&x);
    1116   //  }
    1117   //}
    11181110  lu->n = x;
    11191111  lu->s = FALSE;
    … …  
    11211113  {
    11221114     number luu=(number)lu;
    1123      naNormalize(luu);
     1115     //if (p_IsConstant(lu->n,nacRing)) naCoefNormalize(luu);
     1116     //else   
     1117                         naNormalize(luu);
    11241118     lu=(lnumber)luu;
    11251119  }
    … …  
    11941188    lo=NULL;
    11951189  }
    1196 #if 1
    11971190  else if (lo->n!=NULL)
    11981191  {
    11991192    lo->s = 0;
    12001193    number luu=(number)lo;
    1201     naNormalize(luu);
     1194    // if (p_IsConstant(lo->n,nacRing)) naCoefNormalize(luu);
     1195    // else     
     1196                      naNormalize(luu);
    12021197    lo=(lnumber)luu;
    12031198  }
    1204   else
    1205     lo->s=3;
    1206 #endif
    12071199  naTest((number)lo);
    12081200  return (number)lo;
    … …  
    12231215    return NULL;
    12241216  }
    1225   naNormalize(la);
    12261217  assume(a->z!=NULL && b->z!=NULL);
    12271218  assume(a->n==NULL && b->n==NULL);
    … …  
    12931284    lo->s = 0;
    12941285    number luu=(number)lo;
    1295     naNormalize(luu);
     1286     //if (p_IsConstant(lo->n,nacRing)) naCoefNormalize(luu);
     1287     //else   
     1288                         naNormalize(luu);
    12961289    lo=(lnumber)luu;
    12971290  }
    1298   else
    1299     lo->s=3;
    13001291  naTest((number)lo);
    13011292  return (number)lo;
    … …  
    13541345    lo->z = x;
    13551346    lo->n = NULL;
    1356     lo->s = 2;
    13571347    while (x!=NULL)
    13581348    {
    … …  
    13661356  {
    13671357     number luu=(number)lo;
    1368      naNormalize(luu);
     1358     //if (p_IsConstant(lo->n,nacRing)) naCoefNormalize(luu);
     1359     //else 
     1360                           naNormalize(luu);
    13691361     lo=(lnumber)luu;
    13701362  }
    … …  
    15341526  {
    15351527    phn->s = 0;
    1536     //naNormalize(phn);
    15371528    BOOLEAN has_denom=(ph->n!=NULL);
    15381529    napWrite(ph->z,has_denom/*(ph->n!=NULL)*/);
    … …  
    15971588  a->z = p_ISet(1,nacRing);
    15981589  a->n = NULL;
    1599   a->s = 2;
    16001590  return TRUE;
    16011591#else
    … …  
    17991789    n_Delete(&nz,nacRing);
    18001790  }
    1801 
    1802   if ((p->n!=NULL)
    1803   && (p_IsConstant(p->n,nacRing))
    1804   && (n_IsOne(pGetCoeff(p->n),nacRing)))
    1805   {
    1806     p_Delete(&(p->n), nacRing);
    1807     p->n = NULL;
     1791  if (p->n!=NULL)
     1792  {
     1793    if(!nacGreaterZero(pGetCoeff(p->n)))
     1794    {
     1795      p->z=napNeg(p->z);
     1796      p->n=napNeg(p->n);
     1797    }
     1798
     1799    if (/*(p->n!=NULL) && */
     1800    (p_IsConstant(p->n,nacRing))
     1801    && (n_IsOne(pGetCoeff(p->n),nacRing)))
     1802    {
     1803      p_Delete(&(p->n), nacRing);
     1804      p->n = NULL;
     1805    }
    18081806  }
    18091807}
    … …  
    18141812  lnumber p = (lnumber)pp;
    18151813
    1816   if ((p==NULL) /*|| (p->s==2)*/)
     1814  if (p==NULL)
    18171815    return;
    18181816  naCoefNormalize(pp);
Note: See TracChangeset for help on using the changeset viewer.