43 #ifndef ROOT_TQpResidual 44 #define ROOT_TQpResidual 128 void Project_r3 (Double_t rmin,Double_t rmax);
void CalcResids(TQpDataBase *problem, TQpVar *vars)
Calculate residuals, their norms, and duality complementarity gap, given a problem and variable set...
Bool_t ValidNonZeroPattern()
Check if vector elements as selected through array indices are non-zero.
void Project_r3(Double_t rmin, Double_t rmax)
Perform the projection operation required by Gondzio algorithm: replace each component r3_i of the co...
static void GondzioProjection(TVectorD &v, Double_t rmin, Double_t rmax)
Replace each component r3_i of the complementarity component of the residuals by r3p_i-r3_i, where r3p_i is the projection of r3_i onto the box [rmin, rmax].
TQpResidual & operator=(const TQpResidual &source)
Assignment operator.
Double_t GetResidualNorm()
void Add_r3_xz_alpha(TQpVar *vars, Double_t alpha)
Modify the "complementarity" component of the residuals, by adding the pairwise products of the compl...
void Clear_r3()
set the complementarity component of the residuals to 0.
void Set_r3_xz_alpha(TQpVar *vars, Double_t alpha)
Set the "complementarity" component of the residuals to the pairwise products of the complementary va...
void Clear_r1r2()
set the noncomplementarity components of the residual (the terms arising from the linear equalities i...
TQpResidual()
Constructor.