Remove the "** FIXME: THIS DOES NOT ACTUALLY WORK YET" lines.
[synfig.git] / ETL / trunk / test / spline.cpp
1 /*! ========================================================================
2 ** Extended Template and Library Test Suite
3 ** Spline 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/spline>
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 bspline_basic_test(void)
40 {
41         int ret=0;
42         float f;
43
44         bspline<float> BSpline;
45         etl::clock timer;
46         double t;
47
48         *BSpline.cpoints().insert(BSpline.cpoints().end())=0;
49         *BSpline.cpoints().insert(BSpline.cpoints().end())=-1;
50         *BSpline.cpoints().insert(BSpline.cpoints().end())=0;
51         *BSpline.cpoints().insert(BSpline.cpoints().end())=1;
52         *BSpline.cpoints().insert(BSpline.cpoints().end())=0;
53
54         BSpline.set_m(4);
55         BSpline.reset_knots();
56
57         integral<bspline<float> > inte(BSpline);
58
59
60         /*
61         for(f=0.0;f<1.001;f+=0.05)
62                 fprintf(stderr,"BSpline(%f)= %f\n",f,BSpline(f));
63         */
64
65         fprintf(stderr,"integral of BSpline() on [0,1] = %f\n",inte(0,1.0));
66
67
68         for(f=0.0f,timer.reset();f<1.001f;f+=0.000005f)
69         {
70                 t+=BSpline(f)+BSpline(f+0.1f);
71                 t+=BSpline(f)+BSpline(f+0.1f);
72                 t+=BSpline(f)+BSpline(f+0.1f);
73                 t+=BSpline(f)+BSpline(f+0.1f);
74                 t+=BSpline(f)+BSpline(f+0.1f);
75                 t+=BSpline(f)+BSpline(f+0.1f);
76                 t+=BSpline(f)+BSpline(f+0.1f);
77                 t+=BSpline(f)+BSpline(f+0.1f);
78                 t+=BSpline(f)+BSpline(f+0.1f);
79                 t+=BSpline(f)+BSpline(f+0.1f);
80                 t+=BSpline(f)+BSpline(f+0.1f);
81                 t+=BSpline(f)+BSpline(f+0.1f);
82         }
83         t=timer();
84
85         fprintf(stderr,"BSpline time=%f milliseconds\n",t*1000);
86         return ret;
87 }
88
89 /* === E N T R Y P O I N T ================================================= */
90
91 int main()
92 {
93         int error=0;
94
95         error+=bspline_basic_test();
96
97         return error;
98 }
99