Method incrNodeDisp

SRC/analysis/dof_grp/TransformationDOF_Group.cpp:662–695

//	values in u, components identified by myID;

void
TransformationDOF_Group::incrNodeDisp(const Vector &u)
{
    // call base class method and return if no MP_Constraint
    if (theMP == 0) {
	this->DOF_Group::incrNodeDisp(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;	    
   }    
   
#ifdef TRANSF_INCREMENTAL_MP
   modTotalDisp.addVector(1.0, *modUnbalance, 1.0); // accumulate it for next iteration
#endif // TRANSF_INCREMENTAL_MP

   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->incrTrialDisp(*unbalance);
}


void