Changeset 0d1a36 in git
- Timestamp:
- Mar 24, 2017, 6:38:29 PM (7 years ago)
- Branches:
- (u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
- Children:
- 94d691b55047cb8f22bc2c71361f61736cce019f
- Parents:
- a5db0e98aec0db35009110daba5e79981bcc3fbc
- Files:
-
- 28 edited
Legend:
- Unmodified
- Added
- Removed
-
kernel/GBEngine/janet.h
ra5db0e r0d1a36 16 16 poly root; //poly for parent, NULL for prol 17 17 kBucket_pt root_b; 18 unsignedroot_l;18 int root_l; 19 19 poly history; //parent 20 20 poly lead; //leading monomial for prolongation -
kernel/GBEngine/kInline.h
ra5db0e r0d1a36 568 568 } 569 569 570 KINLINE void sLObject::Tail_Minus_mm_Mult_qq(poly m, poly q, unsignedlq,570 KINLINE void sLObject::Tail_Minus_mm_Mult_qq(poly m, poly q, int lq, 571 571 poly spNoether) 572 572 { … … 581 581 assume(_p != NULL); 582 582 583 unsignedlp=pLength-1;583 int lp=pLength-1; 584 584 pNext(_p) = p_Minus_mm_Mult_qq( pNext(_p), m, q, lp, lq, 585 585 spNoether, tailRing ); -
kernel/GBEngine/kspoly.cc
ra5db0e r0d1a36 793 793 (void) ksCheckCoeff(&lc1, &lc2, currRing->cf); 794 794 795 unsigned l1=0; 796 unsigned l2=0; 795 int l1=0, l2=0; 797 796 798 797 if (p_GetComp(p1, currRing)!=p_GetComp(p2, currRing)) … … 841 840 if (spNoether != NULL) 842 841 { 843 l2 = 0;842 l2 = -1; 844 843 a2 = tailRing->p_Procs->pp_Mult_mm_Noether(a2, m2, spNoether, l2, tailRing); 845 844 assume(l2 == pLength(a2)); -
kernel/GBEngine/kutil.h
ra5db0e r0d1a36 79 79 long FDeg; // pFDeg(p) 80 80 int ecart, 81 length ;// as of pLDeg82 unsigned pLength;// either == 0, or == pLength(p)83 inti_r; // index of TObject in R set, or -1 if not in T81 length, // as of pLDeg 82 pLength, // either == 0, or == pLength(p) 83 i_r; // index of TObject in R set, or -1 if not in T 84 84 BOOLEAN is_normalized; // true, if pNorm was called on p, false otherwise 85 85 // used in incremental sba() with F5C: … … 217 217 KINLINE void SetLmTail(poly lm, poly new_p, int length, 218 218 int use_bucket, ring r); 219 KINLINE void Tail_Minus_mm_Mult_qq(poly m, poly qq, unsignedlq, poly spNoether);219 KINLINE void Tail_Minus_mm_Mult_qq(poly m, poly qq, int lq, poly spNoether); 220 220 KINLINE void Tail_Mult_nn(number n); 221 221 // deletes bucket, makes sure that p and t_p exists -
kernel/GBEngine/syz.cc
ra5db0e r0d1a36 223 223 void syGaussForOne(ideal syz, int elnum, int ModComp,int from,int till) 224 224 { 225 unsignedlu;225 int /*k,j,i,*/lu; 226 226 poly unit1,unit2; 227 227 poly actWith=syz->m[elnum]; -
kernel/GBEngine/syz.h
ra5db0e r0d1a36 26 26 int syzind; 27 27 int order; 28 unsignedlength;28 int length; 29 29 int reference; 30 30 }; -
kernel/GBEngine/syz0.cc
ra5db0e r0d1a36 511 511 ideal result=idInit(16,Fl); 512 512 int i,j,l,k,kkk,/*rkF,*/Sl=0,syComponentOrder=currRing->ComponentOrder; 513 int wend,lini,gencQ=0; 514 unsigned ltR; 513 int /*fstart,*/wend,lini,ltR,gencQ=0; 515 514 intvec *newmodcomp; 516 515 int *Flength; -
kernel/GBEngine/syz1.cc
ra5db0e r0d1a36 855 855 } 856 856 } 857 unsignedlb;857 int lb; 858 858 kBucketClear(syzstr->bucket,&tso.p,&lb); 859 859 } … … 2160 2160 intvec * toStrip) 2161 2161 { 2162 int ii=0,i,tc; 2163 unsigned ltS=0; 2164 unsigned lp; 2162 int ii=0,i,tc,lp,ltS=-1; 2165 2163 poly p,mp=NULL,pp; 2166 2164 SSet sPairs=syzstr->resPairs[index]; … … 2183 2181 tc = pGetComp(sPairs[i].isNotMinimal); 2184 2182 //p = pTakeOutComp1(&tempStripped,tc); 2185 unsignedlu;2183 int lu; 2186 2184 pTakeOutComp(&tempStripped,tc,&p,&lu); 2187 2185 kBucketTakeOutComp(syzstr->bucket,tc,&mp,&lp); … … 2191 2189 p = pNext(mp); 2192 2190 pNext(mp) = NULL; 2193 ltS = 0;2191 ltS = -1; 2194 2192 kBucket_Minus_m_Mult_p(syzstr->bucket,mp,tempStripped,<S); 2195 2193 pDelete(&mp); -
kernel/GBEngine/syz2.cc
ra5db0e r0d1a36 400 400 poly p=pDivide(q,redWith.p); 401 401 pSetCoeff(p,nDiv(pGetCoeff(q),pGetCoeff(redWith.p))); 402 unsigned il=0;402 int il=-1; 403 403 kBucket_Minus_m_Mult_p(syzstr->syz_bucket,p,redWith.syz,&il,NULL); 404 404 pLmDelete(&p); … … 409 409 ideal redWith=syzstr->res[index]; 410 410 if (redWith==NULL) return toRed; 411 unsigned j=IDELEMS(redWith); 412 unsigned i; 411 int j=IDELEMS(redWith),i; 413 412 poly q,result=NULL,resultp; 414 413 … … 505 504 int kres=(*syzstr->Tl)[index]; 506 505 int toGo=0; 507 unsignedil;506 int il; 508 507 SSet redset=syzstr->resPairs[index]; 509 508 poly q; -
kernel/GBEngine/syz3.cc
ra5db0e r0d1a36 472 472 * performs the modification of a single reduction on the syzygy-level 473 473 */ 474 static inline void sySPRedSyz_Kosz(syStrategy syzstr,poly redWith,poly syz,poly q=NULL,unsigned l_syz=0)474 inline void sySPRedSyz_Kosz(syStrategy syzstr,poly redWith,poly syz,poly q=NULL,int l_syz=-1) 475 475 { 476 476 poly p=pDivide(q,redWith); … … 514 514 */ 515 515 static poly syRedTailSyz(poly tored,ideal red,ideal sec_red,int crit_comp,syStrategy syzstr, 516 int * gen_length,int * secgen_length,unsigned * tored_length) 517 { 518 int i=IDELEMS(red)-1,num_mon; 519 unsigned num_tail; 516 int * gen_length,int * secgen_length,int * tored_length) 517 { 518 int i=IDELEMS(red)-1,num_mon,num_tail; 520 519 poly h,hn; 521 520 // BOOLEAN dummy; … … 660 659 int og_idel=IDELEMS(old_generators),ng_place=IDELEMS(new_generators); 661 660 int toReplace=0; 662 int i; 663 unsigned syz_l, j; 661 int i,j,syz_l; 664 662 number /*coefgcd,*/n; 665 663 polyset ogm=old_generators->m; … … 708 706 p = pNeg(p); 709 707 pp = pCopy(old_repr->m[tso.ind2]); 710 unsigned il=0;708 int il=-1; 711 709 while (p!=NULL) 712 710 { … … 718 716 p_Shift(&p,-1,currRing); 719 717 pp = pCopy(old_repr->m[tso.ind1]); 720 il= 0;718 il=-1; 721 719 while (p!=NULL) 722 720 { … … 1028 1026 */ 1029 1027 static poly normalize(poly next_p,ideal add_generators, syStrategy syzstr, 1030 int * g_l, unsigned* p_l,int crit_comp)1028 int * g_l,int * p_l,int crit_comp) 1031 1029 { 1032 1030 int j=0,i=IDELEMS(add_generators); … … 1194 1192 assume (tso.lcm!=NULL); 1195 1193 int ng_place=IDELEMS(new_generators); 1196 unsignedi,j;1194 int i,j; 1197 1195 number n; 1198 1196 poly p; … … 1422 1420 int idel_temp=IDELEMS(temp_generators); 1423 1421 int next_place_add; 1424 int red_deg,l_pairs=IDELEMS(add_generators); 1425 unsigned p_length; 1422 int p_length,red_deg,l_pairs=IDELEMS(add_generators); 1426 1423 poly next_p; 1427 1424 int * gen_length=(int*)omAlloc0(IDELEMS(add_generators)*sizeof(int)); … … 1551 1548 ideal old_generators= syzstr->res[index]; 1552 1549 ideal old_repr= syzstr->orderedRes[index]; 1553 int i,j=0,ii=IDELEMS(old_generators)-1; 1554 unsigned dummy; 1550 int i,j=0,ii=IDELEMS(old_generators)-1,dummy; 1555 1551 poly p; 1556 1552 number n; -
kernel/GBEngine/tgb.cc
ra5db0e r0d1a36 53 53 static BOOLEAN has_t_rep(const int & arg_i, const int & arg_j, slimgb_alg* state); 54 54 static void shorten_tails(slimgb_alg* c, poly monom); 55 static poly redNF2 (poly h,slimgb_alg* c , unsigned&len, number& m,int n=0);55 static poly redNF2 (poly h,slimgb_alg* c , int &len, number& m,int n=0); 56 56 static poly redNFTail (poly h,const int sl,kStrategy strat, int len); 57 57 static int bucket_guess(kBucket* bucket); … … 1835 1835 } 1836 1836 1837 static poly redNF2 (poly h, slimgb_alg * c, unsigned&len, number & m, int n)1837 static poly redNF2 (poly h, slimgb_alg * c, int &len, number & m, int n) 1838 1838 { 1839 1839 m = nInit (1); … … 2791 2791 // now_t_rep(s->j,s->i,c); 2792 2792 number coef; 2793 unsignedmlen = pLength (h);2793 int mlen = pLength (h); 2794 2794 p_Test (h, c->r); 2795 2795 if((!c->nc) & (!(use_noro))) … … 2932 2932 for(j = 0; j < i; j++) 2933 2933 { 2934 unsignedlen;2934 int len; 2935 2935 poly p; 2936 2936 buf[j].flatten (); … … 4703 4703 for(i = 0; i < pn_noro; i++) 4704 4704 { 4705 unsigneddummy_len;4705 int dummy_len; 4706 4706 poly p; 4707 4707 los[i].p = NULL; … … 4866 4866 { 4867 4867 flatten (); 4868 unsignedl;4868 int l; 4869 4869 kBucketClear (bucket, &p, &l); 4870 4870 return l; … … 4934 4934 int rn = erg.reduce_by; 4935 4935 poly red; 4936 unsignedred_len;4936 int red_len; 4937 4937 simple_reducer *pointer; 4938 4938 BOOLEAN work_on_copy = FALSE; -
kernel/GBEngine/tgb_internal.h
ra5db0e r0d1a36 384 384 { 385 385 poly expand; 386 unsignedexpand_length;386 int expand_length; 387 387 int to_reduce_u; 388 388 int to_reduce_l; -
kernel/digitech.cc
ra5db0e r0d1a36 38 38 } 39 39 40 unsigned len;40 int len=0; 41 41 poly erg; 42 42 kBucketClear(erg_bucket,&erg, &len); … … 80 80 } 81 81 omfree(potences); 82 unsignedlen=0;82 int len=0; 83 83 poly erg; 84 84 kBucketClear(erg_bucket,&erg, &len); -
kernel/fast_mult.cc
ra5db0e r0d1a36 665 665 omfree(facult); 666 666 omfree(f_terms); 667 unsignedlen=0;667 int len=0; 668 668 poly erg; 669 669 kBucketClear(erg_bucket,&erg, &len); -
kernel/linear_algebra/MinorProcessor.cc
ra5db0e r0d1a36 1320 1320 or at least an integral domain. */ 1321 1321 void elimOperationBucket(poly &p1, poly &p2, poly &p3, poly &p4, poly &p5, 1322 number &c5, unsignedp5Len)1322 number &c5, int p5Len) 1323 1323 { 1324 1324 #ifdef COUNT_AND_PRINT_OPERATIONS -
kernel/polys.h
ra5db0e r0d1a36 319 319 /// component == comp and *p of all other monoms *lq == pLength(*q) 320 320 /// On return all components pf *q == 0 321 static inline void pTakeOutComp(poly *p, long comp, poly *q, unsigned*lq, const ring R = currRing)321 inline void pTakeOutComp(poly *p, long comp, poly *q, int *lq, const ring R = currRing) 322 322 { 323 323 return p_TakeOutComp(p, comp, q, lq, R); … … 326 326 327 327 /// This is something weird -- Don't use it, unless you know what you are doing 328 static inline polypTakeOutComp(poly * p, int k, const ring R = currRing)328 inline poly pTakeOutComp(poly * p, int k, const ring R = currRing) 329 329 { 330 330 return p_TakeOutComp(p, k, R); … … 343 343 #define pDeleteComp(p,k) p_DeleteComp(p,k,currRing) 344 344 345 static inline void pNorm(poly p, const ring R = currRing){ p_Norm(p, R); } 345 inline void pNorm(poly p, const ring R = currRing){ p_Norm(p, R); } 346 346 347 347 348 #define pSubst(p,n,e) p_Subst(p,n,e,currRing) … … 385 386 /// returns the length of a polynomial (numbers of monomials) 386 387 /// respect syzComp 387 static inline poly pLast(poly a, unsigned&length) { return p_Last (a, length, currRing); }388 static inline poly pLast(poly a) { unsignedl; return pLast(a, l); }388 static inline poly pLast(poly a, int &length) { return p_Last (a, length, currRing); } 389 static inline poly pLast(poly a) { int l; return pLast(a, l); } 389 390 390 391 /*************************************************************** -
libpolys/polys/kbuckets.cc
ra5db0e r0d1a36 366 366 poly p = bucket->buckets[1]; 367 367 poly lm; 368 unsigned pl = bucket->buckets_length[1];368 int pl = bucket->buckets_length[1];//, i; 369 369 int i; 370 370 bucket->buckets[1] = NULL; … … 432 432 } 433 433 434 void kBucketClear(kBucket_pt bucket, poly *p, unsigned*length)434 void kBucketClear(kBucket_pt bucket, poly *p, int *length) 435 435 { 436 436 int i = kBucketCanonicalize(bucket); … … 537 537 /// Bucket number i from bucket is out of length sync, resync 538 538 /// 539 void kBucketAdjust(kBucket_pt bucket, int i) 540 { 539 void kBucketAdjust(kBucket_pt bucket, int i) { 540 541 541 MULTIPLY_BUCKET(bucket,i); 542 542 543 unsignedl1 = bucket->buckets_length[i];543 int l1 = bucket->buckets_length[i]; 544 544 poly p1 = bucket->buckets[i]; 545 545 bucket->buckets[i] = NULL; … … 631 631 assume(*l <= 0 || pLength(q) == *l); 632 632 633 unsignedi, l1;633 int i, l1; 634 634 ring r = bucket->bucket_ring; 635 635 … … 688 688 /// Does not destroy p and m 689 689 /// assume (*l <= 0 || pLength(p) == *l) 690 void kBucket_Minus_m_Mult_p(kBucket_pt bucket, poly m, poly p, unsigned*l,690 void kBucket_Minus_m_Mult_p(kBucket_pt bucket, poly m, poly p, int *l, 691 691 poly spNoether) 692 692 { 693 assume(*l == 0 || pLength(p) == *l);694 unsignedi, l1;693 assume(*l <= 0 || pLength(p) == *l); 694 int i, l1; 695 695 poly p1 = p; 696 696 ring r = bucket->bucket_ring; 697 697 698 if (*l == 0)698 if (*l <= 0) 699 699 { 700 700 l1 = pLength(p1); … … 797 797 assume((!rIsPluralRing(bucket->bucket_ring))||p_IsConstant(m, bucket->bucket_ring)); 798 798 assume(l <= 0 || pLength(p) == l); 799 unsignedi, l1;799 int i, l1; 800 800 poly p1 = p; 801 801 ring r = bucket->bucket_ring; … … 1012 1012 void kBucketTakeOutComp(kBucket_pt bucket, 1013 1013 long comp, 1014 poly *r_p, unsigned*l)1014 poly *r_p, int *l) 1015 1015 { 1016 1016 poly p = NULL, q; 1017 unsignedi, lp = 0, lq;1017 int i, lp = 0, lq; 1018 1018 1019 1019 #ifndef HAVE_PSEUDO_BUCKETS … … 1052 1052 1053 1053 number kBucketPolyRed(kBucket_pt bucket, 1054 poly p1, unsignedl1,1054 poly p1, int l1, 1055 1055 poly spNoether) 1056 1056 { … … 1059 1059 assume(p1 != NULL && 1060 1060 p_DivisibleBy(p1, kBucketGetLm(bucket), r)); 1061 assume(pLength(p1) == l1);1061 assume(pLength(p1) == (int) l1); 1062 1062 1063 1063 poly a1 = pNext(p1), lm = kBucketExtractLm(bucket); -
libpolys/polys/kbuckets.h
ra5db0e r0d1a36 40 40 // Converts Bpoly into a poly and clears bucket 41 41 // i.e., afterwards Bpoly == 0 42 void kBucketClear(kBucket_pt bucket, poly *p, unsigned*length);42 void kBucketClear(kBucket_pt bucket, poly *p, int *length); 43 43 44 44 inline poly kBucketClear(kBucket_pt bucket) 45 45 { 46 unsigneddummy;46 int dummy; 47 47 poly p; 48 48 kBucketClear(bucket, &p, &dummy); … … 81 81 // Returns: Lc(p) / gcd(Lc(p), Lc(q)) 82 82 number kBucketPolyRed(kBucket_pt bucket, 83 poly p, unsignedl,83 poly p, int l, 84 84 poly spNoether); 85 85 … … 95 95 void kBucketTakeOutComp(kBucket_pt bucket, 96 96 long comp, 97 poly *p, unsigned*l);97 poly *p, int *l); 98 98 99 99 ////////////////////////////////////////////////////////////////////////// … … 131 131 /// Does not destroy p and m (TODO: rename into kBucket_Minus_mm_Mult_pp!?) 132 132 /// assume (*l <= 0 || pLength(p) == *l) 133 void kBucket_Minus_m_Mult_p(kBucket_pt bucket, poly m, poly p, unsigned*l,133 void kBucket_Minus_m_Mult_p(kBucket_pt bucket, poly m, poly p, int *l, 134 134 poly spNother = NULL); 135 135 … … 186 186 poly coef[MAX_BUCKET + 1]; // coeff of polys in bucket or NULL : 2..max 187 187 #endif 188 unsignedbuckets_length[MAX_BUCKET + 1]; // length if i-th poly189 unsignedbuckets_used; // max number of used bucket188 int buckets_length[MAX_BUCKET + 1]; // length if i-th poly 189 int buckets_used; // max number of used bucket 190 190 #endif 191 191 ring bucket_ring; -
libpolys/polys/monomials/p_polys.cc
ra5db0e r0d1a36 3435 3435 // Splits *p into two polys: *q which consists of all monoms with 3436 3436 // component == comp and *p of all other monoms *lq == pLength(*q) 3437 void p_TakeOutComp(poly *r_p, long comp, poly *r_q, unsigned*lq, const ring r)3437 void p_TakeOutComp(poly *r_p, long comp, poly *r_q, int *lq, const ring r) 3438 3438 { 3439 3439 spolyrec pp, qq; … … 4471 4471 * respect syzComp 4472 4472 */ 4473 poly p_Last(const poly p, unsigned&l, const ring r)4473 poly p_Last(const poly p, int &l, const ring r) 4474 4474 { 4475 4475 if (p == NULL) -
libpolys/polys/monomials/p_polys.h
ra5db0e r0d1a36 200 200 // returns the length of a polynomial (numbers of monomials) and the last mon. 201 201 // respect syzComp 202 poly p_Last(const poly a, unsigned&l, const ring r);202 poly p_Last(const poly a, int &l, const ring r); 203 203 204 204 /*----------------------------------------------------*/ … … 883 883 884 884 /// like p_Add_q, except that if lp == pLength(lp) lq == pLength(lq) then lp == pLength(p+q) 885 static inline poly p_Add_q(poly p, poly q, unsigned &lp, unsignedlq, const ring r)885 static inline poly p_Add_q(poly p, poly q, int &lp, int lq, const ring r) 886 886 { 887 887 assume( (p != q) || (p == NULL && q == NULL) ); … … 974 974 } 975 975 976 static inline poly p_Minus_mm_Mult_qq(poly p, const poly m, const poly q, 977 unsigned &lp, unsigned lq,const poly spNoether, const ring r)976 static inline poly p_Minus_mm_Mult_qq(poly p, const poly m, const poly q, int &lp, int lq, 977 const poly spNoether, const ring r) 978 978 { 979 979 int shorter; … … 1104 1104 1105 1105 // returns p + m*q destroys p, const: q, m 1106 static inline poly p_Plus_mm_Mult_qq(poly p, poly m, poly q, unsigned &lp, unsignedlq,1106 static inline poly p_Plus_mm_Mult_qq(poly p, poly m, poly q, int &lp, int lq, 1107 1107 const ring r) 1108 1108 { … … 1128 1128 static inline poly p_Plus_mm_Mult_qq(poly p, poly m, poly q, const ring r) 1129 1129 { 1130 unsignedlp = 0, lq = 0;1130 int lp = 0, lq = 0; 1131 1131 return p_Plus_mm_Mult_qq(p, m, q, lp, lq, r); 1132 1132 } … … 1975 1975 // component == comp and *p of all other monoms *lq == pLength(*q) 1976 1976 // On return all components pf *q == 0 1977 void p_TakeOutComp(poly *p, long comp, poly *q, unsigned*lq, const ring r);1977 void p_TakeOutComp(poly *p, long comp, poly *q, int *lq, const ring r); 1978 1978 1979 1979 // This is something weird -- Don't use it, unless you know what you are doing -
libpolys/polys/nc/nc.h
ra5db0e r0d1a36 231 231 // // for p_Plus_mm_Mult_qq in pInline2.h 232 232 // returns p + m*q destroys p, const: q, m 233 poly nc_p_Plus_mm_Mult_qq(poly p, const poly m, const poly q, unsigned&lp,234 const unsigned, const ring r);233 poly nc_p_Plus_mm_Mult_qq(poly p, const poly m, const poly q, int &lp, 234 const int, const ring r); 235 235 236 236 -
libpolys/polys/nc/old.gring.cc
ra5db0e r0d1a36 171 171 172 172 // returns p + m*q destroys p, const: q, m 173 poly nc_p_Plus_mm_Mult_qq(poly p, const poly m, const poly q, unsigned&lp,174 const unsigned, const ring r)173 poly nc_p_Plus_mm_Mult_qq(poly p, const poly m, const poly q, int &lp, 174 const int, const ring r) 175 175 { 176 176 p = p_Add_q(p, nc_mm_Mult_pp( m, q, r ), r); -
libpolys/polys/sbuckets.cc
ra5db0e r0d1a36 174 174 } 175 175 176 void sBucket_Merge_p(sBucket_pt bucket, poly p, unsignedlength)176 void sBucket_Merge_p(sBucket_pt bucket, poly p, int length) 177 177 { 178 178 assume(bucket != NULL); 179 assume(length == 0 || length == pLength(p));179 assume(length <= 0 || length == pLength(p)); 180 180 181 181 if (p == NULL) return; 182 if (length == 0) length = pLength(p);182 if (length <= 0) length = pLength(p); 183 183 184 184 int i = LOG2(length); … … 199 199 } 200 200 201 void sBucket_Add_p(sBucket_pt bucket, poly p, unsignedlength)201 void sBucket_Add_p(sBucket_pt bucket, poly p, int length) 202 202 { 203 203 assume(bucket != NULL); … … 271 271 { 272 272 poly pr = NULL; 273 unsignedlr = 0;274 unsignedi = 0;273 int lr = 0; 274 int i = 0; 275 275 276 276 while (bucket->buckets[i].p == NULL) -
libpolys/polys/sbuckets.h
ra5db0e r0d1a36 82 82 /// Merges p into Spoly: assumes Bpoly and p have no common monoms 83 83 /// destroys p! 84 void sBucket_Merge_p(sBucket_pt bucket, poly p, unsignedlp);84 void sBucket_Merge_p(sBucket_pt bucket, poly p, int lp); 85 85 86 86 /// adds poly p to bucket 87 87 /// destroys p! 88 void sBucket_Add_p(sBucket_pt bucket, poly p, unsignedlp);88 void sBucket_Add_p(sBucket_pt bucket, poly p, int lp); 89 89 90 90 -
libpolys/polys/templates/p_Minus_mm_Mult_qq__T.cc
ra5db0e r0d1a36 140 140 if (spNoether != NULL) 141 141 { 142 unsignedll = 0;142 int ll = 0; 143 143 pNext(a) = r->p_Procs->pp_Mult_mm_Noether(q, m, spNoether, ll, r); 144 144 shorter += ll; -
libpolys/polys/templates/p_Procs.h
ra5db0e r0d1a36 28 28 typedef poly (*pp_Mult_mm_Proc_Ptr)(poly p, const poly m, const ring r); 29 29 typedef poly (*pp_Mult_mm_Noether_Proc_Ptr)(poly p, const poly m, 30 const poly spNoether, unsigned&ll,30 const poly spNoether, int &ll, 31 31 const ring r); 32 32 typedef poly (*p_Add_q_Proc_Ptr)(poly p, poly q, int & shorter, const ring r); -
libpolys/polys/templates/p_kBucketSetLm__T.cc
ra5db0e r0d1a36 38 38 { 39 39 j = 0; 40 for ( unsignedi = 1; i<=bucket->buckets_used; i++)40 for (int i = 1; i<=bucket->buckets_used; i++) 41 41 { 42 42 if (bucket->buckets[i] != NULL) -
libpolys/polys/templates/pp_Mult_mm_Noether__T.cc
ra5db0e r0d1a36 17 17 * 18 18 ***************************************************************/ 19 LINKAGE poly pp_Mult_mm_Noether__T(poly p, const poly m, const poly spNoether, unsigned&ll, const ring ri)19 LINKAGE poly pp_Mult_mm_Noether__T(poly p, const poly m, const poly spNoether, int &ll, const ring ri) 20 20 { 21 21 p_Test(p, ri); … … 38 38 pAssume(!n_IsZero__T(ln,ri->cf)); 39 39 pAssume1(p_GetComp(m, ri) == 0 || p_MaxComp(p, ri) == 0); 40 int l = 0; 40 41 41 42 do … … 63 64 #endif 64 65 { 66 l++; 65 67 q = pNext(q) = r; 66 68 pSetCoeff0(q, n); … … 71 73 } while (p != NULL); 72 74 73 ll = pLength(p); 75 if (ll < 0) 76 ll = l; 77 else 78 ll = pLength(p); 74 79 75 80 pNext(q) = NULL;
Note: See TracChangeset
for help on using the changeset viewer.