Changeset ea125f in git
- Timestamp:
- May 19, 2011, 12:10:13 PM (12 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', '00e2e9c41af3fde1273eb3633f4c0c7c3db2579d')
- Children:
- b4f7ebdb67074ac439d1df267fda83171f69d3b8
- Parents:
- 5f32d50562e470047f7956821b2de0bd34fe3b6b
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/assprimeszerodim.lib
r5f32d50 rea125f 262 262 } 263 263 264 int T = timer; 265 int RT = rtimer; 266 int TT; 267 264 268 def SPR = basering; 269 if(printlevel >= 10) { "========== Start modStd =========="; } 265 270 I = modStd(I,n); 271 if(printlevel >= 10) { "=========== End modStd ==========="; } 272 if(printlevel >= 9) { "modStd takes "+string(rtimer-RT)+" seconds."; } 266 273 int d = vdim(I); 267 274 if(d == -1) { ERROR("Ideal is not zero-dimensional."); } 268 275 if(homog(I) == 1){ return(list(maxideal(1))); } 269 276 poly f = findGen(I); 270 271 int T = timer; 272 int RT = rtimer; 273 int TT; 277 if(printlevel >= 9) { "Coordinate change: "+string(f); } 278 274 279 if(size(f) == nvars(SPR)) 275 280 { 276 281 TT = timer; 277 282 int spT = pTestRad(d,f,I); 278 if(printlevel >=10)279 { 280 "pTestRad(d,f,I) = "+string(spT)+" and takes 281 "+string(timer-TT)+" seconds.";283 if(printlevel >= 9) 284 { 285 "pTestRad(d,f,I) = "+string(spT)+" and takes " 286 +string(timer-TT)+" seconds."; 282 287 } 283 288 if(!spT) … … 289 294 TT = timer; 290 295 I = zeroR(I,n); 291 if(printlevel >=10)296 if(printlevel >= 9) 292 297 { 293 298 "zeroR(I,n) takes "+string(timer-TT)+" seconds."; … … 295 300 TT = timer; 296 301 I = std(I); 297 if(printlevel >=10)302 if(printlevel >= 9) 298 303 { 299 304 "std(I) takes "+string(timer-TT)+" seconds."; … … 307 312 TT = timer; 308 313 I = zeroR(I,n); 309 if(printlevel >=10)314 if(printlevel >= 9) 310 315 { 311 316 "zeroR(I,n) takes "+string(timer-TT)+" seconds."; … … 313 318 TT = timer; 314 319 I = std(I); 315 if(printlevel >=10)320 if(printlevel >= 9) 316 321 { 317 322 "std(I) takes "+string(timer-TT)+" seconds."; … … 322 327 } 323 328 } 324 if(printlevel >=10)329 if(printlevel >= 9) 325 330 { 326 331 "Real-time for radical-check is "+string(rtimer - RT)+" seconds."; … … 341 346 342 347 ring rHelp = 0,T,dp; 343 list CO1,CO2,P ;348 list CO1,CO2,P,H; 344 349 ideal F,G,testF; 345 350 bigint N; 346 351 347 352 list ringL = ringlist(SPR); 348 int i,k,e,int_break ;353 int i,k,e,int_break,s; 349 354 int j = 1; 350 355 int index = 1; … … 437 442 } 438 443 } 439 if(printlevel >=10)444 if(printlevel >= 9) 440 445 { 441 446 "Real-time for computing list in assPrimes is "+string(rtimer - rt)+ … … 452 457 if(F[1]-testF[1]==0) 453 458 { 454 if(printlevel>=10) 455 { 456 "size(L) = "+string(size(L)); 457 } 459 if(printlevel >= 9) { "size(L) = "+string(size(L)); } 460 458 461 F = cleardenom(F[1]); 459 462 … … 461 464 if(e == d) 462 465 { 463 listH = factorize(F[1]);464 465 ints = size(H[1]);466 H = factorize(F[1]); 467 468 s = size(H[1]); 466 469 for(i = 1; i <= s; i++) 467 470 { … … 477 480 map phi = rHelp,var(nvars(SPR)); 478 481 list H = phi(H); 479 if(printlevel >=10)482 if(printlevel >= 9) 480 483 { 481 484 "Real-time without test is "+string(rtimer - RT)+" seconds."; … … 525 528 if(F1 != 0) { int_break = 1; } 526 529 } 527 if(printlevel >=10)530 if(printlevel >= 9) 528 531 { 529 532 "Real-time for test is "+string(rtimer - RT)+" seconds."; … … 579 582 assPrimes(I); 580 583 } 584 581 585 //////////////////////////////////////////////////////////////////////////////// 582 586
Note: See TracChangeset
for help on using the changeset viewer.