Method incrNodeAccel

SRC/analysis/dof_grp/TransformationDOF_Group.cpp:729–756

void
TransformationDOF_Group::incrNodeAccel(const Vector &u)
{
  // call base class method and return if no MP_Constraint
  if (theMP == 0) {
    this->DOF_Group::incrNodeAccel(u);
    return;
  }
  
  const ID &theID = this->getID();
  for (int i=0; i<modNumDOF; i++) {
    int loc = theID(i);
    if (loc >= 0)
      (*modUnbalance)(i) = u(loc);
    else 	
      (*modUnbalance)(i) = 0.0;	    
  }    
  Matrix *T = this->getT();

  // *unbalance = (*T) * (*modUnbalance);
  unbalance->addMatrixVector(0.0, *T, *modUnbalance, 1.0);
  int numDOF = myNode->getNumberDOF();
  for (int i=0; i<numDOF; i++) {
    if (theSPs[i] != 0)
      (*unbalance)(i) = 0.0;
  }
  myNode->incrTrialAccel(*unbalance);
}

const Vector & 
TransformationDOF_Group::getTrialDisp(void)