groebnerCone.cc File Reference

#include <utility>
#include <kernel/GBEngine/kstd1.h>
#include <kernel/ideals.h>
#include <Singular/ipid.h>
#include <polys/monomials/p_polys.h>
#include <polys/monomials/ring.h>
#include <polys/prCopy.h>
#include <gfanlib/gfanlib.h>
#include <gfanlib/gfanlib_matrix.h>
#include <initial.h>
#include <tropicalStrategy.h>
#include <groebnerCone.h>
#include <callgfanlib_conversion.h>
#include <containsMonomial.h>
#include <tropicalCurves.h>
#include <bbcone.h>

Go to the source code of this file.

Functions
static bool	checkPolynomialInput (const ideal I, const ring r)
static bool	checkOrderingAndCone (const ring r, const gfan::ZCone zc)
static bool	checkPolyhedralInput (const gfan::ZCone zc, const gfan::ZVector p)
gfan::ZFan *	toFanStar (groebnerCones setOfCones)
BOOLEAN	flipConeDebug (leftv res, leftv args)
BOOLEAN	groebnerNeighboursDebug (leftv res, leftv args)
BOOLEAN	tropicalNeighboursDebug (leftv res, leftv args)

Function Documentation

static bool checkOrderingAndCone ( const ring  r, const gfan::ZCone  zc  )

static

Definition at line 32 of file groebnerCone.cc.

{
  return true;
  if (r)
  {
    int n = rVar(r); int* w = r->wvhdl[0];
    gfan::ZVector v = wvhdlEntryToZVector(n,w);
    if (r->order[0]==ringorder_ws)
      v = gfan::Integer((long)-1)*v;
    if (!zc.contains(v))
    {
      std::cout << "ERROR: weight of ordering not inside Groebner cone!" << std::endl
                << "cone: " << std::endl
                << toString(&zc)
                << "weight: " << std::endl
                << v << std::endl;
      return false;
    }
    return true;
  }
  return (zc.dimension()==0);
}

static bool checkPolyhedralInput ( const gfan::ZCone  zc, const gfan::ZVector  p  )

static

Definition at line 55 of file groebnerCone.cc.

{
  return zc.containsRelatively(p);
}

static bool checkPolynomialInput ( const ideal  I, const ring  r  )

static

Definition at line 25 of file groebnerCone.cc.

{
  if (r) rTest(r);
  if (I && r) id_Test(I,r);
  return ((!I) || (I && r));
}

BOOLEAN flipConeDebug	(	leftv	res,
		leftv	args
	)

Definition at line 526 of file groebnerCone.cc.

{
  leftv u = args;
  if ((u!=NULL) && (u->Typ()==IDEAL_CMD))
  {
    leftv v = u->next;
    if ((v!=NULL) && (v->Typ()==NUMBER_CMD))
    {
      leftv w = v->next;
      if ((w!=NULL) && (w->Typ()==BIGINTMAT_CMD))
      {
        leftv x = w->next;
        if ((x!=NULL) && (x->Typ()==BIGINTMAT_CMD))
        {
          omUpdateInfo();
          Print("usedBytesBefore=%ld\n",om_Info.UsedBytes);

          ideal I = (ideal) u->CopyD();
          number p = (number) v->CopyD();
          bigintmat* interiorPoint0 = (bigintmat*) w->CopyD();
          bigintmat* facetNormal0 = (bigintmat*) x->CopyD();
          tropicalStrategy debug = tropicalStrategy::debugStrategy(I,p,currRing);

          gfan::ZVector* interiorPoint = bigintmatToZVector(interiorPoint0);
          gfan::ZVector* facetNormal = bigintmatToZVector(facetNormal0);

          groebnerCone sigma(I,currRing,debug);
          groebnerCone theta = sigma.flipCone(*interiorPoint,*facetNormal);

          id_Delete(&I,currRing);
          n_Delete(&p,currRing->cf);
          delete interiorPoint0;
          delete facetNormal0;
          delete interiorPoint;
          delete facetNormal;

          res->rtyp = NONE;
          res->data = NULL;
          return FALSE;
        }
      }
    }
  }
  WerrorS("computeFlipDebug: unexpected parameters");
  return TRUE;
}

BOOLEAN groebnerNeighboursDebug	(	leftv	res,
		leftv	args
	)

Definition at line 573 of file groebnerCone.cc.

{
  leftv u = args;
  if ((u!=NULL) && (u->Typ()==IDEAL_CMD))
  {
    leftv v = u->next;
    if ((v!=NULL) && (v->Typ()==NUMBER_CMD))
    {
      omUpdateInfo();
      Print("usedBytesBefore=%ld\n",om_Info.UsedBytes);

      ideal I = (ideal) u->CopyD();
      number p = (number) v->CopyD();

      tropicalStrategy debug = tropicalStrategy::debugStrategy(I,p,currRing);
      groebnerCone sigma(I,currRing,debug);
      groebnerCones neighbours = sigma.groebnerNeighbours();

      id_Delete(&I,currRing);
      n_Delete(&p,currRing->cf);
      res->rtyp = NONE;
      res->data = NULL;
      return FALSE;
    }
  }
  WerrorS("computeFlipDebug: unexpected parameters");
  return TRUE;
}

gfan::ZFan* toFanStar

(

groebnerCones

setOfCones

)

Definition at line 505 of file groebnerCone.cc.

{
  if (setOfCones.size() > 0)
  {
    groebnerCones::iterator sigma = setOfCones.begin();
    gfan::ZFan* zf = new gfan::ZFan(sigma->getPolyhedralCone().ambientDimension());
    for (; sigma!=setOfCones.end(); sigma++)
    {
      gfan::ZCone zc = sigma->getPolyhedralCone();
      // assume(isCompatible(zf,&zc));
      zf->insert(zc);
    }
    return zf;
  }
  else
    return new gfan::ZFan(gfan::ZFan::fullFan(currRing->N));
}

BOOLEAN tropicalNeighboursDebug	(	leftv	res,
		leftv	args
	)

Definition at line 602 of file groebnerCone.cc.

{
  leftv u = args;
  if ((u!=NULL) && (u->Typ()==IDEAL_CMD))
  {
    leftv v = u->next;
    if ((v!=NULL) && (v->Typ()==NUMBER_CMD))
    {
      omUpdateInfo();
      Print("usedBytesBefore=%ld\n",om_Info.UsedBytes);

      ideal I = (ideal) u->CopyD();
      number p = (number) v->CopyD();

      tropicalStrategy debug = tropicalStrategy::debugStrategy(I,p,currRing);
      groebnerCone sigma(I,currRing,debug);
      groebnerCones neighbours = sigma.groebnerNeighbours();

      id_Delete(&I,currRing);
      n_Delete(&p,currRing->cf);
      res->rtyp = NONE;
      res->data = NULL;
      return FALSE;
    }
  }
  WerrorS("computeFlipDebug: unexpected parameters");
  return TRUE;
}