===================================================================
RCS file: /home/cvs/OpenXM/src/asir-contrib/packages/src/tk_jack.rr,v
retrieving revision 1.1
retrieving revision 1.2
diff -u -p -r1.1 -r1.2
--- OpenXM/src/asir-contrib/packages/src/tk_jack.rr	2012/09/18 11:41:18	1.1
+++ OpenXM/src/asir-contrib/packages/src/tk_jack.rr	2014/03/10 03:04:06	1.2
@@ -1,5 +1,5 @@
 /*
-$OpenXM$
+$OpenXM: OpenXM/src/asir-contrib/packages/src/tk_jack.rr,v 1.1 2012/09/18 11:41:18 takayama Exp $
 Koev-Edelman for higher order derivatives.
 Series solution obtained by an interpolation and Koev-Edelman evaluation.
 Previous version: misc-2011/10/wishart/Prog/mh.rr
@@ -60,8 +60,8 @@ localf checkJack1;
 localf mhdiff;
 localf checkJack2;
 localf listJack;
+localf showStatic;
 
-
 static Alpha$
 Alpha=2$
 static Darray$
@@ -354,8 +354,11 @@ static M_qk$
 M_qk=0$
 def mh_t(A,B,N,M) {
   /* extern M_qk; */
+  printf("P_pmn=%a\n",P_pmn);
   genJack(M,N);
+  printf(" P_pmn=%a\n",P_pmn);
   M_qk = newvect(P_pmn);
+  printf("P_pmn=%a, M_qk=%a\n",P_pmn,M_qk);
   F = 0;
   for (K=0; K<P_pmn; K++) {
      Kap = Parray[K];
@@ -597,6 +600,11 @@ static M_n$  /* digits */
 static M_m$  /* | | <= M_m */
 static M_plist$ /* to save a list. asir only */
 static M_pExec$
+M_kap=0$
+M_n=0$
+M_m=0$
+M_plist=0$
+M_pExec=0$
 
 def pListPartition(M,N) {
   /* initialize */
@@ -651,6 +659,10 @@ static HS_mu$
 static HS_n$
 static HS_plist$
 static HS_hsExec$
+HS_mu=0$
+HS_n=0$
+HS_plist=0$
+HS_hsExec=0$
 def pListHS(Kap,N) {
   HS_n = N;
   HS_mu = newvect(N+1);
@@ -674,7 +686,6 @@ def hsExec_0() {
   printf("hsExec: %a\n",HS_mu);
 }
 
-HS_hsExec = hsExec_0$
 
 /*
   pListHS([0,4,2,1],3);
@@ -685,6 +696,8 @@ HS_hsExec = hsExec_0$
 */
 static P_pki$
 static P_pmn$
+P_pki=0$
+P_pmn=0$
 def pmn(M,N) {
   Min_m_n = (M>N?N:M);
   P_pki=newmat(Min_m_n+1,M+1);
@@ -699,6 +712,7 @@ def pmn(M,N) {
     }
   }
   P_pmn=S;
+  printf("pmn(%a,%a), P_pmn=%a\n",M,N,P_pmn);
   return(S);
 }
 
@@ -738,6 +752,9 @@ def genDarray2(M,N) {
 static M_beta$   /* M_beta[0][*] value of beta_{kappa,mu}, M_beta[1][*] N_mu */
 static M_beta_pt$
 static M_beta_kap$
+M_beta=0$
+M_beta_pt=0$
+M_beta_kap=0$
 def isHStrip(Kap,Nu) {
   N1 = plength(Kap); N2 = plength(Nu);
   if (N2 > N1) return(0);
@@ -800,11 +817,11 @@ def genBeta(Kap) {
   extern M_beta_kap;
   extern HS_hsExec;
   extern P_pmn; */
-  printf("Kappa=%a, M_pmn=%a\n",Kap,P_pmn);
+  printf("Kappa=%a, P_pmn=%a\n",Kap,P_pmn);
   M_beta = newmat(2,P_pmn+1);
   M_beta_pt = 0;
-  for (I=0; I<=M_pmn; I++) for (J=0; J<2; J++) M_beta[J][I] = n;
   N = plength(Kap);
+  for (I=0; I<=P_pmn; I++) for (J=0; J<2; J++) M_beta[J][I] = N;
   HS_hsExec = hsExec_beta;
   M_beta_kap = Kap;
   pListHS(cons(0,Kap),N);
@@ -866,6 +883,9 @@ def psublen(Kap,Mu) {
 static M_jack$ 
 static M_df$ /* Compute differentials? */
 static M_2n$ /* 2^N */
+M_jack=0$
+M_df=0$
+M_2n=0$
 
 def genJack(M,N) {
   /* extern M_jack;
@@ -1023,5 +1043,10 @@ def listJack(M,N) {
   }
   printf("\\end{eqnarray*}\n");
 }
+
+def showStatic() {
+  return [M_qk,M_kap,M_jack,Parray,P_pmn];
+}
+
 endmodule$
 end$
\ No newline at end of file