Spectrum Example

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

        TestAccuracy.h
        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 (TestAccuracy_HEADER_INCLUDED)
#define TestAccuracy_HEADER_INCLUDED

#if defined (_MSC_VER)
        #pragma once
        #pragma warning (4 : 4250) // "Inherits via dominance."
#endif

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

template <class FO>
class TestAccuracy
{

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

public:

        typedef typename FO::DataType   DataType;
        typedef long double     BigFloat;       // To get maximum accuracy during intermediate calculations

   static int           perform_test_single_object (FO &fft);
   static int           perform_test_d (FO &fft, const char *class_name_0);
   static int           perform_test_i (FO &fft, const char *class_name_0);
   static int           perform_test_di (FO &fft, const char *class_name_0);

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

protected:

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

private:

        enum {                  NBR_ACC_TESTS   = 10 * 1000 * 1000      };
   enum {         MAX_NBR_TESTS  = 10000  };

   static void          compute_tf (DataType s [], const DataType x [], long length);
        static void             compute_itf (DataType x [], const DataType s [], long length);
        static int              compare_vect_display (const DataType x_ptr [], const DataType y_ptr [], long len, BigFloat &max_err_rel);
        static BigFloat
                                                compute_power (const DataType x_ptr [], long len);
        static BigFloat
                                                compute_power (const DataType x_ptr [], const DataType y_ptr [], long len);
        static void             compare_vect (const DataType x_ptr [], const DataType y_ptr [], BigFloat &power, long &max_err_pos, long len);

/*\\\ FORBIDDEN MEMBER FUNCTIONS \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/

private:

                                                TestAccuracy ();
                                                ~TestAccuracy ();
                                                TestAccuracy (const TestAccuracy &other);
        TestAccuracy &  operator = (const TestAccuracy &other);
        bool                            operator == (const TestAccuracy &other);
        bool                            operator != (const TestAccuracy &other);

};      // class TestAccuracy

#include        "TestAccuracy.hpp"

#endif  // TestAccuracy_HEADER_INCLUDED

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