janet.h File Reference

#include <kernel/structs.h>

Go to the source code of this file.

Data Structures
struct	Poly
struct	jList

Macros
#define	GCM(sz)   omAlloc((sz))
#define	GCMA(sz)   omAlloc((sz))
#define	GCF(x)   omFree((x))
#define	ListNode   struct LISTNODE
#define	TreeM   struct TREEM
#define	NodeM   struct NODEM

Typedefs
typedef void(*	poly_function )(Poly *)
typedef ListNode *	LCI
typedef ListNode **	LI

Functions
Poly *	FindMinList (jList *)
void	DestroyTree (NodeM *)
NodeM *	create ()
void	ControlProlong (Poly *)
Poly *	is_div_ (TreeM *root, poly item)
void	insert_ (TreeM **tree, Poly *item)
Poly *	NewPoly (poly p=NULL)
void	DestroyPoly ()
void	NFL (Poly *, TreeM *)
void	PNF (Poly *, TreeM *)
void	ClearProl (Poly *x, int i)
void	InitProl (Poly *p)
void	InitHistory (Poly *p)
Poly *	is_present (jList *, poly)
int	GetMult (Poly *, int)
int	GB_length ()
void	InsertInList (jList *, Poly *)
void	ForEachPNF (jList *, int)
void	ClearMult (Poly *, int)
void	ProlVar (Poly *, int)
void	SetMult (Poly *, int)
void	InitLead (Poly *)
void	InsertInCount (jList *, Poly *)
int	GetProl (Poly *, int)
void	SetProl (Poly *, int)
int	ProlCompare (Poly *, Poly *)
int	ValidatePoly (Poly *, TreeM *)
int	ListGreatMoveDegree (jList *, jList *, poly)
int	ListGreatMoveOrder (jList *, jList *, poly)
void	ForEachControlProlong (jList *)
void	NFListQ ()
int	CountList (jList *)
void	DestroyList (jList *)
int	ReducePoly (Poly *x, Poly *y)
int	ReducePolyLead (Poly *x, Poly *y)
void	Define (TreeM **G)
ListNode *	CreateListNode (Poly *x)
void	DestroyListNode (ListNode *x)
void	DestroyFreeNodes ()
int	ComputeBasis (jList *, jList *)
void	Initialization (char *)

Variables
ListNode Poly *	info
ListNode *	next
	NodeM
Poly *	ended
	TreeM

---

Data Structure Documentation

struct Poly

Definition at line 14 of file janet.h.

Data Fields
int	changed
poly	history
poly	lead
char *	mult
int	prolonged
poly	root
kBucket_pt	root_b
int	root_l

struct jList

Definition at line 34 of file janet.h.

Data Fields
ListNode *	root

Macro Definition Documentation

#define GCF

(

x

)

omFree((x))

Definition at line 8 of file janet.h.

#define GCM

(

sz

)

omAlloc((sz))

Definition at line 6 of file janet.h.

#define GCMA

(

sz

)

omAlloc((sz))

Definition at line 7 of file janet.h.

#define ListNode   struct LISTNODE

Definition at line 10 of file janet.h.

#define NodeM   struct NODEM

Definition at line 12 of file janet.h.

#define TreeM   struct TREEM

Definition at line 11 of file janet.h.

Typedef Documentation

typedef ListNode* LCI

Definition at line 48 of file janet.h.

typedef ListNode** LI

Definition at line 51 of file janet.h.

typedef void(* poly_function)(Poly *)

Definition at line 26 of file janet.h.

Function Documentation

void ClearMult	(	Poly *	,
		int
	)

Definition at line 438 of file janet.cc.

{
  x->mult[i/8] &= ~Mask[i%8];
}

void ClearProl	(	Poly *	x,
		int	i
	)

Definition at line 453 of file janet.cc.

{
  (x->mult+offset)[i/8] &= ~Mask[i%8];
}

int ComputeBasis	(	jList *	,
		jList *
	)

Definition at line 1007 of file janet.cc.

{
  // int gb_l,i,ret_value=1;

  T=_lT; Q=_lQ;

//  Debug();

  while((h=FindMinList(Q))!=NULL)
  {
//        PrintS("New element\n");
//  Debug();

        if (!degree_compatible)
        {
          if (!ValidatePoly(h,G))
          {
            DestroyPoly(h);
            continue;
          }

          h->changed=0;

          NFL(h,G);

          if (!h->root)
          {
            DestroyPoly(h);
            continue;
          }
        }

        if (h->root)
        {
          if (pIsConstant(h->root))
          {
            WarnS("Constant in basis\n");
            return 0;
          }

          if (h->changed && ListGreatMove(T,Q,h->root))
          {
//      PrintS("<-\n");
            DestroyTree(G->root);
            G->root=create();
            T2G();
          }
        }

//  PrintS("PNF\n");
        PNF(h,G);
//        Print("{%d}\n",pow_(h->root));
        insert_(&G,h);
        InsertInList(T,h);

//  PrintS("For each PNF\n");
        if (degree_compatible)
            ForEachPNF(T,pow_(h->root));

//  PrintS("Control of prolongations\n");
        if (h->changed)
            ForEachControlProlong(T);
        else
            ControlProlong(h);

//  Debug();

//  PrintS("NFListQ\n");
        if (degree_compatible)
            NFListQ();
//Debug();
    }

//    gb_l=GB_length();

    Print("Length of Janet basis: %d\n",CountList(T));
//    Print("Length of Groebner basis:    %d\n",gb_l);

    DestroyTree(G->root);
    GCF(G);
    DestroyFreeNodes();

    return 1;
}

void ControlProlong

(

Poly *

)

Definition at line 399 of file janet.cc.

{
  for (int i = 0; i< offset; i++)
  {
    (x->mult+offset)[i]&=~((x->mult)[i]);
//    if (!GetMult(x,i) && !GetProl(x,i))
//      ProlVar(x,i);
  }
}

int CountList

(

jList *

)

Definition at line 624 of file janet.cc.

{
  int i=0;
  LCI y=Q->root;

  while(y)
  {
    i++;
    y=y->next;
  }

  return i;
}

NodeM* create

(

)

Definition at line 764 of file janet.cc.

{
  NodeM *y;

  if (FreeNodes == NULL)
  {
    y=(NodeM *)GCM(sizeof(NodeM));
  }
  else
  {
    y=FreeNodes;
    FreeNodes=FreeNodes->left;
  }

  y->left=y->right=NULL;
  y->ended=NULL;
  return y;
}

ListNode* CreateListNode

(

Poly *

x

)

Definition at line 515 of file janet.cc.

{
  ListNode* ret=(ListNode *)GCM(sizeof(ListNode));
  ret->info=x;
  ret->next=NULL;
  return ret;
}

void Define

(

TreeM **

G

)

Definition at line 822 of file janet.cc.

{
  *G=(TreeM *)GCM(sizeof(TreeM));
  (*G)->root=create();
}

void DestroyFreeNodes

(

)

Definition at line 783 of file janet.cc.

{
  NodeM *y;

  while((y=FreeNodes)!=NULL)
  {
    FreeNodes=FreeNodes->left;
    GCF(y);
  }
}

void DestroyList

(

jList *

)

Definition at line 721 of file janet.cc.

{
  LCI y=x->root,z;

  while(y)
  {
    z=y->next;
    DestroyPoly(y->info);
    GCF(y);
    y=z;
  }

  GCF(x);
}

void DestroyListNode

(

ListNode *

x

)

Definition at line 509 of file janet.cc.

{
  DestroyPoly(x->info);
  GCF(x);
}

void DestroyPoly

(

)

void DestroyTree

(

NodeM *

)

Definition at line 811 of file janet.cc.

{
  if (G)
  {
    DestroyTree(G->left);
    DestroyTree(G->right);
    G->left=FreeNodes;
    FreeNodes=G;
  }
}

Poly* FindMinList

(

jList *

)

Definition at line 524 of file janet.cc.

{
  LI min=&(L->root);
  LI l;
  LCI xl;
  Poly *x;

  if (degree_compatible)
  {
    while ((*min) && ((*min)->info->root == NULL))
      min=&((*min)->next);
  }

  if (!(*min)) return NULL;

  l=&((*min)->next);

  while (*l)
  {
    if ((*l)->info->root != NULL)
    {
      if (ProlCompare((*l)->info,(*min)->info))
        min=l;
    }

    l=&((*l)->next);
  }
  x=(*min)->info;
  xl=*min;
  *min=(*min)->next;
  GCF(xl);

  return x;
}

void ForEachControlProlong

(

jList *

)

Definition at line 710 of file janet.cc.

{
  LCI y=x->root;

  while(y)
  {
    ControlProlong(y->info);
    y=y->next;
  }
}

void ForEachPNF	(	jList *	,
		int
	)

Definition at line 699 of file janet.cc.

{
  LCI y=x->root;

  while(y)
  {
    if (pow_(y->info->root) == i) PNF(y->info,G);
    y=y->next;
  }
}

int GB_length

(

)

Definition at line 747 of file janet.cc.

{
  LCI iT=T->root;
  int l=0;

  while(iT)
  {
    if (pow_(iT->info->lead) == pow_(iT->info->history))
      ++l;
    iT=iT->next;
  }

  return l;
}

int GetMult	(	Poly *	,
		int
	)

Definition at line 428 of file janet.cc.

{
  return x->mult[i/8] & Mask[i%8];
}

int GetProl	(	Poly *	,
		int
	)

Definition at line 443 of file janet.cc.

{
  return (x->mult+offset)[i/8] & Mask[i%8];
}

void InitHistory

(

Poly *

p

)

Definition at line 409 of file janet.cc.

{
  if (p->history) pLmFree(&p->history);
  p->history=pLmInit(p->root);
  p->changed=0;
}

void Initialization

(

char *

)

Definition at line 959 of file janet.cc.

{
  offset=(currRing->N % 8 == 0) ? (currRing->N/8)*8 : (currRing->N/8+1)*8;
  if (strstr(Ord,"dp\0") || strstr(Ord,"Dp\0"))
  {
    degree_compatible=1;
    jDeg=p_Deg;
    ListGreatMove=ListGreatMoveDegree;
  }
  else
  {
    degree_compatible=0;
    jDeg=p_Totaldegree;
    ListGreatMove=ListGreatMoveOrder;
  }

  Define(&G);
}

void InitLead

(

Poly *

)

Definition at line 416 of file janet.cc.

{
  if (p->lead) pLmDelete(&p->lead);
  p->lead=pLmInit(p->root);
  p->prolonged=-1;
}

void InitProl

(

Poly *

p

)

Definition at line 423 of file janet.cc.

{
  memset(p->mult+offset,0,sizeof(char)*offset);
}

void insert_	(	TreeM **	tree,
		Poly *	item
	)

Definition at line 919 of file janet.cc.

{
 int power_tmp,i,i_con=currRing->N-1;
 NodeM *curr=(*tree)->root;

 for ( ; (i_con>=0) && !pGetExp(item->root,i_con+1) ; i_con--)
  SetMult(item,i_con);

 for (i = 0; i<= i_con; i++)
 //<=
 {
  power_tmp=pGetExp(item->root,i+1);

  ClearMult(item,i);

  while (power_tmp)
  {
   if (!curr->left)
   {
     SetMult(item,i);
     ClearMultiplicative(curr->right,i);
     curr->left=create();
   };
   curr=curr->left;
   power_tmp--;
  };

  if (i<i_con)
  {
   if (!curr->left) SetMult(item,i);
   if (!curr->right) curr->right=create();
   curr=curr->right;

   ProlVar(item,i);
  }
 }

 curr->ended=item;
}

void InsertInCount	(	jList *	,
		Poly *
	)

Definition at line 578 of file janet.cc.

{
  ListNode *ins;
  LI ix=&(x->root);

  ins=CreateListNode(y);
  ins->next=(ListNode *)(*ix);
  *ix=ins;
  return;
}

void InsertInList	(	jList *	,
		Poly *
	)

Definition at line 559 of file janet.cc.

{
  ListNode *ins;
  LI ix=&(x->root);

  while (*ix)
  {
    if (pLmCmp(y->lead,(*ix)->info->lead) == -1)
      ix=(ListNode **)&((*ix)->next);
    else
      break;
  }

  ins=CreateListNode(y);
  ins->next=(ListNode *)(*ix);
  *ix=ins;
  return;
}

Poly* is_div_	(	TreeM *	root,
		poly	item
	)

Definition at line 860 of file janet.cc.

{
  int power_tmp,i,i_con=currRing->N-1;
  NodeM *curr=tree->root;

  if (!curr) return NULL;
  if (pow_(item) == 0) return NULL;

  for ( ; i_con>=0 && !pGetExp(item,i_con+1) ; i_con--)
    ;

  for (i=0; i <= i_con ; i++)
  {
    power_tmp=pGetExp(item,i+1);

    while (power_tmp)
    {
      if (curr->ended) return curr->ended;

      if (!curr->left)
      {
        if (curr->right)
          return is_div_upper(item,curr->right,i); //??????
        return NULL;
      }

      curr=curr->left;
      power_tmp--;
    }

    if (curr->ended) return curr->ended;

    if (!curr->right) return NULL;

    curr=curr->right;
  }

  if (curr->ended) return curr->ended;
  else return NULL;
}

Poly* is_present	(	jList *	,
		poly
	)

Definition at line 736 of file janet.cc.

{
  LCI iF=F->root;
  while(iF)
    if (pLmCmp(iF->info->root,x) == 0)
      return iF->info;
    else iF=iF->next;

  return NULL;
}

int ListGreatMoveDegree	(	jList *	,
		jList *	,
		poly
	)

Definition at line 606 of file janet.cc.

{
  LCI y=A->root;
  int pow_x=pow_(x);

  if (!y || pow_(y->info->lead) <= pow_x) return 0;

  while(y && pow_(y->info->lead) > pow_x)
  {
    InsertInCount(B,y->info);
    A->root=y->next;
    GCF(y);
    y=A->root;
  }

  return 1;
}

int ListGreatMoveOrder	(	jList *	,
		jList *	,
		poly
	)

Definition at line 589 of file janet.cc.

{
  LCI y=A->root;

  if (!y || pLmCmp(y->info->lead,x) < 0) return 0;

  while(y && pLmCmp(y->info->lead,x) >= 0)
  {
    InsertInCount(B,y->info);
    A->root=y->next;
    GCF(y);
    y=A->root;
  }

  return 1;
}

Poly* NewPoly

(

poly

p = NULL

)

Definition at line 368 of file janet.cc.

{
  Poly *beg=(Poly *)GCM(sizeof(Poly));

  beg->root=p;//(p == NULL ? pInit() : p);
  beg->root_b=NULL;
  beg->root_l=0;
  beg->history=NULL;//pInit();
  beg->lead=NULL;
  beg->mult=(char *)GCMA(sizeof(char)*2*offset);

  for (int i=0; i < currRing->N; i++)
  {
    ClearMult(beg,i);
    ClearProl(beg,i);
  };

  beg->prolonged=-1;

  return beg;
}

void NFL	(	Poly *	,
		TreeM *
	)

Definition at line 195 of file janet.cc.

{
  Poly *f;
  // int g1,f1,gg;

  if ((f=is_div_(F,p->lead))==NULL) return;

  int pX=pow_(p->lead);
  int phX=pow_(p->history);

  if (pX!=phX)
  {
    int phF=pow_(f->history);
    if (pX >= (phX+phF))
    {
      pDelete(&p->root);
      //p->root=NULL;
      return;
    }

/*    poly p2=pInit();
    pLcm(p->history,f->history,p2);
    pSetm(p2);

    if (pLmCmp(p->root,p2) > 0)
    {
      pLmDelete(&p2);
      pDelete(&p->root);
      //p->root=NULL;
      return;
    }

    pLmDelete(&p2);
*/
/*    for(int i=0, gg=0 ; i<currRing->N;i++)
      if ((g1=pGetExp(p->history,i+1)) > (f1=pGetExp(f->history,i+1)))
        gg+=g1;
      else gg+=f1;

    if (pX > gg)
      {
        pDelete(&p->root);
        //x->root=NULL;
        return;
    }
*/
    int pF=pow_(f->lead);

    if ((pX == pF) && (pF == phF))
    {
      pLmFree(&f->history);
      f->history=p_Copy_noCheck(p->history,currRing); /* cf of p->history is NULL */
    }
  }

  //if (TEST_OPT_PROT) { PrintS("R"); mflush(); }
  int /*old_size, */count;
  count=0;
  while(f && p->root)
  {
//    PrintS("R");
//    if (TEST_OPT_PROT) { PrintS("R"); mflush(); }
#if 0
    old_size=nSize(pGetCoeff(p->root));
#endif
    if (ReducePolyLead(p,f) == 0) break;
    if (p->root!=NULL)
    {
      count++;
#if 0
      if ((count>4) && (3<nSize(pGetCoeff(p->root)))
      && (nSize(pGetCoeff(p->root))>old_size))
      {
        pSimpleContent(p->root,old_size);
        count=0;
      }
#else
      if (count>50)
      {
        kBucketClear(p->root_b,&p->root,&p->root_l);
        p_SimpleContent(p->root,2,currRing);
        kBucketInit(p->root_b,p->root,p->root_l);
        count=0;
        //PrintS(".");
      }
#endif
      f=is_div_(F,p->root);
    }
  }
#if 1
  if (p->root_b!=NULL)
  {
    kBucketClear(p->root_b,&p->root,&p->root_l);
    kBucketDestroy(&p->root_b);
    p->root_b=NULL;
  }
#endif

  if (!p->root)
    return;

  InitHistory(p);
  InitProl(p);
  InitLead(p);
  p->changed=1;

  p_Content(p->root,currRing);
  //pSimpleContent(p->root,-1);
  pTest(p->root);
}

void NFListQ

(

)

Definition at line 638 of file janet.cc.

{
  LCI ll;
  int p,p1;
  LI l;

  do
  {
    if (!Q->root) break;

    ll=Q->root;

    p=pow_(Q->root->info->lead);

    while (ll)
    {
      int ploc=pow_(ll->info->lead);
      if (ploc < p) p=ploc;
      ll=ll->next;
    }

    p1=1;

    l=&(Q->root);

    while (*l)
    {
//      PrintS("*");
      int ploc=pow_((*l)->info->lead);

      if (ploc == p)
      {
        if (!ValidatePoly((*l)->info,G))
        {
          ll=(*l);
          *l=(*l)->next;
          DestroyListNode(ll);
          continue;
        };

        (*l)->info->changed=0;
//        PrintS("!");
        NFL((*l)->info,G);
//                                PrintS("$");
        if (!(*l)->info->root)
        {
          ll=(*l);
          *l=(*l)->next;
          DestroyListNode(ll);
          continue;
        };
        p1=0;
      }

      l=&((*l)->next);
    }
  }while(p1);
//  PrintLn();
}

void PNF	(	Poly *	,
		TreeM *
	)

Definition at line 151 of file janet.cc.

{
  if (p->root==NULL) return;

  Poly *f;
  BOOLEAN done=FALSE;
  poly temp=p->root;

//  if (TEST_OPT_PROT) { PrintS("r"); mflush(); }
  int count=0;
  poly pp=p->root;
  int old_size=nSize(pGetCoeff(pp));
  p->root_l=0;
  while(temp->next)
  {
    f=is_div_(F,temp->next);
    if (f)
    {
      if (ReducePoly(p,temp,f)) //temp->next
      {
        count++;
        //if (TEST_OPT_PROT) { PrintS("-"); mflush(); }
        if ((f!=NULL)
        && (count>20)
        && (nSize(pGetCoeff(pp))>old_size)
        )
        {
           //pSimpleContent(pp,2);
           p_Content(pp,currRing);
           count=0;
         //  old_size=nSize(pGetCoeff(pp));
        }
      }
      done=TRUE;
    }
    else
      temp=temp->next;
   }

  if (done) p_Content(p->root,currRing);
  //if (done) pSimpleContent(p->root,-1);
  pTest(p->root);
}

int ProlCompare	(	Poly *	,
		Poly *
	)

Definition at line 470 of file janet.cc.

{
  switch(pLmCmp(item1->lead,item2->lead))
  {
    case -1:
      return 1;

    case 1:
      return 0;

    default:
      if ((item1->root_l<=0)||(item2->root_l<=0))
        return LengthCompare(item1->root,item2->root);
      return item1->root_l<=item2->root_l;
  }
}

void ProlVar	(	Poly *	,
		int
	)

Definition at line 487 of file janet.cc.

{
  Poly *Pr;

  if (!GetProl(temp,i) && !GetMult(temp,i))
  {
    Pr=NewPoly();
    SetProl(temp,i);

    Pr->prolonged=i;
    Pr->history=pLmInit(temp->history);
    Pr->lead=pLmInit(temp->lead);
    pIncrExp(Pr->lead,i+1);
    pSetm(Pr->lead);
     InitProl(temp);

     Pr->changed=0;
//    pTest(Pr->root);
      InsertInCount(Q,Pr);
   }
}

int ReducePoly	(	Poly *	x,
		Poly *	y
	)

int ReducePolyLead	(	Poly *	x,
		Poly *	y
	)

Definition at line 81 of file janet.cc.

{
  if (!x->root || !y->root)
    return 0;

/*  poly b1=pDivide(x->root,y->root);

  number gcd=n_Gcd(pGetCoeff(x->root),pGetCoeff(y->root),currRing->cf);

  number a1=nDiv(pGetCoeff(y->root),gcd);
  pGetCoeff(b1)=nDiv(pGetCoeff(x->root),gcd);

  x->root=pMult_nn(x->root,a1);
  nDelete(&a1);

  x->root=pMinus_mm_Mult_qq(x->root,b1,y->root);

  pDelete(&b1);
*/
#if 1
  if (x->root_b==NULL)
  {
    if (x->root_l<=0) x->root_l=pLength(x->root);
    x->root_b=kBucketCreate(currRing);
    kBucketInit(x->root_b,x->root,x->root_l);
  }
  number coef;
  if (y->root_l<=0) y->root_l=pLength(y->root);
  coef=kBucketPolyRed(x->root_b,y->root,y->root_l,NULL);
  nDelete(&coef);
  x->root=kBucketGetLm(x->root_b);
  if (x->root==NULL)
  {
    kBucketDestroy(&x->root_b);
    x->root_b=NULL;
    x->root_l=0;
  }
#else
  x->root=ksOldSpolyRed(y->root,x->root,NULL);
#endif
//  if (x->root) p_Content(x->root,currRing);
//  if (x->root) pSimpleContent(x->root,5);

  return 1;
}

void SetMult	(	Poly *	,
		int
	)

Definition at line 433 of file janet.cc.

{
  x->mult[i/8] |= Mask[i%8];
}

void SetProl	(	Poly *	,
		int
	)

Definition at line 448 of file janet.cc.

{
  (x->mult+offset)[i/8] |= Mask[i%8];
}

int ValidatePoly	(	Poly *	,
		TreeM *
	)

Definition at line 306 of file janet.cc.

{
  Poly /*f,*/*g;
  // int g1,f1;

  if (x->root) return 1;

  g=is_present(T,x->history); //it's a prolongation - do we have a parent ?

  if (!g)  return 0; //if not - kill him !

  poly lmX=pDivide(x->lead,g->root);
  pGetCoeff(lmX)=nInit(1);

/*  if ((f=is_div_(F,lmX)) != NULL)
  {
    int pX=pow_(lmX);
    int phX=pow_(x->history);

    if (pX!=phX)
    {
      int phF=pow_(f->history);
      if (pX >= (phX+phF))
      {
        pLmDelete(&lmX);
        //x->root=NULL;
        return 0;
      }

      for(int i=0, gg=0 ; i<currRing->N;i++)
        if ((g1=pGetExp(x->history,i+1)) > (f1=pGetExp(f->history,i+1)))
          gg+=g1;
        else
          gg+=f1;

      if (pX > gg)
      {
        pLmDelete(&lmX);
        return 0;
      }
      int pF=pow_(f->root);

      if ((pX == pF) && (pF == phF))
        f->history=x->history;
    }
  }

  pLmDelete(&lmX);

*/
  x->root=pCopy(g->root);
  x->root_l=g->root_l;

  x->root=pMult(x->root,lmX);

  pTest(x->root);

  x->prolonged=-1;

  return 1;
}

Variable Documentation

Poly* ended

Definition at line 42 of file janet.h.

ListNode Poly* info

Definition at line 30 of file janet.h.

ListNode* next

Definition at line 31 of file janet.h.

NodeM

Initial value:

{
 NodeM *left,*right

Definition at line 40 of file janet.h.

TreeM

Initial value:

{
 NodeM *root

Definition at line 46 of file janet.h.

---

Generated on Mon Mar 6 2017 11:17:34 by   doxygen 1.8.8 for Singular debian-1:4.0.3-p3+ds-5~bpo8+1