tropicalStrategy.h File Reference

implementation of the class tropicalStrategy More...

#include <gfanlib/gfanlib_vector.h>
#include <gfanlib/gfanlib_zcone.h>
#include <polys/simpleideals.h>
#include <kernel/ideals.h>
#include <set>
#include <callgfanlib_conversion.h>
#include <containsMonomial.h>
#include <flip.h>
#include <initial.h>
#include <witness.h>
#include <Singular/ipshell.h>
#include <adjustWeights.h>
#include <ppinitialReduction.h>

Go to the source code of this file.

Data Structures
class	tropicalStrategy

Typedefs
typedef gfan::ZVector(*	wAdjAlg1) (gfan::ZVector)
typedef gfan::ZVector(*	wAdjAlg2) (gfan::ZVector, gfan::ZVector)
typedef bool(*	redAlg) (ideal, ring, number)

Functions
BOOLEAN	computeWitnessDebug (leftv res, leftv args)
BOOLEAN	computeFlipDebug (leftv res, leftv args)

Detailed Description

implementation of the class tropicalStrategy

tropicalStrategy is a class that contains information relevant for computing tropical varieties that is dependent on the valuation of the coefficient field. It represents the mutable part of an overall framework that is capable of computing tropical varieties both over coefficient fields without valuation and with valuation (currently: only p-adic valuation over rational numbers)

Definition in file tropicalStrategy.h.

Typedef Documentation

typedef bool(* redAlg) (ideal, ring, number)

Definition at line 35 of file tropicalStrategy.h.

typedef gfan::ZVector(* wAdjAlg1) (gfan::ZVector)

Definition at line 33 of file tropicalStrategy.h.

typedef gfan::ZVector(* wAdjAlg2) (gfan::ZVector, gfan::ZVector)

Definition at line 34 of file tropicalStrategy.h.

Function Documentation

BOOLEAN computeFlipDebug	(	leftv	res,
		leftv	args
	)

Definition at line 946 of file tropicalStrategy.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);
          std::pair<ideal,ring> Js = debug.computeFlip(I,currRing,*interiorPoint,*facetNormal);
          ideal J = Js.first;
          ring s = Js.second;

          id_Delete(&J,s);
          rDelete(s);

          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 computeWitnessDebug	(	leftv	res,
		leftv	args
	)

Definition at line 909 of file tropicalStrategy.cc.

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

          ideal inJ = (ideal) u->CopyD();
          ideal inI = (ideal) v->CopyD();
          ideal I = (ideal) w->CopyD();
          number p = (number) x->CopyD();
          tropicalStrategy debug = tropicalStrategy::debugStrategy(I,p,currRing);
          ideal J = debug.computeWitness(inJ,inI,I,currRing);
          id_Delete(&inJ,currRing);
          id_Delete(&inI,currRing);
          id_Delete(&I,currRing);
          n_Delete(&p,currRing->cf);
          res->rtyp = IDEAL_CMD;
          res->data = (char*) J;
          return FALSE;
        }
      }
    }
  }
  return TRUE;
}