Method writeModelNodes

SRC/recorder/MPCORecorder.cpp:5022–5121

}

int MPCORecorder::writeModelNodes()
{
#ifdef MPCO_TIMING
	mpco::Timer timer("writeModelNodes"); timer.start();
#endif // MPCO_TIMING
	/*
	error flags
	*/
	int retval = 0;
	herr_t status = 0;
	/*
	get all nodes except pressure nodes
	*/
	m_data->nodes.clear();
	ID pressure_node_tags(0, m_data->info.domain->getNumPCs());
	Pressure_ConstraintIter& pressure_constraint_iter = m_data->info.domain->getPCs();
	Pressure_Constraint* pressure_constraint = 0;
	while ((pressure_constraint = pressure_constraint_iter()) != 0) {
		Node* pressure_node = pressure_constraint->getPressureNode();
		if (pressure_node) {
			pressure_node_tags.insert(pressure_node->getTag());
		}
	}
	if (m_data->has_region) { // only nodes in the requested regions
		for (std::vector<int>::const_iterator id_iter = m_data->node_set.begin(); id_iter != m_data->node_set.end(); ++id_iter) {
			int node_id = *id_iter;
			if (pressure_node_tags.getLocationOrdered(node_id) < 0) {
				Node *current_node = m_data->info.domain->getNode(node_id);
				if (current_node != 0) {
					m_data->nodes.push_back(current_node);
				}
			}
		}
	}
	else { // all nodes
		NodeIter& node_iter = m_data->info.domain->getNodes();
		Node* current_node = 0;
		while ((current_node = node_iter()) != 0) {
			int nd = current_node->getTag();
			if (pressure_node_tags.getLocationOrdered(nd) < 0) {
				m_data->nodes.push_back(current_node);
			}
		}
	}

	/*
	quick return
	*/
	size_t num_nodes = m_data->nodes.size();
	if (num_nodes == 0) {
		opserr << "MPCORecorder Error: no nodes to write\n";
		return -1;
	}
	/*
	allocate data for id and coordinates datasets
	*/
	std::vector<int> buffer_id(num_nodes);
	std::vector<double> buffer_coord(num_nodes*m_data->info.num_dimensions);