Method writeModelLocalAxes

SRC/recorder/MPCORecorder.cpp:5223–5353

}

int MPCORecorder::writeModelLocalAxes()
{
#ifdef MPCO_TIMING
	mpco::Timer timer("writeModelLocalAxes"); timer.start();
#endif // MPCO_TIMING
#ifdef MPCO_WRITE_LOC_AX_IS_VERBOSE
	std::stringstream ss;
	ss << "local axes:\n";
#endif // MPCO_WRITE_LOC_AX_IS_VERBOSE
	/*
	error flags
	*/
	int retval = 0;
	herr_t status = 0;
	/*
	allocate vectors for ID and local axes, only for elements that provide a response for local axes
	*/
	std::vector<int> buffer_el_id;
	std::vector<utils::locax::quaternion> buffer_el_data;
	/*
	element request
	*/
	int argc = 1;
	const char **argv = new const char*[argc];
	argv[0] = "localAxes";
	Vector vx(3);
	Vector vy(3);
	Vector vz(3);
	/*
	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;
					/*
					get the element response for the ith section and the ith fiber
					*/
					mpco::element::OutputDescriptor eo_descriptor;
					mpco::element::OutputDescriptorStream eo_stream(&eo_descriptor);
					Response *eo_response = elem->setResponse(argv, argc, eo_stream);