Method determineTangent

SRC/analysis/fe_ele/lagrange/LagrangeEQ_FE.cpp:290–306

}

void  
LagrangeEQ_FE::determineTangent(void)
{
    const Vector &constraint = theEQ->getConstraint();
    int size = constraint.Size();
    
    tang->Zero();    

    (*tang)(size + 1, 0) = -alpha;
    (*tang)(0, size + 1) = -alpha;	
    
    for (int i = 0; i < size; i++) {
        double val = constraint(i) * alpha;
        (*tang)(size + 1, i + 1) = val;
        (*tang)(i + 1, size + 1) = val;
    }
}