From: N. Brouard Date: Thu, 16 May 2019 14:54:33 +0000 (+0000) Subject: Summary: There was some wrong lines added X-Git-Tag: imach-099s7~187 X-Git-Url: https://henry.ined.fr/git/?a=commitdiff_plain;h=daf1d34a8f9a321555f16ae53d7c7ccc93654d83;p=.git Summary: There was some wrong lines added --- diff --git a/src/imach.c b/src/imach.c index 98bf578..aa18cb1 100644 --- a/src/imach.c +++ b/src/imach.c @@ -1,6 +1,9 @@ /* $Id$ $State$ $Log$ + Revision 1.293 2019/05/09 15:17:34 brouard + *** empty log message *** + Revision 1.292 2019/05/09 14:17:20 brouard Summary: Some updates @@ -2966,24 +2969,8 @@ double **pmij(double **ps, double *cov, int ncovmodel, double *x, int nlstate ) 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]); */