Method LagrangeMP_FE

SRC/analysis/fe_ele/lagrange/LagrangeMP_FE.cpp:55–117

//    CDOF = a1*RDOF1 + a2*RDOF2 + ... + an*RDOFn + beta

LagrangeMP_FE::LagrangeMP_FE(int tag, Domain &theDomain, MP_Constraint &TheMP,
			     DOF_Group &theGroup, double Alpha)
:FE_Element(tag, 3,(tag, TheMP.getConstrainedDOFs()).Size()+
	      (TheMP.getRetainedDOFs()).Size() + 
	      (TheMP.getConstrainedDOFs()).Size()), // *see note 1
 alpha(Alpha), theMP(&TheMP), 
 theConstrainedNode(0), theRetainedNode(0),
 theDofGroup(&theGroup), tang(0), resid(0)
{
    const Matrix &constraint = theMP->getConstraint();
    int noRows = constraint.noRows();
    int noCols = constraint.noCols();
    int size = 2*noRows+noCols;
    
    tang = new Matrix(size,size);
    resid = new Vector(size);
    if (tang == 0 || resid == 0 || tang->noCols() == 0 || resid->Size() == 0) {
	opserr << "FATAL LagrangeMP_FE::LagrangeMP_FE() - out of memory\n";
	exit(-1);
    }
    tang->Zero();	
    resid->Zero();

    theRetainedNode = theDomain.getNode(theMP->getNodeRetained());    
    theConstrainedNode = theDomain.getNode(theMP->getNodeConstrained());

    if (theRetainedNode == 0) {
	opserr << "WARNING LagrangeMP_FE::LagrangeMP_FE()";
	opserr << "- no asscoiated Retained Node\n";
	exit(-1);
    }
    
    if (theConstrainedNode == 0) {
	opserr << "WARNING LagrangeMP_FE::LagrangeMP_FE()";
	opserr << "- no asscoiated Constrained Node\n";
	exit(-1);
    }
    
    if (theMP->isTimeVarying() == false) {
	this->determineTangent();
    }
    
    // set the myDOF_Groups tags indicating the attached id's of the
    // DOF_Group objects
    DOF_Group *theConstrainedNodesDOFs = theConstrainedNode->getDOF_GroupPtr();
    if (theConstrainedNodesDOFs == 0) {
	opserr << "WARNING LagrangeMP_FE::LagrangeMP_FE()";
	opserr << " - no DOF_Group with Constrained Node\n";
	exit(-1);	
    }    

    DOF_Group *theRetainedNodesDOFs = theRetainedNode->getDOF_GroupPtr();
    if (theRetainedNodesDOFs == 0) {
	opserr << "WARNING LagrangeMP_FE::LagrangeMP_FE()";
	opserr << " - no DOF_Group with Retained Node\n";
	exit(-1);	
    }