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

Using Dm.id to identify solid instead of Signed distance function

Browse Source
This commit is contained in:
James E McClure 2015-06-26 15:27:23 -04:00
parent 75b1e35564
commit 0136c739af
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
common/TwoPhase.h
View File

@ -468,8 +468,8 @@ void TwoPhase::ComputeLocal(){
//...........................................................................
// Compute volume averages
for (int p=0;p<8;p++){
if ( SDs(i+cube[p][0],j+cube[p][1],k+cube[p][2]) > 0 ){
int np_cube = i+cube[p][0] + (j+cube[p][1])*Nx + (k+cube[p][2])*Nx*Ny;
if ( Dm.id[n] != 0 ){
// 1-D index for this cube corner
n = i+cube[p][0] + (j+cube[p][1])*Nx + (k+cube[p][2])*Nx*Ny;
// compute the norm of the gradient of the phase indicator field
@ -612,8 +612,8 @@ void TwoPhase::ComputeLocalBlob(){
//...........................................................................
// Compute volume averages
for (int p=0;p<8;p++){
if ( SDs(i+cube[p][0],j+cube[p][1],k+cube[p][2]) > 0 ){
int np_cube = i+cube[p][0] + (j+cube[p][1])*Nx + (k+cube[p][2])*Nx*Ny;
if ( Dm.id[n] != 0 ){
// 1-D index for this cube corner
n = i+cube[p][0] + (j+cube[p][1])*Nx + (k+cube[p][2])*Nx*Ny;
// compute the norm of the gradient of the phase indicator field