NDSolver.hh

Go to the documentation of this file.

#pragma once
 
#include "Storage.hh"
#include "Exceptions.hh"
 
namespace cadabra {
 
    
    class NDSolver {
        public:
            NDSolver();
 
            typedef std::vector<double> state_type;
 
            // Integration function.
            void operator()(const state_type &x, state_type &dxdt, const double t);
 
            // Observer function.
            void operator()( const state_type &x , double t );
 
            // Set and check ODE(s).
            void set_ODEs(const Ex&);
            
            // Entry point.
            void integrate();
            
        private:
            // Exception used by the observer function to terminate
            // iteration in case an event happens.
            
            class EventException : public CadabraException {
                public:
                    EventException(std::string="");
            };
            
            std::vector< state_type > states;
            std::vector< double >     times;
 
            Ex ODEs;
    };
    
};