--- imach/src/imach.c	2019/05/09 15:17:34	1.293
+++ imach/src/imach.c	2019/05/16 14:54:33	1.294
@@ -1,6 +1,9 @@
-/* $Id: imach.c,v 1.293 2019/05/09 15:17:34 brouard Exp $
+/* $Id: imach.c,v 1.294 2019/05/16 14:54:33 brouard Exp $
   $State: Exp $
   $Log: imach.c,v $
+  Revision 1.294  2019/05/16 14:54:33  brouard
+  Summary: There was some wrong lines added
+
   Revision 1.293  2019/05/09 15:17:34  brouard
   *** empty log message ***
 
@@ -1087,12 +1090,12 @@ typedef struct {
 #define ODIRSEPARATOR '\\'
 #endif
 
-/* $Id: imach.c,v 1.293 2019/05/09 15:17:34 brouard Exp $ */
+/* $Id: imach.c,v 1.294 2019/05/16 14:54:33 brouard Exp $ */
 /* $State: Exp $ */
 #include "version.h"
 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 fullversion[]="$Revision: 1.293 $ $Date: 2019/05/09 15:17:34 $"; 
+char fullversion[]="$Revision: 1.294 $ $Date: 2019/05/16 14:54:33 $"; 
 char strstart[80];
 char optionfilext[10], optionfilefiname[FILENAMELENGTH];
 int erreur=0, nberr=0, nbwarn=0; /* Error number, number of errors number of warnings  */
@@ -2969,24 +2972,8 @@ double **pmij(double **ps, double *cov,
       ps[ii][ii]=1;
     }
   }
-  /* Added for backcast */ /* Transposed matrix too */
-  for(jj=1; jj<= nlstate+ndeath; jj++){
-    s1=0.;
-    for(ii=1; ii<= nlstate+ndeath; ii++){
-      s1+=ps[ii][jj];
-    }
-    for(ii=1; ii<= nlstate; ii++){
-      ps[ii][jj]=ps[ii][jj]/s1;
-    }
-  }
-  /* Transposition */
-  for(jj=1; jj<= nlstate+ndeath; jj++){
-    for(ii=jj; ii<= nlstate+ndeath; ii++){
-      s1=ps[ii][jj];
-      ps[ii][jj]=ps[jj][ii];
-      ps[jj][ii]=s1;
-    }
-  }
+
+
   /* for(ii=1; ii<= nlstate+ndeath; ii++){ */
   /*   for(jj=1; jj<= nlstate+ndeath; jj++){ */
   /* 	printf(" pmij  ps[%d][%d]=%lf ",ii,jj,ps[ii][jj]); */