From 378e2a8f3f1fd9a12631306454723697c545e6a3 Mon Sep 17 00:00:00 2001 From: daniel-koehn Date: Fri, 12 Aug 2016 03:49:40 +0200 Subject: [PATCH] Fixed some PCG related bugs --- src/PCG.c | 2 +- src/PSV/FWI_PSV.c | 2 ++ 2 files changed, 3 insertions(+), 1 deletion(-) diff --git a/src/PCG.c b/src/PCG.c index 9ad824e..02900a7 100644 --- 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 index 70d7680..90ee5fb 100644 --- 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); }