Method domainChanged

SRC/analysis/integrator/DistributedDisplacementControl.cpp:240–400

int 
DistributedDisplacementControl::domainChanged(void)
{
    // we first create the Vectors needed
    AnalysisModel *theModel = this->getAnalysisModel();
    LinearSOE *theLinSOE = this->getLinearSOE();    
    if (theModel == 0 || theLinSOE == 0) {
	opserr << "WARNING DistributedDisplacementControl::update() ";
	opserr << "No AnalysisModel or LinearSOE has been set\n";
	return -1;
    }    
    
    //int size = theModel->getNumEqn(); // ask model in case N+1 space
    const Vector &b = theLinSOE->getB(); // ask model in case N+1 space
    int size = b.Size();

    // first we determine the id of the nodal dof
    Domain *theDomain = theModel->getDomainPtr();
    if (theDomain == 0) {
      opserr << "BUG WARNING DistributedDisplacementControl::domainChanged() - no Domain associated!!";
      return -1;
    }

    theDofID = -1;
    Node *theNodePtr = theDomain->getNode(theNode);
    if (theNodePtr != 0) {
      DOF_Group *theGroup = theNodePtr->getDOF_GroupPtr();
      if (theGroup == 0) {
	opserr << "BUG DistributedDisplacementControl::domainChanged() - no DOF_Group associated with the node!!\n";
	return -1;
      }
      const ID &theID = theGroup->getID();
      if (theDof < 0 || theDof >= theID.Size()) {
	opserr << "DistributedDisplacementControl::domainChanged() - not a valid dof " << theDof << endln;
	return -1;
      }
      theDofID = theID(theDof);
      if (theDofID < 0) {
	opserr << "DistributedDisplacementControl::domainChanged() - constrained dof not a valid a dof\n";;
	return -1;
      }
    }

    static ID data(1);    

    if (processID != 0) {
      Channel *theChannel = theChannels[0];
      data(0) = theDofID;

      theChannel->sendID(0, 0, data);
      theChannel->recvID(0, 0, data);
      theDofID = data(0);
    } 

    // if main domain, collect all theDofID if not -1 then this is the value & send value out
    else {
      
      for (int j=0; j<numChannels; j++) {