Fixed some PCG related bugs

daniel-koehn · Aug 12, 2016 · 378e2a8 · 378e2a8
1 parent 78fca1b
commit 378e2a8
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/src/PCG.c b/src/PCG.c
@@ -75,7 +75,7 @@ void PCG(float * PCG_new, float * PCG_old, float * PCG_dir, int PCG_class){
 	    if((beta[k]<0.0)&&(PCG_BETA==2)){beta[k] = 0.0;}
 
             if(MYID==0){
-	        printf("\n\n parameter class = %d \t beta = %e\n",MYID,k,beta[k]);
+	        printf("\n\n parameter class = %d \t beta = %e\n",k,beta[k]);
 	    } 
 
         }

diff --git a/src/PSV/FWI_PSV.c b/src/PSV/FWI_PSV.c
@@ -299,6 +299,8 @@ if(GRAD_METHOD==1){
   PCG_vec = PCG_class*NX*NY;  	 /* length of one PCG-parameter class */
 
   PCG_old  =  vector(1,PCG_vec);
+  PCG_new  =  vector(1,PCG_vec);
+  PCG_dir  =  vector(1,PCG_vec);
 
 }