Method addInertiaLoadToUnbalance

SRC/domain/node/Node.cpp:966–996

int
Node::addInertiaLoadToUnbalance(const Vector &accelG, double fact)
{
  // simply return if node has no mass or R matrix
  if (mass == 0 || R == 0)
    return 0;

  // otherwise we must determine MR accelG
  if (accelG.Size() != R->noCols()) {
    opserr << "Node::addInertiaLoadToUnbalance - accelG not of correct dimension";
    return -1;
  }

  // if no load yet create it and assign
  if (unbalLoad == 0) {
      unbalLoad = new Vector(numberDOF); 
      if (unbalLoad == 0 || unbalLoad->Size() != numberDOF) {
	  opserr << "FATAL Node::addunbalLoad - ran out of memory\n";
	  exit(-1);
      }  
  }

  // form - fact * M*R*accelG and add it to the unbalanced load
  //(*unbalLoad) -= ((*mass) * (*R) * accelG)*fact;

  Matrix MR(mass->noRows(), R->noCols());
  MR.addMatrixProduct(0.0, *mass, *R, 1.0);
  unbalLoad->addMatrixVector(1.0, MR, accelG, -fact);

  return 0;
}