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();
}
}
