Changeset f3703e in git
- Timestamp:
- Feb 21, 2014, 6:05:20 PM (10 years ago)
- Branches:
- (u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')
- Children:
- 654726fdd2ca289b1ef3cbf4e931a1502bf145ef
- Parents:
- 2e2f8fd86ee1b100bc0848cbffdbdc7fd2ced7d6
- git-author:
- Max Horn <max@quendi.de>2014-02-21 18:05:20+01:00
- git-committer:
- Max Horn <max@quendi.de>2014-02-21 18:05:34+01:00
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/walk.cc
r2e2f8f rf3703e 4611 4611 clock_t tim; 4612 4612 nstep=0; 4613 int i, k,nwalk,endwalks = 0;4613 int i,nwalk,endwalks = 0; 4614 4614 int nV = currRing->N; 4615 4615 … … 4696 4696 else 4697 4697 { 4698 NORMAL_GW:4699 4698 #ifndef BUCHBERGER_ALG 4700 4699 if(isNolVector(curr_weight) == 0) … … 6790 6789 clock_t tinput = clock(); 6791 6790 #endif 6792 int nsteppert=0, i, k,nV = currRing->N, nwalk=0, npert_tmp=0;6791 int nsteppert=0, i, nV = currRing->N, nwalk=0, npert_tmp=0; 6793 6792 int *npert=(int*)omAlloc(2*nV*sizeof(int)); 6794 6793 ideal Gomega, M,F, G1, Gomega1, Gomega2, M1, F1; -
kernel/khstd.cc
r2e2f8f rf3703e 155 155 #if ADIDEBUG 156 156 PrintS("\nOriginal\n"); 157 int j, l, k;157 int i, j, l, k; 158 158 if (hilb == NULL) 159 159 return; -
libpolys/polys/clapsing.cc
r2e2f8f rf3703e 1464 1464 On(SW_SYMMETRIC_FF); 1465 1465 CFList L; 1466 poly p;1467 1466 if (rField_is_Q(r) || rField_is_Zp(r)) 1468 1467 { -
libpolys/polys/linalg_from_matpol.cc
r2e2f8f rf3703e 7 7 8 8 static void mp_PartClean(matrix, int, int, const ring); 9 static void mp_FinalClean(matrix, const ring);10 9 static int mp_PrepareRow (matrix, int, int, const ring); 11 10 static int mp_PreparePiv (matrix, int, int, const ring); … … 231 230 232 231 232 static poly minuscopy (poly p, const ring R) 233 { 234 poly w; 235 number e; 236 e = n_Init(-1, R); 237 w = p_Copy(p, R); 238 p_Mult_nn(w, e, R); 239 n_Delete(&e, R); 240 return w; 241 } 242 233 243 234 244 /*2 -
libpolys/polys/matpol.cc
r2e2f8f rf3703e 610 610 } 611 611 612 static poly minuscopy (poly p, const ring R)613 {614 poly w;615 number e;616 e = n_Init(-1, R);617 w = p_Copy(p, R);618 p_Mult_nn(w, e, R);619 n_Delete(&e, R);620 return w;621 }622 623 612 /*2 624 613 * insert a monomial into a list, avoid duplicates … … 727 716 for (j=lc-1;j>=0;j--) if(q1[j]) p_Delete(&q1[j], R); 728 717 } 729 }730 731 static void mp_FinalClean(matrix a, const ring)732 {733 omFreeSize((ADDRESS)a->m,a->nrows*a->ncols*sizeof(poly));734 omFreeBin((ADDRESS)a, sip_sideal_bin);735 718 } 736 719 … … 1368 1351 } 1369 1352 1370 static inline BOOLEAN smSmaller(poly a, poly b)1371 {1372 loop1373 {1374 pIter(b);1375 if (b == NULL) return TRUE;1376 pIter(a);1377 if (a == NULL) return FALSE;1378 }1379 }1380 1381 static BOOLEAN sm_IsNegQuot(poly a, const poly b, const poly c, const ring R)1382 {1383 if (p_LmDivisibleByNoComp(c, b, R))1384 {1385 p_ExpVectorDiff(a, b, c, R);1386 // Hmm: here used to be a pSetm(a): but it is unnecessary,1387 // if b and c are correct1388 return FALSE;1389 }1390 else1391 {1392 int i;1393 for (i=rVar(R); i>0; i--)1394 {1395 if(p_GetExp(c,i,R) > p_GetExp(b,i,R))1396 p_SetExp(a,i,p_GetExp(c,i,R)-p_GetExp(b,i,R),R);1397 else1398 p_SetExp(a,i,0,R);1399 }1400 // here we actually might need a pSetm, if a is to be used in1401 // comparisons1402 return TRUE;1403 }1404 }1405 1406 1353 static void mp_ElimBar(matrix a0, matrix re, poly div, int lr, int lc, const ring R) 1407 1354 { -
libpolys/polys/monomials/p_polys.cc
r2e2f8f rf3703e 2814 2814 assume( r != NULL ); assume( r->cf != NULL ); const coeffs C = r->cf; 2815 2815 2816 poly start=ph; 2817 2818 number d, h; 2816 number h; 2819 2817 poly p; 2820 2818 -
libpolys/polys/simpleideals.cc
r2e2f8f rf3703e 34 34 static int idpowerpoint; 35 35 /*index of the actual monomial in idpower*/ 36 static poly * givenideal;37 /*the ideal from which a power is computed*/38 36 39 37 /*2 … … 997 995 } 998 996 999 /*21000 *computes recursively all generators of a certain degree1001 *of the ideal "givenideal"1002 *elms is the number elements in the given ideal1003 *actelm is the actual element to handle1004 *deg is the degree of the power to compute1005 *gendeg is the actual degree of the generator in consideration1006 */1007 static void makepotence(int elms,int actelm,int deg,int gendeg, const ring r)1008 {1009 poly p;1010 int i=0;1011 1012 if ((idpowerpoint == 0) && (actelm ==1))1013 {1014 idpower[idpowerpoint] = p_One(r);1015 gendeg = 0;1016 }1017 while (i<=deg)1018 {1019 if (deg == gendeg)1020 {1021 idpowerpoint++;1022 return;1023 }1024 if (actelm == elms)1025 {1026 p=p_Power(p_Copy(givenideal[actelm-1],r),deg-gendeg,r);1027 idpower[idpowerpoint]=p_Mult_q(idpower[idpowerpoint],p,r);1028 idpowerpoint++;1029 return;1030 }1031 else1032 {1033 p = p_Copy(idpower[idpowerpoint],r);1034 makepotence(elms,actelm+1,deg,gendeg,r);1035 idpower[idpowerpoint] = p;1036 }1037 gendeg++;1038 idpower[idpowerpoint]=p_Mult_q(idpower[idpowerpoint],p_Copy(givenideal[actelm-1],r),r);1039 i++;1040 }1041 }1042 1043 /*21044 *returns the deg-th power of the ideal gid1045 */1046 //ideal idPower(ideal gid,int deg)1047 //{1048 // int i;1049 // ideal id;1050 //1051 // if (deg < 1) deg = 1;1052 // i = binom(IDELEMS(gid)+deg-1,deg);1053 // id=idInit(i,1);1054 // idpower = id->m;1055 // givenideal = gid->m;1056 // idpowerpoint = 0;1057 // makepotence(IDELEMS(gid),1,deg,0);1058 // idpower = NULL;1059 // givenideal = NULL;1060 // idpowerpoint = 0;1061 // return id;1062 //}1063 997 static void id_NextPotence(ideal given, ideal result, 1064 998 int begin, int end, int deg, int restdeg, poly ap, const ring r)
Note: See TracChangeset
for help on using the changeset viewer.