Method addD_ForceSensitivity

SRC/analysis/fe_ele/FE_Element.cpp:1000–1033

}

void  
FE_Element::addD_ForceSensitivity(int gradNumber, const Vector &vect, double fact)
{
  if (myEle != 0) {    
    
    // check for a quick return
    if (fact == 0.0) 
      return;
    if (myEle->isSubdomain() == false) {
      // get the components we need out of the vector
      // and place in a temporary vector
      Vector tmp(numDOF);
      for (int i=0; i<numDOF; i++) {
	int loc = myID(i);
	if (loc >= 0)
	  tmp(i) = vect(loc);
	else
	  tmp(i) = 0.0;		
      }	
      if (theResidual->addMatrixVector(1.0, myEle->getDampSensitivity(gradNumber), tmp, fact) < 0){
	opserr << "WARNING FE_Element::addD_ForceSensitivity() - ";
	opserr << "- addMatrixVector returned error\n";		 
      }		
    }
    else {
      opserr << "WARNING FE_Element::addD_ForceSensitivity() - ";
      opserr << "- this should not be called on a Subdomain!\n";
    }    	    	    				
  }
  else {
    opserr << "WARNING FE_Element::addD_ForceSensitivity() - no Element *given ";
    opserr << "- subclasses must provide implementation\n";
  }    	            
}

void  
FE_Element::addLocalD_ForceSensitivity(int gradNumber, const Vector &accel, double fact)