tropicalDebug.h File Reference

#include "libpolys/polys/simpleideals.h"
#include "gfanlib/gfanlib.h"

Go to the source code of this file.

Functions
bool	checkForNonPositiveEntries (const gfan::ZVector &w)
bool	checkForNonPositiveLaterEntries (const gfan::ZVector &w)
bool	areIdealsEqual (ideal I, ring r, ideal J, ring s)
bool	checkWeightVector (const ideal I, const ring r, const gfan::ZVector &weightVector, bool checkBorder=0)
bool	checkOrderingAndCone (const ring r, const gfan::ZCone zc)
bool	checkPolynomialInput (const ideal I, const ring r)
bool	checkPolyhedralInput (const gfan::ZCone zc, const gfan::ZVector p)

Function Documentation

areIdealsEqual()

bool areIdealsEqual	(	ideal	I,
		ring	r,
		ideal	J,
		ring	s
	)

Definition at line 41 of file tropicalDebug.cc.

{
  ideal Is = idInit(IDELEMS(I));
  nMapFunc identity = n_SetMap(r->cf,s->cf);
  for (int i=0; i<IDELEMS(I); i++)
    Is->m[i] = p_PermPoly(I->m[i],NULL,r,s,identity,NULL,0);
 
  ring origin = currRing;
  if (origin!=s)
    rChangeCurrRing(s);
  ideal stdI = gfanlib_kStd_wrapper(Is,s);
  ideal stdJ = gfanlib_kStd_wrapper(J,s);
  ideal NF1 = kNF(stdI,s->qideal,stdJ);
  ideal NF2 = kNF(stdJ,s->qideal,stdI);
  if (origin!=s)
    rChangeCurrRing(origin);
 
  bool b = true;
  if (NF1 != NULL)
    b = b && idIs0(NF1);
  if (NF2 != NULL)
    b = b && idIs0(NF2);
  if (!b)
    std::cout << "ERROR: input ideals not equal!" << std::endl;
 
  id_Delete(&stdI,s);
  id_Delete(&stdJ,s);
  id_Delete(&NF1,s);
  id_Delete(&NF2,s);
  return (b);
}

checkForNonPositiveEntries()

bool checkForNonPositiveEntries

(

const gfan::ZVector &

w

)

Definition at line 13 of file tropicalDebug.cc.

{
  for (unsigned i=0; i<w.size(); i++)
  {
    if (w[i].sign()<=0)
    {
      std::cout << "ERROR: non-positive weight in weight vector" << std::endl
                << "weight: " << w << std::endl;
      return false;
    }
  }
  return true;
}

checkForNonPositiveLaterEntries()

bool checkForNonPositiveLaterEntries

(

const gfan::ZVector &

w

)

Definition at line 27 of file tropicalDebug.cc.

{
  for (unsigned i=1; i<w.size(); i++)
  {
    if (w[i].sign()<=0)
    {
      std::cout << "ERROR: non-positive weight in weight vector later entries" << std::endl
                << "weight: " << w << std::endl;
      return false;
    }
  }
  return true;
}

checkOrderingAndCone()

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

Definition at line 90 of file tropicalDebug.cc.

{
  if (r)
  {
    if (r->order[0]==ringorder_dp)
      return true;
    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;
      return false;
    }
  }
  return true;
}

checkPolyhedralInput()

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

Definition at line 117 of file tropicalDebug.cc.

{
  return zc.containsRelatively(p);
}

checkPolynomialInput()

bool checkPolynomialInput	(	const ideal	I,
		const ring	r
	)

Definition at line 110 of file tropicalDebug.cc.

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

checkWeightVector()

bool checkWeightVector	(	const ideal	I,
		const ring	r,
		const gfan::ZVector &	weightVector,
		bool	checkBorder = 0
	)

Definition at line 74 of file tropicalDebug.cc.

{
  gfan::ZCone cI = maximalGroebnerCone(I,r);
  if (!cI.contains(weightVector))
  {
    std::cout << "ERROR: weight vector not inside maximal Gröbner cone" << std::endl;
    return false;
  }
  if (checkBorder && cI.containsRelatively(weightVector))
  {
    std::cout << "ERROR: weight vector in the relative interior of maximal Gröbner cone" << std::endl;
    return false;
  }
  return true;
}