Method domainChanged

SRC/analysis/analysis/DomainDecompositionAnalysis.cpp:185–268

}

int
DomainDecompositionAnalysis::domainChanged(void)
{
    // remove existing FE_elements and DOF_Groups from the Analysis
    theModel->clearAll();
    theHandler->clearAll();

    // now we invoke handle() on the constraint handler which
    // causes the creation of FE_Element and DOF_Group objects
    // and their addition to the AnalysisModel.

    numExtEqn = theHandler->handle(&(theSubdomain->getExternalNodes()));

    // we now get a node to number last

    const ID &theExtNodes = theSubdomain->getExternalNodes();
    int idSize = theExtNodes.Size();
    //    int theLastDOF = -1;

    ID theLastDOFs(1);
    int cnt = 0;

    // create an ID containing the tags of the DOF_Groups that are to
    // be numbered last
    for (int i=0; i<idSize; i++) {
	int nodeTag = theExtNodes(i);
	Node *nodePtr = theSubdomain->getNode(nodeTag);
	DOF_Group *dofGrpPtr = nodePtr->getDOF_GroupPtr();
	if (dofGrpPtr != 0) {
	    const ID theID = dofGrpPtr->getID();
	    int size = theID.Size();
	    for (int j=0; j<size; j++)
		if (theID(j) == -3) {
		    theLastDOFs[cnt]  = dofGrpPtr->getTag();
		    cnt++;
		    j = size;
		}
	}
    }

    // we now invoke number() on the numberer which causes
    // equation numbers to be assigned to all the DOFs in the
    // AnalysisModel.    

    theNumberer->numberDOF(theLastDOFs);

    /*************************
    for (int i=0; i<idSize; i++) {
	int nodeTag = theExtNodes(i);
	Node *nodePtr = theSubdomain->getNode(nodeTag);
	DOF_Group *dofPtr = nodePtr->getDOF_GroupPtr();
	if (dofPtr != 0) {
	    const ID theID = dofPtr->getID();
	    int size = theID.Size();
	    for (int j=0; j<size; j++)
		if (theID(j) == -3) {
		    theLastDOF = dofPtr->getTag();
	            i = idSize;