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
|
/*! ========================================================================
** Extended Template and Library Test Suite
** Spline Curve Test
** $Id$
**
** Copyright (c) 2002 Robert B. Quattlebaum Jr.
** Copyright (c) 2010 Nikita Kitaev
**
** This package is free software; you can redistribute it and/or
** modify it under the terms of the GNU General Public License as
** published by the Free Software Foundation; either version 2 of
** the License, or (at your option) any later version.
**
** This package 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.
**
** === N O T E S ===========================================================
**
** ========================================================================= */
/* === H E A D E R S ======================================================= */
#include <ETL/spline>
#include <ETL/angle>
#include <ETL/clock>
#include <ETL/calculus>
#include <stdio.h>
/* === M A C R O S ========================================================= */
using namespace etl;
/* === C L A S S E S ======================================================= */
/* === P R O C E D U R E S ================================================= */
int bspline_basic_test(void)
{
int ret=0;
float f;
bspline<float> BSpline;
etl::clock timer;
double t = 0.0;
BSpline.cpoints().insert(BSpline.cpoints().end(), 0.0);
BSpline.cpoints().insert(BSpline.cpoints().end(), -1.0);
BSpline.cpoints().insert(BSpline.cpoints().end(), 0.0);
BSpline.cpoints().insert(BSpline.cpoints().end(), 1.0);
BSpline.cpoints().insert(BSpline.cpoints().end(), 0.0);
BSpline.set_m(4);
BSpline.reset_knots();
integral<bspline<float> > inte(BSpline);
/*
for(f=0.0;f<1.001;f+=0.05)
fprintf(stderr,"BSpline(%f)= %f\n",f,BSpline(f));
*/
fprintf(stderr,"integral of BSpline() on [0,1] = %f\n",inte(0,1.0));
for(f=0.0f,timer.reset();f<1.001f;f+=0.000005f)
{
t+=BSpline(f)+BSpline(f+0.1f);
t+=BSpline(f)+BSpline(f+0.1f);
t+=BSpline(f)+BSpline(f+0.1f);
t+=BSpline(f)+BSpline(f+0.1f);
t+=BSpline(f)+BSpline(f+0.1f);
t+=BSpline(f)+BSpline(f+0.1f);
t+=BSpline(f)+BSpline(f+0.1f);
t+=BSpline(f)+BSpline(f+0.1f);
t+=BSpline(f)+BSpline(f+0.1f);
t+=BSpline(f)+BSpline(f+0.1f);
t+=BSpline(f)+BSpline(f+0.1f);
t+=BSpline(f)+BSpline(f+0.1f);
}
t=timer();
fprintf(stderr,"BSpline time=%f milliseconds\n",t*1000);
return ret;
}
/* === E N T R Y P O I N T ================================================= */
int main()
{
int error=0;
error+=bspline_basic_test();
return error;
}
|
