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

Fixed bug in analysis.cpp

Browse Source
This commit is contained in:
James E McClure 2015-08-10 11:13:28 -04:00
parent d1fa260c6e
commit 71c0b3cdc3
2 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
analysis/analysis.cpp
View File

@ -515,7 +515,7 @@ static int LocalToGlobalIDs( int nx, int ny, int nz, const RankInfoStruct& rank_
for (size_t i=0; i<LocalIDs.length(); i++) {
if ( LocalIDs(i)>=0 ) {
local.insert(LocalIDs(i));
if ( LocalIDs(i)!=IDs(i) )
if ( LocalIDs(i)!=IDs(i) && IDs(i)>= 0)
map[LocalIDs(i)].remote_ids.insert(IDs(i));
}
}

3
tests/lbpm_segmented_pp.cpp
View File

@ -224,10 +224,11 @@ int main(int argc, char **argv)
Dm.CommunicateMeshHalo(Averages.SDn);
Dm.CommunicateMeshHalo(Averages.SDs);
if (rank==0) printf("computing phase components \n");
int timestep=5;
Averages.Initialize();
if (rank==0) printf("computing phase components \n");
Averages.ComponentAverages();
if (rank==0) printf("sorting phase components \n");
Averages.SortBlobs();
Averages.PrintComponents(timestep);