khstd.h File Reference

#include "kernel/structs.h"
#include "coeffs/bigintmat.h"

Go to the source code of this file.

Functions
void	khCheck (ideal Q, intvec *w, bigintmat *hilb, int &eledeg, int &count, kStrategy strat)
void	khCheckLocInhom (ideal Q, intvec *w, bigintmat *hilb, int &count, kStrategy strat)

Function Documentation

khCheck()

void khCheck	(	ideal	Q,
		intvec *	w,
		bigintmat *	hilb,
		int &	eledeg,
		int &	count,
		kStrategy	strat )

Definition at line 28 of file khstd.cc.

{
  bigintmat *newhilb;
  int deg,l,ln,mw;
  pFDegProc degp;
 
  eledeg--;
  if (eledeg == 0)
  {
    if (strat->ak>0)
    {
      char *used_comp=(char*)omAlloc0(strat->ak+1);
      int i;
      for(i=strat->sl;i>0;i--)
      {
        used_comp[pGetComp(strat->S[i])]='\1';
      }
      for(i=strat->ak;i>0;i--)
      {
        if(used_comp[i]=='\0')
        {
          omFree((ADDRESS)used_comp);
          return;
        }
      }
      omFree((ADDRESS)used_comp);
    }
    degp=currRing->pFDeg;
    // if weights for variables were given to std computations,
    // then pFDeg == degp == kHomModDeg (see kStd)
    if ((degp!=kModDeg) && (degp!=kHomModDeg)) degp=p_Totaldegree;
    // degp = pWDegree;
    l = hilb->cols();
    mw = n_Int(BIMATELEM(*hilb,1,l),coeffs_BIGINT);
    newhilb =hFirstSeries0b(strat->Shdl,Q,strat->kHomW,w,currRing,coeffs_BIGINT);
    ln = newhilb->cols();
    deg = degp(strat->P.p,currRing);
    loop // compare the series in degree deg, try to increase deg -----------
    {
      if (deg < ln) // deg may be out of range
      {
        if (deg < l)
        {
          number s=n_Sub(BIMATELEM(*newhilb,1,deg+1),
                         BIMATELEM(*hilb,1,deg+1),coeffs_BIGINT);
          eledeg = n_Int(s,coeffs_BIGINT);
        }
        else
          eledeg = n_Int(BIMATELEM(*newhilb,1,deg+1),coeffs_BIGINT);
      }
      else
      {
        if (deg < l)
          eledeg = -n_Int(BIMATELEM(*hilb,1,deg+1),coeffs_BIGINT);
        else // we have newhilb = hilb
        {
          while (strat->Ll>=0)
          {
            count++;
            if(TEST_OPT_PROT)
            {
              PrintS("h");
              mflush();
            }
            deleteInL(strat->L,&strat->Ll,strat->Ll,strat);
          }
          delete newhilb;
          return;
        }
      }
      if (eledeg > 0) // elements to delete
        break;
      else if (eledeg <0) // strange....see bug_43
        return;
      deg++;
    } /* loop */
    delete newhilb;
    while ((strat->Ll>=0) && (degp(strat->L[strat->Ll].p,currRing)-mw < deg)) // the essential step
    {
      count++;
      if(TEST_OPT_PROT)
      {
        PrintS("h");
        mflush();
      }
      deleteInL(strat->L,&strat->Ll,strat->Ll,strat);
    }
  }
}

khCheckLocInhom()

void khCheckLocInhom	(	ideal	Q,
		intvec *	w,
		bigintmat *	hilb,
		int &	count,
		kStrategy	strat )

Definition at line 248 of file khstd.cc.

{
  ideal Lm;
  bigintmat *newhilb;
 
  Lm = id_Head(strat->Shdl,currRing);
 
  newhilb =hFirstSeries0b(strat->Shdl,Q,strat->kHomW,w,currRing,coeffs_BIGINT);
 
  if(newhilb->compare(hilb) == 0)
  {
    while (strat->Ll>=0)
    {
      count++;
      if(TEST_OPT_PROT)
      {
        PrintS("h");
        mflush();
      }
      deleteInL(strat->L,&strat->Ll,strat->Ll,strat);
    }
    delete newhilb;
    return;
  }
  id_Delete(&Lm,currRing);
}