Method determineTangent

SRC/analysis/fe_ele/lagrange/LagrangeMP_FE.cpp:323–344

}

void  
LagrangeMP_FE::determineTangent(void)
{
    const Matrix &constraint = theMP->getConstraint();
    int noRows = constraint.noRows();
    int noCols = constraint.noCols();
    int n = noRows+noCols;
    
    tang->Zero();    

    for (int j=0; j<noRows; j++) {
	(*tang)(n+j, j) = -alpha;
	(*tang)(j, n+j) = -alpha;	
    }
    
    for (int i=0; i<noRows; i++)
	for (int j=0; j<noCols; j++) {
	    double val = constraint(i,j) * alpha;
	    (*tang)(n+i, j+noRows) = val;
	    (*tang)(noRows+j, n+i) = val;
	}
}