Changeset e89e56 in git for libfac/charset


Ignore:
Timestamp:
Mar 17, 2008, 6:44:17 PM (16 years ago)
Author:
Hans Schönemann <hannes@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
85194f264939ff85dd7731d524a483b21f35f5ac
Parents:
ca847b0e390392983ef3461c16980b5defee8e0c
Message:
*hannes: fact.tst


git-svn-id: file:///usr/local/Singular/svn/trunk@10621 2c84dea3-7e68-4137-9b89-c4e89433aadc
Location:
libfac/charset
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • libfac/charset/alg_factor.cc

    rca847b re89e56  
    33// emacs edit mode for this file is -*- C++ -*-
    44////////////////////////////////////////////////////////////
    5 static char * rcsid = "$Id: alg_factor.cc,v 1.20 2008-02-22 12:16:02 Singular Exp $";
     5static char * rcsid = "$Id: alg_factor.cc,v 1.21 2008-03-17 17:44:15 Singular Exp $";
    66////////////////////////////////////////////////////////////
    77// FACTORY - Includes
     
    9494  else { l = 1; }
    9595  d= dr-dv+1;
    96   while ( ( dv <= dr  ) && ( !r.isZero()) )
    97   {
     96  while ( ( dv <= dr  ) && ( r != r.genZero()) ){
    9897    test = power(x,dr-dv)*v*LC(r,x);
    9998    if ( dr == 0 ) { r= CanonicalForm(0); }
     
    142141
    143142  h = CanonicalForm(1);
    144   while ( !G.isZero() )
    145   {
     143  while ( G != G.genZero() ) {
    146144     delta= degree(F,v) -degree(G,v);
    147145     beta = power(CanonicalForm(-1), delta+1) * LC(F,v)* power(h, delta);
     
    196194      temp= gcd(R, R.deriv(vf));
    197195      DEBOUTLN(CERR, "sqrf_norm_sub: temp= ", temp);
    198       if (degree(temp,vf) != 0 || temp.isZero() ){ sqfreetest= 0; }
     196      if (degree(temp,vf) != 0 || temp == temp.genZero() ){ sqfreetest= 0; }
    199197      else { sqfreetest= 1; }
    200198      DEBOUTLN(CERR, "sqrf_norm_sub: sqfreetest= ", sqfreetest);
    201199    }
    202200    else{
    203       DEBOUTMSG(CERR, "Starting isSqrFree(R)!");
    204       // Look at isSqrFree!
     201      DEBOUTMSG(CERR, "Starting SqrFreeTest(R)!");
     202      // Look at SqrFreeTest!
    205203      // (z+a^5+w)^4 with z<w<a should not give sqfreetest=1 !
    206204      // for now we use this workaround with Factorize...
     
    221219      for ( i=testlist; i.hasItem(); i++)
    222220        if ( i.getItem().exp() > 1 && degree(i.getItem().factor(), R.mvar()) > 0) { sqfreetest=0; break; }
    223       DEBOUTLN(CERR, "isSqrFree(R)= ", sqfreetest);
     221      DEBOUTLN(CERR, "SqrFreeTest(R)= ", sqfreetest);
    224222    }
    225223    if ( ! sqfreetest ){
     
    261259      temp= gcd(R, R.deriv(vf));
    262260      DEBOUTLN(CERR, "sqrf_norm_sub: temp= ", temp);
    263       if (degree(temp,vf) != 0 || temp.isZero() ){ sqfreetest= 0; }
     261      if (degree(temp,vf) != 0 || temp == temp.genZero() ){ sqfreetest= 0; }
    264262      else { sqfreetest= 1; }
    265263      DEBOUTLN(CERR, "sqrf_norm_sub: sqfreetest= ", sqfreetest);
    266264    }
    267265    else{
    268       DEBOUTMSG(CERR, "Starting isSqrFree(R)!");
    269       // Look at isSqrFree!
     266      DEBOUTMSG(CERR, "Starting SqrFreeTest(R)!");
     267      // Look at SqrFreeTest!
    270268      // (z+a^5+w)^4 with z<w<a should not give sqfreetest=1 !
    271269      // for now we use this workaround with Factorize...
     
    286284      for ( i=testlist; i.hasItem(); i++)
    287285        if ( i.getItem().exp() > 1 && degree(i.getItem().factor(), R.mvar()) > 0) { sqfreetest=0; break; }
    288       DEBOUTLN(CERR, "isSqrFree(R)= ", sqfreetest);
     286      DEBOUTLN(CERR, "SqrFreeTest(R)= ", sqfreetest);
    289287    }
    290288    if ( ! sqfreetest ){
     
    357355
    358356  if ( Astar.length() == 0 ) return 0;
    359   for ( CFListIterator i=Astar; i.hasItem(); i++)
    360   {
     357  for ( CFListIterator i=Astar; i.hasItem(); i++){
    361358    elem= i.getItem();
    362     if ( elem.deriv().isZero() ) return Counter;
     359    if ( elem.deriv() == elem.genZero() ) return Counter;
    363360    else Counter += 1;
    364361  }
     
    757754    Fgcd= pp(Fgcd); Ggcd= pp(Ggcd);
    758755    DEBDECLEVEL(CERR,"newfactoras");
    759     return UnionCFFL(newfactoras(Fgcd,as,success) , newfactoras(Ggcd,as,success));
     756    return myUnion(newfactoras(Fgcd,as,success) , newfactoras(Ggcd,as,success));
    760757  }
    761758  if ( getCharacteristic() > 0 ){
     
    826823    output=newfactoras(i.getItem().factor(),as, success);
    827824    for ( CFFListIterator j=output; j.hasItem(); j++)
    828       Output = appendCFFL(Output,CFFactor(j.getItem().factor(),j.getItem().exp()*i.getItem().exp()));
     825      Output = myappend(Output,CFFactor(j.getItem().factor(),j.getItem().exp()*i.getItem().exp()));
    829826  }
    830827  return Output;
     
    833830/*
    834831$Log: not supported by cvs2svn $
    835 Revision 1.19  2008/01/25 14:19:39  Singular
    836 *hannes: SqrFreeTest -> isSqrFree
    837 
    838 Revision 1.18  2008/01/22 09:51:36  Singular
    839 *hannes: sqrFree/InternalSqrFree -> factory
    840 
    841832Revision 1.17  2007/05/15 14:46:48  Singular
    842833*hannes: factorize in Zp(a)[x...]
  • libfac/charset/algfactor.cc

    rca847b re89e56  
    33// emacs edit mode for this file is -*- C++ -*-
    44////////////////////////////////////////////////////////////
    5 static char * rcsid = "$Id: algfactor.cc,v 1.9 2008-01-22 09:51:37 Singular Exp $";
     5static char * rcsid = "$Id: algfactor.cc,v 1.10 2008-03-17 17:44:16 Singular Exp $";
    66////////////////////////////////////////////////////////////
    77// FACTORY - Includes
     
    429429    success= min(success,csuccess);
    430430    for ( CFFListIterator j=CERR; j.hasItem(); j++)
    431       Output = appendCFFL(Output,CFFactor(j.getItem().factor(),j.getItem().exp()*i.getItem().exp()));
     431      Output = myappend(Output,CFFactor(j.getItem().factor(),j.getItem().exp()*i.getItem().exp()));
    432432  }
    433433  return Output;
Note: See TracChangeset for help on using the changeset viewer.