Ignore:
Timestamp:
Feb 21, 2012, 2:28:22 PM (12 years ago)
Author:
Oleksandr Motsak <motsak@…>
Branches:
(u'spielwiese', '17f1d200f27c5bd38f5dfc6e8a0879242279d1d8')
Children:
74f51f261514079d9bce01fd1a8ffed4de458527
Parents:
5c0183a7fbcc63348c4fbf0ae70d97f6d0fcb203
git-author:
Oleksandr Motsak <motsak@mathematik.uni-kl.de>2012-02-21 14:28:22+01:00
git-committer:
Oleksandr Motsak <motsak@mathematik.uni-kl.de>2012-02-23 20:35:10+01:00
Message:
fix the Schreyer ordering for bigger components...
File:
1 edited

Legend:

Unmodified
Added
Removed
  • libpolys/polys/monomials/p_polys.cc

    r5c0183 re4f491  
    286286          const ideal F = o->data.is.F;
    287287          const int limit = o->data.is.limit;
     288          assume( limit >= 0 );
    288289          const int start = o->data.is.start;
    289            
    290290
    291291          if( F != NULL && c > limit )
     
    298298#endif
    299299#endif
     300//          if( c > limit ) // BUG???
     301            p->exp[start] = 1;
     302//          else
     303//            p->exp[start] = 0;
     304
    300305
    301306            c -= limit;
     
    305310            assume( c < IDELEMS(F) ); // What about others???
    306311
     312            if( c >= IDELEMS(F) )
     313              break;
     314
    307315            const poly pp = F->m[c]; // get reference monomial!!!
     316
     317            assume(pp != NULL);
    308318
    309319#ifndef NDEBUG
     
    313323#endif
    314324#endif
    315 
    316 
    317             assume(pp != NULL);
    318325            if(pp == NULL) break;
    319326
    320327            const int end = o->data.is.end;
    321 
    322328            assume(start <= end);
    323              
    324 //          const int limit = o->data.is.limit;
    325           assume( limit >= 0 );
     329
    326330
    327331//        const int st = o->data.isTemp.start;       
    328 
    329 //          if( c > limit ) // BUG???
    330             p->exp[start] = 1;
    331 //          else
    332 //            p->exp[start] = 0;
    333 
    334332
    335333#ifndef NDEBUG
    336334            Print("p_Setm_General: is(-Temp-) :: c: %d, limit: %d, [st:%d] ===>>> %ld\n", c, limit, start, p->exp[start]);
    337 #endif       
     335#endif
    338336
    339337            // p_ExpVectorAdd(p, pp, r);
Note: See TracChangeset for help on using the changeset viewer.