#ifndef NDEBUG
H.OK();
A.OK();
- assert (active.size() +inactive.size () == opt->cons.size ());
+ assert (active.size() +inactive.size () == opt->cons_.size ());
assert (H.dim() == opt->dim ());
assert (active.size() == A.rows ());
Array<int> allcons;
- for (int i=0; i < opt->cons.size(); i++)
+ for (int i=0; i < opt->cons_.size(); i++)
allcons.push (0);
for (int i=0; i < active.size(); i++)
{
inactive.swap (k,inactive.size()-1);
inactive.pop();
- Vector a (opt->cons[cidx]);
+ Vector a (opt->cons_[cidx]);
// update of matrices
Vector Ha = H*a;
Real aHa = a*Ha;
/*
*/
- Real q = a*opt->quad*a;
+ Real q = a*opt->quad_*a;
Matrix aaq (a,a/q);
H += aaq;
- A -= A*opt->quad*aaq;
+ A -= A*opt->quad_*aaq;
}else
WARN << _("degenerate constraints");
#ifndef NDEBUG
H(op->dim()),
opt (op)
{
- for (int i=0; i < op->cons.size(); i++)
+ for (int i=0; i < op->cons_.size(); i++)
inactive.push (i);
- Choleski_decomposition chol (op->quad);
+ Choleski_decomposition chol (op->quad_);
/*
ugh.