rk.h

//
// Created by milinda on 1/29/17.
//

#ifndef SFCSORTBENCH_RK_H
#define SFCSORTBENCH_RK_H

#include "mesh.h"

namespace ode {
    namespace solver {

        class RK {

        protected :

            ot::Mesh* m_uiMesh;

            unsigned int m_uiOrder;

            unsigned int m_uiNrp;

            MPI_Comm m_uiComm;

            double m_uiTimeBegin;

            double m_uiTimeEnd;

            double m_uiT_h;

            double m_uiCurrentTime;

            double m_uiT_h_prev;

            double m_uiTimeStepCurrent;

            unsigned int m_uiCurrentStep;

            const double RK45_STAGE_1_COEF=16.0/135.0;

            const double RK45_STAGE_1_T=1.0;
            const double RK45_STAGE_1_U=1.0;

            const double RK45_STAGE_2_COEF=0.0;

            const double RK45_STAGE_2_T=(1.0/4.0);
            const double RK45_STAGE_2_U=(1.0/4.0);

            const double RK45_STAGE_3_COEF=6656.0/12825.0;

            const double RK45_STAGE_3_T=(3.0/8.0);
            const double RK45_STAGE_3_U1=(3.0/32.0);
            const double RK45_STAGE_3_U2=(9.0/32.0);


            const double RK45_STAGE_4_COEF=28561.0/56430;
            const double RK45_STAGE_4_T=(12.0/13.0);
            const double RK45_STAGE_4_U1=(1932.0/2197.0);
            const double RK45_STAGE_4_U2=(-7200.0/2197.0);
            const double RK45_STAGE_4_U3=(7296.0/2197.0);

            const double RK45_STAGE_5_COEF=-9.0/50.0;
            const double RK45_STAGE_5_T=1.0;
            const double RK45_STAGE_5_U1=(439.0/216.0);
            const double RK45_STAGE_5_U2=(-8.0);
            const double RK45_STAGE_5_U3=(3680.0/513.0);
            const double RK45_STAGE_5_U4=(-845.0/4104.0);

            const double RK45_STAGE_6_COEF=2.0/55.0;
            const double RK45_STAGE_6_T=0.5;
            const double RK45_STAGE_6_U1=(-8.0/27.0);
            const double RK45_STAGE_6_U2=(2.0);
            const double RK45_STAGE_6_U3=(-3544.0/2565.0);
            const double RK45_STAGE_6_U4=(1859.0/4104.0);
            const double RK45_STAGE_6_U5=(-11.0/40.0);

            std::vector<unsigned int > m_uiRefinedOctIDs;

            std::vector<unsigned int> m_uiCoarsenOctIDs;



        public:
            RK(ot::Mesh * pMesh,double pTBegin,double pTEnd,double pTh);

            ~RK();
/*
            inline void setTimeStepBegin(double tBegin){m_uiTimeBegin=tBegin;}
            *//*
            inline void setTimeStepEnd(double tEnd){m_uiTimeEnd=tEnd;}
            *//*
            inline void setCurrentTimeStep(unsigned int step) {m_uiCurrentStep=step;}
            *//*
            inline void setCurrentTimeValue(double tValue){m_uiCurrentTime=tValue;}
            *//*
            inline void setDT(double dt){m_uiT_h=dt;}*/



            virtual void performSingleIteration(){};

            virtual  void applyInitialConditions(){};

            virtual void applyBoundaryConditions(){};

            virtual void rkSolve(){};
            
            inline void freeMesh(){delete m_uiMesh;};




        };



    } // end of namespace solver

} // end of namespace ode.

#endif //SFCSORTBENCH_RK_H