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

Fixed communication bug in TwoPhase::ComponentAverages

Browse Source
This commit is contained in:
James E McClure 2015-07-13 13:07:52 -04:00
parent 233e36644c
commit 6a7fa484f7
2 changed files with 3 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
analysis/analysis.cpp
View File

@ -22,7 +22,7 @@ int ComputePhaseComponent(IntArray &ComponentLabel,
int Ny = PhaseID.size(1);
int Nz = PhaseID.size(2);
//ComponentLabel.resize(Nx,Ny,Nz);
ComponentLabel.resize(Nx,Ny,Nz);
int cubes_in_blob=0;
int nrecent = 1; // number of nodes added at most recent sweep

5
common/TwoPhase.h
View File

@ -553,8 +553,8 @@ void TwoPhase::ComponentAverages(){
int cube[8][3] = {{0,0,0},{1,0,0},{0,1,0},{1,1,0},{0,0,1},{1,0,1},{0,1,1},{1,1,1}};
NumberComponents_WP = ComputeGlobalPhaseComponent(Nx,Ny,Nz,Dm.rank_info,PhaseID,1,Label_WP);
NumberComponents_NWP = ComputeGlobalPhaseComponent(Nx,Ny,Nz,Dm.rank_info,PhaseID,1,Label_NWP);
NumberComponents_WP = ComputeGlobalPhaseComponent(Nx-2,Ny-2,Nz-2,Dm.rank_info,PhaseID,1,Label_WP);
NumberComponents_NWP = ComputeGlobalPhaseComponent(Nx-2,Ny-2,Nz-2,Dm.rank_info,PhaseID,1,Label_NWP);
ComponentAverages_WP.resize(BLOB_AVG_COUNT,NumberComponents_WP);
ComponentAverages_NWP.resize(BLOB_AVG_COUNT,NumberComponents_NWP);
@ -847,7 +847,6 @@ void TwoPhase::ComponentAverages(){
ComponentAverages_WP(VWNSX,b) = vawns(0);
ComponentAverages_WP(VWNSY,b) = vawns(1);
ComponentAverages_WP(VWNSZ,b) = vawns(2);
}
}
}