1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
|
/* PolyTree.h --
Copyright 2004 Jesus A. De Loera, David Haws, Raymond Hemmecke,
Peter Huggins, Jeremy Tauzer, Ruriko Yoshida
This file is part of LattE.
LattE is free software; you can redistribute it and/or modify it
under the terms of the version 2 of the GNU General Public License
as published by the Free Software Foundation.
LattE is distributed in the hope that it will be useful, but
WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
General Public License for more details.
You should have received a copy of the GNU General Public License
along with LattE; if not, write to the Free Software Foundation,
Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA.
*/
#ifndef POLYTREE_H
#define POLYTREE_H 1
#include <fstream>
#include "latte_ntl.h"
using namespace std;
// Definitions of the types of PolyTree structures
#define POLYTREE_ADD 0
#define POLYTREE_MUL 1
#define POLYTREE_DIV 5
#define POLYTREE_T_NODE 3
#define POLYTREE_S_NODE 2
#define POLYTREE_EXP 4
// Define the types of generators of cones, either (1-rt) or (r-t)
#define ONE_SUB_RT 0
#define R_SUB_T 1
// This is used to hold all the information of Taylor_Expansion function
struct Taylor_Parameters
{
ZZ *Result;
int Degree_of_Expansion;
ZZ *Ten_Power;
};
class PolyTree_Node
{
public:
PolyTree_Node () { Taylor_Expansion_Result_Dirty = 1; };
virtual ~PolyTree_Node () {};
PolyTree_Node **Children;
// 0 mean +, 1 mean *, 2 means S_Node, 3 means T_Node, 4 mean ^ , 5 means /
char Node_Type;
unsigned int Number_of_Children;
virtual int Print(); //Returns 1 if it did print anything, otherwise returns 0
virtual int Check_For_Zero (); //Returns 1 if this node is equal to zero
//This function will return the first Parameters->Degree_of_Expansion number of terms.
//It accepts Taylor_Parameter as input with the Result part of it already
//allocated for enough space for Degree_of_Expansion + 1 terms.
virtual void Taylor_Expansion (Taylor_Parameters *Parameters);
virtual int Print_Rational_Functions_to_File (ofstream &Output_File);
// 1 indicates Taylor_Expansion_Result is invalid
int Taylor_Expansion_Result_Dirty;
ZZ *Taylor_Expansion_Result;
};
//Currently this class is never used
/*class S_Node : public PolyTree_Node
{
public:
ZZ Coefficient;
ZZ Exponent;
int Print ();
};*/
class T_Node : public PolyTree_Node //Also used for ^ type. Exponent holds the exponent.
{
public:
virtual ~T_Node () {};
ZZ Coefficient;
ZZ Exponent;
int Print(); //Returns 1 if it did print anything, otherwise returns 0
int Check_For_Zero (); // Returns 1 if zero, else returns 0
//Same as description in PolyTree_Node
void Taylor_Expansion (Taylor_Parameters *Parameters);
int Print_Rational_Functions_to_File (ofstream &Output_File);
};
// Used to hold information for each generator of a cone.
struct Generator
{
int Form_Type; // 0 means (1-rt) 1 means (r-t)
ZZ R_Exponent;
ZZ T_Exponent;
};
struct Cone_Data
{
int sign; //Holds the sign of this Cone
int order; //Holds the order of this Cone
Generator *Generators_of_Cone; //The denominators
Generator Numerator_Generator;
};
struct PolyTree_Node_List
{
PolyTree_Node *Data;
PolyTree_Node_List *Next;
};
struct T_Node_List
{
T_Node *Data;
T_Node_List *Next;
};
class Node_Controller
{
public:
Node_Controller (int Dimension, int Degree);
~Node_Controller ();
// Returns a pointer to a PolyNode for someone to use
PolyTree_Node *Get_PolyTree_Node ();
// Returns a pointer to a T_Node for someone to use
T_Node *Get_T_Node ();
// Clears all the nodes this Controller controls
void Reset ();
private:
int Dimension_Plus_One;
int Degree_of_Expansion;
PolyTree_Node_List *PolyTree_Node_Head;
PolyTree_Node_List *PolyTree_Node_Unused;
T_Node_List *T_Node_Head;
T_Node_List *T_Node_Unused;
};
#endif
|
