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

fixed memory uct

Browse Source
This commit is contained in:
James E McClure 2018-06-07 11:42:56 -04:00
parent 15cc39709d
commit afc2767416
1 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
tests/lbpm_uCT_pp.cpp
View File

@ -339,7 +339,8 @@ int main(int argc, char **argv)
} }
*/ */
// Write the results to visit {
// Write the results to visit
if (rank==0) printf("Setting up visualization structure \n"); if (rank==0) printf("Setting up visualization structure \n");
// std::vector<IO::MeshDataStruct> meshData(N_levels); // std::vector<IO::MeshDataStruct> meshData(N_levels);
std::vector<IO::MeshDataStruct> meshData(1); std::vector<IO::MeshDataStruct> meshData(1);
@ -442,7 +443,7 @@ int main(int argc, char **argv)
// Write visulization data // Write visulization data
IO::writeData( 0, meshData, comm ); IO::writeData( 0, meshData, comm );
if (rank==0) printf("Finished. \n"); if (rank==0) printf("Finished. \n");
}
// Compute the Minkowski functionals // Compute the Minkowski functionals
MPI_Barrier(comm); MPI_Barrier(comm);
std::shared_ptr<Minkowski> Averages(new Minkowski(Dm[0])); std::shared_ptr<Minkowski> Averages(new Minkowski(Dm[0]));
@ -454,9 +455,9 @@ int main(int argc, char **argv)
double distance=Dist[0](i,j,k); double distance=Dist[0](i,j,k);
Averages->SDn(i,j,k) = distance; Averages->SDn(i,j,k) = distance;
if (distance > 0) if (distance > 0)
Averages->id[n] = 0; Averages->Dm->id[n] = 0;
else else
Averages->id[n] = 1; Averages->Dm->id[n] = 1;
} }
} }
} }