00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028 #ifndef TYSONNOVAKCELLCYCLEMODEL_HPP_
00029 #define TYSONNOVAKCELLCYCLEMODEL_HPP_
00030
00031 #include <boost/serialization/access.hpp>
00032
00033 #include <iostream>
00034
00035 #include "AbstractOdeBasedCellCycleModel.hpp"
00036 #include "TysonNovak2001OdeSystem.hpp"
00037 #include "BackwardEulerIvpOdeSolver.hpp"
00038 #include "CvodeAdaptor.hpp"
00039 #include "Exception.hpp"
00040
00041
00042 #include <boost/serialization/export.hpp>
00043
00056 class TysonNovakCellCycleModel : public AbstractOdeBasedCellCycleModel
00057 {
00058 private:
00059
00060 #ifdef CHASTE_CVODE
00061
00062 static CvodeAdaptor msSolver;
00063 #else
00064
00065 static BackwardEulerIvpOdeSolver msSolver;
00066 #endif //CHASTE_CVODE
00067
00074 TysonNovakCellCycleModel(std::vector<double> parentProteinConcentrations,
00075 double divideTime);
00076
00078 friend class boost::serialization::access;
00085 template<class Archive>
00086 void serialize(Archive & archive, const unsigned int version)
00087 {
00088 archive & boost::serialization::base_object<AbstractOdeBasedCellCycleModel>(*this);
00089 }
00090
00091 public:
00092
00096 TysonNovakCellCycleModel();
00097
00102 void ResetForDivision();
00103
00113 AbstractCellCycleModel* CreateDaughterCellCycleModel();
00114
00121 bool SolveOdeToTime(double currentTime);
00122
00128 double GetOdeStopTime();
00129
00133 double GetSDuration();
00134
00138 double GetG2Duration();
00139
00143 double GetMDuration();
00144
00148 void InitialiseDaughterCell();
00149
00150 };
00151
00152
00153
00154 BOOST_CLASS_EXPORT(TysonNovakCellCycleModel)
00155
00156
00157 #endif