moved component averagign outside of simulation since it was taking a long time to complete for micromodel system

This commit is contained in:
James E McClure
2015-07-14 13:19:57 -04:00
parent 1dec46041e
commit dd005b00bb

View File

@@ -722,11 +722,11 @@ int main(int argc, char **argv)
Averages.ComputeLocal();
Averages.Reduce();
Averages.PrintAll(timestep);
Averages.Initialize();
/* Averages.Initialize();
Averages.ComponentAverages();
Averages.SortBlobs();
Averages.PrintComponents(timestep);
//....................................................................
*/ //....................................................................
}
if (timestep%RESTART_INTERVAL == 0){
@@ -884,6 +884,13 @@ int main(int argc, char **argv)
if (rank==0) printf("Lattice update rate (total)= %f MLUPS \n", MLUPS);
if (rank==0) printf("********************************************************\n");
//************************************************************************/
// Perform component averaging and write tcat averages
Averages.Initialize();
Averages.ComponentAverages();
Averages.SortBlobs();
Averages.PrintComponents(timestep);
//************************************************************************/
DeviceBarrier();
CopyToHost(Averages.Phase.get(),Phi,N*sizeof(double));