#include <lib-algorithms.h>

Public Member Functions
	LEAPs (T &Yin)

T &	operator() (T &Y, double h, AbstFunctor< T > F_space, AbstFunctor< T > F_momentum)

T &	operator() (T &Y, double h, AbstFunctor< T > F_space, AbstFunctor< T > F_momentum, size_t dir)

Private Attributes
T	Y0

T	Yh

Detailed Description

template<class T>
class Algorithms::LEAPs< T >

Definition at line 703 of file lib-algorithms.h.

Constructor & Destructor Documentation

◆ LEAPs()

template<class T >

Algorithms::LEAPs< T >::LEAPs ( T & Yin )

inline

Definition at line 706 of file lib-algorithms.h.

706 : Y0(Yin), Yh(Yin) { }

Algorithms::LEAPs::Y0

T Y0

Definition: lib-algorithms.h:714

Algorithms::LEAPs::Yh

T Yh

Definition: lib-algorithms.h:714

Member Function Documentation

◆ operator()() [1/2]

template<class T >

T & Algorithms::LEAPs< T >::operator()	(	T &	Y,
		double	h,
		AbstFunctor< T > *	F_space,
		AbstFunctor< T > *	F_momentum
	)

Definition at line 718 of file lib-algorithms.h.

                                                                                   {
 //      Take a step using LEAPspace
 
 //      Initialization
         Y0 = Y;
 
         (*F_space)(Y0,Yh);      Yh *= 0.5*h;            //  x  = h * F_space(Y0)
         Y0  += Yh;
         (*F_momentum)(Y0,Yh);   Yh *= h;            //  p  = h * F_momentum(Y0)
         Y0 += Yh;                                   //  Y0 = Y0 + h*Yh
         (*F_space)(Y0,Yh);      Yh *= 0.5*h;            //  x  = h * F_space(Y0)
         Y0  += Yh;
         
         Y = Y0;
 //      ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
         return Y;
     }

◆ operator()() [2/2]

template<class T >

T & Algorithms::LEAPs< T >::operator()	(	T &	Y,
		double	h,
		AbstFunctor< T > *	F_space,
		AbstFunctor< T > *	F_momentum,
		size_t	dir
	)

Definition at line 737 of file lib-algorithms.h.

                                                                                               {
 //      Take a step using LEAPspace
 
 //      Initialization
         Y0 = Y;
 
 //      First step momentum
         (*F_space)(Y0,Yh);      Yh *= 0.5*h;            //  x  = h * F_space(Y0)
         Y0  += Yh;
         (*F_momentum)(Y0,Yh);   Yh *= h;            //  p  = h * F_momentum(Y0)
         Y0 += Yh;                                   //  Y0 = Y0 + h*Yh
         (*F_space)(Y0,Yh);      Yh *= 0.5*h;            //  x  = h * F_space(Y0)
         Y0  += Yh;
 
         Y = Y0;
 //      ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
         return Y;
     }