Homec4science
Diffusion tisue modeling (master) Browse

tismod/src/cellInteraction.cpp
b7305f7b5a8bmaster

Last Change
Raw File

cellInteraction.cpp
View Options

#include "cellInteraction.hpp"
using namespace std;
void CellInteraction::writeState(){
cout<<"state "<<state_dim<<": ";
for(uint i=0;i<state_dim;i++){
cout<<i<<": "<<(*state)[i]<<"; ";
}
cout<<endl;
}
void CellInteraction::initState(){
void *pt = (void *) state;
lattice->forAllCellsAndBoundaries(CellInteraction::initCellState,pt,CellInteraction::initBoundaryState,pt);
}
void CellInteraction::updateDeadCellState(state_type& s){
visitorParams2 vp(&s,this);
lattice->forAllCellsAndBoundaries(CellInteraction::setEmptyCellValue,
(void *)&vp,NULL,NULL);
}
void CellInteraction::initCellState(Cell *ce, void *pt){
visitorParams2 *vp = (visitorParams2 *) pt;
vp->ci->initCell(ce,vp->state);
}
void CellInteraction::setEmptyCellValue(Cell *ce, void *pt){
visitorParams2 *vp = (visitorParams2 *) pt;
if(ce->isEmpty()){
uint idx = ce->stateVectorIndex();
uint dim = vp->ci->cellStateDim();
for(uint i=0;i<dim;i++){
vp->state->at(idx++) = vp->ci->deadCellValue();
}
}
}
void CellInteraction::initBoundaryState(Boundary *bo, void *pt){
visitorParams2 *vp = (visitorParams2 *) pt;
vp->ci->initBoundary(bo,vp->state);
}
void CellInteraction::initCell(Cell *ce, state_type *state){
uint idx = ce->stateVectorIndex();
uint dim = cellStateDim();
for(uint i=0;i<dim;i++){
state->at(idx++) = 0.0;
}
}
void CellInteraction::initBoundary(Boundary *bo, state_type *state){
uint idx = bo->stateVectorIndex();
uint dim = boundaryStateDim();
for(uint i=0;i<dim;i++){
state->at(idx++) = 0.0;
}
}
uint CellInteraction::setStateVectorIndex(uint idx){
uint cnt=idx;
visitorParams vp(&cnt, this);
lattice->forAllCellsAndBoundaries(CellInteraction::setCellIndex,
(void *)&vp,
CellInteraction::setBoundaryIndex,
(void *)&vp);
return cnt;
}
void CellInteraction::setCellIndex(Cell *ce, void *pt){
visitorParams *vp = (visitorParams *) pt;
uint *idx = vp->cnt;
ce->setStateVectorIndex(*idx);
*idx += vp->ci->cellStateDim();
}
void CellInteraction::setBoundaryIndex(Boundary *bo, void *pt){
visitorParams *vp = (visitorParams *) pt;
uint *idx = vp->cnt;
bo->setStateVectorIndex(*idx);
*idx += vp->ci->boundaryStateDim();
}
void CellInteraction::collectCellIdx(Cell *c,void *arg){
vector<uint> *v = (vector<uint> *)arg;
v->push_back(c->stateVectorIndex());
}
c4science · Help