Method mapElements

SRC/recorder/MPCORecorder.cpp:3699–3831

			{}

			void mapElements(Domain *d, bool has_region, const std::vector<int> &subset) {
				/*
				utilties
				*/
				auto lam_get_num_ext_nodes = [](Element* elem) {
					switch (elem->getClassTag()) {
					case ELE_TAG_SFI_MVLEM_3D: return 4;
					default: return elem->getNumExternalNodes();
					}
				};
				/*
				clear previous mappings
				*/
				registered_custom_rules.clear();
				items.clear();
				/*
				quick return
				*/
				if (d == 0) return;
				/*
				auxiliary map to register custom rules
				*/
				std::map<ElementIntegrationRule, int> aux_map_custom_rules;
				/*
				loop over all elements in the domain
				*/
				size_t subset_elem_counter(0);
				ElementIter* element_iter = &(d->getElements());
				Element* current_element = 0;
				//while ((current_element = (*element_iter)()) != 0) {
				while (true) {
					/*
					get next element
					*/
					if (has_region) {
						if (subset_elem_counter == subset.size())
							break;
						current_element = d->getElement(subset[subset_elem_counter++]);
						if (current_element == 0)
							continue; // skip null and go to next iteration
					}
					else {
						current_element = (*element_iter)();
						if (current_element == 0) 
							break; 
					}
					/*
					get class tag, geometry and integration rule type
					*/
					int elem_type = current_element->getClassTag();
					/*
					skip element classes that we don't want to record
					*/
					if (elem_type == ELE_TAG_Subdomain ||
						elem_type == ELE_TAG_ASDEmbeddedNodeElement)
					{
						continue;
					}