void
NodalFloodCount::finalize()
{
// Exchange data in parallel
pack(_packed_data);
_communicator.allgather(_packed_data, false);
unpack(_packed_data);
mergeSets();
// Populate _bubble_maps and _var_index_maps
updateFieldInfo();
// Update the region offsets so we can get unique bubble numbers in multimap mode
updateRegionOffsets();
// Calculate and out output bubble volume data
if (_pars.isParamValid("bubble_volume_file"))
{
calculateBubbleVolumes();
std::vector<Real> data; data.reserve(_all_bubble_volumes.size() + 2);
data.push_back(_fe_problem.timeStep());
data.push_back(_fe_problem.time());
data.insert(data.end(), _all_bubble_volumes.begin(), _all_bubble_volumes.end());
writeCSVFile(getParam<FileName>("bubble_volume_file"), data);
}
// Calculate memory usage
if (_track_memory)
{
_bytes_used += calculateUsage();
_communicator.sum(_bytes_used);
formatBytesUsed();
}
}
void
FeatureFloodCount::finalize()
{
// Exchange data in parallel
pack(_packed_data);
_communicator.allgather(_packed_data, false);
unpack(_packed_data);
mergeSets(true);
// Populate _bubble_maps and _var_index_maps
updateFieldInfo();
// Update the region offsets so we can get unique bubble numbers in multimap mode
updateRegionOffsets();
// Calculate and out output bubble volume data
if (_pars.isParamValid("bubble_volume_file"))
{
calculateBubbleVolumes();
std::vector<Real> data;
data.reserve(_all_bubble_volumes.size() + _total_volume_intersecting_boundary.size() + 2);
// Insert the current timestep and the simulation time into the data vector
data.push_back(_fe_problem.timeStep());
data.push_back(_fe_problem.time());
// Insert the (sorted) bubble volumes into the data vector
data.insert(data.end(), _all_bubble_volumes.begin(), _all_bubble_volumes.end());
// If we are computing the boundary-intersecting volumes, insert
// those numbers into the normalized boundary-intersecting bubble
// volumes into the data vector.
if (_compute_boundary_intersecting_volume)
data.insert(data.end(), _total_volume_intersecting_boundary.begin(), _total_volume_intersecting_boundary.end());
// Finally, write the file
writeCSVFile(getParam<FileName>("bubble_volume_file"), data);
}
// Calculate memory usage
if (_track_memory)
{
_bytes_used += calculateUsage();
_communicator.sum(_bytes_used);
formatBytesUsed();
}
}
void
GrainTracker::finalize()
{
// Don't track grains if the current simulation step is before the specified tracking step
if (_t_step < _tracking_step)
return;
Moose::perf_log.push("finalize()","GrainTracker");
// Exchange data in parallel
pack(_packed_data);
_communicator.allgather(_packed_data, false);
unpack(_packed_data);
mergeSets(false);
Moose::perf_log.push("buildspheres()","GrainTracker");
buildBoundingSpheres(); // Build bounding sphere information
Moose::perf_log.pop("buildspheres()","GrainTracker");
// Now merge sets again but this time we'll add periodic neighbor information
mergeSets(true);
Moose::perf_log.push("trackGrains()","GrainTracker");
trackGrains();
Moose::perf_log.pop("trackGrains()","GrainTracker");
Moose::perf_log.push("remapGrains()","GrainTracker");
if (_remap)
remapGrains();
Moose::perf_log.pop("remapGrains()","GrainTracker");
updateFieldInfo();
Moose::perf_log.pop("finalize()","GrainTracker");
// Calculate and out output bubble volume data
if (_pars.isParamValid("bubble_volume_file"))
{
calculateBubbleVolumes();
std::vector<Real> data; data.reserve(_all_bubble_volumes.size() + 2);
data.push_back(_fe_problem.timeStep());
data.push_back(_fe_problem.time());
data.insert(data.end(), _all_bubble_volumes.begin(), _all_bubble_volumes.end());
writeCSVFile(getParam<FileName>("bubble_volume_file"), data);
}
if (_compute_op_maps)
{
for (std::map<unsigned int, UniqueGrain *>::const_iterator grain_it = _unique_grains.begin();
grain_it != _unique_grains.end(); ++grain_it)
{
if (grain_it->second->status != INACTIVE)
{
std::set<dof_id_type>::const_iterator elem_it_end = grain_it->second->entities_ptr->end();
for (std::set<dof_id_type>::const_iterator elem_it = grain_it->second->entities_ptr->begin(); elem_it != elem_it_end; ++elem_it)
_elemental_data[*elem_it].push_back(std::make_pair(grain_it->first, grain_it->second->variable_idx));
}
}
}
// Calculate memory usage
if (_track_memory)
{
_bytes_used += calculateUsage();
_communicator.sum(_bytes_used);
formatBytesUsed();
}
}