---

Cell.cpp

Go to the documentation of this file.

// Cell.cpp: implementation of the CCell class.
//

#include <iostream.h>           //***this is temporary
#include "Components.h"
#include "Display.h"
#include "DNA.h"
#include "Geometry.h"
#include "Cell.h"
#include <math.h>

// Construction/Destruction

//Constructor
CCell::CCell()
{
        m_pComponents=new CComponents(this);    //Create the components
        m_pDisplay=new CDisplay(this);                  //Create the display
        m_pDNA=new CDNA(this);                                  //Create the DNA
        m_pGeometry = new CGeometry(this);              //Create the geometry
        GenTime=0;
        TotalTime=0;
        NGEN=0;
}

//Destructor
CCell::~CCell()
{
        //Free up memory allocated in the constructor
        delete m_pComponents;
        delete m_pDisplay;                      
        delete m_pDNA;                          
        delete m_pGeometry;
}

//Copy constructor
CCell::CCell(const CCell& rCell)        
{
        m_pComponents=new CComponents(*rCell.m_pComponents,this);       //Create and copy the components
        m_pDisplay=new CDisplay(*rCell.m_pDisplay,this);                        //Create and copy the display
        m_pDNA=new CDNA(*rCell.m_pDNA,this);                                            //Create and copy DNA
        m_pGeometry=new CGeometry(*rCell.m_pGeometry,this);             //Create and copy the geometry
}

//Assignment operator (safe but imperfect)
CCell& CCell::operator =(const CCell & rhs)
{
        if(this == &rhs)                //check addresses, if equal
                return *this;           //return the 1st operand

        delete m_pComponents;
        delete m_pDisplay;              
        delete m_pDNA;          
        delete m_pGeometry;             //delete any existing values
        
        m_pComponents=new CComponents(*rhs.m_pComponents,this); //Create and copy the components
        m_pDisplay=new CDisplay(*rhs.m_pDisplay,this);                  //Create and copy the display
        m_pDNA=new CDNA(*rhs.m_pDNA,this);                                              //Create and copy DNA
        m_pGeometry=new CGeometry(*rhs.m_pGeometry,this);               //Create and copy the geometry
        
        //***DON'T USE YET, m_pCell POINTERS STILL POINT TO SAME PARENT***
        return *this;
}

void CCell::Run()
{
        int NGENMAX=30;
        NGEN = 0;               //Number of generations
        m_pGeometry->CheckVolume();             //Make sure volume and total mass agree
        
        for(;NGEN<=NGENMAX;)                            //Loop for generations
        {
                m_pComponents->StartAverages();
                Generation();
        }
}

void CCell::Division()
{
        ++NGEN;
        m_pComponents->Divide();
        m_pDNA->Divide();
        m_pGeometry->Divide();
}

CGeometry* CCell::GetGeometry()
{
        return m_pGeometry;                             //return pointer to the geometry
}
CComponents* CCell::GetComponents()
{
        return m_pComponents;                   //return pointer to the components
}

CDisplay* CCell::GetDisplay()
{
        return m_pDisplay;                              //return pointer to the display
}

void CCell::Generation()                        //runs one complete generation to division
{
        static double PrintTime=0;
        GenTime=0;                                              //reset GenTime;                                        
        m_pComponents->FirstInt();              //run only on first generation
                                
        do{                                                             //loop until division
                m_pComponents->Integrate();     //integrate one time step
                m_pDNA->AdvanceFork();          //DNA replication
                m_pDNA->EcoInit();
                m_pDNA->Termination();          //does termination occur?
                m_pGeometry->UpdateShape();     //envelope elongation and cross-wall synthesis
                if(TotalTime>PrintTime)
                {
                //      m_pDisplay->PrintComponents();  //Print
                //      PrintTime+=0.02;
                }
        }while(!DivisionCheck());       
        Division();
        m_pDisplay->PrintComponents();          //Print on division
}

bool CCell::DivisionCheck()
{
        double SEP=3.14159*pow(GetGeometry()->GetWidth(),2);                    //Area of hemispherical ends
        //if cell doesn't divide
        if((GetGeometry()->GetSepf()-SEP)<=0.0)         //if pinching off doesn't occur
        {
                GetComponents()->NextY();
                return 0;
        }
        else
                return 1;
}

//DEL double CCell::GetTime()
//DEL {
//DEL   return GenTime;
//DEL }

void CCell::AdvanceTime(double DT)
{
        GenTime+=DT;
        TotalTime+=DT;
}

CDNA* CCell::GetDNA()
{
        return m_pDNA;
}

---