Method writeModelElements

SRC/recorder/MPCORecorder.cpp:5123–5221

}

int MPCORecorder::writeModelElements()
{
#ifdef MPCO_TIMING
	mpco::Timer timer("writeModelElements"); timer.start();
#endif // MPCO_TIMING
	/*
	error flags
	*/
	int retval = 0;
	herr_t status = 0;
	/*
	create element group
	*/
	std::stringstream ss_gp_elements_dir;
	ss_gp_elements_dir << "MODEL_STAGE[" << m_data->info.current_model_stage_id << "]/MODEL/ELEMENTS";
	std::string gp_elements_dir = ss_gp_elements_dir.str();
	hid_t h_gp_elements = h5::group::create(m_data->info.h_file_id, gp_elements_dir.c_str(), H5P_DEFAULT, m_data->info.h_group_proplist, H5P_DEFAULT);
	/*
	map all elements
	*/
	m_data->elements.mapElements(m_data->info.domain, m_data->has_region, m_data->elem_set);
	/*
	loop over all element group-by-tag
	*/
	for (mpco::element::ElementCollection::submap_type::iterator it1 = m_data->elements.items.begin();
		it1 != m_data->elements.items.end(); ++it1) {
		mpco::element::ElementWithSameClassTagCollection &elem_by_tag = it1->second;
		/*
		loop over all element group-by-integration rule
		*/
		for (mpco::element::ElementWithSameClassTagCollection::submap_type::iterator it2 = elem_by_tag.items.begin();
			it2 != elem_by_tag.items.end(); ++it2) {
			mpco::element::ElementWithSameIntRuleCollection &elem_by_rule = it2->second;
			/*
			loop over all element group-by-custom(or default) integration rule
			*/
			for (mpco::element::ElementWithSameIntRuleCollection::submap_type::iterator it3 = elem_by_rule.items.begin();
				it3 != elem_by_rule.items.end(); ++it3) {
				mpco::element::ElementWithSameCustomIntRuleCollection &elem_by_custom_rule = it3->second;
				/*
				create a dataset with all elements in this group
				*/
				if (elem_by_custom_rule.items.size() > 0) {
					/*
					create a name for this dataset using the following format <class_tag>-<class_name>[<integration_rule>:<custom_rule_index>]
					*/
					std::stringstream ss_dset_name;
					ss_dset_name << elem_by_tag.class_tag << "-" << elem_by_tag.class_name
						<< "[" << elem_by_rule.int_rule_type << ":"
						<< elem_by_custom_rule.custom_int_rule_index << "]";
					elem_by_custom_rule.name = ss_dset_name.str();
					/*
					create the dataset and write data to it.
					the first column is the element id, the remaining columns are nodal connectivity
					*/
					std::vector<int> buffer(elem_by_custom_rule.items.size()*(1 + elem_by_tag.num_nodes));
					size_t offset = 0;
					for (mpco::element::ElementWithSameCustomIntRuleCollection::collection_type::iterator it4 = elem_by_custom_rule.items.begin();