BackwardEulerFoxModel2002Modified.hpp

Go to the documentation of this file.

#ifndef _BackwardEulerFoxModel2002Modified_
#define _BackwardEulerFoxModel2002Modified_

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

class BackwardEulerFoxModel2002Modified : public AbstractBackwardEulerCardiacCell<3>
{
private:
    friend class boost::serialization::access;
    template<class Archive>
    void serialize(Archive & archive, const unsigned int version)
    {
        // This calls serialize on the base class.
        archive & boost::serialization::base_object<AbstractBackwardEulerCardiacCell<3> >(*this);
    }
public:
    BackwardEulerFoxModel2002Modified(boost::shared_ptr<AbstractStimulusFunction> pIntracellularStimulus);

    BackwardEulerFoxModel2002Modified(boost::shared_ptr<AbstractIvpOdeSolver> /* unused */,
                                      boost::shared_ptr<AbstractStimulusFunction> pIntracellularStimulus);


    ~BackwardEulerFoxModel2002Modified(void)
    {
    }

    double GetIIonic();

    void ComputeResidual(double var_environment__time, const double rCurrentGuess[3], double rResidual[3]);

    void ComputeJacobian(double var_environment__time, const double rCurrentGuess[3], double rJacobian[3][3]);

protected:
    void UpdateTransmembranePotential(double var_environment__time);

    void ComputeOneStepExceptVoltage(double var_environment__time);
};

// Must be last
#include "SerializationExportWrapper.hpp"
CHASTE_CLASS_EXPORT(BackwardEulerFoxModel2002Modified)

namespace boost
{
namespace serialization
{
template<class Archive>
inline void save_construct_data(
    Archive & ar, const BackwardEulerFoxModel2002Modified * t, const unsigned int file_version)
{
    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, BackwardEulerFoxModel2002Modified * t, const unsigned int file_version)
{
    boost::shared_ptr<AbstractIvpOdeSolver> p_solver;
    boost::shared_ptr<AbstractStimulusFunction> p_stimulus;
    ar >> p_solver;
    ar >> p_stimulus;
    ::new(t)BackwardEulerFoxModel2002Modified(p_solver, p_stimulus);
}
}
} // namespace ...

#endif