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

Debugging labeling error

Browse Source
This commit is contained in:
James E McClure 2015-09-02 12:13:44 -04:00
parent 98c414d1cb
commit 00e2b876db
1 changed files with 6 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
common/TwoPhase.cpp
View File

@ -733,7 +733,12 @@ void TwoPhase::ComponentAverages()
n_nw_pts=n_nw_tris=0;
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);
ComponentAverages_NWP(EULER,LabelNWP) += euler;
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);
}
else ComponentAverages_NWP(EULER,LabelNWP) += euler;
}
}
}