EulerGeodesic.hpp

//
// $Id: EulerGeodesic.hpp,v 1.1 2004/08/12 16:21:33 werner Exp $
//
// $Log: EulerGeodesic.hpp,v $
// Revision 1.1  2004/08/12 16:21:33  werner
// Integration of numerical geodesics now compiles. Working is not yet satisfying.
//
// Revision 1.8  2004/05/12 14:36:22  werner
// Separation of chart definitions on a manifold and the tangential space.
// Introduced convenient coordinate transformations.
//
// Revision 1.7  2004/05/11 18:05:10  werner
//
// Revision 1.6  2004/05/11 16:53:47  werner
// Introduction of coordinate systems for improved type-safety.
// Will support easy coordinate transformations soon.
//
// Revision 1.5  2004/05/06 22:42:16  werner
// Towards a specification of a spacetime via the Acceleration structure.
//
// Revision 1.4  2004/05/05 15:56:52  werner
// Separation of DOP core routines with dynamic size into the ODE library.
//
// Revision 1.3  2004/05/03 13:33:33  werner
// integration improved
//
// Revision 1.2  2004/03/29 11:51:02  werner
// Common interface among simple integrators, DiffMe and Vecal, and preliminiary work on integrating dop853.
//
// Revision 1.1  2004/03/22 11:55:02  werner
// Schwarzschild geodesic integration.
//
// Revision 1.1  2004/02/13 16:36:21  werner
// Initial preliminiary version of the Vector Algebra Library.
//
#ifndef __EulerGeodesic_HPP
#define __EulerGeodesic_HPP "Created 27.02.2001 21:42:27 by werner"

#include <vecal/Matrix.hpp>
#include <ode/Integrator.hpp>

namespace Traum
{

template <class Acceleration>
struct  EulerGeodesic : public IntegratorBase
{
        typedef typename Acceleration::Scalar_t      Scalar_t;
        typedef typename Acceleration::Point_t       Point_t;
        typedef typename Acceleration::Vector_t      Vector_t;
        typedef typename Acceleration::Christoffel_t Christoffel_t;

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

        Point_t  x;
        Vector_t v;

        const Acceleration&Accel;
        Scalar_t ds;

        EulerGeodesic(const Acceleration&A)
        : Accel(A), ds(1)
        {}

        success_code advance(bool backward=false)
        {
                if (backward)
                {
                        x -= v * ds;
                        v += Accel(x, v) * ds;
                }
                else
                {
                        x += v * ds;
                        v -= Accel(x, v) * ds;
                }
                return StepOk;
        }
};


} /* namespace VecAl */ 

#endif /* __EulerGeodesic_HPP */