version 1.297, 2019/05/22 17:56:10
|
version 1.301, 2019/06/04 13:51:20
|
Line 1
|
Line 1
|
/* $Id$ |
/* $Id$ |
$State$ |
$State$ |
$Log$ |
$Log$ |
|
Revision 1.301 2019/06/04 13:51:20 brouard |
|
Summary: Error in 'r'parameter file backcast yearsbproj instead of yearsfproj |
|
|
|
Revision 1.300 2019/05/22 19:09:45 brouard |
|
Summary: version 0.99r19 of May 2019 |
|
|
|
Revision 1.299 2019/05/22 18:37:08 brouard |
|
Summary: Cleaned 0.99r19 |
|
|
|
Revision 1.298 2019/05/22 18:19:56 brouard |
|
*** empty log message *** |
|
|
Revision 1.297 2019/05/22 17:56:10 brouard |
Revision 1.297 2019/05/22 17:56:10 brouard |
Summary: Fix bug by moving date2dmy and nhstepm which gaefin=-1 |
Summary: Fix bug by moving date2dmy and nhstepm which gaefin=-1 |
|
|
Line 1114 typedef struct {
|
Line 1126 typedef struct {
|
/* $State$ */ |
/* $State$ */ |
#include "version.h" |
#include "version.h" |
char version[]=__IMACH_VERSION__; |
char version[]=__IMACH_VERSION__; |
char copyright[]="April 2018,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121), Intel Software 2015-2018"; |
char copyright[]="May 2019,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121), Intel Software 2015-2020"; |
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 5206 void prevalence(double ***probs, double
|
Line 5218 void prevalence(double ***probs, double
|
|
|
void concatwav(int wav[], int **dh, int **bh, int **mw, int **s, double *agedc, double **agev, int firstpass, int lastpass, int imx, int nlstate, int stepm) |
void concatwav(int wav[], int **dh, int **bh, int **mw, int **s, double *agedc, double **agev, int firstpass, int lastpass, int imx, int nlstate, int stepm) |
{ |
{ |
/* Concatenates waves: wav[i] is the number of effective (useful waves) of individual i. |
/* Concatenates waves: wav[i] is the number of effective (useful waves in the sense that a non interview is useless) of individual i. |
Death is a valid wave (if date is known). |
Death is a valid wave (if date is known). |
mw[mi][i] is the mi (mi=1 to wav[i]) effective wave of individual i |
mw[mi][i] is the mi (mi=1 to wav[i]) effective wave of individual i |
dh[m][i] or dh[mw[mi][i]][i] is the delay between two effective waves m=mw[mi][i] |
dh[m][i] or dh[mw[mi][i]][i] is the delay between two effective waves m=mw[mi][i] |
and mw[mi+1][i]. dh depends on stepm. |
and mw[mi+1][i]. dh depends on stepm. s[m][i] exists for any wave from firstpass to lastpass |
*/ |
*/ |
|
|
int i=0, mi=0, m=0, mli=0; |
int i=0, mi=0, m=0, mli=0; |
Line 12227 Please run with mle=-1 to get a correct
|
Line 12239 Please run with mle=-1 to get a correct
|
fputs(line,stdout); |
fputs(line,stdout); |
fputs(line,ficparo); |
fputs(line,ficparo); |
fputs(line,ficlog); |
fputs(line,ficlog); |
|
fputs(line,ficres); |
continue; |
continue; |
}else |
}else |
break; |
break; |
Line 12272 Please run with mle=-1 to get a correct
|
Line 12285 Please run with mle=-1 to get a correct
|
fputs(line,stdout); |
fputs(line,stdout); |
fputs(line,ficparo); |
fputs(line,ficparo); |
fputs(line,ficlog); |
fputs(line,ficlog); |
|
fputs(line,ficres); |
continue; |
continue; |
}else |
}else |
break; |
break; |
Line 12297 Please run with mle=-1 to get a correct
|
Line 12311 Please run with mle=-1 to get a correct
|
fputs(line,stdout); |
fputs(line,stdout); |
fputs(line,ficparo); |
fputs(line,ficparo); |
fputs(line,ficlog); |
fputs(line,ficlog); |
|
fputs(line,ficres); |
continue; |
continue; |
}else |
}else |
break; |
break; |
Line 12330 Please run with mle=-1 to get a correct
|
Line 12345 Please run with mle=-1 to get a correct
|
fputs(line,stdout); |
fputs(line,stdout); |
fputs(line,ficparo); |
fputs(line,ficparo); |
fputs(line,ficlog); |
fputs(line,ficlog); |
|
fputs(line,ficres); |
continue; |
continue; |
}else if(sscanf(line,"prevforecast=%[^\n]\n",modeltemp)) |
}else if(sscanf(line,"prevforecast=%[^\n]\n",modeltemp)) |
parameterline=11; |
parameterline=11; |
Line 12353 Please run with mle=-1 to get a correct
|
Line 12369 Please run with mle=-1 to get a correct
|
prvforecast = 1; |
prvforecast = 1; |
} |
} |
else if((num_filled=sscanf(line,"prevforecast=%d yearsfproj=%lf mobil_average=%d\n",&prevfcast,&yrfproj,&mobilavproj)) !=EOF){/* && (num_filled == 3))*/ |
else if((num_filled=sscanf(line,"prevforecast=%d yearsfproj=%lf mobil_average=%d\n",&prevfcast,&yrfproj,&mobilavproj)) !=EOF){/* && (num_filled == 3))*/ |
printf(" Num_filled=%d, yearsfproj=%lf, mobil_average=%d\n",prevfcast,yrfproj,mobilavproj); |
printf("prevforecast=%d yearsfproj=%lf mobil_average=%d\n",prevfcast,yrfproj,mobilavproj); |
|
fprintf(ficlog,"prevforecast=%d yearsfproj=%lf mobil_average=%d\n",prevfcast,yrfproj,mobilavproj); |
|
fprintf(ficres,"prevforecast=%d yearsfproj=%lf mobil_average=%d\n",prevfcast,yrfproj,mobilavproj); |
prvforecast = 2; |
prvforecast = 2; |
} |
} |
else { |
else { |
Line 12374 Please run with mle=-1 to get a correct
|
Line 12392 Please run with mle=-1 to get a correct
|
prvbackcast = 1; |
prvbackcast = 1; |
} |
} |
else if((num_filled=sscanf(line,"prevbackcast=%d yearsbproj=%lf mobil_average=%d\n",&prevbcast,&yrbproj,&mobilavproj)) ==3){/* && (num_filled == 3))*/ |
else if((num_filled=sscanf(line,"prevbackcast=%d yearsbproj=%lf mobil_average=%d\n",&prevbcast,&yrbproj,&mobilavproj)) ==3){/* && (num_filled == 3))*/ |
printf(" Num_filled=%d, yearsbproj=%lf, mobil_average=%d\n",prevbcast,yrbproj,mobilavproj); |
printf("prevbackcast=%d yearsbproj=%lf mobil_average=%d\n",prevbcast,yrbproj,mobilavproj); |
|
fprintf(ficlog,"prevbackcast=%d yearsbproj=%lf mobil_average=%d\n",prevbcast,yrbproj,mobilavproj); |
|
fprintf(ficres,"prevbackcast=%d yearsbproj=%lf mobil_average=%d\n",prevbcast,yrbproj,mobilavproj); |
prvbackcast = 2; |
prvbackcast = 2; |
} |
} |
else { |
else { |
Line 12383 Please run with mle=-1 to get a correct
|
Line 12403 Please run with mle=-1 to get a correct
|
goto end; |
goto end; |
} |
} |
break; |
break; |
/* /\*fscanf(ficpar,"backcast=%d starting-back-date=%lf/%lf/%lf final-back-date=%lf/%lf/%lf mobil_average=%d\n",&backcast,&jback1,&mback1,&anback1,&jback2,&mback2,&anback2,&mobilavproj);*\/ */ |
|
/* if((num_filled=sscanf(line,"backcast=%d starting-back-date=%lf/%lf/%lf final-back-date=%lf/%lf/%lf mobil_average=%d\n",&backcast,&jback1,&mback1,&anback1,&jback2,&mback2,&anback2,&mobilavproj)) !=EOF){ */ |
|
/* if (num_filled != 8) { */ |
|
/* printf("Error: Not 8 (data)parameters in line but %d, for example:backcast=1 starting-back-date=1/1/1990 final-back-date=1/1/1970 mobil_average=1\n, your line=%s . Probably you are running an older format.\n",num_filled,line); */ |
|
/* fprintf(ficlog,"Error: Not 8 (data)parameters in line but %d, for example:backcast=1 starting-back-date=1/1/1990 final-back-date=1/1/1970 mobil_average=1\n, your line=%s . Probably you are running an older format.\n",num_filled,line); */ |
|
/* goto end; */ |
|
/* } */ |
|
/* printf("backcast=%d starting-back-date=%.lf/%.lf/%.lf final-back-date=%.lf/%.lf/%.lf mobil_average=%d\n",backcast,jback1,mback1,anback1,jback2,mback2,anback2,mobilavproj); */ |
|
/* fprintf(ficparo,"backcast=%d starting-back-date=%.lf/%.lf/%.lf final-back-date=%.lf/%.lf/%.lf mobil_average=%d\n",backcast,jback1,mback1,anback1,jback2,mback2,anback2,mobilavproj); */ |
|
/* fprintf(ficlog,"backcast=%d starting-back-date=%.lf/%.lf/%.lf final-back-date=%.lf/%.lf/%.lf mobil_average=%d\n",backcast,jback1,mback1,anback1,jback2,mback2,anback2,mobilavproj); */ |
|
/* fprintf(ficres,"backcast=%d starting-back-date=%.lf/%.lf/%.lf final-back-date=%.lf/%.lf/%.lf mobil_average=%d\n",backcast,jback1,mback1,anback1,jback2,mback2,anback2,mobilavproj); */ |
|
/* /\* day and month of proj2 are not used but only year anproj2.*\/ */ |
|
/* dateback1=anback1+(mback1-1)/12.+(jback1-1)/365.; */ |
|
/* dateback2=anback2+(mback2-1)/12.+(jback2-1)/365.; */ |
|
/* } */ |
|
/* break; */ |
|
case 13: |
case 13: |
if((num_filled=sscanf(line,"result:%[^\n]\n",resultline)) !=EOF){ |
if((num_filled=sscanf(line,"result:%[^\n]\n",resultline)) !=EOF){ |
if (num_filled == 0){ |
if (num_filled == 0){ |