1419 int i=0,
j,next_new_el;
1420 int idel_temp=
IDELEMS(temp_generators);
1422 int p_length,red_deg,l_pairs=
IDELEMS(add_generators);
1435 next_place_add =
IDELEMS(add_generators);
1436 while ((next_place_add>0) && (add_generators->m[next_place_add-1]==
NULL))
1439 next_new_el = next_place_add;
1441 while ((
i<idel_temp) && (temp_generators->m[
i]!=
NULL))
1444 #ifdef USE_REGULARITY 1449 while ((
i<idel_temp) && (temp_generators->m[
i]!=
NULL))
1458 while ((
i<idel_temp) && (
p_FDeg(temp_generators->m[
i],
currRing)==next_deg))
1460 next_p = temp_generators->m[
i];
1461 temp_generators->m[
i] =
NULL;
1462 next_p =
normalize(next_p,add_generators,syzstr,gen_length,&p_length,
1474 gen_length,secgen_length,&p_length);
1477 if (next_place_add>=
IDELEMS(add_generators))
1482 (
IDELEMS(add_generators)+16)*
sizeof(
int));
1483 IDELEMS(add_generators) += 16;
1486 add_generators->m[next_place_add] = next_p;
1503 #ifdef WITH_SCHREYER_ORD 1504 add_repr->m[next_place_add] =
pHead(add_generators->m[next_place_add]);
1506 add_repr->m[next_place_add] =
pOne();
1511 gen_length[next_place_add] = p_length;
1522 if ((next_place_add>next_new_el) || (next_deg<0))
1526 while (pairs_left && ((next_deg<0) || (red_deg<= next_deg)))
1529 add_repr,new_generators,new_repr,crit_comp,&next_new_el);
1531 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 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)
static unsigned pLength(poly a)
#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)
ring currRing
Widely used global variable which specifies the current polynomial ring for Singular interpreter and ...
BOOLEAN idIs0(ideal h)
returns true if h is the zero ideal
#define pCopy(p)
return a copy of the poly