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))
1452 pDelete(&temp_generators->m[i]);
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;
1493 (*syzstr->
Tl)[index]+16-
IDELEMS(totake[index]));
1494 for (
j=
IDELEMS(totake[index]);
j<(*syzstr->
Tl)[index]+16;
j++)
1495 totake[index]->
m[
j] =
NULL;
1496 IDELEMS(totake[index]) = (*syzstr->
Tl)[index]+16;
1499 totake[
index]->m[(*syzstr->
Tl)[index]] =
pCopy(next_p);
1501 totake[
index]->m[(*syzstr->
Tl)[index]] =
pHead(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();
1509 ((*syzstr->
Tl)[index])++;
1510 pSetComp(add_repr->m[next_place_add],(*syzstr->
Tl)[index]);
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,
#define pSetmComp(p)
TODO:
#define omRealloc0Size(addr, o_size, size)
static poly normalize(poly next_p, ideal add_generators, syStrategy syzstr, int *g_l, int *p_l, int crit_comp)
static void updatePairsHIndex(SSet *resPairs, int *l_pairs, syStrategy, int index, ideal add_generators, ideal, ideal, ideal, int, int *first_new)
#define omFreeSize(addr, size)
static number & pGetCoeff(poly p)
return an alias to the leading coefficient of p assumes that p != NULL NOTE: not copy ...
#define pGetComp(p)
Component.
static int pLength(poly a)
ring currRing
Widely used global variable which specifies the current polynomial ring for Singular interpreter and ...
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 long p_FDeg(const poly p, const ring r)
#define pHead(p)
returns newly allocated copy of Lm(p), coef is copied, next=NULL, p might be NULL ...
static int index(p_Length length, p_Ord ord)
static poly syRedTailSyz(poly tored, ideal red, ideal sec_red, int crit_comp, syStrategy syzstr, int *gen_length, int *secgen_length, int *tored_length)
ideal idInit(int idsize, int rank)
initialise an ideal / module
void pEnlargeSet(poly **p, int l, int increment)
void pNorm(poly p, const ring R=currRing)
BOOLEAN idIs0(ideal h)
returns true if h is the zero ideal
#define pCopy(p)
return a copy of the poly