Method update

SRC/analysis/integrator/CentralDifferenceNoDamping.cpp:208–252

int
CentralDifferenceNoDamping::update(const Vector &X)
{
  updateCount++;
  if (updateCount > 1) {
    opserr << "ERROR CentralDifferenceNoDamping::update() - called more than once -";
    opserr << " Central Difference integraion schemes require a LINEAR solution algorithm\n";
    return -1;
  }
  
  AnalysisModel *theModel = this->getAnalysisModel();

  if (theModel == 0) {
    opserr << "ERROR CentralDifferenceNoDamping::update() - no AnalysisModel set\n";
    return -2;
  }	
  
  // check domainChanged() has been called, i.e. Ut will not be zero
  if (U == 0) {
    opserr << "WARNING CentralDifferenceNoDamping::update() - domainChange() failed or not called\n";
    return -2;
  }	

  // check deltaU is of correct size
  if (X.Size() != U->Size()) {
    opserr << "WARNING CentralDifferenceNoDamping::update() - Vectors of incompatible size ";
    opserr << " expecting " << U->Size() << " obtained " << X.Size() << endln;
    return -3;
  }

  //  determine the acceleration at time t 
  (*Udotdot) = X;

  //  determine the vel at t+ 0.5 * delta t 
  Udot->addVector(1.0, X, deltaT);
  
  //  determine the displacement at t+delta t 
  U->addVector(1.0, *Udot, deltaT);

  // update the disp & responses at the DOFs
  theModel->setDisp(*U);
  theModel->updateDomain();

  return 0;
}    

int
CentralDifferenceNoDamping::commit(void)