/*! ========================================================================
** Extended Template and Library Test Suite
** Hermite Curve Test
** $Id$
**
** Copyright (c) 2002 Robert B. Quattlebaum Jr.
**
** 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/clock>
#include <ETL/hermite>
#include <ETL/angle>
#include <ETL/fixed>
#include <ETL/surface>
#include <ETL/gaussian>
#include <ETL/calculus>
#include <stdio.h>
/* === M A C R O S ========================================================= */
using namespace etl;
#define HERMITE_TEST_ITERATIONS (100000)
/* === C L A S S E S ======================================================= */
/* === P R O C E D U R E S ================================================= */
template <class Angle>
void angle_cos_speed_test(void)
{
Angle a,b,c,d;
float tmp,tmp2;
for(tmp=-1;tmp<1;tmp+=0.000002f)
{
a=(typename Angle::cos)(tmp);
b=(typename Angle::cos)(tmp);
c=(typename Angle::cos)(tmp);
d=(typename Angle::cos)(tmp);
tmp2=((typename Angle::cos)(a)).get();
tmp2=((typename Angle::cos)(b)).get();
tmp2=((typename Angle::cos)(c)).get();
tmp2=((typename Angle::cos)(d)).get();
}
if (tmp2 == 0) return; // disable unused warning
}
template <class Angle>
void angle_sin_speed_test(void)
{
Angle a,b,c,d;
float tmp,tmp2;
for(tmp=-1.0;tmp<1.0;tmp+=0.000002f)
{
a=(typename Angle::sin)(tmp);
b=(typename Angle::sin)(tmp);
c=(typename Angle::sin)(tmp);
d=(typename Angle::sin)(tmp);
tmp2=((typename Angle::sin)(a)).get();
tmp2=((typename Angle::sin)(b)).get();
tmp2=((typename Angle::sin)(c)).get();
tmp2=((typename Angle::sin)(d)).get();
}
if (tmp2 == 0) return; // disable unused warning
}
template <class Angle>
void angle_tan_speed_test(void)
{
Angle a,b,c,d;
float tmp,tmp2;
for(tmp=-1;tmp<1;tmp+=0.000002f)
{
a=(typename Angle::tan)(tmp);
b=(typename Angle::tan)(tmp);
c=(typename Angle::tan)(tmp);
d=(typename Angle::tan)(tmp);
tmp2=((typename Angle::tan)(a)).get();
tmp2=((typename Angle::tan)(b)).get();
tmp2=((typename Angle::tan)(c)).get();
tmp2=((typename Angle::tan)(d)).get();
}
if (tmp2 == 0) return; // disable unused warning
}
template <class Angle, class mytan>
void angle_atan2_speed_test(void)
{
Angle a,b,c;
float x,y;
for(y=-10;y<10;y+=0.05f)
for(x=-10;x<10;x+=0.05f)
{
a=mytan(y,x);
a=mytan(x,y);
b=mytan(y,x);
b=mytan(x,y);
c=mytan(y,x);
c=mytan(x,y);
a=mytan(y,x);
a=mytan(x,y);
b=mytan(y,x);
b=mytan(x,y);
c=mytan(y,x);
c=mytan(x,y);
}
}
int surface_and_gaussian_blur_test()
{
int ret=0;
etl::clock MyTimer;
float endtime;
{
surface<float> my_surface(1000,1000);
MyTimer.reset();
gaussian_blur(my_surface.begin(),my_surface.end(),30,30);
endtime=MyTimer();
printf("surface_and_gaussian_blur_test<float>: %f seconds\n",endtime);
}
{
surface<double> my_surface(1000,1000);
MyTimer.reset();
gaussian_blur(my_surface.begin(),my_surface.end(),30,30);
endtime=MyTimer();
printf("surface_and_gaussian_blur_test<double>: %f seconds\n",endtime);
}
{
surface<fixed> my_surface(1000,1000);
MyTimer.reset();
gaussian_blur(my_surface.begin(),my_surface.end(),30,30);
endtime=MyTimer();
printf("surface_and_gaussian_blur_test<fixed>: %f seconds\n",endtime);
}
return ret;
}
int hermite_int_test()
{
int ret=0;
hermite<int> Hermie;
hermite<int>::time_type f;
int i;
etl::clock timer;
etl::clock::value_type t;
Hermie.p1()=0;
Hermie.t1()=40000;
Hermie.p2()=0;
Hermie.t2()=40000;
Hermie.sync();
{float t;
for(f=0.0f,i=0,timer.reset();i<HERMITE_TEST_ITERATIONS;i++,f+=0.000005f)
{
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
}
}
t=timer();
printf("hermite<int>:time=%f milliseconds\n",t*1000);
return ret;
}
int hermite_float_test(void)
{
int ret=0;
float f; int i;
hermite<float> Hermie;
etl::clock timer;
double t;
Hermie.p1()=0;
Hermie.t1()=1;
Hermie.p2()=0;
Hermie.t2()=1;
Hermie.sync();
{float t;
for(f=0.0f,i=0,timer.reset();i<HERMITE_TEST_ITERATIONS;i++,f+=0.000005f)
{
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
}
}
t=timer();
printf("hermite<float>:time=%f milliseconds\n",t*1000);
return ret;
}
int hermite_double_test(void)
{
int ret=0,i;
float f;
hermite<double> Hermie;
etl::clock timer;
double t;
Hermie.p1()=0;
Hermie.t1()=1;
Hermie.p2()=0;
Hermie.t2()=1;
Hermie.sync();
for(f=0.0f,i=0,timer.reset();i<HERMITE_TEST_ITERATIONS;i++,f+=0.000005f)
{
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
t+=Hermie(f)+Hermie(f+0.1f);
}
t=timer();
printf("hermite<double>:time=%f milliseconds\n",t*1000);
return ret;
}
int hermite_fixed_test(void)
{
int ret=0;
int i;
hermite<fixed> Hermie;
hermite<fixed>::time_type f;
hermite<fixed>::time_type inc(0.0005f), inc2(1.10);
fixed sum(0);
etl::clock timer;
double t;
Hermie.p1()=0;
Hermie.t1()=1;
Hermie.p2()=0;
Hermie.t2()=1;
Hermie.sync();
{fixed t;
for(i=0,f=0,timer.reset();i<HERMITE_TEST_ITERATIONS;i++,f+=inc)
{
sum+=Hermie(f)+Hermie(f+inc2);
sum+=Hermie(f)+Hermie(f+inc2);
sum+=Hermie(f)+Hermie(f+inc2);
sum+=Hermie(f)+Hermie(f+inc2);
sum+=Hermie(f)+Hermie(f+inc2);
sum+=Hermie(f)+Hermie(f+inc2);
sum+=Hermie(f)+Hermie(f+inc2);
sum+=Hermie(f)+Hermie(f+inc2);
sum+=Hermie(f)+Hermie(f+inc2);
sum+=Hermie(f)+Hermie(f+inc2);
sum+=Hermie(f)+Hermie(f+inc2);
sum+=Hermie(f)+Hermie(f+inc2);
}
}
t=timer();
printf("hermite<fixed>:time=%f milliseconds\n",t*1000);
return ret;
}
int hermite_angle_test(void)
{
int ret=0,i;
float f;
hermite<angle> Hermie;
etl::clock timer;
angle tmp;
double t;
Hermie.p1()=angle::degrees(0);
Hermie.t1()=angle::degrees(45);
Hermie.p2()=angle::degrees(-45);
Hermie.t2()=angle::degrees(180);
Hermie.sync();
for(f=0.0f,i=0,timer.reset();i<HERMITE_TEST_ITERATIONS;i++,f+=0.000005f)
{
tmp+=Hermie(f)+Hermie(f+0.1f);
tmp+=Hermie(f)+Hermie(f+0.1f);
tmp+=Hermie(f)+Hermie(f+0.1f);
tmp+=Hermie(f)+Hermie(f+0.1f);
tmp+=Hermie(f)+Hermie(f+0.1f);
tmp+=Hermie(f)+Hermie(f+0.1f);
tmp+=Hermie(f)+Hermie(f+0.1f);
tmp+=Hermie(f)+Hermie(f+0.1f);
tmp+=Hermie(f)+Hermie(f+0.1f);
tmp+=Hermie(f)+Hermie(f+0.1f);
tmp+=Hermie(f)+Hermie(f+0.1f);
tmp+=Hermie(f)+Hermie(f+0.1f);
}
t=timer();
fprintf(stderr,"hermite<angle>:time=%f milliseconds\n",t*1000);
return ret;
}
/* === E N T R Y P O I N T ================================================= */
int main()
{
int error=0;
error+=surface_and_gaussian_blur_test();
error+=hermite_float_test();
error+=hermite_double_test();
error+=hermite_int_test();
error+=hermite_fixed_test();
error+=hermite_angle_test();
return error;
}