Method writeSections

SRC/recorder/MPCORecorder.cpp:5355–5853

}

int MPCORecorder::writeSections()
{
#ifdef MPCO_TIMING
	mpco::Timer timer("writeSections"); timer.start();
#endif // MPCO_TIMING
	/*
	error flags
	*/
	int retval = 0;
	herr_t status = 0;
	/*
	auxiliary variables
	*/
	std::map<int, mpco::element::SectionAssignment> sec_assignments;
	std::map<mpco::element::FiberSectionData, mpco::element::SectionAssignment> aux_sec_assignments;
	std::map<int, std::string> sec_id_and_names;
	std::map<int, std::string> mat_id_and_names;
	/*
	an auxiliary map to store, for each element the number of gauss points (with sections) and the number of fibers.
	key = element tag, value = vector with size equal to the number of gauss points, and each item is the number of fibers
	for that gauss point
	*/
	m_data->elem_ngauss_nfiber_info.clear();
	/*
	create section assignments group
	*/
	std::stringstream ss_gp_section_assignments_dir;
	ss_gp_section_assignments_dir << "MODEL_STAGE[" << m_data->info.current_model_stage_id << "]/MODEL/SECTION_ASSIGNMENTS";
	std::string gp_section_assignments_dir = ss_gp_section_assignments_dir.str();
	hid_t h_gp_section_assignments = h5::group::create(m_data->info.h_file_id, gp_section_assignments_dir.c_str(), H5P_DEFAULT, m_data->info.h_group_proplist, H5P_DEFAULT);
	/*
	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;
				/*
				loop over all alements 
				*/
				for (std::vector<Element*>::iterator it4 = elem_by_custom_rule.items.begin(); 
					it4 != elem_by_custom_rule.items.end(); ++it4) {
					Element *elem = *it4;
					/*
					prepare the request
					*/
					std::string request1 = utils::shell::isShellElementTag(elem->getClassTag()) ? "material" : "section";