Method setNodeDisp

SRC/analysis/dof_grp/TransformationDOF_Group.cpp:479–564

//	values in u, components identified by myID;

void
TransformationDOF_Group::setNodeDisp(const Vector &u)
{
#ifdef TRANSF_INCREMENTAL_MP
    // save the previous mod trial here
    static Vector modTrialDispOld;
    modTrialDispOld = modTotalDisp; // at previous iteration
#endif // TRANSF_INCREMENTAL_MP

  // call base class method and return if no MP_Constraint
  if (theMP == 0) {
    this->DOF_Group::setNodeDisp(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;
  }


  if (needRetainedData == 0) {
    int retainedNode = theMP->getNodeRetained();
    Domain *theDomain = myNode->getDomain();
    Node *retainedNodePtr = theDomain->getNode(retainedNode);
    const Vector &responseR = retainedNodePtr->getTrialDisp();
    const ID &retainedDOF = theMP->getRetainedDOFs();

    for (int i=numConstrainedNodeRetainedDOF, j=0; i<modNumDOF; i++, j++) {
      int loc = theID(i);
      if (loc < 0)
	(*modUnbalance)(i) = responseR(retainedDOF(j));
    }
  }

#ifdef TRANSF_INCREMENTAL_MP
  modTotalDisp = *modUnbalance; // save it for next iteration
#endif // TRANSF_INCREMENTAL_MP

  // at this point the modUnbalance contains the reduced total displacement.
  // remove the trial one to obtain the increment, so that we transform only the increment
#ifdef TRANSF_INCREMENTAL_MP
  modUnbalance->addVector(1.0, modTrialDispOld, -1.0);
#ifdef TRANSF_INCREMENTAL_MP_DEBUG
  opserr << " N = " << myNode->getTag() << "\n";
  opserr << " solut: " << u;
  opserr << " oldtr: " << modTrialDispOld;
  opserr << " trial: " << myNode->getTrialDisp();
  opserr << " commi: " << myNode->getDisp();
  opserr << " incre: " << myNode->getIncrDisp();
  opserr << " delta: " << myNode->getIncrDeltaDisp();
  Domain* dom = myNode->getDomain();
  Node* ret = dom->getNode(theMP->getNodeRetained());
  opserr << " R VEL: " << ret->getTrialVel();