Functions | Variables
tropicalVariety.cc File Reference
#include <callgfanlib_conversion.h>
#include <std_wrapper.h>
#include <bbfan.h>
#include <groebnerCone.h>
#include <tropicalVarietyOfPolynomials.h>
#include <tropicalVarietyOfIdeals.h>
#include <coeffs/numbers.h>
#include <misc/options.h>
#include <kernel/structs.h>
#include <iostream>
#include <gfanlib/gfanlib_zfan.h>

Go to the source code of this file.

Functions

static void setOptionRedSB ()
 
static void undoSetOptionRedSB ()
 
static gfan::ZFan * toZFan (std::set< gfan::ZCone > maxCones)
 
BOOLEAN tropicalVariety (leftv res, leftv args)
 

Variables

BITSET bitsetSave1
 
BITSET bitsetSave2
 

Function Documentation

static void setOptionRedSB ( )
static

Definition at line 20 of file tropicalVariety.cc.

21 {
22  SI_SAVE_OPT(bitsetSave1,bitsetSave2);
23  si_opt_1|=Sy_bit(OPT_REDSB);
24 }
OPT_REDSB
#define OPT_REDSB
Definition: options.h:71
si_opt_1
unsigned si_opt_1
Definition: options.c:5
Sy_bit
#define Sy_bit(x)
Definition: options.h:30
bitsetSave2
BITSET bitsetSave2
Definition: tropicalVariety.cc:14
bitsetSave1
BITSET bitsetSave1
Definition: tropicalVariety.cc:14
SI_SAVE_OPT
#define SI_SAVE_OPT(A, B)
Definition: options.h:19
static gfan::ZFan* toZFan ( std::set< gfan::ZCone >  maxCones)
static

Definition at line 35 of file tropicalVariety.cc.

36 {
37  std::set<gfan::ZCone>::iterator sigma = maxCones.begin();
38  gfan::ZFan* zf = new gfan::ZFan(sigma->ambientDimension());
39  for (; sigma!=maxCones.end(); sigma++)
40  zf->insert(*sigma);
41  return zf;
42 }
BOOLEAN tropicalVariety ( leftv  res,
leftv  args 
)

Definition at line 45 of file tropicalVariety.cc.

46 {
47  leftv u = args;
48  if ((u!=NULL) && (u->Typ()==POLY_CMD))
49  {
50  poly g = (poly) u->Data();
51  leftv v = u->next;
52  if (v==NULL)
53  {
54  try
55  {
56  ideal I = idInit(1);
57  I->m[0] = g;
58  tropicalStrategy currentStrategy(I,currRing);
59  std::set<gfan::ZCone> maxCones = tropicalVariety(g,currRing,&currentStrategy);
60  res->rtyp = fanID;
61  res->data = (char*) toZFan(maxCones);
62  I->m[0] = NULL;
63  id_Delete(&I,currRing);
64  return FALSE;
65  }
66  catch (const std::exception& ex)
67  {
68  Werror("ERROR: %s",ex.what());
69  return TRUE;
70  }
71  }
72  if ((v!=NULL) && (v->Typ()==NUMBER_CMD))
73  {
74  try
75  {
76  ideal I = idInit(1);
77  I->m[0] = g;
78  number p = (number) v->Data();
79  tropicalStrategy currentStrategy(I,p,currRing);
80  ideal startingIdeal = currentStrategy.getStartingIdeal();
81  ring startingRing = currentStrategy.getStartingRing();
82  poly gStart = startingIdeal->m[0];
83  std::set<gfan::ZCone> maxCones = tropicalVariety(gStart,startingRing,&currentStrategy);
84  res->rtyp = fanID;
85  res->data = (char*) toZFan(maxCones);
86  I->m[0] = NULL;
87  id_Delete(&I,currRing);
88  return FALSE;
89  }
90  catch (const std::exception& ex)
91  {
92  Werror("ERROR: %s",ex.what());
93  return TRUE;
94  }
95 
96  }
97  }
98  if ((u!=NULL) && (u->Typ()==IDEAL_CMD))
99  {
100  ideal I = (ideal) u->Data();
101  leftv v = u->next;
102 
103  if (idSize(I)==1)
104  {
105  poly g = I->m[0];
106  if (v==NULL)
107  {
108  try
109  {
110  tropicalStrategy currentStrategy(I,currRing);
111  std::set<gfan::ZCone> maxCones = tropicalVariety(g,currRing,&currentStrategy);
112  res->rtyp = fanID;
113  res->data = (char*) toZFan(maxCones);
114  return FALSE;
115  }
116  catch (const std::exception& ex)
117  {
118  Werror("ERROR: %s",ex.what());
119  return TRUE;
120  }
121  }
122  if ((v!=NULL) && (v->Typ()==NUMBER_CMD))
123  {
124  try
125  {
126  number p = (number) v->Data();
127  tropicalStrategy currentStrategy(I,p,currRing);
128  ideal startingIdeal = currentStrategy.getStartingIdeal();
129  ring startingRing = currentStrategy.getStartingRing();
130  poly gStart = startingIdeal->m[0];
131  std::set<gfan::ZCone> maxCones = tropicalVariety(gStart,startingRing,&currentStrategy);
132  res->rtyp = fanID;
133  res->data = (char*) toZFan(maxCones);
134  return FALSE;
135  }
136  catch (const std::exception& ex)
137  {
138  Werror("ERROR: %s",ex.what());
139  return TRUE;
140  }
141  }
142  }
143 
144  if (v==NULL)
145  {
146  try
147  {
148  setOptionRedSB();
149  ideal stdI;
150  if (!hasFlag(u,FLAG_STD))
151  stdI = gfanlib_kStd_wrapper(I,currRing);
152  else
153  stdI = id_Copy(I,currRing);
154  tropicalStrategy currentStrategy(stdI,currRing);
155  gfan::ZFan* tropI = tropicalVariety(currentStrategy);
156  res->rtyp = fanID;
157  res->data = (char*) tropI;
158  undoSetOptionRedSB();
159  id_Delete(&stdI,currRing);
160  return FALSE;
161  }
162  catch (const std::exception& ex)
163  {
164  Werror("ERROR: %s",ex.what());
165  return TRUE;
166  }
167  }
168  if ((v!=NULL) && (v->Typ()==NUMBER_CMD))
169  {
170  try
171  {
172  number p = (number) v->Data();
173  ideal stdI;
174  if (!hasFlag(u,FLAG_STD))
175  stdI = gfanlib_kStd_wrapper(I,currRing);
176  else
177  stdI = id_Copy(I,currRing);
178  tropicalStrategy currentStrategy(stdI,p,currRing);
179  gfan::ZFan* tropI = tropicalVariety(currentStrategy);
180  res->rtyp = fanID;
181  res->data = (char*) tropI;
182  id_Delete(&stdI,currRing);
183  return FALSE;
184  }
185  catch (const std::exception& ex)
186  {
187  Werror("ERROR: %s",ex.what());
188  return TRUE;
189  }
190  }
191  return FALSE;
192  }
193  WerrorS("tropicalVariety: unexpected parameters");
194  return TRUE;
195 }
sleftv
Class used for (list of) interpreter objects.
Definition: subexpr.h:83
FALSE
#define FALSE
Definition: auxiliary.h:140
setOptionRedSB
static void setOptionRedSB()
Definition: tropicalVariety.cc:20
p
return P p
Definition: myNF.cc:203
id_Copy
ideal id_Copy(ideal h1, const ring r)
Definition: simpleideals.cc:398
ideal
const ideal
Definition: gb_hack.h:42
id_Delete
void id_Delete(ideal *h, ring r)
Definition: simpleideals.cc:119
TRUE
#define TRUE
Definition: auxiliary.h:144
toZFan
static gfan::ZFan * toZFan(std::set< gfan::ZCone > maxCones)
Definition: tropicalVariety.cc:35
g
g
Definition: cfModGcd.cc:4031
WerrorS
void WerrorS(const char *s)
Definition: feFopen.cc:23
fanID
int fanID
Definition: bbfan.cc:18
sleftv::Typ
int Typ()
Definition: subexpr.cc:949
sleftv::data
void * data
Definition: subexpr.h:89
currRing
ring currRing
Widely used global variable which specifies the current polynomial ring for Singular interpreter and ...
Definition: polys.cc:12
gfanlib_kStd_wrapper
ideal gfanlib_kStd_wrapper(ideal I, ring r, tHomog h=testHomog)
Definition: std_wrapper.cc:5
poly
polyrec * poly
Definition: hilb.h:10
tropicalVariety
BOOLEAN tropicalVariety(leftv res, leftv args)
Definition: tropicalVariety.cc:45
FLAG_STD
#define FLAG_STD
Definition: ipid.h:108
sleftv::next
leftv next
Definition: subexpr.h:87
idInit
ideal idInit(int idsize, int rank)
Definition: simpleideals.cc:40
v
const Variable & v
< [in] a sqrfree bivariate poly
Definition: facBivar.h:37
NULL
#define NULL
Definition: omList.c:10
hasFlag
#define hasFlag(A, F)
Definition: ipid.h:111
sleftv::rtyp
int rtyp
Definition: subexpr.h:92
sleftv::Data
void * Data()
Definition: subexpr.cc:1091
idSize
static int idSize(const ideal id)
Count the effective size of an ideal (without the trailing allocated zero-elements) ...
Definition: ideals.h:46
undoSetOptionRedSB
static void undoSetOptionRedSB()
Definition: tropicalVariety.cc:30
Werror
void Werror(const char *fmt,...)
Definition: reporter.cc:199
static void undoSetOptionRedSB ( )
static

Definition at line 30 of file tropicalVariety.cc.

31 {
32  SI_RESTORE_OPT(bitsetSave1,bitsetSave2);
33 }
SI_RESTORE_OPT
#define SI_RESTORE_OPT(A, B)
Definition: options.h:22
bitsetSave2
BITSET bitsetSave2
Definition: tropicalVariety.cc:14
bitsetSave1
BITSET bitsetSave1
Definition: tropicalVariety.cc:14

Variable Documentation

BITSET bitsetSave1

Definition at line 14 of file tropicalVariety.cc.

BITSET bitsetSave2

Definition at line 14 of file tropicalVariety.cc.

Generated on Mon Aug 3 2015 12:41:10 by   doxygen 1.8.9.1 for Singular a9b8ed0|M