Changeset f2e5d3 in git
- Timestamp:
- Sep 11, 2004, 4:31:57 PM (20 years ago)
- Branches:
- (u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
- Children:
- 496ccbf85092197727a74e5175eb81f14d1d85e2
- Parents:
- ac386536495e0f0b17da0d31bac81682cc5f37b7
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/tgb.cc
rac38653 rf2e5d3 2014 2014 2015 2015 found_in_row=row; 2016 assume(pn< mat->get_rows());2016 assume(pn<=mat->get_rows()); 2017 2017 for(i=row+1;i<pn;i++){ 2018 2018 int first=mat->min_col_not_zero_in_row(i); … … 2087 2087 v=done_index-1-v; 2088 2088 t->set(i,v,nCopy(p_i->coef)); 2089 p_i=p_i->next; 2089 2090 } 2090 2091 } … … 2275 2276 //Print("Done_index:%i",done_index); 2276 2277 if(done_index>0) 2278 { 2277 2279 pTest(done[0]); 2280 } 2278 2281 for(i=1;i<done_index;i++) 2279 2282 { … … 2287 2290 #ifdef TGB_DEBUG 2288 2291 for(i=0;i<done_index;i++) 2292 { 2289 2293 pTest(done[i]); 2294 } 2290 2295 #endif 2291 2296 poly* m; … … 2296 2301 m_size+=pLength(p[i]); 2297 2302 } 2303 m=(poly*) omalloc(m_size*sizeof(poly)); 2298 2304 //q=(poly*) omalloc(m_size*sizeof(poly)); 2299 2305 #ifdef TGB_DEBUG 2306 2307 for(i=0;i<done_index;i++) 2308 { 2309 2310 pTest(done[i]); 2311 } 2312 #endif 2313 2314 assume(done_index==0); 2300 2315 for(i=0;i<p_index;i++) 2301 2316 { 2317 assume(done_index==0); 2302 2318 poly p_i=p[i]; 2319 assume(done_index==0); 2320 pTest(p[i]); 2321 assume(done_index==0); 2303 2322 while(p_i) 2304 2323 { 2324 assume(done_index==0); 2305 2325 m[m_index]=pLmInit(p_i); 2326 assume(done_index==0); 2306 2327 pSetCoeff(m[m_index],nInit(1)); 2307 p_i=p_i->next; 2328 assume(done_index==0); 2329 p_i=p_i->next; 2330 assume(done_index==0); 2308 2331 m_index++; 2309 } 2310 } 2332 assume(done_index==0); 2333 } 2334 } 2335 #ifdef TGB_DEBUG 2336 2337 for(i=0;i<done_index;i++) 2338 { 2339 2340 pTest(done[i]); 2341 } 2342 #endif 2311 2343 int q_size=m_index; 2312 2344 poly* q=(poly*) omalloc(q_size*sizeof(poly)); 2313 2345 int q_index=0; 2314 2346 //next Step additional reductors 2347 #ifdef TGB_DEBUG 2348 2349 for(i=0;i<done_index;i++) 2350 { 2351 2352 pTest(done[i]); 2353 } 2354 #endif 2315 2355 while(m_index>0) 2316 2356 { 2357 #ifdef TGB_DEBUG 2358 2359 for(i=0;i<done_index;i++) 2360 { 2361 2362 pTest(done[i]); 2363 } 2364 #endif 2317 2365 qsort(m, m_index,sizeof(poly),pLmCmp_func); 2318 2366 2319 2367 2320 2368 pos=0; 2369 #ifdef TGB_DEBUG 2370 2371 for(i=0;i<done_index;i++) 2372 { 2373 2374 pTest(done[i]); 2375 } 2376 #endif 2321 2377 for(i=1;i<m_index;i++) 2322 2378 { 2379 pTest(m[i]); 2380 pTest(m[pos]); 2323 2381 if((!(pLmEqual(m[i],m[pos])))) 2324 2382 m[++pos]=m[i]; … … 2344 2402 q=(poly*) omalloc(q_size*sizeof(poly)); 2345 2403 } 2346 for(i=0;i<m_index;i++)2347 {2348 BOOLEAN in_done=FALSE;2349 pTest(m[i]);2350 2404 #ifdef TGB_DEBUG 2351 2405 int j; 2352 2406 for(j=0;j<done_index;j++) 2353 2407 { 2354 Print("%i %i\n",j,i); 2408 2409 pTest(done[j]); 2410 } 2411 #endif 2412 for(i=0;i<m_index;i++) 2413 { 2414 BOOLEAN in_done=FALSE; 2415 pTest(m[i]); 2416 #ifdef TGB_DEBUG 2417 2418 for(j=0;j<done_index;j++) 2419 { 2420 //Print("%i %i\n",j,i); 2355 2421 pTest(done[j]); 2356 2422 } … … 2682 2748 while(c->to_destroy) 2683 2749 { 2684 pDelete(& c->to_destroy->p);2750 pDelete(&(c->to_destroy->p)); 2685 2751 poly_list_node* old=c->to_destroy; 2686 2752 c->to_destroy=c->to_destroy->next;
Note: See TracChangeset
for help on using the changeset viewer.