Method addM_Force

SRC/analysis/fe_ele/FE_Element.cpp:765–799

}

void  
FE_Element::addM_Force(const Vector &accel, double fact)
{
    if (myEle != 0) {    

	// check for a quick return
	if (fact == 0.0 || !myEle->isActive()) 
	    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) = accel(loc);
		else
		    tmp(i) = 0.0;		
	    }	 
		
	    if (theResidual->addMatrixVector(1.0, myEle->getMass(), tmp, fact) < 0){
		opserr << "WARNING FE_Element::addM_Force() - ";
		opserr << "- addMatrixVector returned error\n";		 
	    }		
	}
	else {
	    opserr << "WARNING FE_Element::addM_Force() - ";
	    opserr << "- this should not be called on a Subdomain!\n";
	}    	    	    				
    }
    else {
	opserr << "WARNING FE_Element::addM_Force() - no Element *given ";
	opserr << "- subclasses must provide implementation\n";
    }    	            
}

void  
FE_Element::addD_Force(const Vector &accel, double fact)