From ea8845e70aeb417bbdc70612f9a8f13b9d460b8d Mon Sep 17 00:00:00 2001 From: fred Date: Sun, 24 Mar 2002 20:07:56 +0000 Subject: [PATCH] lilypond-0.1.57 --- lily/include/ineq-constrained-qp.hh | 14 ++++++-------- lily/include/qlpsolve.hh | 4 ++-- lily/qlpsolve.cc | 14 +++++++------- 3 files changed, 15 insertions(+), 17 deletions(-) diff --git a/lily/include/ineq-constrained-qp.hh b/lily/include/ineq-constrained-qp.hh index a7bd31693d..5cb5c43eb9 100644 --- a/lily/include/ineq-constrained-qp.hh +++ b/lily/include/ineq-constrained-qp.hh @@ -18,12 +18,12 @@ class Ineq_constrained_qp { friend class Active_constraints; - Array cons; - Array consrhs; + Array cons_; + Array consrhs_; public: - Matrix quad; - Vector lin; - Real const_term; + Matrix quad_; + Vector lin_; + Real const_term_; /** @@ -40,9 +40,7 @@ public: /** @return the number of variables in the problem */ - int dim() const{ - return lin.dim(); - } + int dim() const; /** add a constraint diff --git a/lily/include/qlpsolve.hh b/lily/include/qlpsolve.hh index 6c018c5a9f..db985df6c5 100644 --- a/lily/include/qlpsolve.hh +++ b/lily/include/qlpsolve.hh @@ -37,8 +37,8 @@ class Active_constraints { public: String status() const; - Vector vec (int k) const { return opt->cons[k]; } - Real rhs (int k) const { return opt->consrhs[k]; } + Vector vec (int k) const { return opt->cons_[k]; } + Real rhs (int k) const { return opt->consrhs_[k]; } /** drop constraint. drop constraint k from the active set. k is the index of the diff --git a/lily/qlpsolve.cc b/lily/qlpsolve.cc index 7baa0a1a4a..3d45207eab 100644 --- a/lily/qlpsolve.cc +++ b/lily/qlpsolve.cc @@ -41,12 +41,12 @@ Active_constraints::OK() #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 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++) { @@ -79,7 +79,7 @@ Active_constraints::add (int k) 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; @@ -123,10 +123,10 @@ Active_constraints::drop (int k) /* */ - 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 @@ -143,9 +143,9 @@ Active_constraints::Active_constraints (Ineq_constrained_qp const *op) 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. -- 2.39.5