41#define WITH_SCHREYER_ORD
66 for (
i=old_hilbs->
length()+deg-1;
i>=deg;
i--)
67 (*shifted_old_hilbs)[
i] = (*old_hilbs)[
i-deg];
69 (*old_hilb1)[
i] = (*old_hilbs)[
i]-(*shifted_old_hilbs)[
i];
71 (*old_hilb1)[
i] = -(*shifted_old_hilbs)[
i];
73 (*new_hilb1)[
i] = (*new_hilbs)[
i];
75 while ((
i<biggest_length) && isRegular)
77 isRegular = isRegular && ((*old_hilb1)[
i] == (*new_hilb1)[
i]);
84 delete shifted_old_hilbs;
98 Print(
"aktueller hoechster index ist: %d\n",(*syzstr->
Tl)[
i]);
99 Print(
"der %d-te modul ist:\n",
i);
101 PrintS(
"Seine Darstellung:\n");
112 ideal old_repr,
int old_tl, poly next_generator,
resolvente totake)
115 poly gen=
pCopy(next_generator),
p;
116 poly neg_gen=
pCopy(next_generator);
117 ideal current_ideal,current_repr;
119 poly w_gen=
pHead(next_generator);
124 neg_gen =
pNeg(neg_gen);
138 current_ideal = syzstr->
res[
index-1];
140 current_tl = (*syzstr->
Tl)[
index-1];
144 current_ideal = old_ideal;
145 current_repr = old_repr;
148 if (!
idIs0(current_ideal))
153 current_ideal->rank+current_tl);
155 current_ideal->rank);
161 while ((start>0) && (syzstr->
res[
index]->m[start-1]==
NULL)) start--;
175 current_ideal->rank+current_tl);
181 while ((start_ttk>0) && (totake[
index]->
m[start_ttk-1]==
NULL)) start_ttk--;
194 if (current_ideal->m[
i]!=
NULL)
214 totake[
index]->m[
i+start_ttk] =
224 (*syzstr->
Tl)[
index] += current_tl;
247 if (resPairs[
i].p1!=
NULL)
248 assume(resPairs[
i].p1==old_generators->m[resPairs[
i].ind1]);
249 if (resPairs[
i].p2!=
NULL)
250 assume(resPairs[
i].p2==old_generators->m[resPairs[
i].ind2]);
268 while ((syzIndex!=0) && (
res[syzIndex]==
NULL)) syzIndex--;
280 p =
res[syzIndex]->m[
i];
292 PrintS(
"error in the resolvent\n");
308 int index,ideal new_generators,ideal new_repr,
int crit_comp)
310 if (
idIs0(new_generators))
return;
311 ideal old_generators=syzstr->
res[
index];
313 int i=0,
j,
k,kk,og_elem=0,og_idel=
IDELEMS(old_generators),
l=*l_pairs,jj,ll,j1;
321 while ((og_elem<og_idel) && (old_generators->m[og_elem]!=
NULL))
323 if ((
index>0) && (
pGetComp(old_generators->m[og_elem])<=crit_comp))
327 while ((
l>0) && ((*resPairs)[
l-1].
lcm==
NULL))
l--;
328 while ((
i<
IDELEMS(new_generators)) && (new_generators->m[
i]!=
NULL))
331 if (
IDELEMS(old_generators)==og_elem)
341 while ((
j<og_elem) && (old_generators->m[
j]!=
NULL) &&
343 while ((
j<og_elem) && (old_generators->m[
j]!=
NULL) &&
345 for (jj=og_elem;jj>
j;jj--)
347 old_generators->m[jj] = old_generators->m[jj-1];
348 old_repr->m[jj] = old_repr->m[jj-1];
350 old_generators->m[
j] = new_generators->m[
i];
351 new_generators->m[
i] =
NULL;
352 old_repr->m[
j] = new_repr->m[
i];
353 new_repr->m[
i] =
NULL;
355 for (jj=0;jj<*l_pairs;jj++)
359 if ((*resPairs)[jj].ind1>=
j) (*resPairs)[jj].ind1++;
360 if ((*resPairs)[jj].ind2>=
j) (*resPairs)[jj].ind2++;
364 for (jj=og_ini;jj<og_elem;jj++)
369 pLcm(old_generators->m[jj],old_generators->m[
j],
p);
392 p2 =
pMDivide(prs[j1],old_generators->m[j1]);
397 while ((ti<
l) && (((*resPairs)[ti].ind1!=j1)|| ((*resPairs)[ti].ind2!=jj))) ti++;
417 for (jj=og_ini;jj<og_elem;jj++)
424 for (ll=0;ll<*l_pairs;ll++)
426 temp[ll].p = (*resPairs)[ll].p;
427 temp[ll].p1 = (*resPairs)[ll].p1;
428 temp[ll].p2 = (*resPairs)[ll].p2;
429 temp[ll].syz = (*resPairs)[ll].syz;
430 temp[ll].lcm = (*resPairs)[ll].lcm;
431 temp[ll].ind1 = (*resPairs)[ll].ind1;
432 temp[ll].ind2 = (*resPairs)[ll].ind2;
433 temp[ll].syzind = (*resPairs)[ll].syzind;
434 temp[ll].order = (*resPairs)[ll].order;
435 temp[ll].isNotMinimal = (*resPairs)[ll].isNotMinimal;
444 tso.p1 = old_generators->m[jj];
445 tso.p2 = old_generators->m[
j];
449 tso.isNotMinimal =
NULL;
454Print(
"erzeuge Paar im Modul %d,%d mit: \n",
index,tso.order);
515 int * gen_length,
int * secgen_length,
int * tored_length)
517 int i=
IDELEMS(red)-1,num_mon,num_tail;
521 while ((
i>0) && (red->m[
i-1]==
NULL))
i--;
528 num_tail = *tored_length-1;
564 *tored_length = num_mon+num_tail;
587 if ((
j>=
i) || (
p==
NULL))
break;
610 int * gen_length,
int* orp_l,
int * tored_l,
int * syzred_l)
612 int num_mon,num_tail,syz_l;
622 num_tail = *tored_l-1;
626 dummy = syRedSyzPair(syzstr,
index,gen_length,orp_l);
639 *tored_l = num_mon+num_tail;
652 ideal new_repr,
int * ogm_l,
int * orp_l)
654 SObject tso = resPairs[itso];
656 ideal old_generators=syzstr->
res[
index];
658 int og_idel=
IDELEMS(old_generators),ng_place=
IDELEMS(new_generators);
668 int syz_place=
IDELEMS(syzygies);
670 while ((syz_place>0) && (syzygies->m[syz_place-1]==
NULL)) syz_place--;
671 while ((ng_place>0) && (new_generators->m[ng_place-1]==
NULL)) ng_place--;
672 while ((og_idel>0) && (old_generators->m[og_idel-1]==
NULL)) og_idel--;
673 assume (tso.ind1<og_idel);
674 assume (tso.ind2<og_idel);
675 assume (tso.ind1!=tso.ind2);
676 assume (tso.p1 == old_generators->m[tso.ind1]);
677 assume (tso.p2 == old_generators->m[tso.ind2]);
678 tso.p1 = old_generators->m[tso.ind1];
679 tso.p2 = old_generators->m[tso.ind2];
680 if ((tso.p1!=
NULL) && (tso.p2!=
NULL))
688 toReplace = tso.ind1+1;
690 toReplace = tso.ind2+1;
694 if ((
index==0) && (old_generators->rank==1) &&
706 pp =
pCopy(old_repr->m[tso.ind2]);
716 pp =
pCopy(old_repr->m[tso.ind1]);
739 tso.syz =
pCopy(old_repr->m[tso.ind1]);
748 p =
pCopy(old_repr->m[tso.ind2]);
755 tso.syz =
pAdd(
p,tso.syz);
767Print(
"reduziere Paar im Module %d mit: \n",
index);
775 if ((tso.p!=
NULL) && (!deleteP))
787 if (
j>=ng_place)
break;
844 PrintS(
"Shit happens!\n");
847Print(
"erhalte Paar im Module %d mit: \n",
index);
855 pDelete(&old_generators->m[toReplace-1]);
856 pDelete(&old_repr->m[toReplace-1]);
857 for (
i=toReplace-1;
i<og_idel-1;
i++)
859 old_generators->m[
i] = old_generators->m[
i+1];
860 old_repr->m[
i] = old_repr->m[
i+1];
862 old_generators->m[og_idel-1] =
NULL;
863 old_repr->m[og_idel-1] =
NULL;
864 for (
i=itso+1;
i<
l;
i++)
868 if ((resPairs[
i].ind1==toReplace-1)||(resPairs[
i].ind2==toReplace-1))
872 if (resPairs[
i].ind1>=toReplace)
873 (resPairs[
i].ind1)--;
874 if (resPairs[
i].ind2>=toReplace)
875 (resPairs[
i].ind2)--;
885 if (ng_place>=
IDELEMS(new_generators))
899 new_generators->m[ng_place] = tso.p;
901 new_repr->m[ng_place] = tso.syz;
915 else if (
pGetComp(tso.syz)<=crit_comp)
921 if (syz_place>=
IDELEMS(syzygies))
926 syzygies->m[syz_place] = tso.syz;
928 pNorm(syzygies->m[syz_place]);
931 resPairs[itso] = tso;
940 ideal new_generators,ideal new_repr,
int crit_comp,
syStrategy syzstr,
944 int i,
j,actdeg=resPairs[0].order;
964 while (resPairs[
i].order==actdeg)
968 new_generators, new_repr,ogm_l,orp_l);
975 if (!
idIs0(new_generators))
977 else if (resPairs[0].
lcm==
NULL)
984 actdeg = resPairs[0].order;
1002 int index,
int next_comp)
1005 int l_pairs=2*og_idel;
1007 if ((
idIs0(new_generators)) || (new_generators->m[0]==
NULL))
1009 WerrorS(
"Hier ist was faul!\n");
1016 new_generators,new_repr,next_comp);
1017 if (
redPairs(resPairs,l_pairs,syzygies, new_generators,new_repr,
1018 next_comp,syzstr,
index))
break;
1028 int * g_l,
int * p_l,
int crit_comp)
1057 int index,ideal add_generators,ideal ,ideal ,
1058 ideal ,
int ,
int* first_new)
1060 int i=*first_new,
l=*l_pairs,
j,ll,j1,add_idel=
IDELEMS(add_generators);
1061 ideal
pairs=
idInit(add_idel,add_generators->rank);
1067 while ((
l>0) && ((*resPairs)[
l-1].
lcm==
NULL))
l--;
1068 while ((
i<add_idel) && (add_generators->m[
i]!=
NULL))
1075 pLcm(add_generators->m[
j],add_generators->m[
i],
p);
1098 p2 =
pMDivide(prs[j1],add_generators->m[j1]);
1103 while ((ti<
l) && (((*resPairs)[ti].ind1!=j1)|| ((*resPairs)[ti].ind2!=
j))) ti++;
1130 for (ll=0;ll<*l_pairs;ll++)
1132 temp[ll].p = (*resPairs)[ll].p;
1133 temp[ll].p1 = (*resPairs)[ll].p1;
1134 temp[ll].p2 = (*resPairs)[ll].p2;
1135 temp[ll].syz = (*resPairs)[ll].syz;
1136 temp[ll].lcm = (*resPairs)[ll].lcm;
1137 temp[ll].ind1 = (*resPairs)[ll].ind1;
1138 temp[ll].ind2 = (*resPairs)[ll].ind2;
1139 temp[ll].syzind = (*resPairs)[ll].syzind;
1140 temp[ll].order = (*resPairs)[ll].order;
1141 temp[ll].isNotMinimal = (*resPairs)[ll].isNotMinimal;
1150 tso.p1 = add_generators->m[
j];
1151 tso.p2 = add_generators->m[
i];
1155 tso.isNotMinimal =
NULL;
1160Print(
"erzeuge Paar im Modul %d,%d mit: \n",
index,tso.order);
1183 ideal new_generators, ideal new_repr,
int * next_place_add,
int ** g_l,
1187 syStrategy syzstr,
int ,ideal add_generators, ideal add_repr,
1188 ideal new_generators, ideal new_repr,
int * next_place_add,
int ** g_l,
1192 SObject tso = resPairs[itso];
1194 int ng_place=
IDELEMS(new_generators);
1202 assume (tso.ind1<*next_place_add);
1203 assume (tso.ind2<*next_place_add);
1204 assume (tso.ind1!=tso.ind2);
1205 assume (tso.p1 == add_generators->m[tso.ind1]);
1206 assume (tso.p2 == add_generators->m[tso.ind2]);
1207 tso.p1 = add_generators->m[tso.ind1];
1208 tso.p2 = add_generators->m[tso.ind2];
1209 if ((tso.p1!=
NULL) && (tso.p2!=
NULL))
1219 poly p1=
pCopy(tso.p1);
1220 poly p2=
pCopy(tso.p2);
1221 poly pp1,pp2,tp1,tp2;
1222 poly sp1=
pCopy(add_repr->m[tso.ind1]),sp2=
pCopy(add_repr->m[tso.ind2]);
1253 sp1 =
pMult(p2,sp1);
1254 sp2 =
pMult(p1,sp2);
1255 tso.p =
pSub(tp1,tp2);
1256 tso.syz =
pSub(sp1,sp2);
1264 tso.syz =
pCopy(add_repr->m[tso.ind1]);
1265 poly tt =
pMDivide(tso.lcm,tso.p1);
1273 p =
pCopy(add_repr->m[tso.ind2]);
1280 tso.syz =
pAdd(
p,tso.syz);
1284Print(
"reduziere Paar im Module %d mit: \n",
index);
1299 if (
j>=*next_place_add)
break;
1319 PrintS(
"Shit happens!\n");
1322Print(
"erhalte Paar im Module %d mit: \n",
index);
1340 if (*next_place_add>=
IDELEMS(add_generators))
1345 (
IDELEMS(add_generators)+16)*
sizeof(
int));
1346 IDELEMS(add_generators) += 16;
1350 add_generators->m[*next_place_add] = tso.p;
1351 add_repr->m[*next_place_add] = tso.syz;
1352 (*g_l)[*next_place_add] = tso.length;
1353 (*next_place_add)++;
1357 while ((ng_place>0) && (new_generators->m[ng_place-1]==
NULL) &&
1358 (new_repr->m[ng_place-1]==
NULL)) ng_place--;
1359 if (ng_place>=
IDELEMS(new_generators))
1362 IDELEMS(new_generators) += 16;
1366 new_generators->m[ng_place] = tso.p;
1367 new_repr->m[ng_place] = tso.syz;
1371 resPairs[itso] = tso;
1379 int index,ideal add_generators,ideal add_repr,ideal new_generators,
1380 ideal new_repr,
int crit_comp,
int * red_deg,
int * next_place_add,
int **g_l,
1389 while ((
i<l_pairs) && (resPairs[
i].order==*red_deg))
1396 new_generators, new_repr,next_place_add,g_l,deg_soc);
1403 *red_deg = resPairs[0].order;
1416 ideal new_generators, ideal new_repr, ideal add_generators,
1420 int i=0,
j,next_new_el;
1421 int idel_temp=
IDELEMS(temp_generators);
1423 int p_length,red_deg,l_pairs=
IDELEMS(add_generators);
1436 next_place_add =
IDELEMS(add_generators);
1437 while ((next_place_add>0) && (add_generators->m[next_place_add-1]==
NULL))
1440 next_new_el = next_place_add;
1442 while ((
i<idel_temp) && (temp_generators->m[
i]!=
NULL))
1445#ifdef USE_REGULARITY
1450 while ((
i<idel_temp) && (temp_generators->m[
i]!=
NULL))
1459 while ((
i<idel_temp) && (
p_FDeg(temp_generators->m[
i],
currRing)==next_deg))
1461 next_p = temp_generators->m[
i];
1462 temp_generators->m[
i] =
NULL;
1463 next_p =
normalize(next_p,add_generators,syzstr,gen_length,&p_length,
1475 gen_length,secgen_length,&p_length);
1478 if (next_place_add>=
IDELEMS(add_generators))
1483 (
IDELEMS(add_generators)+16)*
sizeof(
int));
1484 IDELEMS(add_generators) += 16;
1487 add_generators->m[next_place_add] = next_p;
1504#ifdef WITH_SCHREYER_ORD
1505 add_repr->m[next_place_add] =
pHead(add_generators->m[next_place_add]);
1507 add_repr->m[next_place_add] =
pOne();
1512 gen_length[next_place_add] = p_length;
1523 if ((next_place_add>next_new_el) || (next_deg<0))
1527 while (pairs_left && ((next_deg<0) || (red_deg<= next_deg)))
1530 add_repr,new_generators,new_repr,crit_comp,&next_new_el);
1532 add_repr,new_generators,new_repr,crit_comp,&red_deg,&next_place_add,&gen_length,
1549 ideal old_generators= syzstr->
res[
index];
1551 int i,
j=0,ii=
IDELEMS(old_generators)-1,dummy;
1558 if (old_generators->m[
i]!=
NULL)
1565 if (new_generators->m[
i]!=
NULL)
1574 if ((
j>=ii) || (
p==
NULL))
break;
1575 if ((old_generators->m[
j]!=
NULL) &&
1599 new_repr->m[
j] =
NULL;
1611 ideal temp_generators =
idInit(
IDELEMS(new_generators),new_generators->rank);
1613 ideal add_generators =
idInit(
IDELEMS(new_generators),new_generators->rank);
1616 int j,jj,
k,deg_p,idel_temp=
IDELEMS(temp_generators);
1621 if (new_generators->m[
j]!=
NULL)
1623 p = new_generators->m[
j];
1624 new_generators->m[
j] =
NULL;
1632 if (deg_p<min_deg) min_deg = deg_p;
1635 while ((
k<idel_temp) && (temp_generators->m[
k]!=
NULL) &&
1637 for (jj=idel_temp-1;jj>
k;jj--)
1639 temp_generators->m[jj] = temp_generators->m[jj-1];
1641 temp_generators->m[
k] =
p;
1646 add_generators,add_repr,syzstr,
index,next_comp,totake);
1650 while ((jj>0) && (add_generators->m[jj-1]==
NULL)) jj--;
1662 syzstr->
res[
index]->m[
j+
k] = add_generators->m[
k];
1664 add_generators->m[
k] =
NULL;
1665 add_repr->m[
k] =
NULL;
1701 int ng_idel=
IDELEMS(new_generators);
1702 ideal new_repr =
idInit(ng_idel, crit_comp+ng_idel);
1712 for (
i=0;
i<ng_idel;
i++)
1714 if (new_generators->m[
i]!=
NULL)
1734#ifdef WITH_SCHREYER_ORD
1735 new_repr->m[
i] =
pHead(new_generators->m[
i]);
1737 new_repr->m[
i] =
pOne();
1746PrintS(
"Add new generators:\n");
1748PrintS(
"with representations:\n");
1769 int crit_comp,t_comp,next_deg,old_tl;
1770 ideal temp=
NULL,old_ideal,old_repr;
1785 if (origR!=syzstr->
syRing)
1795 if (temp->m[0]!=
NULL)
1803 if (jj>maxdeg) maxdeg = jj;
1805 while (!
idIs0(temp))
1810 if (temp->m[
j]!=
NULL)
1820 temp1->m[
k] = temp->m[
j];
1831#ifdef USE_REGULARITY
1832 int last_generator=
IDELEMS(temp)-1;
1833 while ((last_generator>=0) && (temp->m[last_generator]==
NULL))
1843 ideal temp_gens,old_std;
1845 if (rk_arg>1) syzstr->
res[0]->rank = rk_arg;
1851 if (temp->m[
i]!=
NULL)
1853 new_generators->m[0] =
kNF(syzstr->
res[0],
currRing->qideal,temp->m[
i]);
1855 pNorm(new_generators->m[0]);
1857 next_gen =
pCopy(new_generators->m[0]);
1859 if (!
idIs0(new_generators))
1868 old_tl = (*syzstr->
Tl)[0];
1871 t_comp = (*syzstr->
Tl)[
index];
1872 if (
index==0) crit_comp = t_comp;
1882#ifdef USE_REGULARITY
1883 if ((
index==0) && (!isRegular) && (
i==last_generator))
1915 PrintS(
"\n not regular\n");
1926 new_generators = temp_gens;
1930 if (
idIs0(new_generators))
break;
1933 if (!
idIs0(new_generators))
1937 if (new_generators->m[
j]!=
NULL)
1940 new_generators->m[
j] =
NULL;
1958 Print(
"The %d-th syzygy-module is now:\n",
index);
2000 PrintS(
"Da ist was faul!!!\n");
2001 Print(
"Aber: Regularitaet %d, Grad %ld\n",
2021 if (origR!=syzstr->
syRing)
static int si_max(const int a, const int b)
static FORCE_INLINE number n_Gcd(number a, number b, const coeffs r)
in Z: return the gcd of 'a' and 'b' in Z/nZ, Z/2^kZ: computed as in the case Z in Z/pZ,...
void WerrorS(const char *s)
intvec * hFirstSeries(ideal A, intvec *module_w, ideal Q, intvec *wdegree)
#define idDelete(H)
delete an ideal
BOOLEAN idIs0(ideal h)
returns true if h is the zero ideal
poly initial(const poly p, const ring r, const gfan::ZVector &w)
Returns the initial form of p with respect to w.
static BOOLEAN length(leftv result, leftv arg)
KINLINE poly ksOldCreateSpoly(poly p1, poly p2, poly spNoether, ring r)
void kBucketClear(kBucket_pt bucket, poly *p, int *length)
void kBucket_Minus_m_Mult_p(kBucket_pt bucket, poly m, poly p, int *l, poly spNoether)
Bpoly == Bpoly - m*p; where m is a monom Does not destroy p and m assume (*l <= 0 || pLength(p) == *l...
void kBucketDestroy(kBucket_pt *bucket_pt)
void kBucketInit(kBucket_pt bucket, poly lm, int length)
kBucket_pt kBucketCreate(const ring bucket_ring)
Creation/Destruction of buckets.
number kBucketPolyRed(kBucket_pt bucket, poly p1, int l1, poly spNoether)
const poly kBucketGetLm(kBucket_pt bucket)
poly kNF(ideal F, ideal Q, poly p, int syzComp, int lazyReduce)
int lcm(unsigned long *l, unsigned long *a, unsigned long *b, unsigned long p, int dega, int degb)
static number & pGetCoeff(poly p)
return an alias to the leading coefficient of p assumes that p != NULL NOTE: not copy
#define omFreeSize(addr, size)
#define omRealloc0Size(addr, o_size, size)
static int index(p_Length length, p_Ord ord)
void p_Shift(poly *p, int i, const ring r)
shifts components of the vector p by i
void pEnlargeSet(poly **p, int l, int increment)
static int pLength(poly a)
static long p_FDeg(const poly p, const ring r)
#define __p_Mult_nn(p, n, r)
void rChangeCurrRing(ring r)
VAR ring currRing
Widely used global variable which specifies the current polynomial ring for Singular interpreter and ...
Compatibility layer for legacy polynomial operations (over currRing)
#define pHead(p)
returns newly allocated copy of Lm(p), coef is copied, next=NULL, p might be NULL
#define pGetComp(p)
Component.
#define pSetCoeff(p, n)
deletes old coeff before setting the new one
#define pLmDelete(p)
assume p != NULL, deletes Lm(p)->coef and Lm(p)
#define pGetExp(p, i)
Exponent.
#define pSetmComp(p)
TODO:
#define pDivisibleBy(a, b)
returns TRUE, if leading monom of a divides leading monom of b i.e., if there exists a expvector c > ...
#define pCopy(p)
return a copy of the poly
#define pLmDivisibleByNoComp(a, b)
like pLmDivisibleBy, does not check components
ideal idrMoveR(ideal &id, ring src_r, ring dest_r)
ideal idrCopyR(ideal id, ring src_r, ring dest_r)
ideal idrMoveR_NoSort(ideal &id, ring src_r, ring dest_r)
void PrintS(const char *s)
ring rAssure_C_dp(const ring r)
void rDelete(ring r)
unconditionally deletes fields in r
ring rAssure_dp_C(const ring r)
ideal idInit(int idsize, int rank)
initialise an ideal / module
ideal id_Head(ideal h, const ring r)
returns the ideals of initial terms
long id_RankFreeModule(ideal s, ring lmRing, ring tailRing)
return the maximal component number found in any polynomial in s
void idShow(const ideal id, const ring lmRing, const ring tailRing, const int debugPrint)
void idSkipZeroes(ideal ide)
gives an ideal/module the minimal possible size
void id_Shift(ideal M, int s, const ring r)
static poly normalize(poly next_p, ideal add_generators, syStrategy syzstr, int *g_l, int *p_l, int crit_comp)
static ideal syAppendSyz(ideal new_generators, syStrategy syzstr, int index, int crit_comp, resolvente totake)
static void syCreateRegularExtension(syStrategy syzstr, ideal old_ideal, ideal old_repr, int old_tl, poly next_generator, resolvente totake)
void syReorder_Kosz(syStrategy syzstr)
static BOOLEAN syRedSyz(kBucket_pt bucket, ideal red, int crit_comp, int *g_l)
static ideal normalizeOldPart(ideal new_generators, ideal new_repr, syStrategy syzstr, int index, int)
static void redOnePairHIndex(SSet resPairs, int itso, int crit_comp, syStrategy syzstr, int, ideal add_generators, ideal add_repr, ideal new_generators, ideal new_repr, int *next_place_add, int **g_l, poly deg_soc)
syStrategy syKosz(ideal arg, int *length)
static void syTestPairs(SSet resPairs, int length, ideal old_generators)
static void redOnePair(SSet resPairs, int itso, int l, ideal syzygies, int crit_comp, syStrategy syzstr, int index, ideal new_generators, ideal new_repr, int *ogm_l, int *orp_l)
static BOOLEAN reducePairsHIndex(SSet resPairs, int l_pairs, syStrategy syzstr, int index, ideal add_generators, ideal add_repr, ideal new_generators, ideal new_repr, int crit_comp, int *red_deg, int *next_place_add, int **g_l, resolvente totake)
static BOOLEAN redPairs(SSet resPairs, int l_pairs, ideal syzygies, ideal new_generators, ideal new_repr, int crit_comp, syStrategy syzstr, int index)
static void updatePairs(SSet *resPairs, int *l_pairs, syStrategy syzstr, int index, ideal new_generators, ideal new_repr, int crit_comp)
static void updatePairsHIndex(SSet *resPairs, int *l_pairs, syStrategy, int index, ideal add_generators, ideal, ideal, ideal, int, int *first_new)
static void procedeNextGenerators(ideal temp_generators, ideal, ideal new_generators, ideal new_repr, ideal add_generators, ideal add_repr, syStrategy syzstr, int index, int crit_comp, resolvente totake)
void sySPRedSyz_Kosz(syStrategy syzstr, poly redWith, poly syz, poly q=NULL, int l_syz=-1)
static ideal kosz_std(ideal new_generators, ideal new_repr, syStrategy syzstr, int index, int next_comp)
static poly syRedTailSyz(poly tored, ideal red, ideal sec_red, int crit_comp, syStrategy syzstr, int *gen_length, int *secgen_length, int *tored_length)
static BOOLEAN syIsRegular(ideal old_ideal, ideal new_ideal, int deg)
static ideal kosz_ext(ideal new_generators, ideal new_repr, syStrategy syzstr, int index, int next_comp, resolvente totake)
intvec * syBetti(resolvente res, int length, int *regularity, intvec *weights, BOOLEAN tomin, int *row_shift)
resolvente sySchreyerResolvente(ideal arg, int maxlength, int *length, BOOLEAN isMonomial=FALSE, BOOLEAN notReplace=FALSE)
void syCompactifyPairSet(SSet sPairs, int sPlength, int first)
void syEnterPair(syStrategy syzstr, SObject *so, int *sPlength, int index)
void syInitializePair(SObject *so)
void syDeletePair(SObject *so)