version 1.153, 2014/06/20 16:45:46
|
version 1.159, 2014/09/01 10:34:10
|
Line 1
|
Line 1
|
/* $Id$ |
/* $Id$ |
$State$ |
$State$ |
$Log$ |
$Log$ |
|
Revision 1.159 2014/09/01 10:34:10 brouard |
|
Summary: WIN32 |
|
Author: Brouard |
|
|
|
Revision 1.158 2014/08/27 17:11:51 brouard |
|
*** empty log message *** |
|
|
|
Revision 1.157 2014/08/27 16:26:55 brouard |
|
Summary: Preparing windows Visual studio version |
|
Author: Brouard |
|
|
|
In order to compile on Visual studio, time.h is now correct and time_t |
|
and tm struct should be used. difftime should be used but sometimes I |
|
just make the differences in raw time format (time(&now). |
|
Trying to suppress #ifdef LINUX |
|
Add xdg-open for __linux in order to open default browser. |
|
|
|
Revision 1.156 2014/08/25 20:10:10 brouard |
|
*** empty log message *** |
|
|
|
Revision 1.155 2014/08/25 18:32:34 brouard |
|
Summary: New compile, minor changes |
|
Author: Brouard |
|
|
|
Revision 1.154 2014/06/20 17:32:08 brouard |
|
Summary: Outputs now all graphs of convergence to period prevalence |
|
|
Revision 1.153 2014/06/20 16:45:46 brouard |
Revision 1.153 2014/06/20 16:45:46 brouard |
Summary: If 3 live state, convergence to period prevalence on same graph |
Summary: If 3 live state, convergence to period prevalence on same graph |
Author: Brouard |
Author: Brouard |
Line 455
|
Line 482
|
#include <stdio.h> |
#include <stdio.h> |
#include <stdlib.h> |
#include <stdlib.h> |
#include <string.h> |
#include <string.h> |
|
|
|
#ifdef _WIN32 |
|
#include <io.h> |
|
#else |
#include <unistd.h> |
#include <unistd.h> |
|
#endif |
|
|
#include <limits.h> |
#include <limits.h> |
#include <sys/types.h> |
#include <sys/types.h> |
#include <sys/stat.h> |
#include <sys/stat.h> |
#include <errno.h> |
#include <errno.h> |
extern int errno; |
/* extern int errno; */ |
|
|
|
/* #ifdef LINUX */ |
|
/* #include <time.h> */ |
|
/* #include "timeval.h" */ |
|
/* #else */ |
|
/* #include <sys/time.h> */ |
|
/* #endif */ |
|
|
#ifdef LINUX |
|
#include <time.h> |
#include <time.h> |
#include "timeval.h" |
|
#else |
|
#include <sys/time.h> |
|
#endif |
|
|
|
#ifdef GSL |
#ifdef GSL |
#include <gsl/gsl_errno.h> |
#include <gsl/gsl_errno.h> |
Line 500 extern int errno;
|
Line 534 extern int errno;
|
#define AGESUP 130 |
#define AGESUP 130 |
#define AGEBASE 40 |
#define AGEBASE 40 |
#define AGEGOMP 10. /**< Minimal age for Gompertz adjustment */ |
#define AGEGOMP 10. /**< Minimal age for Gompertz adjustment */ |
#ifdef UNIX |
#ifdef _WIN32 |
#define DIRSEPARATOR '/' |
|
#define CHARSEPARATOR "/" |
|
#define ODIRSEPARATOR '\\' |
|
#else |
|
#define DIRSEPARATOR '\\' |
#define DIRSEPARATOR '\\' |
#define CHARSEPARATOR "\\" |
#define CHARSEPARATOR "\\" |
#define ODIRSEPARATOR '/' |
#define ODIRSEPARATOR '/' |
|
#else |
|
#define DIRSEPARATOR '/' |
|
#define CHARSEPARATOR "/" |
|
#define ODIRSEPARATOR '\\' |
#endif |
#endif |
|
|
/* $Id$ */ |
/* $Id$ */ |
/* $State$ */ |
/* $State$ */ |
|
|
char version[]="Imach version 0.98nT, January 2014,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121)"; |
char version[]="Imach version 0.98nX, August 2014,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121)"; |
char fullversion[]="$Revision$ $Date$"; |
char fullversion[]="$Revision$ $Date$"; |
char strstart[80]; |
char strstart[80]; |
char optionfilext[10], optionfilefiname[FILENAMELENGTH]; |
char optionfilext[10], optionfilefiname[FILENAMELENGTH]; |
Line 587 char popfile[FILENAMELENGTH];
|
Line 621 char popfile[FILENAMELENGTH];
|
|
|
char optionfilegnuplot[FILENAMELENGTH], optionfilehtm[FILENAMELENGTH], optionfilehtmcov[FILENAMELENGTH] ; |
char optionfilegnuplot[FILENAMELENGTH], optionfilehtm[FILENAMELENGTH], optionfilehtmcov[FILENAMELENGTH] ; |
|
|
struct timeval start_time, end_time, curr_time, last_time, forecast_time; |
/* struct timeval start_time, end_time, curr_time, last_time, forecast_time; */ |
struct timezone tzp; |
/* struct timezone tzp; */ |
extern int gettimeofday(); |
/* extern int gettimeofday(); */ |
struct tm tmg, tm, tmf, *gmtime(), *localtime(); |
struct tm tml, *gmtime(), *localtime(); |
long time_value; |
|
extern long time(); |
extern time_t time(); |
|
|
|
struct tm start_time, end_time, curr_time, last_time, forecast_time; |
|
time_t rstart_time, rend_time, rcurr_time, rlast_time, rforecast_time; /* raw time */ |
|
struct tm tm; |
|
|
char strcurr[80], strfor[80]; |
char strcurr[80], strfor[80]; |
|
|
char *endptr; |
char *endptr; |
Line 1053 double brent(double ax, double bx, doubl
|
Line 1092 double brent(double ax, double bx, doubl
|
{ |
{ |
int iter; |
int iter; |
double a,b,d,etemp; |
double a,b,d,etemp; |
double fu,fv,fw,fx; |
double fu=0,fv,fw,fx; |
double ftemp; |
double ftemp; |
double p,q,r,tol1,tol2,u,v,w,x,xm; |
double p,q,r,tol1,tol2,u,v,w,x,xm; |
double e=0.0; |
double e=0.0; |
Line 1236 void powell(double p[], double **xi, int
|
Line 1275 void powell(double p[], double **xi, int
|
xits=vector(1,n); |
xits=vector(1,n); |
*fret=(*func)(p); |
*fret=(*func)(p); |
for (j=1;j<=n;j++) pt[j]=p[j]; |
for (j=1;j<=n;j++) pt[j]=p[j]; |
|
rcurr_time = time(NULL); |
for (*iter=1;;++(*iter)) { |
for (*iter=1;;++(*iter)) { |
fp=(*fret); |
fp=(*fret); |
ibig=0; |
ibig=0; |
del=0.0; |
del=0.0; |
last_time=curr_time; |
rlast_time=rcurr_time; |
(void) gettimeofday(&curr_time,&tzp); |
/* (void) gettimeofday(&curr_time,&tzp); */ |
printf("\nPowell iter=%d -2*LL=%.12f %ld sec. %ld sec.",*iter,*fret, curr_time.tv_sec-last_time.tv_sec, curr_time.tv_sec-start_time.tv_sec);fflush(stdout); |
rcurr_time = time(NULL); |
fprintf(ficlog,"\nPowell iter=%d -2*LL=%.12f %ld sec. %ld sec.",*iter,*fret, curr_time.tv_sec-last_time.tv_sec, curr_time.tv_sec-start_time.tv_sec); fflush(ficlog); |
curr_time = *localtime(&rcurr_time); |
/* fprintf(ficrespow,"%d %.12f %ld",*iter,*fret,curr_time.tv_sec-start_time.tv_sec); */ |
printf("\nPowell iter=%d -2*LL=%.12f %ld sec. %ld sec.",*iter,*fret, rcurr_time-rlast_time, rcurr_time-rstart_time);fflush(stdout); |
|
fprintf(ficlog,"\nPowell iter=%d -2*LL=%.12f %ld sec. %ld sec.",*iter,*fret,rcurr_time-rlast_time, rcurr_time-rstart_time); fflush(ficlog); |
|
/* fprintf(ficrespow,"%d %.12f %ld",*iter,*fret,curr_time.tm_sec-start_time.tm_sec); */ |
for (i=1;i<=n;i++) { |
for (i=1;i<=n;i++) { |
printf(" %d %.12f",i, p[i]); |
printf(" %d %.12f",i, p[i]); |
fprintf(ficlog," %d %.12lf",i, p[i]); |
fprintf(ficlog," %d %.12lf",i, p[i]); |
Line 1254 void powell(double p[], double **xi, int
|
Line 1296 void powell(double p[], double **xi, int
|
fprintf(ficlog,"\n"); |
fprintf(ficlog,"\n"); |
fprintf(ficrespow,"\n");fflush(ficrespow); |
fprintf(ficrespow,"\n");fflush(ficrespow); |
if(*iter <=3){ |
if(*iter <=3){ |
tm = *localtime(&curr_time.tv_sec); |
tml = *localtime(&rcurr_time); |
strcpy(strcurr,asctime(&tm)); |
strcpy(strcurr,asctime(&tml)); |
/* asctime_r(&tm,strcurr); */ |
/* asctime_r(&tm,strcurr); */ |
forecast_time=curr_time; |
rforecast_time=rcurr_time; |
itmp = strlen(strcurr); |
itmp = strlen(strcurr); |
if(strcurr[itmp-1]=='\n') /* Windows outputs with a new line */ |
if(strcurr[itmp-1]=='\n') /* Windows outputs with a new line */ |
strcurr[itmp-1]='\0'; |
strcurr[itmp-1]='\0'; |
printf("\nConsidering the time needed for this last iteration #%d: %ld seconds,\n",*iter,curr_time.tv_sec-last_time.tv_sec); |
printf("\nConsidering the time needed for this last iteration #%d: %ld seconds,\n",*iter,rcurr_time-rlast_time); |
fprintf(ficlog,"\nConsidering the time needed for this last iteration #%d: %ld seconds,\n",*iter,curr_time.tv_sec-last_time.tv_sec); |
fprintf(ficlog,"\nConsidering the time needed for this last iteration #%d: %ld seconds,\n",*iter,rcurr_time-rlast_time); |
for(niterf=10;niterf<=30;niterf+=10){ |
for(niterf=10;niterf<=30;niterf+=10){ |
forecast_time.tv_sec=curr_time.tv_sec+(niterf-*iter)*(curr_time.tv_sec-last_time.tv_sec); |
rforecast_time=rcurr_time+(niterf-*iter)*(rcurr_time-rlast_time); |
tmf = *localtime(&forecast_time.tv_sec); |
forecast_time = *localtime(&rforecast_time); |
/* asctime_r(&tmf,strfor); */ |
/* asctime_r(&tmf,strfor); */ |
strcpy(strfor,asctime(&tmf)); |
strcpy(strfor,asctime(&forecast_time)); |
itmp = strlen(strfor); |
itmp = strlen(strfor); |
if(strfor[itmp-1]=='\n') |
if(strfor[itmp-1]=='\n') |
strfor[itmp-1]='\0'; |
strfor[itmp-1]='\0'; |
printf(" - if your program needs %d iterations to converge, convergence will be \n reached in %s i.e.\n on %s (current time is %s);\n",niterf, asc_diff_time(forecast_time.tv_sec-curr_time.tv_sec,tmpout),strfor,strcurr); |
printf(" - if your program needs %d iterations to converge, convergence will be \n reached in %s i.e.\n on %s (current time is %s);\n",niterf, asc_diff_time(rforecast_time-rcurr_time,tmpout),strfor,strcurr); |
fprintf(ficlog," - if your program needs %d iterations to converge, convergence will be \n reached in %s i.e.\n on %s (current time is %s);\n",niterf, asc_diff_time(forecast_time.tv_sec-curr_time.tv_sec,tmpout),strfor,strcurr); |
fprintf(ficlog," - if your program needs %d iterations to converge, convergence will be \n reached in %s i.e.\n on %s (current time is %s);\n",niterf, asc_diff_time(rforecast_time-rcurr_time,tmpout),strfor,strcurr); |
} |
} |
} |
} |
for (i=1;i<=n;i++) { |
for (i=1;i<=n;i++) { |
Line 3892 fprintf(fichtm," \n<ul><li><b>Graphs</b>
|
Line 3934 fprintf(fichtm," \n<ul><li><b>Graphs</b>
|
before but expressed in per year i.e. quasi incidences if stepm is small and probabilities too: <a href=\"%s%d_2.png\">%s%d_2.png</a><br> \ |
before but expressed in per year i.e. quasi incidences if stepm is small and probabilities too: <a href=\"%s%d_2.png\">%s%d_2.png</a><br> \ |
<img src=\"%s%d_2.png\">",stepm,subdirf2(optionfilefiname,"pe"),jj1,subdirf2(optionfilefiname,"pe"),jj1,subdirf2(optionfilefiname,"pe"),jj1); |
<img src=\"%s%d_2.png\">",stepm,subdirf2(optionfilefiname,"pe"),jj1,subdirf2(optionfilefiname,"pe"),jj1,subdirf2(optionfilefiname,"pe"),jj1); |
/* Period (stable) prevalence in each health state */ |
/* Period (stable) prevalence in each health state */ |
for(cpt=1; cpt<nlstate;cpt++){ |
for(cpt=1; cpt<=nlstate;cpt++){ |
fprintf(fichtm,"<br>- Period (stable) prevalence in each health state : <a href=\"%s%d_%d.png\">%s%d_%d.png</a><br> \ |
fprintf(fichtm,"<br>- Convergence from each state (1 to %d) to period (stable) prevalence in state %d <a href=\"%s%d_%d.png\">%s%d_%d.png</a><br> \ |
<img src=\"%s%d_%d.png\">",subdirf2(optionfilefiname,"p"),cpt,jj1,subdirf2(optionfilefiname,"p"),cpt,jj1,subdirf2(optionfilefiname,"p"),cpt,jj1); |
<img src=\"%s%d_%d.png\">",nlstate, cpt, subdirf2(optionfilefiname,"p"),cpt,jj1,subdirf2(optionfilefiname,"p"),cpt,jj1,subdirf2(optionfilefiname,"p"),cpt,jj1); |
} |
} |
for(cpt=1; cpt<=nlstate;cpt++) { |
for(cpt=1; cpt<=nlstate;cpt++) { |
fprintf(fichtm,"\n<br>- Life expectancy by health state (%d) at initial age and its decomposition into health expectancies : <a href=\"%s%d%d.png\">%s%d%d.png</a> <br> \ |
fprintf(fichtm,"\n<br>- Life expectancy by health state (%d) at initial age and its decomposition into health expectancies in each alive state (1 to %d) : <a href=\"%s%d%d.png\">%s%d%d.png</a> <br> \ |
<img src=\"%s%d%d.png\">",cpt,subdirf2(optionfilefiname,"exp"),cpt,jj1,subdirf2(optionfilefiname,"exp"),cpt,jj1,subdirf2(optionfilefiname,"exp"),cpt,jj1); |
<img src=\"%s%d%d.png\">",cpt,nlstate,subdirf2(optionfilefiname,"exp"),cpt,jj1,subdirf2(optionfilefiname,"exp"),cpt,jj1,subdirf2(optionfilefiname,"exp"),cpt,jj1); |
} |
} |
} /* end i1 */ |
} /* end i1 */ |
}/* End k1 */ |
}/* End k1 */ |
Line 4084 plot [%.f:%.f] \"%s\" every :::%d::%d u
|
Line 4126 plot [%.f:%.f] \"%s\" every :::%d::%d u
|
k=3; |
k=3; |
fprintf(ficgp,"\n#\n#\n#CV preval stable (period): 'pij' files, cov=%d state=%d",k1, cpt); |
fprintf(ficgp,"\n#\n#\n#CV preval stable (period): 'pij' files, cov=%d state=%d",k1, cpt); |
fprintf(ficgp,"\nset out \"%s%d_%d.png\" \n",subdirf2(optionfilefiname,"p"),cpt,k1); |
fprintf(ficgp,"\nset out \"%s%d_%d.png\" \n",subdirf2(optionfilefiname,"p"),cpt,k1); |
l=(nlstate+ndeath)*(cpt-1)+1; |
|
fprintf(ficgp,"set xlabel \"Age\" \nset ylabel \"Probability\" \n\ |
fprintf(ficgp,"set xlabel \"Age\" \nset ylabel \"Probability\" \n\ |
set ter png small size 320, 240\n\ |
set ter png small size 320, 240\n\ |
unset log y\n\ |
unset log y\n\ |
Line 4094 plot [%.f:%.f] ", ageminpar, agemaxpar)
|
Line 4135 plot [%.f:%.f] ", ageminpar, agemaxpar)
|
fprintf(ficgp,"\"%s\"",subdirf2(fileres,"pij")); |
fprintf(ficgp,"\"%s\"",subdirf2(fileres,"pij")); |
else |
else |
fprintf(ficgp,", '' "); |
fprintf(ficgp,", '' "); |
fprintf(ficgp," u ($1==%d ? ($3):1/0):($%d/($%d",k1,k+l,k+l); |
l=(nlstate+ndeath)*(i-1)+1; |
|
fprintf(ficgp," u ($1==%d ? ($3):1/0):($%d/($%d",k1,k+l+(cpt-1),k+l); |
for (j=1; j<= (nlstate-1) ; j ++) |
for (j=1; j<= (nlstate-1) ; j ++) |
fprintf(ficgp,"+$%d",k+l+j); |
fprintf(ficgp,"+$%d",k+l+j); |
fprintf(ficgp,")) t \"prev(%d,%d)\" w l",i,cpt); |
fprintf(ficgp,")) t \"prev(%d,%d)\" w l",i,cpt); |
Line 5172 calandcheckages(int imx, int maxwav, dou
|
Line 5214 calandcheckages(int imx, int maxwav, dou
|
} |
} |
else if(agev[m][i] >*agemax){ |
else if(agev[m][i] >*agemax){ |
*agemax=agev[m][i]; |
*agemax=agev[m][i]; |
printf(" Max anint[%d][%d]=%.0f annais[%d]=%.0f, agemax=%.2f\n",m,i,anint[m][i], i,annais[i], *agemax); |
/* printf(" Max anint[%d][%d]=%.0f annais[%d]=%.0f, agemax=%.2f\n",m,i,anint[m][i], i,annais[i], *agemax);*/ |
} |
} |
/*agev[m][i]=anint[m][i]-annais[i];*/ |
/*agev[m][i]=anint[m][i]-annais[i];*/ |
/* agev[m][i] = age[i]+2*m;*/ |
/* agev[m][i] = age[i]+2*m;*/ |
Line 5301 int main(int argc, char *argv[])
|
Line 5343 int main(int argc, char *argv[])
|
/* setlocale (LC_MESSAGES, ""); */ |
/* setlocale (LC_MESSAGES, ""); */ |
|
|
/* gettimeofday(&start_time, (struct timezone*)0); */ /* at first time */ |
/* gettimeofday(&start_time, (struct timezone*)0); */ /* at first time */ |
(void) gettimeofday(&start_time,&tzp); |
rstart_time = time(NULL); |
|
/* (void) gettimeofday(&start_time,&tzp);*/ |
|
start_time = *localtime(&rstart_time); |
curr_time=start_time; |
curr_time=start_time; |
tm = *localtime(&start_time.tv_sec); |
/*tml = *localtime(&start_time.tm_sec);*/ |
tmg = *gmtime(&start_time.tv_sec); |
/* strcpy(strstart,asctime(&tml)); */ |
strcpy(strstart,asctime(&tm)); |
strcpy(strstart,asctime(&start_time)); |
|
|
/* printf("Localtime (at start)=%s",strstart); */ |
/* printf("Localtime (at start)=%s",strstart); */ |
/* tp.tv_sec = tp.tv_sec +86400; */ |
/* tp.tm_sec = tp.tm_sec +86400; */ |
/* tm = *localtime(&start_time.tv_sec); */ |
/* tm = *localtime(&start_time.tm_sec); */ |
/* tmg.tm_year=tmg.tm_year +dsign*dyear; */ |
/* tmg.tm_year=tmg.tm_year +dsign*dyear; */ |
/* tmg.tm_mon=tmg.tm_mon +dsign*dmonth; */ |
/* tmg.tm_mon=tmg.tm_mon +dsign*dmonth; */ |
/* tmg.tm_hour=tmg.tm_hour + 1; */ |
/* tmg.tm_hour=tmg.tm_hour + 1; */ |
/* tp.tv_sec = mktime(&tmg); */ |
/* tp.tm_sec = mktime(&tmg); */ |
/* strt=asctime(&tmg); */ |
/* strt=asctime(&tmg); */ |
/* printf("Time(after) =%s",strstart); */ |
/* printf("Time(after) =%s",strstart); */ |
/* (void) time (&time_value); |
/* (void) time (&time_value); |
Line 5334 int main(int argc, char *argv[])
|
Line 5378 int main(int argc, char *argv[])
|
i=strlen(pathr); |
i=strlen(pathr); |
if(pathr[i-1]=='\n') |
if(pathr[i-1]=='\n') |
pathr[i-1]='\0'; |
pathr[i-1]='\0'; |
|
i=strlen(pathr); |
|
if(pathr[i-1]==' ') /* This may happen when dragging on oS/X! */ |
|
pathr[i-1]='\0'; |
for (tok = pathr; tok != NULL; ){ |
for (tok = pathr; tok != NULL; ){ |
printf("Pathr |%s|\n",pathr); |
printf("Pathr |%s|\n",pathr); |
while ((val = strsep(&tok, "\"" )) != NULL && *val == '\0'); |
while ((val = strsep(&tok, "\"" )) != NULL && *val == '\0'); |
Line 5389 int main(int argc, char *argv[])
|
Line 5436 int main(int argc, char *argv[])
|
path=%s \n\ |
path=%s \n\ |
optionfile=%s\n\ |
optionfile=%s\n\ |
optionfilext=%s\n\ |
optionfilext=%s\n\ |
optionfilefiname=%s\n",pathimach,pathtot,path,optionfile,optionfilext,optionfilefiname); |
optionfilefiname='%s'\n",pathimach,pathtot,path,optionfile,optionfilext,optionfilefiname); |
|
|
printf("Local time (at start):%s",strstart); |
printf("Local time (at start):%s",strstart); |
fprintf(ficlog,"Local time (at start): %s",strstart); |
fprintf(ficlog,"Local time (at start): %s",strstart); |
fflush(ficlog); |
fflush(ficlog); |
/* (void) gettimeofday(&curr_time,&tzp); */ |
/* (void) gettimeofday(&curr_time,&tzp); */ |
/* printf("Elapsed time %d\n", asc_diff_time(curr_time.tv_sec-start_time.tv_sec,tmpout)); */ |
/* printf("Elapsed time %d\n", asc_diff_time(curr_time.tm_sec-start_time.tm_sec,tmpout)); */ |
|
|
/* */ |
/* */ |
strcpy(fileres,"r"); |
strcpy(fileres,"r"); |
Line 5405 int main(int argc, char *argv[])
|
Line 5452 int main(int argc, char *argv[])
|
/*---------arguments file --------*/ |
/*---------arguments file --------*/ |
|
|
if((ficpar=fopen(optionfile,"r"))==NULL) { |
if((ficpar=fopen(optionfile,"r"))==NULL) { |
printf("Problem with optionfile %s with errno=%s\n",optionfile,strerror(errno)); |
printf("Problem with optionfile '%s' with errno='%s'\n",optionfile,strerror(errno)); |
fprintf(ficlog,"Problem with optionfile %s with errno=%s\n",optionfile,strerror(errno)); |
fprintf(ficlog,"Problem with optionfile '%s' with errno='%s'\n",optionfile,strerror(errno)); |
fflush(ficlog); |
fflush(ficlog); |
/* goto end; */ |
/* goto end; */ |
exit(70); |
exit(70); |
Line 6627 Interval (in months) between two waves:
|
Line 6674 Interval (in months) between two waves:
|
} |
} |
printf("See log file on %s\n",filelog); |
printf("See log file on %s\n",filelog); |
/* gettimeofday(&end_time, (struct timezone*)0);*/ /* after time */ |
/* gettimeofday(&end_time, (struct timezone*)0);*/ /* after time */ |
(void) gettimeofday(&end_time,&tzp); |
/*(void) gettimeofday(&end_time,&tzp);*/ |
tm = *localtime(&end_time.tv_sec); |
rend_time = time(NULL); |
tmg = *gmtime(&end_time.tv_sec); |
end_time = *localtime(&rend_time); |
strcpy(strtend,asctime(&tm)); |
/* tml = *localtime(&end_time.tm_sec); */ |
|
strcpy(strtend,asctime(&end_time)); |
printf("Local time at start %s\nLocal time at end %s",strstart, strtend); |
printf("Local time at start %s\nLocal time at end %s",strstart, strtend); |
fprintf(ficlog,"Local time at start %s\nLocal time at end %s\n",strstart, strtend); |
fprintf(ficlog,"Local time at start %s\nLocal time at end %s\n",strstart, strtend); |
printf("Total time used %s\n", asc_diff_time(end_time.tv_sec -start_time.tv_sec,tmpout)); |
printf("Total time used %s\n", asc_diff_time(rend_time -rstart_time,tmpout)); |
|
|
printf("Total time was %ld Sec.\n", end_time.tv_sec -start_time.tv_sec); |
printf("Total time was %.0lf Sec.\n", difftime(rend_time,rstart_time)); |
fprintf(ficlog,"Total time used %s\n", asc_diff_time(end_time.tv_sec -start_time.tv_sec,tmpout)); |
fprintf(ficlog,"Total time used %s\n", asc_diff_time(rend_time -rstart_time,tmpout)); |
fprintf(ficlog,"Total time was %ld Sec.\n", end_time.tv_sec -start_time.tv_sec); |
fprintf(ficlog,"Total time was %.0lf Sec.\n", difftime(rend_time,rstart_time)); |
/* printf("Total time was %d uSec.\n", total_usecs);*/ |
/* printf("Total time was %d uSec.\n", total_usecs);*/ |
/* if(fileappend(fichtm,optionfilehtm)){ */ |
/* if(fileappend(fichtm,optionfilehtm)){ */ |
fprintf(fichtm,"<br>Local time at start %s<br>Local time at end %s<br>\n</body></html>",strstart, strtend); |
fprintf(fichtm,"<br>Local time at start %s<br>Local time at end %s<br>\n</body></html>",strstart, strtend); |
Line 6656 Interval (in months) between two waves:
|
Line 6704 Interval (in months) between two waves:
|
printf("Current directory %s!\n",pathcd); |
printf("Current directory %s!\n",pathcd); |
/*strcat(plotcmd,CHARSEPARATOR);*/ |
/*strcat(plotcmd,CHARSEPARATOR);*/ |
sprintf(plotcmd,"gnuplot"); |
sprintf(plotcmd,"gnuplot"); |
#ifndef UNIX |
#ifdef _WIN32 |
sprintf(plotcmd,"\"%sgnuplot.exe\"",pathimach); |
sprintf(plotcmd,"\"%sgnuplot.exe\"",pathimach); |
#endif |
#endif |
if(!stat(plotcmd,&info)){ |
if(!stat(plotcmd,&info)){ |
printf("Error gnuplot program not found: %s\n",plotcmd);fflush(stdout); |
printf("Error or gnuplot program not found: '%s'\n",plotcmd);fflush(stdout); |
if(!stat(getenv("GNUPLOTBIN"),&info)){ |
if(!stat(getenv("GNUPLOTBIN"),&info)){ |
printf("Error gnuplot program not found: %s Environment GNUPLOTBIN not set.\n",plotcmd);fflush(stdout); |
printf("Error or gnuplot program not found: '%s' Environment GNUPLOTBIN not set.\n",plotcmd);fflush(stdout); |
}else |
}else |
strcpy(pplotcmd,plotcmd); |
strcpy(pplotcmd,plotcmd); |
#ifdef UNIX |
#ifdef __unix |
strcpy(plotcmd,GNUPLOTPROGRAM); |
strcpy(plotcmd,GNUPLOTPROGRAM); |
if(!stat(plotcmd,&info)){ |
if(!stat(plotcmd,&info)){ |
printf("Error gnuplot program not found: %s\n",plotcmd);fflush(stdout); |
printf("Error gnuplot program not found: '%s'\n",plotcmd);fflush(stdout); |
}else |
}else |
strcpy(pplotcmd,plotcmd); |
strcpy(pplotcmd,plotcmd); |
#endif |
#endif |
Line 6676 Interval (in months) between two waves:
|
Line 6724 Interval (in months) between two waves:
|
strcpy(pplotcmd,plotcmd); |
strcpy(pplotcmd,plotcmd); |
|
|
sprintf(plotcmd,"%s %s",pplotcmd, optionfilegnuplot); |
sprintf(plotcmd,"%s %s",pplotcmd, optionfilegnuplot); |
printf("Starting graphs with: %s\n",plotcmd);fflush(stdout); |
printf("Starting graphs with: '%s'\n",plotcmd);fflush(stdout); |
|
|
if((outcmd=system(plotcmd)) != 0){ |
if((outcmd=system(plotcmd)) != 0){ |
printf("\n Problem with gnuplot command %s, err=%d\n", plotcmd, outcmd); |
printf("gnuplot command might not be in your path: '%s', err=%d\n", plotcmd, outcmd); |
printf("\n Trying on same directory\n"); |
printf("\n Trying if gnuplot resides on the same directory that IMaCh\n"); |
sprintf(plotcmd,"%sgnuplot %s", pathimach, optionfilegnuplot); |
sprintf(plotcmd,"%sgnuplot %s", pathimach, optionfilegnuplot); |
if((outcmd=system(plotcmd)) != 0) |
if((outcmd=system(plotcmd)) != 0) |
printf("\n Still a problem with gnuplot command %s, err=%d\n", plotcmd, outcmd); |
printf("\n Still a problem with gnuplot command %s, err=%d\n", plotcmd, outcmd); |
} |
} |
printf(" Wait..."); |
printf(" Successful, please wait..."); |
while (z[0] != 'q') { |
while (z[0] != 'q') { |
/* chdir(path); */ |
/* chdir(path); */ |
printf("\nType e to edit output files, g to graph again and q for exiting: "); |
printf("\nType e to edit results with your browser, g to graph again and q for exit: "); |
scanf("%s",z); |
scanf("%s",z); |
/* if (z[0] == 'c') system("./imach"); */ |
/* if (z[0] == 'c') system("./imach"); */ |
if (z[0] == 'e') { |
if (z[0] == 'e') { |
#ifdef OSX |
#ifdef __APPLE__ |
sprintf(pplotcmd, "open %s", optionfilehtm); |
sprintf(pplotcmd, "open %s", optionfilehtm); |
|
#elif __linux |
|
sprintf(pplotcmd, "xdg-open %s", optionfilehtm); |
#else |
#else |
sprintf(pplotcmd, "%s", optionfilehtm); |
sprintf(pplotcmd, "%s", optionfilehtm); |
#endif |
#endif |
Line 6709 Interval (in months) between two waves:
|
Line 6759 Interval (in months) between two waves:
|
scanf("%s",z); |
scanf("%s",z); |
} |
} |
} |
} |
|
|
|
|
|
|