TODO

   1 Modify luksan code (e.g. luksan/plis.c) to use full NLopt stopping
   2 criteria, and in particular the absolute tolerances and maxtime.
   3
   4 Python interface.
   5
   6 SLSQP routine from scipy.
   7
   8 Add nonlinearly constrained multistart algorithm.
   9
  10 More object-based interface, e.g. nlopt_create(dimensions, function,
  11 algorithm), then nlopt_add_constraint(....) etcetera?  This would be
  12 more scalable and backwards compatible, so that I could add new
  13 features re-entrantly without breaking backwards compatibility or
  14 making the arargument list increasingly unwieldy.