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

Debugging common/TwoPhase.cpp for labeling error

Browse Source
This commit is contained in:
James E McClure 2015-09-02 12:56:35 -04:00
parent 00e2b876db
commit 53fe2b72be
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
common/TwoPhase.cpp
View File

@ -734,9 +734,9 @@ void TwoPhase::ComponentAverages()
geomavg_MarchingCubes(SDn,fluid_isovalue,i,j,k,nw_pts,n_nw_pts,nw_tris,n_nw_tris);
euler = geomavg_EulerCharacteristic(nw_pts,nw_tris,n_nw_pts,n_nw_tris,i,j,k);
if (LabelNWP < 0){
printf("Error: rank=%i \n"Dm.rank);
printf("mesh =%i,%i,%i \n",i+iprocx*(Nx-2),j+jprocx*(Ny-2),k+kprocx*(Nz-2));
printf("Label=%i \n"LabelNWP);
printf("Error: rank=%i \n",Dm.rank);
printf("mesh =%i,%i,%i \n",i+Dm.iproc*(Nx-2),j+Dm.jproc*(Ny-2),k+Dm.kproc*(Nz-2));
printf("Label=%i \n",LabelNWP);
}
else ComponentAverages_NWP(EULER,LabelNWP) += euler;
}