Method handle

SRC/analysis/handler/AutoConstraintHandler.cpp:310–522

}

int
AutoConstraintHandler::handle(const ID* nodesLast)
{
	// first check links exist to a Domain and an AnalysisModel object
	Domain* theDomain = this->getDomainPtr();
	AnalysisModel* theModel = this->getAnalysisModelPtr();
	Integrator* theIntegrator = this->getIntegratorPtr();
	if ((theDomain == 0) || (theModel == 0) || (theIntegrator == 0)) {
		opserr << "WARNING AutoConstraintHandler::handle() - ";
		opserr << " setLinks() has not been called\n";
		return -1;
	}

	// create a multimap (key = NodeID, value = SP_Constraint).
	// multimap allows for duplicate keys because we may have multiple SP
	// constraints on the same node.
	// in this way we have :
	// a) a fast access to the key (to tell whether a node is constrained or not)
	// b) and a range of SP_Constraints acting on that node
	std::unordered_multimap<int, SP_Constraint*> sp_map;
	{
		SP_ConstraintIter& theSPs = theDomain->getDomainAndLoadPatternSPs();
		SP_Constraint* theSP;
		while ((theSP = theSPs()) != 0) 
			sp_map.emplace(std::make_pair(theSP->getNodeTag(), theSP));
	}

	// create a DOF_Group for each Node and add it to the AnalysisModel
	// create a TrasformationDOF_Group for nodes constrained in SP_Constraints
	int countDOF = 0;
	{
		NodeIter& theNodes = theDomain->getNodes();
		Node* nodPtr;
		int numDofGrp = 0;
		while ((nodPtr = theNodes()) != 0) {

			// process this node
			DOF_Group* dofPtr = 0;
			int nodeTag = nodPtr->getTag();

			// if the node is constrained in an SP constraint
			// handle it with the transformation method
			auto it_range = sp_map.equal_range(nodeTag);
			if (it_range.first != it_range.second) {

				// this node is constrained in 1 or more SP_Constraints.
				// handle it with a TransformationDOF_Group.
				// note: the TransformationConstraintHandler pointer is useless in TDOF...
				TransformationDOF_Group* tDofPtr =
					new TransformationDOF_Group(numDofGrp++, nodPtr, nullptr);
				dofPtr = tDofPtr;

				// count and add all SP_Constraints.
				// if verbose, keep track of conflicting constraints
				int numSPs = 0;
				for (auto it = it_range.first; it != it_range.second; it++) {
					tDofPtr->addSP_Constraint(*(it->second));
					numSPs++;