/* $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
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]); */