Method domainChanged

SRC/analysis/integrator/DisplacementControl.cpp:362–556

int 
DisplacementControl::domainChanged(void)
{
   // we first create the Vectors needed
  //opserr<<" this is the domain change function"<<endln;//Abbas
   AnalysisModel *theModel = this->getAnalysisModel();
   LinearSOE *theLinSOE = this->getLinearSOE(); 
    if (theModel == 0 || theLinSOE == 0) {
      opserr << "WARNING DisplacementControl::update() ";
      opserr << "No AnalysisModel or LinearSOE has been set\n";
      return -1;
   }    
   int size = theModel->getNumEqn(); // ask model in case N+1 space

   if (deltaUhat == 0 || deltaUhat->Size() != size) { // create new Vector
      if (deltaUhat != 0)
	 delete deltaUhat;   // delete the old
      deltaUhat = new Vector(size);
      if (deltaUhat == 0 || deltaUhat->Size() != size) { // check got it
	 opserr << "FATAL DisplacementControl::domainChanged() - ran out of memory for";
	 opserr << " deltaUhat Vector of size " << size << endln;
	 exit(-1);
      }
   }
  

   if (deltaUbar == 0 || deltaUbar->Size() != size) { // create new Vector
      if (deltaUbar != 0)
	 delete deltaUbar;   // delete the old
      deltaUbar = new Vector(size);
      if (deltaUbar == 0 || deltaUbar->Size() != size) { // check got it
	opserr << "FATAL DisplacementControl::domainChanged() - ran out of memory for";
	opserr << " deltaUbar Vector of size " << size << endln;
	exit(-1);
      }
   }

   if (deltaU == 0 || deltaU->Size() != size) { // create new Vector
     if (deltaU != 0)
       delete deltaU;   // delete the old
     deltaU = new Vector(size);
     if (deltaU == 0 || deltaU->Size() != size) { // check got it
       opserr << "FATAL DisplacementControl::domainChanged() - ran out of memory for";
       opserr << " deltaU Vector of size " << size << endln;
       exit(-1);
      }
   }

   if (deltaUstep == 0 || deltaUstep->Size() != size) { 
      if (deltaUstep != 0)
	delete deltaUstep;  
      deltaUstep = new Vector(size);
      if (deltaUstep == 0 || deltaUstep->Size() != size) { 
	 opserr << "FATAL DisplacementControl::domainChanged() - ran out of memory for";
	 opserr << " deltaUstep Vector of size " << size << endln;
	 exit(-1);
      }
   }