/*****************************************************************************

        TestSpeed.hpp
        Copyright (c) 2005 Laurent de Soras

--- Legal stuff ---

This library is free software; you can redistribute it and/or
modify it under the terms of the GNU Lesser General Public
License as published by the Free Software Foundation; either
version 2.1 of the License, or (at your option) any later version.

This library 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
Lesser General Public License for more details.

You should have received a copy of the GNU Lesser General Public
License along with this library; if not, write to the Free Software
Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

*Tab=3***********************************************************************/



#if defined (TestSpeed_CURRENT_CODEHEADER)
	#error Recursive inclusion of TestSpeed code header.
#endif
#define	TestSpeed_CURRENT_CODEHEADER

#if ! defined (TestSpeed_CODEHEADER_INCLUDED)
#define	TestSpeed_CODEHEADER_INCLUDED



/*\\\ INCLUDE FILES \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/

#include	"test_fnc.h"
#include	"stopwatch/StopWatch.h"
#include	"TestWhiteNoiseGen.h"

#include	<typeinfo>

#include	<cstdio>



/*\\\ PUBLIC \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/



template <class FO>
int	TestSpeed <FO>::perform_test_single_object (FO &fft)
{
	assert (&fft != 0);

   int            ret_val = 0;

	const std::type_info &	ti = typeid (fft);
	const char *	class_name_0 = ti.name ();

   if (ret_val == 0)
   {
	   perform_test_d (fft, class_name_0);
   }
   if (ret_val == 0)
   {
	   perform_test_i (fft, class_name_0);
   }
	if (ret_val == 0)
   {
      perform_test_di (fft, class_name_0);
   }

   if (ret_val == 0)
   {
      printf ("\n");
   }

   return (ret_val);
}



template <class FO>
int	TestSpeed <FO>::perform_test_d (FO &fft, const char *class_name_0)
{
	assert (&fft != 0);
   assert (class_name_0 != 0);

	const long		len = fft.get_length ();
   const long     nbr_tests = limit (
      static_cast <long> (NBR_SPD_TESTS / len / len),
      1L,
      static_cast <long> (MAX_NBR_TESTS)
   );

	TestWhiteNoiseGen <DataType>	noise;
	std::vector <DataType>	x (len, 0);
	std::vector <DataType>	s (len);
	noise.generate (&x [0], len);

   printf (
		"%s::do_fft () speed test [%ld samples]... ",
		class_name_0,
		len
	);
	fflush (stdout);

	stopwatch::StopWatch	chrono;
	chrono.start ();
	for (long test = 0; test < nbr_tests; ++ test)
	{
		fft.do_fft (&s [0], &x [0]);
		chrono.stop_lap ();
	}

	printf ("%.1f clocks/sample\n", chrono.get_time_best_lap (len));

	return (0);
}



template <class FO>
int	TestSpeed <FO>::perform_test_i (FO &fft, const char *class_name_0)
{
	assert (&fft != 0);
   assert (class_name_0 != 0);

	const long		len = fft.get_length ();
   const long     nbr_tests = limit (
      static_cast <long> (NBR_SPD_TESTS / len / len),
      1L,
      static_cast <long> (MAX_NBR_TESTS)
   );

	TestWhiteNoiseGen <DataType>	noise;
	std::vector <DataType>	x (len);
	std::vector <DataType>	s (len, 0);
	noise.generate (&s [0], len);

   printf (
		"%s::do_ifft () speed test [%ld samples]... ",
		class_name_0,
		len
	);
	fflush (stdout);

	stopwatch::StopWatch	chrono;
	chrono.start ();
	for (long test = 0; test < nbr_tests; ++ test)
	{
		fft.do_ifft (&s [0], &x [0]);
		chrono.stop_lap ();
	}

	printf ("%.1f clocks/sample\n", chrono.get_time_best_lap (len));

	return (0);
}



template <class FO>
int	TestSpeed <FO>::perform_test_di (FO &fft, const char *class_name_0)
{
	assert (&fft != 0);
   assert (class_name_0 != 0);

	const long		len = fft.get_length ();
   const long     nbr_tests = limit (
      static_cast <long> (NBR_SPD_TESTS / len / len),
      1L,
      static_cast <long> (MAX_NBR_TESTS)
   );

	TestWhiteNoiseGen <DataType>	noise;
	std::vector <DataType>	x (len, 0);
	std::vector <DataType>	s (len);
	std::vector <DataType>	y (len);
	noise.generate (&x [0], len);

   printf (
		"%s::do_fft () / do_ifft () / rescale () speed test [%ld samples]... ",
		class_name_0,
		len
	);
	fflush (stdout);

	stopwatch::StopWatch	chrono;

	chrono.start ();
	for (long test = 0; test < nbr_tests; ++ test)
	{
		fft.do_fft (&s [0], &x [0]);
		fft.do_ifft (&s [0], &y [0]);
		fft.rescale (&y [0]);
		chrono.stop_lap ();
	}

	printf ("%.1f clocks/sample\n", chrono.get_time_best_lap (len));

	return (0);
}



/*\\\ PROTECTED \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/



/*\\\ PRIVATE \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/



#endif	// TestSpeed_CODEHEADER_INCLUDED

#undef TestSpeed_CURRENT_CODEHEADER



/*\\\ EOF \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/