bug fix in cdirect -- typo caused slope to be underestimated in some cases, and will...

darcs-hash:20110329225957-c8de0-96a76c698ce0a934f7f0c2f5334aa259f073b0ef.gz

diff --git a/cdirect/cdirect.c b/cdirect/cdirect.c
index d37e8913e73e37fa1bba20d731635420368bd0fe..b9b909a53e2fa55a468f826217a368503e3f65d4 100644 (file)
--- a/cdirect/cdirect.c
+++ b/cdirect/cdirect.c
@@ -414,7 +414,7 @@ static nlopt_result divide_good_rects(params *p)
          if (im >= 0)
               K1 = (hull[i][1] - hull[im][1]) / (hull[i][0] - hull[im][0]);
          if (ip < nhull)
-              K1 = (hull[i][1] - hull[ip][1]) / (hull[i][0] - hull[ip][0]);
+              K2 = (hull[i][1] - hull[ip][1]) / (hull[i][0] - hull[ip][0]);
          K = MAX(K1, K2);
          if (hull[i][1] - K * hull[i][0]
              <= p->minf - magic_eps * fabs(p->minf) || ip == nhull) {