Acceleration.hpp

//
// $Id: Acceleration.hpp,v 1.2 2004/05/11 18:05:10 werner Exp $
//
// $Log: Acceleration.hpp,v $
// Revision 1.2  2004/05/11 18:05:10  werner
// *** empty log message ***
//
// Revision 1.1  2004/05/06 22:42:16  werner
// Towards a specification of a spacetime via the Acceleration structure.
//
//
#ifndef __ACCELERATION_HPP
#define __ACCELERATION_HPP "Created 03.05.2004 21:44:45 by werner benger"

#include "Geodesic.hpp"
#include <ode/Integrate853.hpp>

namespace Traum
{
using namespace VecAl;

template <class SpaceTime> 
        class AffineAcceleration
{
        const SpaceTime&ST;
public:
        typedef typename SpaceTime::Scalar_t     Scalar_t;
        typedef typename SpaceTime::Point_t      Point_t;
        typedef typename SpaceTime::Vector_t     Vector_t;
        typedef typename SpaceTime::Christoffel_t Christoffel_t;

        typedef typename Point_t::Vector_t PointComponents_t;
        enum { Dims = PointComponents_t::SIZE };

        AffineAcceleration(const SpaceTime&R)
        : ST(R)
        {}

        Vector_t operator()(const Point_t&q, const Vector_t&dot_q) const
        {
        Christoffel_t Gamma;
                ST.getChristoffel(Gamma, q);
                return Vector_t( Gamma(dot_q, dot_q) );
        }
};

template <class SpaceTime> 
class TimeAcceleration : public AffineAcceleration<SpaceTime>
{
public:
        typedef typename AffineAcceleration<SpaceTime>::Point_t  Point_t;
        typedef typename AffineAcceleration<SpaceTime>::Vector_t Vector_t;

        Vector_t operator()(const Point_t&q, const Vector_t&dot_q) const
        {
        Vector_t v = AffineAcceleration<SpaceTime>::operator()(q, dot_q);
                return -v + dot_q*v[0];
        }
};


} // namespace Vecal

#endif // __ACCELERATION_HPP