cellAdhesion.cpp
b7305f7b5a8budi

cellAdhesion.cpp
View Options

			#include "cellAdhesion.hpp"

			void CellAdhesion::operator() (const state_type &y , state_type &dydt , const double t){
			uint nbcells = lattice-> getNbCells();
			for (uint c=0;c<nbcells;c++){
			Cell& ce = lattice->getCell(c);
			uint sind = ce.stateVectorIndex();
			double no = allnotch[sind];
			// double a = y[sind];
			BoundaryIterator it(&ce);
			double da = 0;//putative complex
			for(Boundary *b=it.first();!it.isLast();b=it.next()){
			if(b->opposite()!=NULL){
			uint ind = b->neighbor()->stateVectorIndex();
			double len = b->scaledLength(lattice->getScale());
			// da += kadnalldelta[ind]noleny[sind];
			da += alldelta[ind]*len;
			}
			}
			dydt[sind] = getAdhesion(day[sind])-gay[sind];
			// dydt[sind] = (da-ga)*y[sind];

			}
			}


			void CellAdhesion::initState(double *s){
			int n= lattice->getNbCells();
			// allnotch.resize(n);
			// alldelta.resize(n);
			//resize allnotch, alldelta & state ???
			for (int c=0;c<n;c++){
			uint idx = lattice->getCell(c).stateVectorIndex();
			alldelta[idx] = s[c*3];
			allnotch[idx] = s[c*3+1];
			state[idx] = s[c*3+2];
			}
			}