Fix 1682016: default to 700x510 instead; big enough to be useful, not too big to...
[synfig.git] / ETL / trunk / test / hermite.cpp
1 /*! ========================================================================
2 ** Extended Template and Library Test Suite
3 ** Hermite Curve Test
4 ** $Id$
5 **
6 ** Copyright (c) 2002 Robert B. Quattlebaum Jr.
7 **
8 ** This package is free software; you can redistribute it and/or
9 ** modify it under the terms of the GNU General Public License as
10 ** published by the Free Software Foundation; either version 2 of
11 ** the License, or (at your option) any later version.
12 **
13 ** This package is distributed in the hope that it will be useful,
14 ** but WITHOUT ANY WARRANTY; without even the implied warranty of
15 ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
16 ** General Public License for more details.
17 **
18 ** === N O T E S ===========================================================
19 **
20 ** ========================================================================= */
21
22 /* === H E A D E R S ======================================================= */
23
24 #include <ETL/hermite>
25 #include <ETL/angle>
26 #include <ETL/clock>
27 #include <ETL/calculus>
28 #include <stdio.h>
29
30 /* === M A C R O S ========================================================= */
31
32 using namespace etl;
33
34 /* === C L A S S E S ======================================================= */
35
36
37 /* === P R O C E D U R E S ================================================= */
38
39 int basic_test(void)
40 {
41         int ret=0;
42         float f;
43
44         hermite<float> Hermie;
45         etl::clock timer;
46         double t;
47
48         Hermie.p1()=0;
49         Hermie.t1()=1;
50         Hermie.p2()=0;
51         Hermie.t2()=1;
52
53         Hermie.sync();
54
55         integral<hermite<float> > inte(Hermie);
56
57
58         fprintf(stderr,"integral of curve() on [0,1] = %f\n",inte(0,1.0));
59         fprintf(stderr,"integral of curve() on [-1,3] = %f\n",inte(-1.0,3.0));
60         Hermie.set_rs(-1.0,7.0);
61         inte=integral<hermite<float> >(Hermie);
62         fprintf(stderr,"integral of curve()[%f,%f] on [-1,7] = %f\n",Hermie.get_r(),Hermie.get_s(),inte(-1.0,7.0));
63         fprintf(stderr,"integral of curve()[%f,%f] on [0,1] = %f\n",Hermie.get_r(),Hermie.get_s(),inte(0,1.0));
64         Hermie.set_rs(0.0,1.0);
65
66
67         for(f=0.0f,timer.reset();f<1.001f;f+=0.000005f)
68         {
69                 t+=Hermie(f)+Hermie(f+0.1f);
70                 t+=Hermie(f)+Hermie(f+0.1f);
71                 t+=Hermie(f)+Hermie(f+0.1f);
72                 t+=Hermie(f)+Hermie(f+0.1f);
73                 t+=Hermie(f)+Hermie(f+0.1f);
74                 t+=Hermie(f)+Hermie(f+0.1f);
75                 t+=Hermie(f)+Hermie(f+0.1f);
76                 t+=Hermie(f)+Hermie(f+0.1f);
77                 t+=Hermie(f)+Hermie(f+0.1f);
78                 t+=Hermie(f)+Hermie(f+0.1f);
79                 t+=Hermie(f)+Hermie(f+0.1f);
80                 t+=Hermie(f)+Hermie(f+0.1f);
81         }
82         t=timer();
83
84         fprintf(stderr,"time=%f milliseconds\n",t*1000);
85         return ret;
86 }
87
88 int angle_test(void)
89 {
90         int ret=0;
91         float f;
92
93         hermite<angle> Hermie;
94         etl::clock timer;
95         angle tmp;
96         double t;
97
98         Hermie.p1()=angle::degrees(0);
99         Hermie.t1()=angle::degrees(45);
100
101         Hermie.p2()=angle::degrees(-45);
102         Hermie.t2()=angle::degrees(180);
103
104         Hermie.sync();
105
106
107         for(f=0.0f,timer.reset();f<1.001f;f+=0.000005f)
108         {
109                 tmp+=Hermie(f)+Hermie(f+0.1f);
110                 tmp+=Hermie(f)+Hermie(f+0.1f);
111                 tmp+=Hermie(f)+Hermie(f+0.1f);
112                 tmp+=Hermie(f)+Hermie(f+0.1f);
113                 tmp+=Hermie(f)+Hermie(f+0.1f);
114                 tmp+=Hermie(f)+Hermie(f+0.1f);
115                 tmp+=Hermie(f)+Hermie(f+0.1f);
116                 tmp+=Hermie(f)+Hermie(f+0.1f);
117                 tmp+=Hermie(f)+Hermie(f+0.1f);
118                 tmp+=Hermie(f)+Hermie(f+0.1f);
119                 tmp+=Hermie(f)+Hermie(f+0.1f);
120                 tmp+=Hermie(f)+Hermie(f+0.1f);
121         }
122         t=timer();
123
124         fprintf(stderr,"angle time=%f milliseconds\n",t*1000);
125
126         return ret;
127 }
128
129 int fixed_test(void)
130 {
131         int ret=0;
132         float f;
133
134         hermite<float,fixed> Hermie;
135         etl::clock timer;
136         double t;
137
138         Hermie.p1()=0;
139         Hermie.t1()=1;
140         Hermie.p2()=0;
141         Hermie.t2()=1;
142
143         Hermie.sync();
144
145
146
147         for(f=0.0f,timer.reset();f<1.001f;f+=0.005f)
148         {
149                 t+=Hermie(f)+Hermie(f+0.1f);
150                 t+=Hermie(f)+Hermie(f+0.1f);
151                 t+=Hermie(f)+Hermie(f+0.1f);
152                 t+=Hermie(f)+Hermie(f+0.1f);
153                 t+=Hermie(f)+Hermie(f+0.1f);
154                 t+=Hermie(f)+Hermie(f+0.1f);
155                 t+=Hermie(f)+Hermie(f+0.1f);
156                 t+=Hermie(f)+Hermie(f+0.1f);
157                 t+=Hermie(f)+Hermie(f+0.1f);
158                 t+=Hermie(f)+Hermie(f+0.1f);
159                 t+=Hermie(f)+Hermie(f+0.1f);
160                 t+=Hermie(f)+Hermie(f+0.1f);
161         }
162         t=timer();
163
164         fprintf(stderr,"time=%f milliseconds\n",t*1000);
165         return ret;
166 }
167
168
169 int anglefixed_test(void)
170 {
171         int ret=0;
172         float f;
173
174         hermite<angle,fixed> Hermie;
175         etl::clock timer;
176         angle tmp;
177         double t;
178
179         Hermie.p1()=angle::degrees(0);
180         Hermie.t1()=angle::degrees(45);
181
182         Hermie.p2()=angle::degrees(-45);
183         Hermie.t2()=angle::degrees(180);
184
185         Hermie.sync();
186
187
188         for(f=0.0f,timer.reset();f<1.001f;f+=0.0005f)
189         {
190                 tmp+=Hermie(f)+Hermie(f+0.1f);
191                 tmp+=Hermie(f)+Hermie(f+0.1f);
192                 tmp+=Hermie(f)+Hermie(f+0.1f);
193                 tmp+=Hermie(f)+Hermie(f+0.1f);
194                 tmp+=Hermie(f)+Hermie(f+0.1f);
195                 tmp+=Hermie(f)+Hermie(f+0.1f);
196                 tmp+=Hermie(f)+Hermie(f+0.1f);
197                 tmp+=Hermie(f)+Hermie(f+0.1f);
198                 tmp+=Hermie(f)+Hermie(f+0.1f);
199                 tmp+=Hermie(f)+Hermie(f+0.1f);
200                 tmp+=Hermie(f)+Hermie(f+0.1f);
201                 tmp+=Hermie(f)+Hermie(f+0.1f);
202         }
203         t=timer();
204
205         fprintf(stderr,"angle fixed time=%f milliseconds\n",t*1000);
206
207         return ret;
208 }
209
210 int float_intersection_test()
211 {
212         int ret=0;
213
214         hermite<float> curve1(0,1,0,1);
215         hermite<float> curve2(-1,2,-1,-2);
216         double t1,t2;
217         float d;
218
219         t1=curve1.intersect(curve2);
220         t2=curve2.intersect(curve1);
221
222         d=curve1(t1)-curve2(t2);
223
224         fprintf(stderr,"float:Intersection difference: %f (t1=%f, t2=%f)\n",d,t1,t2);
225
226         if(d>0.01)
227         {
228                 fprintf(stderr,"float:FAILED INTERSECTION TEST.\n");
229                 ret++;
230         }
231
232         return ret;
233 }
234
235 /* === E N T R Y P O I N T ================================================= */
236
237 int main()
238 {
239         int error=0;
240
241         error+=basic_test();
242         error+=angle_test();
243         error+=fixed_test();
244         error+=float_intersection_test();
245         return error;
246 }