/*! ========================================================================

** 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=""></etl>

#include <etl hermite=""></etl>

#include <etl angle=""></etl>

#include <etl fixed=""></etl>

#include <etl surface=""></etl>

#include <etl gaussian=""></etl>

#include <etl calculus=""></etl>

#include <stdio.h></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=""></class>

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=""></class>

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=""></class>

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=""></class>

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);</float>

		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);</float>

	}

	{

		surface<double> my_surface(1000,1000);</double>

		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);</double>

	}

	{

		surface<fixed> my_surface(1000,1000);</fixed>

		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);</fixed>

	}

	return ret;

}

int hermite_int_test()

{

	int ret=0;

	hermite<int> Hermie;</int>

	hermite<int>::time_type f;</int>

	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

	{

		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);</int>

	return ret;

}

int hermite_float_test(void)

{

	int ret=0;

	float f; int i;

	hermite<float> Hermie;</float>

	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

	{

		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);</float>

	return ret;

}

int hermite_double_test(void)

{

	int ret=0,i;

	float f;

	hermite<double> Hermie;</double>

	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

	{

		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);</double>

	return ret;

}

int hermite_fixed_test(void)

{

	int ret=0;

    int i;

	hermite<fixed> Hermie;</fixed>

	hermite<fixed>::time_type f;</fixed>

	hermite<fixed>::time_type inc(0.0005f), inc2(1.10);</fixed>

	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

	{

		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);</fixed>

	return ret;

}

int hermite_angle_test(void)

{

	int ret=0,i;

	float f;

	hermite<angle> Hermie;</angle>

	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

	{

		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);</angle>

	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;

}