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

tweak labels

Browse Source
This commit is contained in:
James E McClure 2018-12-10 09:02:01 -05:00
parent c3a1c30b4c
commit d4a8a3e039
2 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
models/ColorModel.cpp
View File

@ -208,7 +208,7 @@ void ScaLBL_ColorModel::AssignComponentLabels(double *phase)
VALUE=LabelList[idx]; VALUE=LabelList[idx];
AFFINITY=AffinityList[idx]; AFFINITY=AffinityList[idx];
double volume_fraction = double(label_count_global[idx])/double((Nx-2)*(Ny-2)*(Nz-2)*nprocz); double volume_fraction = double(label_count_global[idx])/double((Nx-2)*(Ny-2)*(Nz-2)*nprocz);
printf(" label=%i, affinity=%f, volume fraction==%f\n",int(VALUE),AFFINITY,volume_fraction); printf(" label=%i, affinity=%f, volume fraction=%f\n",int(VALUE),AFFINITY,volume_fraction);
} }
} }

5
tests/lbpm_morphopen_pp.cpp
View File

@ -167,9 +167,10 @@ int main(int argc, char **argv)
for (int j=0;j<ny;j++){ for (int j=0;j<ny;j++){
for (int i=0;i<nx;i++){ for (int i=0;i<nx;i++){
int n = k*nx*ny+j*nx+i; int n = k*nx*ny+j*nx+i;
char VALUE=id[n]; char LOCVAL = id[n];
for (unsigned int idx=0; idx < NLABELS; idx++){ for (unsigned int idx=0; idx < NLABELS; idx++){
if (VALUE == LabelList[idx]){ char VALUE=LabelList[idx];
if (LOCVAL == VALUE){
idx = NLABELS; idx = NLABELS;
char NEWVALUE=HistoryLabels[idx]; char NEWVALUE=HistoryLabels[idx];
if (SignDist(i,j,k) < 1.0) id[n] = NEWVALUE; if (SignDist(i,j,k) < 1.0) id[n] = NEWVALUE;