EI/VirtualLeaf Changeset - e730ce1f2d31 · Centrum Wiskunde & Informatica (CWI)

Changeset - e730ce1f2d31

Parent rev.

Child rev.

[Not reviewed]

rd fick laws default

0 1 0

Michael Guravage - 14 years ago 2011-08-26 17:04:13
michael.guravage@cwi.nl

Added second reaction-diffusion protocol: Fick's laws of diffusion.

--
user: Michael Guravage <michael.guravage@cwi.nl>
branch 'default'
changed src/protocols/MyDiffusionModel/mydiffusionmodel.cpp

1 file changed with 9 insertions and 1 deletions:

src/protocols/MyDiffusionModel/mydiffusionmodel.cpp

0 comments (0 inline, 0 general)

src/protocols/MyDiffusionModel/mydiffusionmodel.cpp

➞

Show inline comments

@@ @@ -22,25 +22,25 @@ @@
 #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 QString( "Diffusion dependent coloriing" );
+  return QString( "Fick's Laws of diffusion" );
+}
 // 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) {
@@ @@ -56,24 +56,32 @@ void Mydiffusionmodel::CellHouseKeeping( @@
   if (c->Chemical(0) < 0.5){
     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;
     double phi = w->Length() * ( par->D[c] ) * ( w->C2()->Chemical(c) - w->C1()->Chemical(c) );
     dchem_c1[c] += phi;
     dchem_c2[c] -= phi;
+  }
+}
 void Mydiffusionmodel::WallDynamics(Wall *w, double *dw1, double *dw2) {
   // add biochemical networks for reactions occuring at walls here
+}
 void Mydiffusionmodel::CellDynamics(CellBase *c, double *dchem) {
   // add biochemical networks for intracellular reactions here
   double Y = c->Chemical(0),
          A = c->Chemical(1),
          H = c->Chemical(2),
          S = c->Chemical(3);

0 comments (0 inline, 0 general)