testAxialVector.cc

Go to the documentation of this file.

#include <utl/CoordinateSystem.h>
#include <utl/AxialVector.h>
#include <utl/Vector.h>
#include <utl/TransformationMatrix.h>
#include <utl/AugerUnits.h>
#include <tst/Verify.h>
#include <utl/Triple.h>

#include <cppunit/extensions/HelperMacros.h>

using namespace utl;
using namespace tst;


class testAxialVector : public CppUnit::TestFixture {

  CPPUNIT_TEST_SUITE(testAxialVector);
  CPPUNIT_TEST(testAxialVectorTransformation);
  CPPUNIT_TEST_SUITE_END();

public:
  void setUp() { }

  void tearDown() { }

  void
  testAxialVectorTransformation()
  {
    const CoordinateSystemPtr CS(CoordinateSystem::GetRootCoordinateSystem());
    const AxialVector p(0,0,0, CS);

    const Vector t(1,0,0, CS);
    const CoordinateSystemPtr CTrans(CoordinateSystem::Translation(t, CS));

    p.TransformTo(CTrans);
    CPPUNIT_ASSERT(p.GetCoordinateSystem() == CTrans);
    CPPUNIT_ASSERT(Verify<CloseTo>(
      p.GetCoordinates(CTrans), Triple(0,0,0)));
    CPPUNIT_ASSERT(p.GetCoordinateSystem() == CTrans);
    CPPUNIT_ASSERT(Verify<CloseTo>(
      p.GetCoordinates(CS), Triple(0,0,0)));
    CPPUNIT_ASSERT(p.GetCoordinateSystem() == CS);

    const AxialVector horizontal(1,0,0, CS);
    const CoordinateSystemPtr CRot(CoordinateSystem::RotationZ(90*deg, CS));
    horizontal.TransformTo(CRot);
    CPPUNIT_ASSERT(horizontal.GetCoordinateSystem() == CRot);
    CPPUNIT_ASSERT(Verify<CloseTo>(
      horizontal.GetCoordinates(CRot), Triple(0,-1,0)));
    CPPUNIT_ASSERT(horizontal.GetCoordinateSystem() == CRot);
    CPPUNIT_ASSERT(Verify<CloseTo>(
      horizontal.GetCoordinates(CS), Triple(1,0,0)));
    CPPUNIT_ASSERT(horizontal.GetCoordinateSystem() == CS);
  }

};


CPPUNIT_TEST_SUITE_REGISTRATION(testAxialVector);


// Configure (x)emacs for this file ...
// Local Variables:
// mode:c++
// compile-command: "make -C .. -k run-testGeometry"
// End: