rk3adaptive.h

#ifndef zh2_rk3_adaptive
#define zh2_rk3_adaptive

#include "integrator.h"

namespace h2comb {

    
    class AdaptiveRungeKutta3Integrator : public Integrator
    {
    public: 
        AdaptiveRungeKutta3Integrator(Integrable *s, double acc) : Integrator(s), accuracy(acc) {}
        void integrate(double T, double deltaTguess);
        
    private:
        void step();
        
        double accuracy;
        
        double t, deltaT;
        
        int dim;
        
        double *state, *k1, *k2, *k3, *intermediateState;
        double infnormIntState();
    };

}

#endif