int main(void)
{
halInit();
chSysInit();
chMtxInit(&printMtx);
#if DEBUG_TIME_MEASURING == TRUE
rtcRetrieve(&RTC_DRIVER, &timingData);
PRINT("Started. Minute: %d Seconds: %d", timingData.time.minute,
timingData.time.seconds);
#endif
initialiseDebugHw();
initialiseCC3000();
initialiseControl(&controlInfo);
initialiseSensorHw();
clarityTransportInformation tcp;
memset(&tcp, 0, sizeof(tcp));
tcp.type = CLARITY_TRANSPORT_TCP;
tcp.addr.type = CLARITY_ADDRESS_URL;
strncpy(tcp.addr.addr.url, SERVER_URL, CLARITY_MAX_URL_LENGTH);
tcp.port = SERVER_PORT;
PRINT("Starting...", NULL);
if (clarityInit(&cc3000ApiMutex, cc3000Unresponsive, &ap, debugPrint) != CLARITY_SUCCESS)
{
PRINT_ERROR();
}
clarityHttpPersistant persistant;
memset(&persistant,0,sizeof(persistant));
persistant.closeOnComplete = false;
clarityTimeDate time;
rtcRetrieve(&RTC_DRIVER, &time);
clarityRegisterProcessStarted();
if (time.date.year < 14)
{
PRINT("Time needs updated.", NULL);
if (updateRtcWithSntp() != 0)
{
PRINT_ERROR();
}
}
else
{
PRINT("Time doesn't need updated.", NULL);
}
if (eepromWasLastShutdownOk() != EEPROM_ERROR_OK)
{
PRINT("Last shutdown was not OK.", NULL);
if (httpPostShutdownError(&tcp, &persistant) == 0)
{
if (eepromAcknowledgeLastShutdownError() != EEPROM_ERROR_OK)
{
PRINT_ERROR();
}
}
else
{
PRINT_ERROR();
}
}
else
{
PRINT("Last shutdown was OK.", NULL);
}
PRINT("Posting Lux.", NULL);
if (httpPostLux(&tcp, &persistant) != CLARITY_SUCCESS)
{
PRINT_ERROR();
}
PRINT("Posting Temperature.", NULL);
if (httpPostTemperature(&tcp, &persistant) != CLARITY_SUCCESS)
{
PRINT_ERROR();
}
persistant.closeOnComplete = true;
PRINT("Posting Pressure.", NULL);
if (httpPostPressure(&tcp, &persistant) != CLARITY_SUCCESS)
{
PRINT_ERROR();
}
#if 1
if (clarityHttpServerStart(&controlInfo) != CLARITY_SUCCESS)
{
PRINT_ERROR();
}
while(1)
{
if (palReadPad(BUTTON_PORT, BUTTON_PAD))
{
break;
}
PRINT("main sleeping", NULL);
chThdSleep(MS2ST(100));
}
PRINT("Shutting down...", NULL);
if (clarityHttpServerStop() != CLARITY_SUCCESS)
{
PRINT("clarityHttpServerStop() failed", NULL);
}
#endif
clarityRegisterProcessFinished();
PRINT("Done.", NULL);
if (clarityShutdown() != CLARITY_SUCCESS)
{
PRINT("clarityShutdown() failed", NULL);
}
PRINT("Shut down.", NULL);
deinitialiseCC3000();
#if DEBUG_TIME_MEASURING == TRUE
rtcRetrieve(&RTC_DRIVER, &timingData);
PRINT("Started. Minute: %d Seconds: %d", timingData.time.minute,
timingData.time.seconds);
#endif
configureRtcAlarmAndStandby(&RTC_DRIVER, STANDBY_TIME_S);
return 0;
}
int main (int argc, char* argv[])
{
int id; // process rank
int p; // number of processes
//int num;
int nMax;
MPI_Init (&argc, &argv);
MPI_Comm_rank (MPI_COMM_WORLD, &id);
MPI_Comm_size (MPI_COMM_WORLD, &p);
//FILE *fin;
controlStruct control;
acfStruct acfStructure;
noiseStruct noiseStructure;
char fname[1024]; // read in parameter file
char Tname[1024]; // read in scintillation timescale
int nt;
double *tdiss;
//char oname[1024]; // output file name
//char pname[1024]; // file to plot
char dname[1024]; // graphics device
int i;
int n = 1; // number of dynamic spectrum to simulate
double flux0;
double flux1;
int plotMode = 0; // plot existing dynamic spectrum, off by default
control.noplot = 1; // don't show dynamic spectrum while simulationg, on by default
// read options
for (i=0;i<argc;i++)
{
if (strcmp(argv[i],"-f") == 0)
{
strcpy(fname,argv[++i]);
strcpy(Tname,argv[++i]);
printf ("Parameters are in %s\n", fname);
}
//else if (strcmp(argv[i],"-o")==0)
//{
// strcpy(oname,argv[++i]);
// //printf ("Dynamic spectrum is output into %s\n", oname);
//}
else if (strcmp(argv[i],"-n")==0)
{
n = atoi (argv[++i]);
printf ("Number of dynamic spectrum to simulate %d\n", n);
}
//else if (strcmp(argv[i],"-p")==0) // just plot
//{
// strcpy(pname,argv[++i]);
// plotMode = 1;
// printf ("Plotting exiting dynamic spectrum.\n");
//}
//else if (strcmp(argv[i],"-noplot")==0)
//{
// control.noplot = 1; // Don't show dynamic spetrum while simulationg
//}
//else if (strcmp(argv[i],"-dev")==0)
//{
// strcpy(dname,argv[++i]);
// printf ("Graphic device: %s\n", dname);
//}
}
nt = readDissNum (Tname);
tdiss = (double *)malloc(sizeof(double)*nt);
readDiss (Tname, tdiss);
sprintf(dname, "%s", "1/xs");
/*
if ((fin=fopen(oname, "w"))==NULL)
{
printf ("Can't open file...\n");
exit(1);
}
*/
if (plotMode==0)
{
// Simulate dynamic spectrum
// read parameters
initialiseControl(&control);
readParams (fname,dname,n,&control);
//readParams (fname,oname,dname,n,&control);
printf ("Finished reading parameters.\n");
//preAllocateMemory (&acfStructure, &control);
//allocateMemory (&acfStructure);
//num = (int)((control.scint_ts1-control.scint_ts0)/control.scint_ts_step);
//for (tdiff=control.scint_ts0; tdiff<control.scint_ts1; tdiff+=control.scint_ts_step)
//for (i=id; i<=num; i+=p)
for (i=id; i<nt; i+=p)
{
control.nsub = tdiss[i];
control.tsub = control.T/(double)(control.nsub);
control.scint_ts = control.tsub/2.0;
control.whiteLevel = control.whiteLevel0*sqrt(control.nsub*control.nchan); // 0.1 gives 1 to a 10*10 dynamic spectrum
calNoise (&noiseStructure, &control);
//printf ("%lf\n", noiseStructure.detection);
//printf ("tdiff fdiff: %lf %lf\n", tdiff, fdiff);
flux0 = control.cFlux0;
flux1 = control.cFlux1;
acfStructure.probability = 1.0;
calculateScintScale (&acfStructure, &control);
//printf ("calculateScintScal\n");
nMax = 0;
while (fabs(acfStructure.probability-0.8)>=PRECISION && nMax <= 100)
{
control.cFlux = flux0+(flux1-flux0)/2.0;
//printf ("%lf %lf %.8lf %.3f\n", tdiff, fdiff, control.cFlux, acfStructure.probability);
// simulate dynamic spectra
calculateNDynSpec (&acfStructure, &control, &noiseStructure);
//if (control.noplot==0 && n == 1)
//{
// // plot while simulating
// heatMap (&acfStructure, dname);
//}
// merged into calculateNDynSpec
//qualifyVar (&acfStructure, &noiseStructure, &control);
if (acfStructure.probability>0.8)
{
flux1 = control.cFlux;
}
else
{
flux0 = control.cFlux;
}
nMax++;
//printf ("%lf %f %d\n", control.cFlux, acfStructure.probability, nMax);
}
printf ("%d %lf %lf %f %lf %d\n", control.nsub, control.whiteLevel, noiseStructure.detection, control.cFlux, acfStructure.probability, nMax);
deallocateMemory (&acfStructure);
}
fflush (stdout);
MPI_Finalize ();
printf ("Threshold: %f \n", noiseStructure.detection);
// deallocate memory
}
//else
//{
// plotDynSpec(pname, dname);
//}
/*
if (fclose(fin))
{
printf ("Can't close file...\n");
exit(1);
}
*/
free(tdiss);
return 0;
}
