Changeset 6f801b in git
- Timestamp:
- Feb 22, 2008, 1:16:03 PM (15 years ago)
- Branches:
- (u'spielwiese', '0d6b7fcd9813a1ca1ed4220cfa2b104b97a0a003')
- Children:
- 9667879f6c3f848cdfbb5c6c0e9a14ea6e130e81
- Parents:
- 8078c8a05a777673a09af83fb19f5f3e9c1f24ff
- Location:
- libfac
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
libfac/charset/alg_factor.cc
r8078c8 r6f801b 3 3 // emacs edit mode for this file is -*- C++ -*- 4 4 //////////////////////////////////////////////////////////// 5 static char * rcsid = "$Id: alg_factor.cc,v 1. 19 2008-01-25 14:19:39Singular Exp $";5 static char * rcsid = "$Id: alg_factor.cc,v 1.20 2008-02-22 12:16:02 Singular Exp $"; 6 6 //////////////////////////////////////////////////////////// 7 7 // FACTORY - Includes … … 94 94 else { l = 1; } 95 95 d= dr-dv+1; 96 while ( ( dv <= dr ) && ( r != r.genZero()) ){ 96 while ( ( dv <= dr ) && ( !r.isZero()) ) 97 { 97 98 test = power(x,dr-dv)*v*LC(r,x); 98 99 if ( dr == 0 ) { r= CanonicalForm(0); } … … 141 142 142 143 h = CanonicalForm(1); 143 while ( G != G.genZero() ) { 144 while ( !G.isZero() ) 145 { 144 146 delta= degree(F,v) -degree(G,v); 145 147 beta = power(CanonicalForm(-1), delta+1) * LC(F,v)* power(h, delta); … … 194 196 temp= gcd(R, R.deriv(vf)); 195 197 DEBOUTLN(CERR, "sqrf_norm_sub: temp= ", temp); 196 if (degree(temp,vf) != 0 || temp == temp.genZero() ){ sqfreetest= 0; }198 if (degree(temp,vf) != 0 || temp.isZero() ){ sqfreetest= 0; } 197 199 else { sqfreetest= 1; } 198 200 DEBOUTLN(CERR, "sqrf_norm_sub: sqfreetest= ", sqfreetest); … … 259 261 temp= gcd(R, R.deriv(vf)); 260 262 DEBOUTLN(CERR, "sqrf_norm_sub: temp= ", temp); 261 if (degree(temp,vf) != 0 || temp == temp.genZero() ){ sqfreetest= 0; }263 if (degree(temp,vf) != 0 || temp.isZero() ){ sqfreetest= 0; } 262 264 else { sqfreetest= 1; } 263 265 DEBOUTLN(CERR, "sqrf_norm_sub: sqfreetest= ", sqfreetest); … … 355 357 356 358 if ( Astar.length() == 0 ) return 0; 357 for ( CFListIterator i=Astar; i.hasItem(); i++){ 359 for ( CFListIterator i=Astar; i.hasItem(); i++) 360 { 358 361 elem= i.getItem(); 359 if ( elem.deriv() == elem.genZero() ) return Counter;362 if ( elem.deriv().isZero() ) return Counter; 360 363 else Counter += 1; 361 364 } … … 830 833 /* 831 834 $Log: not supported by cvs2svn $ 835 Revision 1.19 2008/01/25 14:19:39 Singular 836 *hannes: SqrFreeTest -> isSqrFree 837 832 838 Revision 1.18 2008/01/22 09:51:36 Singular 833 839 *hannes: sqrFree/InternalSqrFree -> factory -
libfac/charset/csutil.cc
r8078c8 r6f801b 2 2 //////////////////////////////////////////////////////////// 3 3 // emacs edit mode for this file is -*- C++ -*- 4 static char * rcsid = "$Id: csutil.cc,v 1.1 5 2007-10-16 15:49:49Singular Exp $";4 static char * rcsid = "$Id: csutil.cc,v 1.16 2008-02-22 12:16:02 Singular Exp $"; 5 5 ///////////////////////////////////////////////////////////// 6 6 // FACTORY - Includes … … 138 138 if (dg <= df) {l=LC(gg); gg = gg -LC(gg)*power(v,dg);} 139 139 else { l = 1; } 140 while ( ( dg <= df ) && ( ff != ff.genZero()) )140 while ( ( dg <= df ) && ( !ff.isZero()) ) 141 141 { 142 142 // CERR << "Start gcd..." << "\n"; … … 195 195 else { l = 1; } 196 196 n= 0; 197 while ( ( dg <= df ) && ( ff != ff.genZero()) )197 while ( ( dg <= df ) && ( !ff.isZero()) ) 198 198 { 199 199 test= power(v,df-dg) * gg * LC(ff); … … 372 372 { 373 373 elem= premasb(i.getItem(),as); 374 if ( elem != elem.genZero() ) output.append(elem);374 if ( !elem.isZero() ) output.append(elem); 375 375 } 376 376 return output; … … 457 457 { 458 458 test = mydivremt(r,testelem,a,b); 459 if ( test && b == r.genZero() ) r = a;459 if ( test && b.isZero() ) r = a; 460 460 else break; 461 461 } … … 472 472 { 473 473 test = divremt(r,testelem,a,b); 474 if ( test && b == r.genZero() )474 if ( test && b.isZero() ) 475 475 { 476 476 Remembern.FS2= Union(Remembern.FS2, CFList(testelem)); … … 988 988 /* 989 989 $Log: not supported by cvs2svn $ 990 Revision 1.15 2007/10/16 15:49:49 Singular 991 SAGE: remove alg content in myfitting 992 990 993 Revision 1.14 2007/05/15 14:46:48 Singular 991 994 *hannes: factorize in Zp(a)[x...] -
libfac/factor/Factor.cc
r8078c8 r6f801b 1 1 /* Copyright 1996 Michael Messollen. All rights reserved. */ 2 2 /////////////////////////////////////////////////////////////////////////////// 3 static char * rcsid = "$Id: Factor.cc,v 1.4 0 2008-01-25 14:19:40Singular Exp $ ";3 static char * rcsid = "$Id: Factor.cc,v 1.41 2008-02-22 12:16:03 Singular Exp $ "; 4 4 static char * errmsg = "\nYou found a bug!\nPlease inform (Michael Messollen) michael@math.uni-sb.de \nPlease include above information and your input (the ideal/polynomial and characteristic) in your bug-report.\nThank you."; 5 5 /////////////////////////////////////////////////////////////////////////////// … … 190 190 191 191 retvalue= mydivremt(num,denum,a,b); 192 if ( retvalue && b == num.genZero() ) // num/denum from FFdomain192 if ( retvalue && b.isZero() ) // num/denum from FFdomain 193 193 return a; 194 194 else // num/denum is rational … … 257 257 numerator= numerator(x*lt,levelF); // x <- x*lt 258 258 denumerator= power(lt,degree(F,levelF)-1); // == lt^(1-degree(F,x) 259 while ( numerator.genZero() == is_rational(numerator, denumerator))259 while (is_rational(numerator, denumerator).isZero()) 260 260 numerator*= lt; 261 261 intermediate= is_rational(numerator,denumerator); … … 270 270 intermediate*= j.getItem().factor(); 271 271 test1= mydivremt( intermediate,F,a,b); 272 if ( test1 && b == intermediate.genZero() ) // Yupp!272 if ( test1 && b.isZero() ) // Yupp! 273 273 { 274 274 IntermediateList.append(CFFactor(1/a,1)); … … 288 288 test= content(numerator,x); 289 289 test1= mydivremt(denumerator,test,a,b); 290 if ( test1 && b == numerator.genZero() ) // Yupp!290 if ( test1 && b.isZero() ) // Yupp! 291 291 { 292 292 save_denumerator*= a; … … 308 308 // save_denumerator should be a multiple of the leading coeff 309 309 test1= mydivremt(save_denumerator,lt,a,b); 310 if ( test1 && b == save_denumerator.genZero() ) // Yupp!310 if ( test1 && b.isZero() ) // Yupp! 311 311 // We have to multiply one of the factors with 312 312 // the multiplicity of the save_denumerator <-> lc … … 440 440 return find_irreducible( degree_of_Extension, gen, Variable(1) ); 441 441 } 442 442 443 443 444 /////////////////////////////////////////////////////////////// … … 1329 1330 /* 1330 1331 $Log: not supported by cvs2svn $ 1332 Revision 1.40 2008/01/25 14:19:40 Singular 1333 *hannes: SqrFreeTest -> isSqrFree 1334 1331 1335 Revision 1.39 2008/01/22 09:51:37 Singular 1332 1336 *hannes: sqrFree/InternalSqrFree -> factory -
libfac/factor/MVMultiHensel.cc
r8078c8 r6f801b 2 2 /////////////////////////////////////////////////////////////////////////////// 3 3 // emacs edit mode for this file is -*- C++ -*- 4 // static char * rcsid = "$Id: MVMultiHensel.cc,v 1.1 5 2008-01-22 09:51:37Singular Exp $";4 // static char * rcsid = "$Id: MVMultiHensel.cc,v 1.16 2008-02-22 12:16:03 Singular Exp $"; 5 5 /////////////////////////////////////////////////////////////////////////////// 6 6 // FACTORY - Includes … … 502 502 CanonicalForm a,b; 503 503 // we add a division test now for intermediat.One and intermediat.Two 504 if ( mydivremt (mF,intermediat.One,a,b) && b == mF.genZero() )504 if ( mydivremt (mF,intermediat.One,a,b) && b.isZero() ) 505 505 Retlistinter.append(CFFactor(intermediat.One,1) ); 506 if ( mydivremt (mF,intermediat.Two,a,b) && b == mF.genZero() )506 if ( mydivremt (mF,intermediat.Two,a,b) && b.isZero() ) 507 507 Retlistinter.append(CFFactor(intermediat.Two,1) ); 508 508 … … 520 520 /* 521 521 $Log: not supported by cvs2svn $ 522 Revision 1.15 2008/01/22 09:51:37 Singular 523 *hannes: sqrFree/InternalSqrFree -> factory 524 522 525 Revision 1.14 2007/06/14 14:16:35 Singular 523 526 *hannes: Factorize2 etc. -
libfac/factor/Truefactor.cc
r8078c8 r6f801b 2 2 /////////////////////////////////////////////////////////////////////////////// 3 3 // emacs edit mode for this file is -*- C++ -*- 4 //static char * rcsid = "@(#) $Id: Truefactor.cc,v 1.1 2 2008-01-22 09:51:37Singular Exp $";4 //static char * rcsid = "@(#) $Id: Truefactor.cc,v 1.13 2008-02-22 12:16:03 Singular Exp $"; 5 5 /////////////////////////////////////////////////////////////////////////////// 6 6 // Factory - Includes … … 241 241 //CERR << " 1:" << change_poly(factor.factor(),SubstitutionList,1) <<"\n"; 242 242 c= mydivremt(U,test_f,a,b); 243 if ( c && b == U.genZero() && !hasAlgVar(test_f))243 if ( c && b.isZero() && !hasAlgVar(test_f)) 244 244 // factor.getFactor() divides U 245 245 { … … 450 450 /* 451 451 $Log: not supported by cvs2svn $ 452 Revision 1.12 2008/01/22 09:51:37 Singular 453 *hannes: sqrFree/InternalSqrFree -> factory 454 452 455 Revision 1.11 2008/01/07 13:34:56 Singular 453 456 *hannes: omse optiomzations(isOne)
Note: See TracChangeset
for help on using the changeset viewer.