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.