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
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
|
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html><head><meta http-equiv="Content-Type" content="text/html;charset=UTF-8">
<title>polylib: eval_ehrhart.c Source File</title>
<link href="doxygen.css" rel="stylesheet" type="text/css">
<link href="tabs.css" rel="stylesheet" type="text/css">
</head><body>
<!-- Generated by Doxygen 1.5.6 -->
<div class="navigation" id="top">
<div class="tabs">
<ul>
<li><a href="main.html"><span>Main Page</span></a></li>
<li><a href="annotated.html"><span>Classes</span></a></li>
<li class="current"><a href="files.html"><span>Files</span></a></li>
</ul>
</div>
<h1>eval_ehrhart.c</h1><a href="eval__ehrhart_8c.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/************************************************/</span>
<a name="l00002"></a>00002 <span class="comment">/* eval_ehrhart.c */</span>
<a name="l00003"></a>00003 <span class="comment">/* functions to evaluate an Ehrhart polynomial. */</span>
<a name="l00004"></a>00004 <span class="comment">/* written by Emmanuel Jeannot (c) 1997. */</span>
<a name="l00005"></a>00005 <span class="comment">/* Emmanuel.Jeannot@ens-lyon.fr */</span>
<a name="l00006"></a>00006 <span class="comment">/* http://www.ens-lyon.fr/~ejeannot */</span>
<a name="l00007"></a>00007 <span class="comment">/* */</span>
<a name="l00008"></a>00008 <span class="comment">/* modified 1998, 2000, Vincent Loechner */</span>
<a name="l00009"></a>00009 <span class="comment">/* (ArithmetiqueLib, Param_Names) */</span>
<a name="l00010"></a>00010 <span class="comment">/************************************************/</span>
<a name="l00011"></a>00011
<a name="l00012"></a>00012 <span class="preprocessor">#include <stdio.h></span>
<a name="l00013"></a>00013 <span class="preprocessor">#include <math.h></span>
<a name="l00014"></a>00014 <span class="preprocessor">#include <<a class="code" href="assert_8h.html">assert.h</a>></span>
<a name="l00015"></a>00015 <span class="preprocessor">#include <stdlib.h></span>
<a name="l00016"></a>00016
<a name="l00017"></a>00017 <span class="preprocessor">#include <<a class="code" href="polylib_8h.html">polylib/polylib.h</a>></span>
<a name="l00018"></a>00018
<a name="l00019"></a>00019 <span class="comment">/* #define EVAL_EHRHART_DEBUG */</span>
<a name="l00020"></a>00020
<a name="l00021"></a>00021 <span class="comment">/********************************************************/</span>
<a name="l00022"></a>00022 <span class="comment">/* function in domain */</span>
<a name="l00023"></a>00023 <span class="comment">/* check if the parameters in list_args */</span>
<a name="l00024"></a>00024 <span class="comment">/* verifies the constraints of Domain P */</span>
<a name="l00025"></a>00025 <span class="comment">/********************************************************/</span>
<a name="l00026"></a><a class="code" href="eval__ehrhart_8h.html#eaa7a17d905454d792bc8bb5409f941d">00026</a> <span class="keywordtype">int</span> <a class="code" href="eval__ehrhart_8c.html#eaa7a17d905454d792bc8bb5409f941d">in_domain</a>(<a class="code" href="structpolyhedron.html">Polyhedron</a> *P, Value *list_args) {
<a name="l00027"></a>00027
<a name="l00028"></a>00028 <span class="keywordtype">int</span> col,row;
<a name="l00029"></a>00029 Value v; <span class="comment">/* value of the constraint of a row when</span>
<a name="l00030"></a>00030 <span class="comment"> parameters are instanciated*/</span>
<a name="l00031"></a>00031
<a name="l00032"></a>00032 <span class="keywordflow">if</span>( !P )
<a name="l00033"></a>00033 <span class="keywordflow">return</span>( 0 );
<a name="l00034"></a>00034
<a name="l00035"></a>00035 <a class="code" href="polyhedron_8h.html#54a764faf683e794f4db092c8d050104">POL_ENSURE_INEQUALITIES</a>(P);
<a name="l00036"></a>00036
<a name="l00037"></a>00037 <a class="code" href="source_2arith_2arithmetique_8h.html#f71a2ca0294a19cff0cdcbdcc052ee27">value_init</a>(v);
<a name="l00038"></a>00038
<a name="l00039"></a>00039 <span class="comment">/* P->Constraint constraint matrice of polyhedron P */</span>
<a name="l00040"></a>00040 <span class="keywordflow">for</span>(row=0;row<P-><a class="code" href="structpolyhedron.html#2db05293e731089b9198d45cf027bd1c">NbConstraints</a>;row++) {
<a name="l00041"></a>00041 <a class="code" href="source_2arith_2arithmetique_8h.html#864613888dc46f15679aa4f63e468f89">value_assign</a>(v,P-><a class="code" href="structpolyhedron.html#b2691b8eb5d2dc43a509ac2094d2bc73">Constraint</a>[row][P-><a class="code" href="structpolyhedron.html#2a02cea8b7ba3dde415041b8b2373bc8">Dimension</a>+1]); <span class="comment">/*constant part*/</span>
<a name="l00042"></a>00042 <span class="keywordflow">for</span>(col=1;col<P-><a class="code" href="structpolyhedron.html#2a02cea8b7ba3dde415041b8b2373bc8">Dimension</a>+1;col++) {
<a name="l00043"></a>00043 <a class="code" href="source_2arith_2arithmetique_8h.html#9900fbd029b36f5887e587642a064a52">value_addmul</a>(v, P-><a class="code" href="structpolyhedron.html#b2691b8eb5d2dc43a509ac2094d2bc73">Constraint</a>[row][col], list_args[col-1]);
<a name="l00044"></a>00044 }
<a name="l00045"></a>00045 <span class="keywordflow">if</span> (<a class="code" href="source_2arith_2arithmetique_8h.html#47d32925340d2dc99ef2d4215080a60d">value_notzero_p</a>(P-><a class="code" href="structpolyhedron.html#b2691b8eb5d2dc43a509ac2094d2bc73">Constraint</a>[row][0])) {
<a name="l00046"></a>00046
<a name="l00047"></a>00047 <span class="comment">/*if v is not >=0 then this constraint is not respected */</span>
<a name="l00048"></a>00048 <span class="keywordflow">if</span> (<a class="code" href="source_2arith_2arithmetique_8h.html#f031da52c8a160fd1ee6e4c793f8c78a">value_neg_p</a>(v)) {
<a name="l00049"></a>00049 <a class="code" href="source_2arith_2arithmetique_8h.html#b9b282921e85a0527d462d331533d619">value_clear</a>(v);
<a name="l00050"></a>00050 <span class="keywordflow">return</span>( <a class="code" href="eval__ehrhart_8c.html#eaa7a17d905454d792bc8bb5409f941d">in_domain</a>(P-><a class="code" href="structpolyhedron.html#ab97126cb9e56431d40eaf4e68b3953d">next</a>, list_args) );
<a name="l00051"></a>00051 }
<a name="l00052"></a>00052 }
<a name="l00053"></a>00053 <span class="keywordflow">else</span> {
<a name="l00054"></a>00054
<a name="l00055"></a>00055 <span class="comment">/*if v is not = 0 then this constraint is not respected */</span>
<a name="l00056"></a>00056 <span class="keywordflow">if</span> (<a class="code" href="source_2arith_2arithmetique_8h.html#47d32925340d2dc99ef2d4215080a60d">value_notzero_p</a>(v)) {
<a name="l00057"></a>00057 <a class="code" href="source_2arith_2arithmetique_8h.html#b9b282921e85a0527d462d331533d619">value_clear</a>(v);
<a name="l00058"></a>00058 <span class="keywordflow">return</span>( <a class="code" href="eval__ehrhart_8c.html#eaa7a17d905454d792bc8bb5409f941d">in_domain</a>(P-><a class="code" href="structpolyhedron.html#ab97126cb9e56431d40eaf4e68b3953d">next</a>, list_args) );
<a name="l00059"></a>00059 }
<a name="l00060"></a>00060 }
<a name="l00061"></a>00061 }
<a name="l00062"></a>00062
<a name="l00063"></a>00063 <span class="comment">/* if not return before this point => all the constraints are respected */</span>
<a name="l00064"></a>00064 <a class="code" href="source_2arith_2arithmetique_8h.html#b9b282921e85a0527d462d331533d619">value_clear</a>(v);
<a name="l00065"></a>00065 <span class="keywordflow">return</span> 1;
<a name="l00066"></a>00066 } <span class="comment">/* in_domain */</span>
<a name="l00067"></a>00067
<a name="l00068"></a>00068 <span class="comment">/****************************************************/</span>
<a name="l00069"></a>00069 <span class="comment">/* function compute enode */</span>
<a name="l00070"></a>00070 <span class="comment">/* compute the value of enode p with parameters */</span>
<a name="l00071"></a>00071 <span class="comment">/* list "list_args */</span>
<a name="l00072"></a>00072 <span class="comment">/* compute the polynomial or the periodic */</span>
<a name="l00073"></a>00073 <span class="comment">/****************************************************/</span>
<a name="l00074"></a>00074
<a name="l00075"></a><a class="code" href="eval__ehrhart_8c.html#f8ffd38faa2d17abcb0ed78b6d3f53b0">00075</a> <span class="keyword">static</span> <span class="keywordtype">double</span> <a class="code" href="eval__ehrhart_8c.html#f8ffd38faa2d17abcb0ed78b6d3f53b0">compute_enode</a>(<a class="code" href="struct__enode.html">enode</a> *<a class="code" href="vector_8c.html#a45b2e3dcf291527c5aedc420819adfc">p</a>, Value *list_args) {
<a name="l00076"></a>00076
<a name="l00077"></a>00077 <span class="keywordtype">int</span> i;
<a name="l00078"></a>00078 Value <a class="code" href="polyparam_8c.html#742204794ea328ba293fe59cec79b990">m</a>, param;
<a name="l00079"></a>00079 <span class="keywordtype">double</span> res=0.0;
<a name="l00080"></a>00080
<a name="l00081"></a>00081 <span class="keywordflow">if</span> (!p)
<a name="l00082"></a>00082 <span class="keywordflow">return</span>(0.);
<a name="l00083"></a>00083
<a name="l00084"></a>00084 <a class="code" href="source_2arith_2arithmetique_8h.html#f71a2ca0294a19cff0cdcbdcc052ee27">value_init</a>(m);
<a name="l00085"></a>00085 <a class="code" href="source_2arith_2arithmetique_8h.html#f71a2ca0294a19cff0cdcbdcc052ee27">value_init</a>(param);
<a name="l00086"></a>00086
<a name="l00087"></a>00087 <span class="keywordflow">if</span> (p-><a class="code" href="struct__enode.html#36787a943f0f24530a8934c48a69d282">type</a> == <a class="code" href="types_8h.html#01e5f42033d0f34c1a182a25a7409feda486649798da33c51389df14eaccd3ce">polynomial</a>) {
<a name="l00088"></a>00088 <span class="keywordflow">if</span> (p-><a class="code" href="struct__enode.html#ba9a89a7893638f2ad9d77c9d0247797">size</a> > 1)
<a name="l00089"></a>00089 <a class="code" href="source_2arith_2arithmetique_8h.html#864613888dc46f15679aa4f63e468f89">value_assign</a>(param,list_args[p-><a class="code" href="struct__enode.html#1230adfbd7f0ca1f8628a03be4e6ef35">pos</a>-1]);
<a name="l00090"></a>00090
<a name="l00091"></a>00091 <span class="comment">/* Compute the polynomial using Horner's rule */</span>
<a name="l00092"></a>00092 <span class="keywordflow">for</span> (i=p-><a class="code" href="struct__enode.html#ba9a89a7893638f2ad9d77c9d0247797">size</a>-1;i>0;i--) {
<a name="l00093"></a>00093 res +=<a class="code" href="eval__ehrhart_8c.html#3f1631cacb65ff7aa29df190b4d30d00">compute_evalue</a>(&p-><a class="code" href="struct__enode.html#215794f8b8ab87573bd46f14a8e916d4">arr</a>[i],list_args);
<a name="l00094"></a>00094 res *=VALUE_TO_DOUBLE(param);
<a name="l00095"></a>00095 }
<a name="l00096"></a>00096 res +=<a class="code" href="eval__ehrhart_8c.html#3f1631cacb65ff7aa29df190b4d30d00">compute_evalue</a>(&p-><a class="code" href="struct__enode.html#215794f8b8ab87573bd46f14a8e916d4">arr</a>[0],list_args);
<a name="l00097"></a>00097 }
<a name="l00098"></a>00098 <span class="keywordflow">else</span> <span class="keywordflow">if</span> (p-><a class="code" href="struct__enode.html#36787a943f0f24530a8934c48a69d282">type</a> == <a class="code" href="types_8h.html#01e5f42033d0f34c1a182a25a7409fed3cb930eb7594d2236ab4fc6dd622c207">periodic</a>) {
<a name="l00099"></a>00099 <a class="code" href="source_2arith_2arithmetique_8h.html#864613888dc46f15679aa4f63e468f89">value_assign</a>(m,list_args[p-><a class="code" href="struct__enode.html#1230adfbd7f0ca1f8628a03be4e6ef35">pos</a>-1]);
<a name="l00100"></a>00100
<a name="l00101"></a>00101 <span class="comment">/* Choose the right element of the periodic */</span>
<a name="l00102"></a>00102 <a class="code" href="source_2arith_2arithmetique_8h.html#8cc56567a4a29271559ac0fd5f6c5bfa">value_set_si</a>(param,p-><a class="code" href="struct__enode.html#ba9a89a7893638f2ad9d77c9d0247797">size</a>);
<a name="l00103"></a>00103 <a class="code" href="source_2arith_2arithmetique_8h.html#60a1720f53f9951170f929062609ddcf">value_pmodulus</a>(m,m,param);
<a name="l00104"></a>00104 res = <a class="code" href="eval__ehrhart_8c.html#3f1631cacb65ff7aa29df190b4d30d00">compute_evalue</a>(&p-><a class="code" href="struct__enode.html#215794f8b8ab87573bd46f14a8e916d4">arr</a>[VALUE_TO_INT(m)],list_args);
<a name="l00105"></a>00105 }
<a name="l00106"></a>00106 <a class="code" href="source_2arith_2arithmetique_8h.html#b9b282921e85a0527d462d331533d619">value_clear</a>(m);
<a name="l00107"></a>00107 <a class="code" href="source_2arith_2arithmetique_8h.html#b9b282921e85a0527d462d331533d619">value_clear</a>(param);
<a name="l00108"></a>00108 <span class="keywordflow">return</span> res;
<a name="l00109"></a>00109 } <span class="comment">/* compute_enode */</span>
<a name="l00110"></a>00110
<a name="l00111"></a>00111 <span class="comment">/*************************************************/</span>
<a name="l00112"></a>00112 <span class="comment">/* return the value of Ehrhart Polynomial */</span>
<a name="l00113"></a>00113 <span class="comment">/* It returns a double, because since it is */</span>
<a name="l00114"></a>00114 <span class="comment">/* a recursive function, some intermediate value */</span>
<a name="l00115"></a>00115 <span class="comment">/* might not be integral */</span>
<a name="l00116"></a>00116 <span class="comment">/*************************************************/</span>
<a name="l00117"></a>00117
<a name="l00118"></a><a class="code" href="eval__ehrhart_8h.html#3f1631cacb65ff7aa29df190b4d30d00">00118</a> <span class="keywordtype">double</span> <a class="code" href="eval__ehrhart_8c.html#3f1631cacb65ff7aa29df190b4d30d00">compute_evalue</a>(<a class="code" href="struct__evalue.html">evalue</a> *e,Value *list_args) {
<a name="l00119"></a>00119
<a name="l00120"></a>00120 <span class="keywordtype">double</span> res;
<a name="l00121"></a>00121
<a name="l00122"></a>00122 <span class="keywordflow">if</span> (<a class="code" href="source_2arith_2arithmetique_8h.html#47d32925340d2dc99ef2d4215080a60d">value_notzero_p</a>(e-><a class="code" href="struct__evalue.html#e4aada7b06d44dcbf9f4d2daf1834e29">d</a>)) {
<a name="l00123"></a>00123 <span class="keywordflow">if</span> (<a class="code" href="source_2arith_2arithmetique_8h.html#5f258141f9e07b66da01ff7f4fe3d56d">value_notone_p</a>(e-><a class="code" href="struct__evalue.html#e4aada7b06d44dcbf9f4d2daf1834e29">d</a>))
<a name="l00124"></a>00124 res = VALUE_TO_DOUBLE(e->x.n) / VALUE_TO_DOUBLE(e-><a class="code" href="struct__evalue.html#e4aada7b06d44dcbf9f4d2daf1834e29">d</a>);
<a name="l00125"></a>00125 <span class="keywordflow">else</span>
<a name="l00126"></a>00126 res = VALUE_TO_DOUBLE(e->x.n);
<a name="l00127"></a>00127 }
<a name="l00128"></a>00128 <span class="keywordflow">else</span>
<a name="l00129"></a>00129 res = <a class="code" href="eval__ehrhart_8c.html#f8ffd38faa2d17abcb0ed78b6d3f53b0">compute_enode</a>(e->x.p,list_args);
<a name="l00130"></a>00130 <span class="keywordflow">return</span> res;
<a name="l00131"></a>00131 } <span class="comment">/* compute_evalue */</span>
<a name="l00132"></a>00132
<a name="l00133"></a>00133
<a name="l00134"></a>00134 <span class="comment">/****************************************************/</span>
<a name="l00135"></a>00135 <span class="comment">/* function compute_poly : */</span>
<a name="l00136"></a>00136 <span class="comment">/* Check for the good validity domain */</span>
<a name="l00137"></a>00137 <span class="comment">/* return the number of point in the Polyhedron */</span>
<a name="l00138"></a>00138 <span class="comment">/* in allocated memory */</span>
<a name="l00139"></a>00139 <span class="comment">/* Using the Ehrhart pseudo-polynomial */</span>
<a name="l00140"></a>00140 <span class="comment">/****************************************************/</span>
<a name="l00141"></a><a class="code" href="eval__ehrhart_8h.html#f970e7f69cbf548188a8f3754001a150">00141</a> Value *<a class="code" href="eval__ehrhart_8c.html#f970e7f69cbf548188a8f3754001a150">compute_poly</a>(<a class="code" href="struct__enumeration.html">Enumeration</a> *en,Value *list_args) {
<a name="l00142"></a>00142
<a name="l00143"></a>00143 Value *tmp;
<a name="l00144"></a>00144 <span class="comment">/* double d; int i; */</span>
<a name="l00145"></a>00145
<a name="l00146"></a>00146 tmp = (Value *) malloc (<span class="keyword">sizeof</span>(Value));
<a name="l00147"></a>00147 <a class="code" href="assert_8h.html#07d17d6d5d1074c0969bc5d3c3d1d84a">assert</a>(tmp != NULL);
<a name="l00148"></a>00148 <a class="code" href="source_2arith_2arithmetique_8h.html#f71a2ca0294a19cff0cdcbdcc052ee27">value_init</a>(*tmp);
<a name="l00149"></a>00149 <a class="code" href="source_2arith_2arithmetique_8h.html#8cc56567a4a29271559ac0fd5f6c5bfa">value_set_si</a>(*tmp,0);
<a name="l00150"></a>00150
<a name="l00151"></a>00151 <span class="keywordflow">if</span>(!en)
<a name="l00152"></a>00152 <span class="keywordflow">return</span>(tmp); <span class="comment">/* no ehrhart polynomial */</span>
<a name="l00153"></a>00153 <span class="keywordflow">if</span>(en-><a class="code" href="struct__enumeration.html#57e35419d63ab12898542ef58d8decac">ValidityDomain</a>) {
<a name="l00154"></a>00154 <span class="keywordflow">if</span>(!en-><a class="code" href="struct__enumeration.html#57e35419d63ab12898542ef58d8decac">ValidityDomain</a>-><a class="code" href="structpolyhedron.html#2a02cea8b7ba3dde415041b8b2373bc8">Dimension</a>) { <span class="comment">/* no parameters */</span>
<a name="l00155"></a>00155 <a class="code" href="source_2arith_2arithmetique_8h.html#3a9e8def0399e848b9f715238e2378fa">value_set_double</a>(*tmp,<a class="code" href="eval__ehrhart_8c.html#3f1631cacb65ff7aa29df190b4d30d00">compute_evalue</a>(&en-><a class="code" href="struct__enumeration.html#a499b062a10ce16356e9b52d25cae0d9">EP</a>,list_args)+.25);
<a name="l00156"></a>00156 <span class="keywordflow">return</span>(tmp);
<a name="l00157"></a>00157 }
<a name="l00158"></a>00158 }
<a name="l00159"></a>00159 <span class="keywordflow">else</span>
<a name="l00160"></a>00160 <span class="keywordflow">return</span>(tmp); <span class="comment">/* no Validity Domain */</span>
<a name="l00161"></a>00161 <span class="keywordflow">while</span>(en) {
<a name="l00162"></a>00162 <span class="keywordflow">if</span>(<a class="code" href="eval__ehrhart_8c.html#eaa7a17d905454d792bc8bb5409f941d">in_domain</a>(en-><a class="code" href="struct__enumeration.html#57e35419d63ab12898542ef58d8decac">ValidityDomain</a>,list_args)) {
<a name="l00163"></a>00163
<a name="l00164"></a>00164 <span class="preprocessor">#ifdef EVAL_EHRHART_DEBUG</span>
<a name="l00165"></a>00165 <span class="preprocessor"></span> <a class="code" href="polyparam_8c.html#fa315bc856959f9a781359f69b74d7a4">Print_Domain</a>(stdout,en-><a class="code" href="struct__enumeration.html#57e35419d63ab12898542ef58d8decac">ValidityDomain</a>,NULL);
<a name="l00166"></a>00166 <a class="code" href="ehrhart_8c.html#99718adf05e065181f35f5ca847f94b0">print_evalue</a>(stdout,&en-><a class="code" href="struct__enumeration.html#a499b062a10ce16356e9b52d25cae0d9">EP</a>,NULL);
<a name="l00167"></a>00167 <span class="preprocessor">#endif</span>
<a name="l00168"></a>00168 <span class="preprocessor"></span>
<a name="l00169"></a>00169 <span class="comment">/* d = compute_evalue(&en->EP,list_args);</span>
<a name="l00170"></a>00170 <span class="comment"> i = d;</span>
<a name="l00171"></a>00171 <span class="comment"> printf("(double)%lf = %d\n", d, i ); */</span>
<a name="l00172"></a>00172 <a class="code" href="source_2arith_2arithmetique_8h.html#3a9e8def0399e848b9f715238e2378fa">value_set_double</a>(*tmp,<a class="code" href="eval__ehrhart_8c.html#3f1631cacb65ff7aa29df190b4d30d00">compute_evalue</a>(&en-><a class="code" href="struct__enumeration.html#a499b062a10ce16356e9b52d25cae0d9">EP</a>,list_args)+.25);
<a name="l00173"></a>00173 <span class="keywordflow">return</span>(tmp);
<a name="l00174"></a>00174 }
<a name="l00175"></a>00175 <span class="keywordflow">else</span>
<a name="l00176"></a>00176 en=en-><a class="code" href="struct__enumeration.html#d2fd42474e015c9abc6aea0e2bea08b8">next</a>;
<a name="l00177"></a>00177 }
<a name="l00178"></a>00178 <a class="code" href="source_2arith_2arithmetique_8h.html#8cc56567a4a29271559ac0fd5f6c5bfa">value_set_si</a>(*tmp,0);
<a name="l00179"></a>00179 <span class="keywordflow">return</span>(tmp); <span class="comment">/* no compatible domain with the arguments */</span>
<a name="l00180"></a>00180 } <span class="comment">/* compute_poly */</span>
<a name="l00181"></a>00181
<a name="l00182"></a>00182
<a name="l00183"></a>00183
<a name="l00184"></a>00184
</pre></div></div>
<hr size="1"><address style="text-align: right;"><small>Generated on Tue Sep 15 18:33:59 2009 for polylib by
<a href="http://www.doxygen.org/index.html">
<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.6 </small></address>
</body>
</html>
|
