Changeset 500

Timestamp:

08/17/07 14:09:25 (1 year ago)

Author:

dmitrey.kroshko

Message:

add some examples

Files:

• trunk/openopt/scikits/openopt/examples/milp_1.py (modified) (1 diff)
• trunk/openopt/scikits/openopt/examples/nlp_1.py (modified) (1 diff)
• trunk/openopt/scikits/openopt/examples/nlp_2.py (added)

trunk/openopt/scikits/openopt/examples/milp_1.py

@@ +1 @@
+__docformat__ = "restructuredtext en"
+
-from numpy import *
+from scikits.openopt import MILP
+
-f = [1, 2, 3, 4, 5, 4, 2, 1]
-intVars = [5, 8]

trunk/openopt/scikits/openopt/examples/nlp_1.py

@@ +1 @@
+"""
+Example:
+(x0-5)^2 + (x1-5)^2 + ... +(x149-5)^2 -> min
+subjected to 
+lb<= x <= ub:
+x4 <= 4
+8 <= x5 <= 15
+
+# Ax <= b
+x0+...+x149 >= 800
+x10+x11 <= 9
+
+# Aeq x = beq
+x100+x101 = 8
+
+# c(x) <= 0
+2*x0^4-32 <= 0
+x1^2+x2^2-8 <= 0
+
+# h(x) = 0
+1e6*(x[149]-1)**4 = 0
+(x[148]-1.5)**4 = 0
+"""
+
+
-from scikits.openopt import NLP
-
-
-
-
-
+
+
+
+
+
+
+
+
-p.c = lambda x: [2* x[0] **4-32, x[1]**2+x[2]**2 - 8]
+#c(x) gradients (optional):
+
+def DC(x):
+    r = zeros((p.n, len(p.c)))
+    r[0,0] = 2 * 4 * x[0]**3
+    r[1,1] = 2 * x[1]
+    r[2,1] = 2 * x[2]
+    return r    
+p.dc = DC
+
-h1 = lambda x: 1e6*(x[-1]-1)**4
-h2 = lambda x: (x[-2]-1.5)**4
-
+
+
+
+
+
+
+
+
+
-p.lb = -6*ones(p.n)
-p.ub = 6*ones(p.n)
-
-
-
+
+
+
+
+
+
+
+
+
+
+
-
-r = p.solve('lincher')
+# r.xf and r.ff are optim point and optim objFun value
+# r.ff should be something like 98.6134737452