imach/src/imach.c - diff

Return to imach.c CVS log

Up to [Local Repository] / imach / src

-version 1.293, 2019/05/09 15:17:34
+version 1.294, 2019/05/16 14:54:33
  Line 1
  /* $Id$
    $State$
    $Log$
+   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 ***
- Line 2969  double **pmij(double **ps, double *cov,
+ Line 2972  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]); */