Maleckar2009OdeSystem.hpp

Go to the documentation of this file.

#ifndef _MALECKAR2009ODESYSTEM_
#define _MALECKAR2009ODESYSTEM_


#include "ChasteSerialization.hpp"
#include <boost/serialization/base_object.hpp>
#include "AbstractCardiacCell.hpp"
#include "AbstractStimulusFunction.hpp"
#include "AbstractIvpOdeSolver.hpp"

class Maleckar2009OdeSystem : public AbstractCardiacCell
{

private:
    double mScaleFactorGks;
    double mScaleFactorIto;
    double mScaleFactorGkr;
    double mScaleFactorGna;
    double mScaleFactorAch;
    double mScaleFactorGNaK;
    double mScaleFactorGNaCa;
    double mScaleFactorGCaL;
    double mScaleFactorGKur;
    double mScaleFactorGK1;
    double mScaleFactorAZD;

    void VerifyStateVariables();

    friend class boost::serialization::access;

    template<class Archive>
    void serialize(Archive & archive, const unsigned int version)
    {
        archive & boost::serialization::base_object<AbstractCardiacCell>(*this);
    }

public:

    Maleckar2009OdeSystem(boost::shared_ptr<AbstractIvpOdeSolver> pSolver,
                 boost::shared_ptr<AbstractStimulusFunction> pIntracellularStimulus);

    ~Maleckar2009OdeSystem(void);

    void SetScaleFactorGks(double sfgks);

    void SetScaleFactorIto(double sfito);

    void SetScaleFactorGkr(double sfgkr);

    void SetScaleFactorGna(double sfgna);

    void SetScaleFactorAch(double sfach);

    void SetScaleFactorGNaK(double sfgnak);

    void SetScaleFactorGNaCa(double sfgnaca);

    void SetScaleFactorGCaL(double sfgcal);

    void SetScaleFactorGKur(double sfgkur);

    void SetScaleFactorGK1(double sfgk1);

    void SetScaleFactorAZD(double sfazd);

    double GetIIonic();

    void EvaluateYDerivatives(
            double var_environment__time,
            const std::vector<double> &rY,
            std::vector<double> &rDY);

};

// Needs to be included last
#include "SerializationExportWrapper.hpp"
CHASTE_CLASS_EXPORT(Maleckar2009OdeSystem)

namespace boost
{
    namespace serialization
    {
        template<class Archive>
        inline void save_construct_data(
            Archive & ar, const Maleckar2009OdeSystem * t, const unsigned int fileVersion)
        {
            const boost::shared_ptr<AbstractIvpOdeSolver> p_solver = t->GetSolver();
            const boost::shared_ptr<AbstractStimulusFunction> p_stimulus = t->GetStimulusFunction();
            ar << p_solver;
            ar << p_stimulus;
        }

        template<class Archive>
        inline void load_construct_data(
            Archive & ar, Maleckar2009OdeSystem * t, const unsigned int fileVersion)
        {
            boost::shared_ptr<AbstractIvpOdeSolver> p_solver;
            boost::shared_ptr<AbstractStimulusFunction> p_stimulus;
            ar >> p_solver;
            ar >> p_stimulus;
            ::new(t)Maleckar2009OdeSystem(p_solver, p_stimulus);
        }

    }

}

#endif