RDE_lib2/NonLinearGroupElement.h

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

Copyright 2010 Terry Lyons, Stephen Buckley, Djalil Chafai, 
Greg Gyurkó and Arend Janssen. 

Distributed under the terms of the GNU General Public License, 
Version 3. (See accompanying file License.txt)

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



#pragma once

#ifndef __NONLINEARGROUPELEMENT__
#define __NONLINEARGROUPELEMENT__


#include <map>
#include <vector>
#include "POLYLIE.h"
#include "lietovectorfield.h"
#include "AbstractSolutionPoint.h"


namespace SPB
{

        template <typename my_alg_type_IN, typename my_alg_type_OUT>
        class NonLinearGroupElement
        {
                typedef typename my_alg_type_IN::LIE LIE_IN;

        public:
                NonLinearGroupElement(void)
                {
                };

                NonLinearGroupElement(const LIE_IN & liearg, lietovectorfield< my_alg_type_IN, POLYLIE<my_alg_type_OUT> > & lie2vf)
                {
                        equation = lie2vf.lie2vectorfield(liearg);
                }; 

                ~NonLinearGroupElement(void)
                {
                };


                AbstractSolutionPoint<my_alg_type_OUT> evaluate(AbstractSolutionPoint<my_alg_type_OUT> & InitialPoint)
                {
                        OdeSolver<my_alg_type_OUT> temp;
                        return temp.solve(InitialPoint, equation);
                };

        private:
                POLYLIE<my_alg_type_OUT> equation;
        };
}
#endif // __NONLINEARGROUPELEMENT__