void Cell::output_func(adevs::Bag<CellEvent>& yb)
{
CellEvent e;
// Assume we are dying
e.value = Dead;
// Check in case this in not true
if (check_born_rule())
{
e.value = Alive;
}
// Set the initial visualization value
if (vis_phase != NULL) *vis_phase = e.value;
// Generate an event for each neighbor
for (long int dx = -1; dx <= 1; dx++)
{
for (long int dy = -1; dy <= 1; dy++)
{
e.x = (x+dx)%w;
e.y = (y+dy)%h;
if (e.x < 0) e.x = w-1;
if (e.y < 0) e.y = h-1;
// Don't send to self
if (e.x != x || e.y != y)
{
yb.insert(e);
}
}
}
}
void FourZoneBuildingExt::output_func(const double *q, const bool* state_event,
adevs::Bag<OMC_ADEVS_IO_TYPE>& yb)
{
FourZoneBuilding::output_func(q,state_event,yb);
update_vars(q,false);
PortValue<BuildingEvent*> pv;
pv.port = BuildingModelInterface::tempData;
pv.value = new TemperatureEvent(OUTDOOR_THERMOMETER,0,get_outdoor_d1());
yb.insert(pv);
pv.value = new TemperatureEvent(THERMOSTAT_THERMOMETER,0,get_z1_t1());
yb.insert(pv);
pv.value = new TemperatureEvent(THERMOSTAT_THERMOMETER,1,get_z2_t1());
yb.insert(pv);
pv.value = new TemperatureEvent(THERMOSTAT_THERMOMETER,2,get_z3_t1());
yb.insert(pv);
pv.value = new TemperatureEvent(THERMOSTAT_THERMOMETER,3,get_z4_t1());
yb.insert(pv);
pv.port = extraTempData;
pv.value = new TemperatureEvent(THERMOSTAT_THERMOMETER,4,get_znoise_t1());
yb.insert(pv);
}
void SolutionQueue::output_func(adevs::Bag<QCell>& yb){
yb.insert(solutionlist.front());
}
void Cell::output_func(adevs::Bag<vec>& yb){
/* EXPERIMENTAL BEHAVIOR CODE */
// // If the cell hasn't attached
// if(getTimeAttached() == 0.0){
// // check to see if the cell could attach
// vec dist = p + ta()*v;
//
// if(dist[2] < BioChem::getInstance()->getFilterTopPos()){
// // see what would happen if the cell did attach
// dist[2] = BioChem::getInstance()->getFilterTopPos();
//
// // look at all the other cells
// CellPopulation* pop = CellPopulation::getInstance();
// CellPopulation::iterator iter;
// bool settled = false;
// for(iter = pop->begin(); iter != pop->end(); iter++){
// vec position = (*iter)->getPos();
// do{
// // could this cell collide?
// if(getDistance(position, dist) < 1.0){
// // avoid collision by stacking
// vec stackem;
// stackem[0] = r.uniform(-1.0,1.0);
// stackem[1] = r.uniform(-1.0,1.0);
// stackem[2] = 0.0;
// dist = position + stackem;
// }
// } while (!settled);
// }
//
// p += ta()*v;
// p[2] = BioChem::getInstance()->getFilterTopPos();
// } else {
// // set the position to the new position
// p = dist;
// }
// }
/* Move in current direction */
if (getTime() >= getTimeAttached()) p += ta()*v;
// if (getTime() >= 3.0) p += ta()*v;
// Don't go below 0.0
if (p[0] < (double)0.3)
p[0] = (double)0.3;
if (p[1] < (double)0.3)
p[1] = (double)0.3;
// Don't go above Width or Depth
if (p[0] > BioChem::getInstance()->getDomain()[0]-(double)0.3)
p[0] = BioChem::getInstance()->getDomain()[0]-(double)0.3;
if (p[1] > BioChem::getInstance()->getDomain()[1]-(double)0.3)
p[1] = BioChem::getInstance()->getDomain()[1]-(double)0.3;
// Don't go above or below the filter
if (p[2] < BioChem::getInstance()->getFilterBottomPos())
p[2] = BioChem::getInstance()->getFilterBottomPos();
if (p[2] > BioChem::getInstance()->getFilterTopPos())
p[2] = BioChem::getInstance()->getFilterTopPos();
/* Output new position */
yb.insert(p);
}