Method writeSets

SRC/recorder/MPCORecorder.cpp:5855–5969

}

int MPCORecorder::writeSets()
{
#ifdef MPCO_TIMING
	mpco::Timer timer("writeSets"); timer.start();
#endif // MPCO_TIMING
	/*
	error flags
	*/
	int retval = 0;
	herr_t status = 0;
	/*
	get all regions.
	note if the user is saving results not on the whole model, but only on a subset,
	the of each reagion we write only the nodes/elements contained in the subset
	and if the region (so modified) is empty, we skip it.
	*/
	std::set<int> nset;
	std::set<int> eset;
	if (m_data->has_region) {
		// create ordered sets, search will be faster
		for (size_t i = 0; i < m_data->node_set.size(); i++)
			nset.insert(m_data->node_set[i]);
		for (size_t i = 0; i < m_data->elem_set.size(); i++)
			eset.insert(m_data->elem_set[i]);
	}
	ID tags;
	m_data->info.domain->getRegionTags(tags);
	if (tags.Size() > 0) {
		/*
		write the sets group
		*/
		std::stringstream ss_gp_sets;
		ss_gp_sets << "MODEL_STAGE[" << m_data->info.current_model_stage_id << "]/MODEL/SETS";
		std::string gp_sets = ss_gp_sets.str();
		hid_t h_gp_sets = h5::group::create(m_data->info.h_file_id, gp_sets.c_str(), H5P_DEFAULT, m_data->info.h_group_proplist, H5P_DEFAULT);
		/*
		for each region
		*/
		for (int i = 0; i < tags.Size(); i++) {
			/*
			get current region
			*/
			MeshRegion *region = m_data->info.domain->getRegion(tags(i));
			if (region == NULL) continue;
			/*
			write the current set group
			*/
			std::stringstream ss_gp_curr_set;
			ss_gp_curr_set << "SET_" << tags(i);
			std::string gp_curr_set = ss_gp_curr_set.str();
			hid_t h_gp_curr_set = h5::group::create(h_gp_sets, gp_curr_set.c_str(), H5P_DEFAULT, m_data->info.h_group_proplist, H5P_DEFAULT);
			/*
			write the current set NODES dataset
			*/
			{
				const ID &nodes = region->getNodes();
				if (nodes.Size() > 0) {
					std::vector<int> aux;