OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

debugging

Browse Source
This commit is contained in:
James E McClure
2015-06-03 08:41:52 -04:00
parent 8ba62c09e7
commit 34ea196569
1 changed files with 4 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
tests/BlobAnalyzeParallel.cpp
View File

@@ -230,9 +230,9 @@ int main(int argc, char **argv)
vz = f5-f6+f11-f12-f13+f14+f15-f16-f17+f18;
// Assign array components needed for averaging
Averages.SDs(i,j,k)=SignDist(i,j,k);
Averages.Phase(i,j,k)=(da-db)/(da+db);
Averages.Phase_tplus(i,j,k)=(da-db)/(da+db);
Averages.Phase_tminus(i,j,k)=(da-db)/(da+db);
Averages.Phase(i,j,k)=Phase(i,j,k);//(da-db)/(da+db);
Averages.Phase_tplus(i,j,k)=Phase(i,j,k);//(da-db)/(da+db);
Averages.Phase_tminus(i,j,k)=Phase(i,j,k);//(da-db)/(da+db);
Averages.Press(i,j,k)=press;
Averages.Vel_x(i,j,k)=vx;
Averages.Vel_y(i,j,k)=vy;
@@ -255,6 +255,7 @@ int main(int argc, char **argv)
nblobs = ComputeGlobalBlobIDs(nx,ny,nz,rank_info,
Averages.Phase,Averages.SDs,vF,vS,Averages.BlobLabel);
if ( rank==0 ) { printf("Identified %i blobs\n",nblobs); }
for (int k=0; k<nz+2; k++){