int Model::GreedyAlgo()
{
_start = clock();
initSystem();
for (t = 0; t < T; t++){
startAssignment();
bool assignComplete = false; //flag to indicate whether all requests are satisfied or pending for next time period
while(1)
{
reviewSystem();
//assignment: assign the highest arc delay/demand to a candidate vehicle
Arc* pArc;
Vehicle* pVeh;
bool isVehtoArcAssigned = assignVehtoArc(pArc, pVeh);
updateSystem(pArc, pVeh, isVehtoArcAssigned);
assignComplete = true;
k = 0;
for (vector<Arc>::iterator iArc = arc.begin(); iArc != arc.end(); iArc++, k++){
if(iArc->from == iArc->to) continue;
if(!arcIsAssigned[&(*iArc)]){
assignComplete = false; //assign is complete only if all arcs are assigned
break;
}
}
if(assignComplete)
break;
//if node load maximum is reached -- cancel vehicles toward this station and/or cancel previous route
} //done while loop
k = 0;
for (vector<Arc>::iterator iArc = arc.begin(); iArc != arc.end(); iArc++, k++){
if (iArc->from != iArc->to) continue;
v = 0;
for (vector<Vehicle>::iterator iVeh = sys.vehicle.begin(); iVeh != sys.vehicle.end(); ++iVeh, v++){
if (!(iVeh->assigned) && iVeh->to == iArc->to && iVeh->time == t){
sys.route[k][v][t] = 1; //stay in the station if the vehicle is not assigned toward any other stations
sys.NodeLoad[iArc->to][t]++; //update node load
iVeh->time = t + 1; //update time
iVeh->powLvl += maxC; //update battery level
}
}
}
}
// compute the total cost
computeTotalCost();
_end = clock();
cmp_time = (double)(_end - _start) / CLOCKS_PER_SEC;
sprintf(path, "%soutput/GreedyAlgo_%dV_%dN_%dArc_%dT.txt", directoryPath.c_str(), nVeh, N, nArc, T);
heuristicSol();
return 1;
}
void SoSpring::inputChanged(SoField * whichField)
{
if (whichField == &trigger)
{
// Start the thing!
updateSystem();
timer->setInterval(timeStep.getValue());
timer->schedule();
evaluate();
}
}
void Engine::updateSystems()
{
// restart the game clock
m_clock.restart();
// update the time
m_taskManager.addTask(
[this]
{
this->m_time = this->m_clock.getElapsedTime();
}, true, true);
// perform per-system update
// launches repeating tasks for the registered systems
for (auto& system : m_systems)
{
updateSystem(system.get(), true, false);
}
// launch task manager
m_taskManager.start();
}
void performAction(int fd,char *buf)
/*--------------------------------------------
* esegue operazione richiesta dal client (web o telnet)
* int fd - descrittore del socket
* char *buf - stringa con il comando
* attualmente implementate:
* checkBuonGiorno
* checkBuonaNotte
* presenze_attiva
* presenze_disattiva
* irrigazione_attiva id_circuito
* irrigazione_disattiva id_circuito
* bg_attiva
* bn_attiva
* bg_disattiva
* bn_disattiva
* get
* channels (restituisce i descrittori dei canali)
* values (restituisce i valori dei canali)
* latest_values (restituisce gli ultimi valori eng dei canali)
* outchannels
* outputs
* mm_names
* mm_values
* knx_values
* value adId channel_n
* panel panel_id
* presenze_attivo
* irrigazione_attivo irrigazione_n
* buongiorno_attivo
* buonanotte_attivo
* irrigazione_rt irrigazione_n (real time table)
* irrigazione irrigazione_n (mysql tabel)
* pioggia
* effemeridi
* sottostato inizio_estate inizio_inverno (restituisce antiintrusione,pioggia,giorno,estate)
* set
* knx_value id_knx value
* output system_n channel_n value
* output ch_id value
* irrigazione id_irrigazione circuito ora_start ripetitivita tempo_off durata validita
* toggle system_n channel_n
* reset system_n (sends reset signal)
* status
* reload (spegne e riaccende tutti i sistemi)
* system (ricarica tabella system)
* analog (ricarica tabella analog)
* digital (ricarica tabella digital)
* out (ricarica tabella digitalOut)
* presenze (ricarica scenari_presenze)
* bgbn (ricarica scenari bgbn)
* dig system_n channel_n (restituisce corrispondente canale digitale)
* ana system_n channel_n (restituisce corrispondente canale analogico)
* out system_n channel_n (restituisce corrispondente uscita digitale)
* update
* analog channel_id
* digital channel_id
* out channel_id
* bg
* bn
* presenze
* irrigazione
* system
* update system_n (ricarica scheda system_n)
* disable system_n (disabilita scheda system_n)
* enable system_n (abilita scheda system_n)
* chn
* disable system_n adId channel_n (disabilita canale channel_n tipo adId su scheda system_n)
* enable system_n adId channel_n (abilita canale channel_n tipo adId su scheda system_n)
*-----------------------------------------*/
{
char buffer[255];
int value;
int i,j,k;
int state;
MYSQL *connection=NULL, mysql;
int systemId;
int chnId;
int adId;
char query[255];
char *c,*d;
struct timeval tim;
double t1;
int rem_time;
if(strstr(buf,"checkBuonGiorno")==buf)
{
*buonGiornoAttivo=checkBuonGiornoAttivo();
sprintf(buffer,"%d %s",*buonGiornoAttivo,TERMINATOR);
if (send(fd, buffer,strlen(buffer), 0) == -1)
perror("send");
return;
}
if(strstr(buf,"checkBuonaNotte")==buf)
{
*buonaNotteAttivo=checkBuonaNotteAttivo();
sprintf(buffer,"%d %s",*buonaNotteAttivo,TERMINATOR);
if (send(fd, buffer,strlen(buffer), 0) == -1)
perror("send");
return;
}
if(strstr(buf,"presenze_attiva")==buf)
{
*scenariPresenzeAttivo=1;
sprintf(buffer,"%d %s",*scenariPresenzeAttivo,TERMINATOR);
if (send(fd, buffer,strlen(buffer), 0) == -1)
perror("send");
return;
}
if(strstr(buf,"irrigazione_attiva ")==buf)
{
i=atoi(&buf[19]);
j=id_irrigazioneToId(i);
if(j>=0)
{
if(circuiti_attivi(j)>0)
{
irrigazioneTable[j].attivo=1;
irrigazione_attiva(j,irrigazioneTable[j].msg_avvio_manuale);
}
sprintf(buffer,"%d %s",irrigazioneTable[j].attivo,TERMINATOR);
if (send(fd, buffer,strlen(buffer), 0) == -1)
perror("send");
}
else
{
sprintf(buffer,"-1 %s",TERMINATOR);
if (send(fd, buffer,strlen(buffer), 0) == -1)
perror("send");
}
return;
}
if(strstr(buf,"bg_attiva")==buf)
{
*buonaNotteAttivo=0;
for(i=0;i<SCENARIBGBNCOUNT;i++)
{
if((scenariBgBnTable[i].bg==0)&&(scenariBgBnTable[i].attivo))
{
j=id_digitalOutToId(scenariBgBnTable[i].id_digital_out);
digitalOutTable[j].start_time=0;
digitalOutTable[j].req_value=digitalOutTable[j].value;
}
}
*buonGiornoAttivo=1;
for(i=0;i<SCENARIBGBNCOUNT;i++)
{
if((scenariBgBnTable[i].bg)&&(scenariBgBnTable[i].attivo))
{
j=id_digitalOutToId(scenariBgBnTable[i].id_digital_out);
digitalOutTable[j].start_time=time(NULL);
digitalOutTable[j].req_value=digitalOutTable[j].on_value;
digitalOutTable[j].po_delay=scenariBgBnTable[i].ritardo;
}
}
sprintf(buffer,"%d %s",*buonGiornoAttivo,TERMINATOR);
if (send(fd, buffer,strlen(buffer), 0) == -1)
perror("send");
return;
}
if(strstr(buf,"bn_attiva")==buf)
{
*buonGiornoAttivo=0;
for(i=0;i<SCENARIBGBNCOUNT;i++)
{
if((scenariBgBnTable[i].bg)&&(scenariBgBnTable[i].attivo))
{
j=id_digitalOutToId(scenariBgBnTable[i].id_digital_out);
digitalOutTable[j].start_time=0;
digitalOutTable[j].req_value=digitalOutTable[j].value;
}
}
*buonaNotteAttivo=1;
for(i=0;i<SCENARIBGBNCOUNT;i++)
{
if((scenariBgBnTable[i].bg==0)&&(scenariBgBnTable[i].attivo))
{
j=id_digitalOutToId(scenariBgBnTable[i].id_digital_out);
digitalOutTable[j].start_time=time(NULL);
digitalOutTable[j].req_value=1-digitalOutTable[j].on_value;
digitalOutTable[j].po_delay=scenariBgBnTable[i].ritardo;
}
}
sprintf(buffer,"%d %s",*buonaNotteAttivo,TERMINATOR);
if (send(fd, buffer,strlen(buffer), 0) == -1)
perror("send");
return;
}
if(strstr(buf,"bg_disattiva")==buf)
{
*buonGiornoAttivo=0;
for(i=0;i<SCENARIBGBNCOUNT;i++)
{
if((scenariBgBnTable[i].bg)&&(scenariBgBnTable[i].attivo))
{
j=id_digitalOutToId(scenariBgBnTable[i].id_digital_out);
digitalOutTable[j].start_time=0;
digitalOutTable[j].req_value=digitalOutTable[j].value;
}
}
sprintf(buffer,"%d %s",*buonGiornoAttivo,TERMINATOR);
if (send(fd, buffer,strlen(buffer), 0) == -1)
perror("send");
return;
}
if(strstr(buf,"bn_disattiva")==buf)
{
*buonaNotteAttivo=0;
for(i=0;i<SCENARIBGBNCOUNT;i++)
{
if((scenariBgBnTable[i].bg==0)&&(scenariBgBnTable[i].attivo))
{
j=id_digitalOutToId(scenariBgBnTable[i].id_digital_out);
digitalOutTable[j].start_time=0;
digitalOutTable[j].req_value=digitalOutTable[j].value;
}
}
sprintf(buffer,"%d %s",*buonaNotteAttivo,TERMINATOR);
if (send(fd, buffer,strlen(buffer), 0) == -1)
perror("send");
return;
}
if(strstr(buf,"presenze_disattiva")==buf)
{
*scenariPresenzeAttivo=0;
sprintf(buffer,"%d %s",*scenariPresenzeAttivo,TERMINATOR);
if (send(fd, buffer,strlen(buffer), 0) == -1)
perror("send");
return;
}
if(strstr(buf,"irrigazione_disattiva ")==buf)
{
i=atoi(&buf[21]);
i=id_irrigazioneToId(i);
if(i>=0)
{
irrigazioneTable[i].attivo=0;
irrigazione_disattiva(i,irrigazioneTable[i].msg_arresto_manuale);
sprintf(buffer,"%d %s",irrigazioneTable[i].attivo,TERMINATOR);
if (send(fd, buffer,strlen(buffer), 0) == -1)
perror("send");
}
else
{
sprintf(buffer,"-1 %s",TERMINATOR);
if (send(fd, buffer,strlen(buffer), 0) == -1)
perror("send");
}
return;
}
if(strstr(buf,"set")==buf)
{
if(strstr(&buf[3]," knx_value ")==&buf[3])
{
if(c=strchr(&buf[14],' '))
{
*c='\0';
c++;
setKnx(atoi(&buf[14]),atof(c));
}
}
if(strstr(&buf[3]," output ")==&buf[3])
{
if(c=strchr(&buf[11],' '))
{
*c='\0';
c++;
if(d=strchr(c,' '))
{
*d='\0';
d++;
systemId=systemNumToId(atoi(&buf[11]),NUMSYSTEMS);
chnId=atoi(c)-1;
printf("%d %d\n",digitalOutTable[deviceToid[2][systemId][chnId]].id_digital_out, atoi(d));
if(((systemId<NUMSYSTEMS)&&(chnId<systems[systemId].out_ch_d))
&&(digitalOutTable[deviceToid[2][systemId][chnId]].id_digital_out!=-1))
setOutput(digitalOutTable[deviceToid[2][systemId][chnId]].id_digital_out,atoi(d));
}
else
setOutput(atoi(&buf[11]),atoi(c));
return;
}
}
if(strstr(&buf[3]," irrigazione")==&buf[3])
{
set_irrigazione(fd,&buf[15]);
return;
}
return;
}
if(strstr(buf,"toggle ")==buf)
{
if(c=strchr(&buf[7],' '))
{
*c='\0';
c++;
systemId=systemNumToId(atoi(&buf[7]),NUMSYSTEMS);
chnId=atoi(c)-1;
if(((systemId<NUMSYSTEMS)&&(chnId<systems[systemId].out_ch_d))
&&(digitalOutTable[deviceToid[2][systemId][chnId]].id_digital_out!=-1)
&&(digitalOutTable[deviceToid[2][systemId][chnId]].value!=-1))
setOutput(digitalOutTable[deviceToid[2][systemId][chnId]].id_digital_out,
(1-digitalOutTable[deviceToid[2][systemId][chnId]].value));
}
return;
}
if(strstr(buf,"reset ")==buf)
{
sendResetSignal(systemNumToId(atoi(&buf[6]),NUMSYSTEMS));
return;
}
if(strstr(buf,"status")==buf)
{
for(i=0;i<NUMSYSTEMS;i++)
buffer[i]=systems[i].status;
for(j=0;j<strlen(TERMINATOR);j++)
buffer[j+i]=TERMINATOR[j];
if (send(fd, buffer, NUMSYSTEMS+strlen(TERMINATOR), 0) == -1)
perror("send");
return;
}
if(strstr(buf,"reload ")==buf)
{
if(strstr(&buf[7]," system")==&buf[7])
{
loadSystemTable(1);
return;
}
if(strstr(&buf[7]," analog")==&buf[7])
{
loadAnalogTable(1);
return;
}
if(strstr(&buf[7]," digital")==&buf[7])
{
loadDigitalTable(1);
return;
}
if(strstr(&buf[7]," out")==&buf[7])
{
loadDigitalOutTable(1);
return;
}
if(strstr(&buf[7]," presenze")==&buf[7])
{
loadScenariPresenzeTable(1);
return;
}
if(strstr(&buf[7]," bgbn")==&buf[7])
{
loadScenariBgBnTable(1);
return;
}
return;
}
if(strstr(buf,"reload")==buf)
{
*reloadCommand=1;
return;
}
if(strstr(buf,"dig ")==buf)
{
if(c=strchr(&buf[4],' '))
{
*c='\0';
c++;
i=deviceToid[1][systemNumToId(atoi(&buf[4]),NUMSYSTEMS)][atoi(c)-1];
printf("id_digital\tdescription\tdevice_num\tch_num\tenabled\n");
printf("%d\t%s\t%d\t%d\t%d\n",digitalTable[i].id_digital,
digitalTable[i].description,digitalTable[i].device_num,
digitalTable[i].ch_num,digitalTable[i].enabled);
}
return;
}
if(strstr(buf,"ana ")==buf)
{
if(c=strchr(&buf[4],' '))
{
*c='\0';
c++;
i=deviceToid[0][systemNumToId(atoi(&buf[4]),NUMSYSTEMS)][atoi(c)-1];
printf("id_analog\tdescription\tdevice_num\tch_num\tenabled\n");
printf("%d\t%s\t%d\t%d\t%d\n",analogTable[i].id_analog,
analogTable[i].description,analogTable[i].device_num,
analogTable[i].ch_num,analogTable[i].enabled);
}
return;
}
if(strstr(buf,"out ")==buf)
{
if(c=strchr(&buf[4],' '))
{
*c='\0';
c++;
i=deviceToid[2][systemNumToId(atoi(&buf[4]),NUMSYSTEMS)][atoi(c)-1];
printf("id_digital_out\tdescription\tdevice_num\tch_num\tvalue\tpo_delay\tstart\n");
printf("%d\t%s\t%d\t%d\t%d\t%d\t%d\n",digitalOutTable[i].id_digital_out,
digitalOutTable[i].description,digitalOutTable[i].device_num,
digitalOutTable[i].ch_num,digitalOutTable[i].value,
digitalOutTable[i].po_delay,(int)digitalOutTable[i].start_time);
}
return;
}
if(strstr(buf,"update")==buf)
{
if(strstr(&buf[6]," analog ")==&buf[6])
{
updateAnalogChannel(atoi(&buf[14]));
return;
}
if(strstr(&buf[6]," digital ")==&buf[6])
{
updateDigitalChannel(atoi(&buf[15]));
return;
}
if(strstr(&buf[6]," out ")==&buf[6])
{
updateDigitalOutChannel(fd,atoi(&buf[11]));
return;
}
if(strstr(&buf[6]," bg")==&buf[6])
{
updateDayNight(fd);
return;
}
if(strstr(&buf[6]," bn")==&buf[6])
{
updateDayNight(fd);
return;
}
if(strstr(&buf[6]," presenze")==&buf[6])
{
updatePresenze(fd);
return;
}
if(strstr(&buf[6]," irrigazione")==&buf[6])
{
int id_irrigazione=atoi(&buf[18]);
if(id_irrigazione<1)
id_irrigazione=-1;
updateIrrigazioneCircuiti(fd,id_irrigazione);
return;
}
return;
}
if(strstr(buf,"get")==buf)
{
if(strstr(&buf[3]," channels")==&buf[3])
{
for(i=0;i<NUMSYSTEMS;i++)
{
if(systems[i].enabled)
{
for(j=0;j<systems[i].in_ch_an;j++)
{
if(analogTable[deviceToid[0][i][j]].enabled)
{
k=deviceToid[0][i][j];
sprintf(buffer,"0 %d %d %s`",analogTable[k].device_num,
analogTable[k].ch_num,analogTable[k].description);
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
}
for(j=0;j<systems[i].in_ch_d;j++)
{
if(digitalTable[deviceToid[1][i][j]].enabled)
{
k=deviceToid[1][i][j];
sprintf(buffer,"1 %d %d %s`",digitalTable[k].device_num,
digitalTable[k].ch_num,digitalTable[k].description);
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
}
}
}
if (send(fd, TERMINATOR, strlen(TERMINATOR), 0) == -1)
perror("send");
return;
}
if(strstr(&buf[3]," values")==&buf[3])
{
for(i=0;i<NUMSYSTEMS;i++)
{
if(systems[i].enabled)
{
for(j=0;j<systems[i].in_ch_an;j++)
{
if(analogTable[deviceToid[0][i][j]].enabled)
{
k=deviceToid[0][i][j];
sprintf(buffer,"0 %d %d %f %f %f %f %f %d %d`",
analogTable[k].device_num,
analogTable[k].ch_num,
analogTable[k].value_eng,
analogTable[k].value_eng1,
analogTable[k].value_eng2,
analogTable[k].value_eng3,
analogTable[k].value_eng4,
analogTable[k].value,
(analogTable[k].enabled && systems[i].enabled));
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
}
for(j=0;j<systems[i].in_ch_d;j++)
{
if(digitalTable[deviceToid[1][i][j]].enabled)
{
k=deviceToid[1][i][j];
sprintf(buffer,"1 %d %d %d %d %d %d %d %d %d`",
digitalTable[k].device_num,
digitalTable[k].ch_num,
digitalTable[k].value,
digitalTable[k].value1,
digitalTable[k].value2,
digitalTable[k].value3,
digitalTable[k].value4,
digitalTable[k].value,
(digitalTable[k].enabled && systems[i].enabled));
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
}
}
}
for(i=0;i<NUMMULTIMETERS;i++)
{
if(multimeters[i].enabled)
{
for(j=0;j<multimeters[i].out_ch_1+multimeters[i].out_ch_2;j++)
{
if(readingTable[deviceToid[3][i][j]].enabled)
{
k=deviceToid[3][i][j];
sprintf(buffer,"3 %d %d %d %d`",
readingTable[k].multimeter_num,
readingTable[k].ch_num,
readingTable[k].value,
(readingTable[k].enabled && multimeters[i].enabled));
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
}
}
}
if (send(fd, TERMINATOR, strlen(TERMINATOR), 0) == -1)
perror("send");
return;
}
if(strstr(&buf[3]," latest_values")==&buf[3])
{
for(i=0;i<NUMSYSTEMS;i++)
{
if(systems[i].enabled)
{
for(j=0;j<systems[i].in_ch_an;j++)
{
if(analogTable[deviceToid[0][i][j]].enabled)
{
k=deviceToid[0][i][j];
sprintf(buffer,"0 %d %d %f %d`",
analogTable[k].device_num,
analogTable[k].ch_num,
analogTable[k].value_eng,
(analogTable[k].enabled && systems[i].enabled));
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
}
for(j=0;j<systems[i].in_ch_d;j++)
{
if(digitalTable[deviceToid[1][i][j]].enabled)
{
k=deviceToid[1][i][j];
sprintf(buffer,"1 %d %d %d %d`",
digitalTable[k].device_num,
digitalTable[k].ch_num,
digitalTable[k].value,
(digitalTable[k].enabled && systems[i].enabled));
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
}
for(j=0;j<systems[i].out_ch_d;j++)
{
k=deviceToid[2][i][j];
sprintf(buffer,"2 %d %d %d %d`",
digitalOutTable[k].device_num,
digitalOutTable[k].ch_num,
digitalOutTable[k].value,
systems[i].enabled);
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
}
}
for(i=0;i<NUMMULTIMETERS;i++)
{
if(multimeters[i].enabled)
{
for(j=0;j<multimeters[i].out_ch_1+multimeters[i].out_ch_2;j++)
{
if(readingTable[deviceToid[3][i][j]].enabled)
{
k=deviceToid[3][i][j];
sprintf(buffer,"3 %d %d %d %d`",
readingTable[k].multimeter_num,
readingTable[k].ch_num,
readingTable[k].value,
(readingTable[k].enabled && multimeters[i].enabled));
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
}
}
}
if (send(fd, TERMINATOR, strlen(TERMINATOR), 0) == -1)
perror("send");
return;
}
if(strstr(&buf[3]," mm_names")==&buf[3])
{
for(i=0;i<NUMMULTIMETERS;i++)
{
if(multimeters[i].enabled)
{
for(j=0;j<multimeters[i].out_ch_1+multimeters[i].out_ch_2;j++)
{
if(readingTable[deviceToid[3][i][j]].enabled)
{
k=deviceToid[3][i][j];
sprintf(buffer,"%d %d %s`",readingTable[k].multimeter_num,
readingTable[k].ch_num,readingTable[k].description);
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
}
}
}
if (send(fd, TERMINATOR, strlen(TERMINATOR), 0) == -1)
perror("send");
return;
}
if(strstr(&buf[3]," mm_values")==&buf[3])
{
for(i=0;i<NUMMULTIMETERS;i++)
{
if(multimeters[i].enabled)
{
for(j=0;j<multimeters[i].out_ch_1+multimeters[i].out_ch_2;j++)
{
if(readingTable[deviceToid[3][i][j]].enabled)
{
k=deviceToid[3][i][j];
sprintf(buffer,"%d %d %d %d`",
readingTable[k].multimeter_num,
readingTable[k].ch_num,
readingTable[k].value,
(readingTable[k].enabled && multimeters[i].enabled));
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
}
}
}
if (send(fd, TERMINATOR, strlen(TERMINATOR), 0) == -1)
perror("send");
return;
}
if(strstr(&buf[3]," knx_values")==&buf[3])
{
for(i=0;i<NUMKNXCHANNELS;i++)
{
j=id_knx_gatewayToId(knxTable[i].id_knx_gateway);
if((knxTable[i].input_output==1)&&
(knxGateways[j].enabled))
{
if(knxTable[i].data_type=='F')
sprintf(buffer,"%d %d %f %d`",
knxTable[i].id_knx_gateway,
knxTable[i].id_knx_line,
knxTable[i].value_eng,
knxTable[i].enabled);
else
sprintf(buffer,"%d %d %d %d`",
knxTable[i].id_knx_gateway,
knxTable[i].id_knx_line,
knxTable[i].value,
knxTable[i].enabled);
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
}
if (send(fd, TERMINATOR, strlen(TERMINATOR), 0) == -1)
perror("send");
return;
}
if(strstr(&buf[3]," outchannels")==&buf[3])
{
for(i=0;i<NUMSYSTEMS;i++)
{
if(systems[i].enabled)
{
for(j=0;j<systems[i].out_ch_d;j++)
{
if(digitalOutTable[deviceToid[2][i][j]].id_digital_out!=-1)
{
k=deviceToid[2][i][j];
sprintf(buffer,"%d %d %s`",digitalOutTable[k].device_num,
digitalOutTable[k].ch_num,digitalOutTable[k].description);
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
}
}
}
if (send(fd, TERMINATOR, strlen(TERMINATOR), 0) == -1)
perror("send");
return;
}
if(strstr(&buf[3]," outputs")==&buf[3])
{
for(i=0;i<NUMSYSTEMS;i++)
{
if(systems[i].enabled)
{
for(j=0;j<systems[i].out_ch_d;j++)
{
if(digitalOutTable[deviceToid[2][i][j]].id_digital_out!=-1)
{
rem_time=digitalOutTable[deviceToid[2][i][j]].po_delay -
(time(NULL)-digitalOutTable[deviceToid[2][i][j]].start_time);
if(rem_time<0)
rem_time=0;
k=deviceToid[2][i][j];
sprintf(buffer,"%d %d %d %d`",
digitalOutTable[k].device_num,
digitalOutTable[k].ch_num,
digitalOutTable[k].value,
rem_time);
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
}
}
}
if (send(fd, TERMINATOR, strlen(TERMINATOR), 0) == -1)
perror("send");
return;
}
if(strstr(&buf[3]," value ")==&buf[3])
{
if(strlen(&buf[10])>2)
{
buf[11]='\0';
adId=atoi(&buf[10]);
if(adId==0)
{
i=id_analogToId(atoi(&buf[12]));
sprintf(buffer,"%f",analogTable[i].value_eng);
}
else
{
i=id_digitalToId(atoi(&buf[12]));
sprintf(buffer,"%d",digitalTable[i].value);
}
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
if (send(fd, TERMINATOR, strlen(TERMINATOR), 0) == -1)
perror("send");
}
return;
}
if(strstr(&buf[3]," panel ")==&buf[3])
{
if((strlen(&buf[10])>=0)&&(atoi(&buf[10])<NUMPANELS))
{
for(i=0;i<16;i++)
{
sprintf(buffer,"%d %f`",panels[atoi(&buf[10])][i],
analogTable[id_analogToId(panels[atoi(&buf[10])][i])].value_eng);
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
if (send(fd, TERMINATOR, strlen(TERMINATOR), 0) == -1)
perror("send");
}
return;
}
if(strstr(&buf[3]," presenze_attivo")==&buf[3])
{
sprintf(buffer,"%d %s",*scenariPresenzeAttivo,TERMINATOR);
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
return;
}
if(strstr(&buf[3]," irrigazione_attivo ")==&buf[3])
{
i=atoi(&buf[23]);
i=id_irrigazioneToId(i);
if(i>=0)
{
sprintf(buffer,"%d %s",(irrigazioneTable[i].current_circuito>0),TERMINATOR);
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
}
else
{
sprintf(buffer,"-1 %s",TERMINATOR);
if (send(fd, buffer,strlen(buffer), 0) == -1)
perror("send");
}
return;
}
if(strstr(&buf[3]," buongiorno_attivo")==&buf[3])
{
*buonGiornoAttivo=checkBuonGiornoAttivo();
sprintf(buffer,"%d %s",*buonGiornoAttivo,TERMINATOR);
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
return;
}
if(strstr(&buf[3]," buonanotte_attivo")==&buf[3])
{
*buonaNotteAttivo=checkBuonaNotteAttivo();
sprintf(buffer,"%d %s",*buonaNotteAttivo,TERMINATOR);
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
return;
}
if(strstr(&buf[3]," irrigazione_rt ")==&buf[3])
{
get_irrigazione_rt(fd,&buf[18]);
return;
}
if(strstr(&buf[3]," irrigazione ")==&buf[3])
{
get_irrigazione(fd,&buf[15]);
return;
}
if(strstr(&buf[3]," pioggia")==&buf[3])
{
int pioggia=-1;
if(*id_digital_pioggia!=-1)
{
if(digitalTable[*id_digital_pioggia].value!=-1)
pioggia=(digitalTable[*id_digital_pioggia].value==1);
}
sprintf(buffer,"%d %s",pioggia,TERMINATOR);
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
return;
}
if(strstr(&buf[3]," effemeridi")==&buf[3])
{
time_t timer;
struct tm now;
char dawn[6];
char sunset[6];
timer=time(NULL);
now=*localtime(&timer);
i2h(dawn,effemeridiTable[now.tm_mon][now.tm_mday - 1].dawn);
i2h(sunset,effemeridiTable[now.tm_mon][now.tm_mday - 1].sunset);
sprintf(buffer,"%s %s %s",dawn,sunset,TERMINATOR);
if (send(fd, buffer, strlen(buffer), 0) == -1)
perror("send");
return;
}
if(strstr(&buf[3]," sottostato ")==&buf[3])
{
get_sottostato(fd,&buf[15]);
return;
}
return;
}
if(strstr(buf,"system")==buf)
{
if(strstr(&buf[6]," update ")==&buf[6])
{
strcpy(buf,&buf[14]);
while(c=strchr(buf,'|'))
{
*c='\0';
i=atoi(buf);
updateSystem(i);
strcpy(buf,c+1);
}
}
if(strstr(&buf[6]," disable")==&buf[6])
{
if(strlen(&buf[14]))
{
systemId=atoi(&buf[14]);
mysql_init(&mysql);
connection = mysql_real_connect(&mysql,"localhost","root", "minair","socket", 0, NULL,0);
if( connection == NULL )
return;
sprintf(query,"UPDATE system SET enabled=0 WHERE device_num=%d",systems[systemId].device_num);
state = mysql_query(connection, query);
if(state==0)
{
systems[systemId].enabled=0;
systems[systemId].status='d';
}
mysql_close(connection);
return;
}
}
if(strstr(&buf[6]," enable")==&buf[6])
{
if(strlen(&buf[13]))
{
systemId=atoi(&buf[13]);
mysql_init(&mysql);
connection = mysql_real_connect(&mysql,"localhost","root", "minair","socket", 0, NULL,0);
if( connection == NULL )
return;
sprintf(query,"UPDATE system SET enabled=1 WHERE device_num=%d",systems[systemId].device_num);
state = mysql_query(connection, query);
if(state==0)
systems[systemId].enabled=1;
mysql_close(connection);
return;
}
}
}
if(strstr(buf,"chn")==buf)
{
if(strstr(&buf[3]," disable ")==&buf[3])
{
if(strlen(&buf[12])&&(c=strchr(&buf[12],' '))&&(d=strchr(c+1,' ')))
{
*c='\0';
*d='\0';
systemId=atoi(&buf[12]);
adId=atoi(c+1);
chnId=atoi(d+1);
mysql_init(&mysql);
connection = mysql_real_connect(&mysql,"localhost","root", "minair","socket", 0, NULL,0);
if( connection == NULL )
return;
if(adId==0)
sprintf(query,"UPDATE analog SET enabled=0 WHERE device_num=%d AND ch_num=%d",analogTable[deviceToid[0][systemId][chnId]].device_num,analogTable[deviceToid[0][systemId][chnId]].ch_num);
else
sprintf(query,"UPDATE digital SET enabled=0 WHERE device_num=%d AND ch_num=%d",digitalTable[deviceToid[1][systemId][chnId]].device_num,digitalTable[deviceToid[1][systemId][chnId]].ch_num);
state = mysql_query(connection, query);
printf("%s\n",query);
if(state==0)
{
if(adId==0)
analogTable[deviceToid[adId][systemId][chnId]].enabled=0;
else
digitalTable[deviceToid[adId][systemId][chnId]].enabled=0;
}
return;
}
}
if(strstr(&buf[3]," enable ")==&buf[3])
{
if(strlen(&buf[11])&&(c=strchr(&buf[11],' '))&&(d=strchr(c+1,' ')))
{
*c='\0';
*d='\0';
systemId=atoi(&buf[11]);
adId=atoi(c+1);
chnId=atoi(d+1);
mysql_init(&mysql);
connection = mysql_real_connect(&mysql,"localhost","root", "minair","socket", 0, NULL,0);
if( connection == NULL )
return;
if(adId==0)
sprintf(query,"UPDATE analog SET enabled=1 WHERE device_num=%d AND ch_num=%d",analogTable[deviceToid[0][systemId][chnId]].device_num,analogTable[deviceToid[0][systemId][chnId]].ch_num);
else
sprintf(query,"UPDATE digital SET enabled=1 WHERE device_num=%d AND ch_num=%d",digitalTable[deviceToid[1][systemId][chnId]].device_num,digitalTable[deviceToid[1][systemId][chnId]].ch_num);
state = mysql_query(connection, query);
printf("%s\n",query);
if(state==0)
{
if(adId==0)
analogTable[deviceToid[adId][systemId][chnId]].enabled=1;
else
digitalTable[deviceToid[adId][systemId][chnId]].enabled=1;
}
return;
}
}
}
return;
}
/**
* Simulates the system over the next interval of time. The timestep
* is given as argument.
*
* @param ds A C_FLOAT64 specifying the timestep
*/
C_FLOAT64 CTauLeapMethod::doSingleStep(C_FLOAT64 ds)
{
C_FLOAT64 Lambda, Tmp, Tau, Tau1, Tau2;
updatePropensities();
mAvgDX = 0.0;
mSigDX = 0.0;
CMathReaction * pReaction = mReactions.array();
const C_FLOAT64 * pAmu = mAmu.array();
const C_FLOAT64 * pAmuEnd = pAmu + mNumReactions;
const C_FLOAT64 * pFirstSpecies = mContainerState.array() + mFirstReactionSpeciesIndex;
for (; pAmu != pAmuEnd; ++pAmu, ++pReaction)
{
const CMathReaction::Balance & Balance = pReaction->getNumberBalance();
const CMathReaction::SpeciesBalance * it = Balance.array();
const CMathReaction::SpeciesBalance * end = it + Balance.size();
for (; it != end; ++it)
{
mAvgDX[it->first - pFirstSpecies] += it->second **pAmu;
mSigDX[it->first - pFirstSpecies] += it->second * it->second **pAmu;
}
}
Tau1 = Tau2 = std::numeric_limits< C_FLOAT64 >::infinity();
const C_FLOAT64 * pSpecies = mContainerState.array() + mFirstReactionSpeciesIndex;
const C_FLOAT64 * pSpeciesEnd = pSpecies + mNumReactionSpecies;
C_FLOAT64 * pAvgDX = mAvgDX.array();
C_FLOAT64 * pSigDX = mSigDX.array();
for (; pSpecies != pSpeciesEnd; ++pSpecies, ++pAvgDX, ++pSigDX)
{
if ((Tmp = mEpsilon * fabs(*pSpecies)) < 1.0)
Tmp = 1.0;
*pAvgDX = Tmp / fabs(*pAvgDX);
*pSigDX = (Tmp * Tmp) / fabs(*pSigDX);
if (Tau1 > *pAvgDX)
Tau1 = *pAvgDX;
if (Tau2 > *pSigDX)
Tau2 = *pSigDX;
}
Tau = std::min(Tau1, Tau2);
if (ds < Tau)
Tau = ds;
pAmu = mAmu.array();
C_FLOAT64 * pK = mK.array();
C_FLOAT64 * pKEnd = pK + mNumReactions;
for (; pAmu != pAmuEnd; ++pAmu, ++pK)
{
Lambda = *pAmu * Tau;
if (Lambda < 0.0)
CCopasiMessage(CCopasiMessage::EXCEPTION, MCTrajectoryMethod + 10);
else if (Lambda > 2.0e9)
CCopasiMessage(CCopasiMessage::EXCEPTION, MCTrajectoryMethod + 26);
*pK = mpRandomGenerator->getRandomPoisson(Lambda);
}
while (!updateSystem())
{
Tau *= 0.5;
pK = mK.array();
for (; pK != pKEnd; ++pK)
{
*pK *= 0.5;
if (*pK < floor(*pK + 0.75))
{
*pK += mpRandomGenerator->getRandomCC() < 0.5 ? - 0.5 : 0.5;
}
}
}
return Tau;
}
/**
* Simulates the system over the next interval of time. The timestep
* is given as argument.
*
* @param ds A C_FLOAT64 specifying the timestep
*/
C_FLOAT64 CTauLeapMethod::doSingleStep(C_FLOAT64 ds)
{
C_FLOAT64 Lambda, Tmp, Tau, Tau1, Tau2;
updatePropensities();
mAvgDX = 0.0;
mSigDX = 0.0;
std::vector< CReactionDependencies >::const_iterator itReaction = mReactionDependencies.begin();
const C_FLOAT64 * pAmu = mAmu.array();
const C_FLOAT64 * pAmuEnd = pAmu + mNumReactions;
for (; pAmu != pAmuEnd; ++pAmu, ++itReaction)
{
const C_FLOAT64 * pMultiplicity = itReaction->mSpeciesMultiplier.array();
const C_FLOAT64 * pMultiplicityEnd = pMultiplicity + itReaction->mSpeciesMultiplier.size();
const size_t * pIndex = itReaction->mMethodSpeciesIndex.array();
for (; pMultiplicity != pMultiplicityEnd; ++pMultiplicity, ++pIndex)
{
mAvgDX[*pIndex] += *pMultiplicity * *pAmu;
mSigDX[*pIndex] += *pMultiplicity * *pMultiplicity * *pAmu;
}
}
Tau1 = Tau2 = std::numeric_limits< C_FLOAT64 >::infinity();
const C_FLOAT64 * pNumber = mMethodState.beginIndependent() + mFirstReactionSpeciesIndex - 1;
const C_FLOAT64 * pNumberEnd = pNumber + mNumReactionSpecies;
C_FLOAT64 * pAvgDX = mAvgDX.array();
C_FLOAT64 * pSigDX = mSigDX.array();
for (; pNumber != pNumberEnd; ++pNumber, ++pAvgDX, ++pSigDX)
{
if ((Tmp = mEpsilon * fabs(*pNumber)) < 1.0)
Tmp = 1.0;
*pAvgDX = Tmp / fabs(*pAvgDX);
*pSigDX = (Tmp * Tmp) / fabs(*pSigDX);
if (Tau1 > *pAvgDX)
Tau1 = *pAvgDX;
if (Tau2 > *pSigDX)
Tau2 = *pSigDX;
}
Tau = std::min(Tau1, Tau2);
if (ds < Tau)
Tau = ds;
pAmu = mAmu.array();
C_FLOAT64 * pK = mK.array();
C_FLOAT64 * pKEnd = pK + mNumReactions;
for (; pAmu != pAmuEnd; ++pAmu, ++pK)
{
Lambda = *pAmu * Tau;
if (Lambda < 0.0)
CCopasiMessage(CCopasiMessage::EXCEPTION, MCTrajectoryMethod + 10);
else if (Lambda > 2.0e9)
CCopasiMessage(CCopasiMessage::EXCEPTION, MCTrajectoryMethod + 26);
*pK = mpRandomGenerator->getRandomPoisson(Lambda);
}
while (!updateSystem())
{
Tau *= 0.5;
pK = mK.array();
for (; pK != pKEnd; ++pK)
{
*pK *= 0.5;
if (*pK < floor(*pK + 0.75))
{
*pK += mpRandomGenerator->getRandomCC() < 0.5 ? - 0.5 : 0.5;
}
}
}
return Tau;
}
