int main(int argc, char *argv[])
{
/* Command line arguments
* -s simulator
* -n no of simulations
*/
void (*forwardSimulate)(st_part_at* prop_part, double tstep);
forwardSimulate = getForwardSimulator(argc, argv);
int n = getIters(argc, argv, 0, 1);
char *output_file = getOutputFile(argc, argv);
char *dir = getDir(argc, argv);
char *path = createPath(dir, output_file);
printf("Writing to %s\n", path);
FILE *out = fopen(path, "w");
int i, k;
r = gsl_rng_alloc(gsl_rng_mt19937);
/* For timings */
struct timeb t_start, t_end;
double t_diff;
/* Values for scaling parameter*/
int N_C = 41;
double sim_time;/*Simulation time in secs, for a single sim*/
double sc;
/*The following three vectors were generated in R and copied in */
double sc_vec[41] = {1,1.3,1.6,2,2.5,3.2,4,5,6.3,7.9,10,
12.6,15.8,20,25.1,31.6,39.8,50.1, 63.1,79.4,100,
125.9,158.5,199.5,251.2,316.2,398.1,501.2,631,794.3,
1000,1258.9,1584.9,1995.3,2511.9,3162.3,3981.1,5011.9,6309.6,7943.3,
10000};
double X1[41] = {1,1,1,2,2,3,4,5,6,7,9,11,12,14,16,17,19,20,
21,21,22,23,23,23,24,24,24, 24,24,24,24,24,24,24,24,24,24,24,24,24,24};
double X2[41] = {20,20,20,20,20,20,20,20,20,20,21,21,22,
24,26,28,32,36,43,50,60,73,89,110,135,168,209,260,325,407,
510,639,802,1007,1265,1591,2000,2515,3164,3981,5010};
/* Initialise LNA */
pSIv = SetUpSIRI();
pODEIv = pODEInfoAlloc(pgetvoidRI(pSIv));
/* Initialise parameters
* Other paramters and IC are changed in the
* loop
*/
st_part_at *part;
part = timingsInitPrior(0, 0);
VSET(part->params, 0, 2.0);
VSET(part->params, 2, 1.0/50.0);
VSET(part->params, 3, 1.0);
for(k=0; k<N_C; k++) {
sc = sc_vec[k];
VSET(part->params, 1,sc); VSET(part->params, 4, 1.0/(50*sc));
ftime(&t_start);
for(i=0; i<n; i++){
VSET(part->sps, 0, X1[k]); VSET(part->sps, 1, X2[k]);
//VSET(part->sps, 0, 0); VSET(part->sps, 1, 0);
updateRes(part->res);
forwardSimulate(part, 100);
}
ftime(&t_end);
t_diff = toMilliseconds(t_start, t_end);
sim_time = (t_diff/1000.0)/n;
printf("%f, %f\n", sc, sim_time);
fprintf(out, "%f, %f\n", sc, sim_time);
}
ODEInfoDeAlloc(pODEIv);
DestroySIRI(pSIv);
gsl_rng_free(r);
return(GSL_SUCCESS);
}
QPair<int, int> QMailHeartbeatTimer::interval() const
{
const Q_D(QMailHeartbeatTimer);
return qMakePair(toMilliseconds(d->timer->minimumInterval()), toMilliseconds(d->timer->maximumInterval()));
}
int64 milliseconds() const {
return toMilliseconds();
}
