testTabulatedFunctionComplexLgAmpPhase.cc

Go to the documentation of this file.

#include <iostream>

#include <tst/Verify.h>

#include <utl/ComplexLgAmpPhase.h>
#include <utl/TabulatedFunctionComplexLgAmpPhase.h>
#include <utl/AugerException.h>
#include <utl/AugerUnits.h>

#include <cppunit/extensions/HelperMacros.h>

using namespace std;
using namespace utl;
using namespace tst;

class testTabulatedFunctionComplexLgAmpPhase : public CppUnit::TestFixture {

  CPPUNIT_TEST_SUITE(testTabulatedFunctionComplexLgAmpPhase);
  CPPUNIT_TEST(testInterpolation);
  CPPUNIT_TEST(testMultiplication);
  CPPUNIT_TEST_SUITE_END();

private:

public:
  void
  setUp()
  {
  }

  void
  tearDown()
  {
  }

  void
  testInterpolation()
  {
    // configure the tests
    const unsigned int numpoints = 10;
    const unsigned int numcheckpoints = 100;
    const double phasestepsize = 50.*degree;
    
    // create our working TabulatedFunctionComplexLgAmpPhase
    TabulatedFunctionComplexLgAmpPhase tfc1;
    
    // fill in some test data
    for (unsigned int i=0; i<=numpoints; ++i)
    tfc1.PushBack(double(i),ComplexLgAmpPhase(double(numpoints*i),double(phasestepsize*i)));
    
    // test number of entries
    CPPUNIT_ASSERT(Verify<Equal>(tfc1.GetNPoints(),numpoints+1));

    // now test if interpolated values are ok
    for (unsigned int i=0; i<numcheckpoints; ++i) {
      double point = double(i*numpoints/double(numcheckpoints));
      CPPUNIT_ASSERT(Verify<CloseTo>(tfc1.Y(point).GetLgAmplitude(), double(numpoints*point))); // amplitude
      CPPUNIT_ASSERT(Verify<CloseTo>(tfc1.Y(point).GetPhase(), double(phasestepsize*point))); // phase
    }  
  }

  void
  testMultiplication()
  {
    // configure the tests
    const unsigned int numpoints = 10;
    const unsigned int numcheckpoints = 100;
    const double phasestepsize = 50.*degree;
    const double offset = 2.0;
    const unsigned int densityfactor = 2;
    
    // create our working TabulatedFunctionComplexLgAmpPhasees
    TabulatedFunctionComplexLgAmpPhase tfc1;
    TabulatedFunctionComplexLgAmpPhase tfc2;
    
    // fill in some test data
    for (unsigned int i=0; i<=numpoints; ++i)
      tfc1.PushBack(double(i),ComplexLgAmpPhase(double(numpoints*i),double(phasestepsize*i)));
    for (unsigned int i=0; i<=numpoints*densityfactor; ++i)
      tfc2.PushBack(double(offset+i/double(densityfactor)),ComplexLgAmpPhase(double(numpoints*i),double(phasestepsize*i)));
    
    TabulatedFunctionComplexLgAmpPhase tfc3=tfc1;
    tfc3*=tfc2;
    
    TabulatedFunctionComplexLgAmpPhase tfc4=tfc2;
    tfc4*=tfc1;
    
    // screen output for debugging purposes
    cout << "TFC1:\n";
    for (TabulatedFunctionComplexLgAmpPhase::Iterator it = tfc1.Begin(); it != tfc1.End(); ++it)
      cout << it->X() << "\t" << it->Y().GetLgAmplitude() << "\t" << it->Y().GetPhase()/degree << endl;
    cout << "\n" << endl;
    
    cout << "TFC2:\n";
    for (TabulatedFunctionComplexLgAmpPhase::Iterator it = tfc2.Begin(); it != tfc2.End(); ++it)
      cout << it->X() << "\t" << it->Y().GetLgAmplitude() << "\t" << it->Y().GetPhase()/degree << endl;
    cout << "\n" << endl;

    cout << "TFC3:\n";
    for (TabulatedFunctionComplexLgAmpPhase::Iterator it = tfc3.Begin(); it != tfc3.End(); ++it)
      cout << it->X() << "\t" << it->Y().GetLgAmplitude() << "\t" << it->Y().GetPhase()/degree << endl;
    cout << "\n" << endl;

    cout << "TFC4:\n";
    for (TabulatedFunctionComplexLgAmpPhase::Iterator it = tfc4.Begin(); it != tfc4.End(); ++it)
      cout << it->X() << "\t" << it->Y().GetLgAmplitude() << "\t" << it->Y().GetPhase()/degree << endl;
    cout << "\n" << endl;

    // now test if interpolated values are ok
    cout << "TFC3 versus manual multiplication:" << endl;
    for (unsigned int i=0; i<numcheckpoints; ++i) {
      double point = double(offset+i*(numpoints-offset)/double(numcheckpoints));
      CPPUNIT_ASSERT(Verify<CloseTo>(tfc3.Y(point).GetLgAmplitude(),tfc4.Y(point).GetLgAmplitude()));  // compare amplitudes of the two multiplied TFCs
      CPPUNIT_ASSERT(Verify<CloseTo>(tfc3.Y(point).GetPhase(),tfc4.Y(point).GetPhase()));          // compare amplitudes of the two multiplied TFCs
      cout << point << "\t" << tfc3.Y(point).GetLgAmplitude() << "\t" << tfc3.Y(point).GetPhase()/degree << "\tversus\t";
      cout << point << "\t" << (tfc1.Y(point)*tfc2.Y(point)).GetLgAmplitude() << "\t" << (tfc1.Y(point)*tfc2.Y(point)).GetPhase()/degree << "\n";
      CPPUNIT_ASSERT(Verify<CloseTo>(tfc3.Y(point).GetLgAmplitude(),(tfc1.Y(point)*tfc2.Y(point)).GetLgAmplitude())); // compare amplitudes of tfc3 with manual multiplication
      CPPUNIT_ASSERT(Verify<CloseTo>(tfc3.Y(point).GetPhase(),(tfc1.Y(point)*tfc2.Y(point)).GetPhase())); // compare amplitudes of tfc3 with manual multiplication
    }  
  }

/*
  void
  testBinException()
  {
  }
*/

};

CPPUNIT_TEST_SUITE_REGISTRATION(testTabulatedFunctionComplexLgAmpPhase);

// Configure (x)emacs for this file ...
// Local Variables:
// mode: c++
// End: