--- imach/src/imach.c	2017/06/30 15:48:31	1.276
+++ imach/src/imach.c	2017/07/19 14:09:02	1.278
@@ -1,6 +1,12 @@
-/* $Id: imach.c,v 1.276 2017/06/30 15:48:31 brouard Exp $
+/* $Id: imach.c,v 1.278 2017/07/19 14:09:02 brouard Exp $
   $State: Exp $
   $Log: imach.c,v $
+  Revision 1.278  2017/07/19 14:09:02  brouard
+  Summary: Bug for mobil_average=0 and prevforecast fixed(?)
+
+  Revision 1.277  2017/07/17 08:53:49  brouard
+  Summary: BOM files can be read now
+
   Revision 1.276  2017/06/30 15:48:31  brouard
   Summary: Graphs improvements
 
@@ -1031,12 +1037,12 @@ typedef struct {
 #define ODIRSEPARATOR '\\'
 #endif
 
-/* $Id: imach.c,v 1.276 2017/06/30 15:48:31 brouard Exp $ */
+/* $Id: imach.c,v 1.278 2017/07/19 14:09:02 brouard Exp $ */
 /* $State: Exp $ */
 #include "version.h"
 char version[]=__IMACH_VERSION__;
 char copyright[]="February 2016,INED-EUROREVES-Institut de longevite-Japan Society for the Promotion of Science (Grant-in-Aid for Scientific Research 25293121), Intel Software 2015-2018";
-char fullversion[]="$Revision: 1.276 $ $Date: 2017/06/30 15:48:31 $"; 
+char fullversion[]="$Revision: 1.278 $ $Date: 2017/07/19 14:09:02 $"; 
 char strstart[80];
 char optionfilext[10], optionfilefiname[FILENAMELENGTH];
 int erreur=0, nberr=0, nbwarn=0; /* Error number, number of errors number of warnings  */
@@ -7942,8 +7948,8 @@ set ter svg size 640, 480\nunset log y\n
    double *agemingoodr, *agemaxgoodr; 
   
   
-   /* modcovmax=2*cptcoveff;/\* Max number of modalities. We suppose  */
-   /* 		   a covariate has 2 modalities, should be equal to ncovcombmax  *\/ */
+   /* modcovmax=2*cptcoveff;  Max number of modalities. We suppose  */
+   /* 		   a covariate has 2 modalities, should be equal to ncovcombmax   */
 
    sumnewp = vector(1,ncovcombmax);
    sumnewm = vector(1,ncovcombmax);
@@ -8267,11 +8273,11 @@ set ter svg size 640, 480\nunset log y\n
 	for(j=1; j<=nlstate+ndeath;j++) {
 	  ppij=0.;
 	  for(i=1; i<=nlstate;i++) {
-	    /* if (mobilav>=1)  */
-	    ppij=ppij+p3mat[i][j][h]*prev[(int)agec][i][k];
-	    /* else { */ /* even if mobilav==-1 we use mobaverage */
-	    /* 	ppij=ppij+p3mat[i][j][h]*probs[(int)(agec)][i][k]; */
-	    /* } */
+	    if (mobilav>=1)
+	     ppij=ppij+p3mat[i][j][h]*prev[(int)agec][i][k];
+	    else { /* even if mobilav==-1 we use mobaverage, probs may not sums to 1 */
+	    	ppij=ppij+p3mat[i][j][h]*probs[(int)(agec)][i][k];
+	    }
 	    fprintf(ficresf," %.3f", p3mat[i][j][h]);
 	  } /* end i */
 	  fprintf(ficresf," %.3f", ppij);
@@ -10613,7 +10619,7 @@ int main(int argc, char *argv[])
   int vpopbased=0;
   int nres=0;
   int endishere=0;
-
+  int noffset=0;
   int ncurrv=0; /* Temporary variable */
   
   char ca[32], cb[32];
@@ -10839,17 +10845,52 @@ int main(int argc, char *argv[])
     fflush(ficlog);
     goto end;
   }
+      /*-------- Rewriting parameter file ----------*/
+  strcpy(rfileres,"r");    /* "Rparameterfile */
+  strcat(rfileres,optionfilefiname);    /* Parameter file first name */
+  strcat(rfileres,".");    /* */
+  strcat(rfileres,optionfilext);    /* Other files have txt extension */
+  if((ficres =fopen(rfileres,"w"))==NULL) {
+    printf("Problem writing new parameter file: %s\n", rfileres);goto end;
+    fprintf(ficlog,"Problem writing new parameter file: %s\n", rfileres);goto end;
+    fflush(ficlog);
+    goto end;
+  }
+  fprintf(ficres,"#IMaCh %s\n",version);
 
+				      
   /* Reads comments: lines beginning with '#' */
   numlinepar=0;
-
-    /* First parameter line */
+  /* Is it a BOM UTF-8 Windows file? */
+  /* First parameter line */
   while(fgets(line, MAXLINE, ficpar)) {
+    noffset=0;
+    if( line[0] == (char)0xEF && line[1] == (char)0xBB) /* EF BB BF */
+    {
+      noffset=noffset+3;
+      printf("# File is an UTF8 Bom.\n"); // 0xBF
+    }
+    else if( line[0] == (char)0xFE && line[1] == (char)0xFF)
+    {
+      noffset=noffset+2;
+      printf("# File is an UTF16BE BOM file\n");
+    }
+    else if( line[0] == 0 && line[1] == 0)
+    {
+      if( line[2] == (char)0xFE && line[3] == (char)0xFF){
+	noffset=noffset+4;
+	printf("# File is an UTF16BE BOM file\n");
+      }
+    } else{
+      ;/*printf(" Not a BOM file\n");*/
+    }
+  
     /* If line starts with a # it is a comment */
-    if (line[0] == '#') {
+    if (line[noffset] == '#') {
       numlinepar++;
       fputs(line,stdout);
       fputs(line,ficparo);
+      fputs(line,ficres);
       fputs(line,ficlog);
       continue;
     }else
@@ -10870,6 +10911,7 @@ int main(int argc, char *argv[])
       numlinepar++;
       fputs(line,stdout);
       fputs(line,ficparo);
+      fputs(line,ficres);
       fputs(line,ficlog);
       continue;
     }else
@@ -10892,6 +10934,7 @@ int main(int argc, char *argv[])
       numlinepar++;
       fputs(line,stdout);
       fputs(line,ficparo);
+      fputs(line,ficres);
       fputs(line,ficlog);
       continue;
     }else
@@ -11153,16 +11196,6 @@ Please run with mle=-1 to get a correct
     
     fflush(ficlog);
     
-    /*-------- Rewriting parameter file ----------*/
-    strcpy(rfileres,"r");    /* "Rparameterfile */
-    strcat(rfileres,optionfilefiname);    /* Parameter file first name*/
-    strcat(rfileres,".");    /* */
-    strcat(rfileres,optionfilext);    /* Other files have txt extension */
-    if((ficres =fopen(rfileres,"w"))==NULL) {
-      printf("Problem writing new parameter file: %s\n", rfileres);goto end;
-      fprintf(ficlog,"Problem writing new parameter file: %s\n", rfileres);goto end;
-    }
-    fprintf(ficres,"#%s\n",version);
   }    /* End of mle != -3 */
   
   /*  Main data