#include <QObject>

#include <QtGui>

#include <QDebug>

#include "simplugin.h"

#include "parameter.h"

#include "wallbase.h"

#include "cellbase.h"

#include "mydiffusionmodel.h"

QString Mydiffusionmodel::ModelID(void) {

  // specify the name of your model here

}

// return the number of chemicals your model uses

int Mydiffusionmodel::NChem(void) { return 4; }

// To be executed after cell division

void Mydiffusionmodel::OnDivide(ParentInfo *parent_info, CellBase *daughter1, CellBase *daughter2) {

  // rules to be executed after cell division go here

  // (e.g., cell differentiation rules)

}

void Mydiffusionmodel::SetCellColor(CellBase *c, QColor *color) { 

  // add cell coloring rules here

  double red=c->Chemical(1)/(1.+c->Chemical(1));

  double green=c->Chemical(0)/(1.+c->Chemical(0));

  double blue=c->Chemical(3)/(1.+c->Chemical(3));

  color->setRgbF(red,green,blue);

}

void Mydiffusionmodel::CellHouseKeeping(CellBase *c) {

  // add cell behavioral rules here

  c->EnlargeTargetArea(par->cell_expansion_rate);

  if(c->Area() > 2*c->BaseArea()){

    c->Divide();

  }

}

void Mydiffusionmodel::CelltoCellTransport(Wall *w, double *dchem_c1, double *dchem_c2) {

  // add biochemical transport rules here

  // Passive fluxes (Fick's law)

  for (int c=0;c<NChem();c++) {

    if (w->C1()->BoundaryPolP() || w->C2()->BoundaryPolP()) return;