RDE_lib2/SolutionPath.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 __SOLUTIONPATH__
#define __SOLUTIONPATH__


#include "Path.h"
#include "GroupElement.h"
#include "LieMatrix.h"
#include "dyadic_interval.h"
#include <vector>



template <typename my_alg_type_IN, typename my_alg_type_OUT>
class SolutionPath
{
private:
        Path<my_alg_type_IN> & theControl;

        lietovectorfield< my_alg_type_IN, LieMatrix<my_alg_type_OUT> > l2vf;


public:
        SolutionPath(Path<my_alg_type_IN> & theControlIn, std::vector< LieMatrix<my_alg_type_OUT> > theEquationIn)
                : theControl(theControlIn), l2vf(theEquationIn)
        {
        };

        ~SolutionPath(void)
        {
        };


        GroupElement<my_alg_type_OUT> DescribePath(dyadic_interval increment)
        {
                LieMatrix<my_alg_type_OUT> temp;
                temp = l2vf.lie2vectorfield(theControl.DescribePath(increment, increment.n));
                return temp.exp();
        };

        GroupElement<my_alg_type_OUT> DescribePath(dyadic_interval increment, const int accuracy)
        {
                LieMatrix<my_alg_type_OUT> temp;
                temp = l2vf.lie2vectorfield(theControl.DescribePath(increment, accuracy));
                return temp.exp();
        };

};// END CLASS DEFINITION SolutionPath

#endif // __SOLUTIONPATH__