From ba0c857ee73e809a1f7e87d4eeeb64b94ea997d7 Mon Sep 17 00:00:00 2001
From: "N. Brouard" <brouard@ined.fr>
Date: Tue, 2 Sep 2014 11:19:57 +0000
Subject: [PATCH] *** empty log message ***

---
 src/hpijx.h | 55 +++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 55 insertions(+)
 create mode 100644 src/hpijx.h

diff --git a/src/hpijx.h b/src/hpijx.h
new file mode 100644
index 0000000..e9560bd
--- /dev/null
+++ b/src/hpijx.h
@@ -0,0 +1,55 @@
+    /*------------- h Pij x at various ages ------------*/
+  
+    strcpy(filerespij,"pij");  strcat(filerespij,fileres);
+    if((ficrespij=fopen(filerespij,"w"))==NULL) {
+      printf("Problem with Pij resultfile: %s\n", filerespij);goto end;
+      fprintf(ficlog,"Problem with Pij resultfile: %s\n", filerespij);goto end;
+    }
+    printf("Computing pij: result on file '%s' \n", filerespij);
+    fprintf(ficlog,"Computing pij: result on file '%s' \n", filerespij);
+  
+    stepsize=(int) (stepm+YEARM-1)/YEARM;
+    /*if (stepm<=24) stepsize=2;*/
+
+    agelim=AGESUP;
+    hstepm=stepsize*YEARM; /* Every year of age */
+    hstepm=hstepm/stepm; /* Typically 2 years, = 2/6 months = 4 */ 
+
+    /* hstepm=1;   aff par mois*/
+    pstamp(ficrespij);
+    fprintf(ficrespij,"#****** h Pij x Probability to be in state j at age x+h being in i at x ");
+    i1=pow(2,cptcoveff);
+    for(cptcov=1,k=0;cptcov<=i1;cptcov++){
+      /*for(cptcod=1;cptcod<=ncodemax[cptcov];cptcod++){*/
+	k=k+1;
+	fprintf(ficrespij,"\n#****** ");
+	for(j=1;j<=cptcoveff;j++) 
+	  fprintf(ficrespij,"V%d=%d ",Tvaraff[j],nbcode[Tvaraff[j]][codtab[k][j]]);
+	fprintf(ficrespij,"******\n");
+	
+	for (agedeb=fage; agedeb>=bage; agedeb--){ /* If stepm=6 months */
+	  nhstepm=(int) rint((agelim-agedeb)*YEARM/stepm); /* Typically 20 years = 20*12/6=40 */ 
+	  nhstepm = nhstepm/hstepm; /* Typically 40/4=10 */
+
+	  /*	  nhstepm=nhstepm*YEARM; aff par mois*/
+
+	  p3mat=ma3x(1,nlstate+ndeath,1, nlstate+ndeath, 0,nhstepm);
+	  oldm=oldms;savm=savms;
+	  hpxij(p3mat,nhstepm,agedeb,hstepm,p,nlstate,stepm,oldm,savm, k);  
+	  fprintf(ficrespij,"# Cov Agex agex+h hpijx with i,j=");
+	  for(i=1; i<=nlstate;i++)
+	    for(j=1; j<=nlstate+ndeath;j++)
+	      fprintf(ficrespij," %1d-%1d",i,j);
+	  fprintf(ficrespij,"\n");
+	  for (h=0; h<=nhstepm; h++){
+	    fprintf(ficrespij,"%d %3.f %3.f",k,agedeb, agedeb+ h*hstepm/YEARM*stepm );
+	    for(i=1; i<=nlstate;i++)
+	      for(j=1; j<=nlstate+ndeath;j++)
+		fprintf(ficrespij," %.5f", p3mat[i][j][h]);
+	    fprintf(ficrespij,"\n");
+	  }
+	  free_ma3x(p3mat,1,nlstate+ndeath,1, nlstate+ndeath, 0,nhstepm);
+	  fprintf(ficrespij,"\n");
+	}
+      /*}*/
+    }
-- 
2.43.0