Changeset 569 for PythonPackages/OpenOpt/openopt/doc/restricted_dom.py
- Timestamp:
- 06/13/2010 12:24:08 PM (2 years ago)
- Files:
Legend:
- Unmodified
- Added
- Removed
- Modified
- Copied
- Moved
PythonPackages/OpenOpt/openopt/doc/restricted_dom.py
r358 r569 24 24 from openopt import NLP 25 25 26 n = 10 26 n = 100 27 27 an = arange(n) # array [0, 1, 2, ..., n-1] 28 28 x0 = n+15*(1+cos(an)) 29 29 30 # from all OO-connected NLP solvers 31 # only ralg can handle x0 out of dom objFunc: 32 # x0 = n+15*(cos(an)) 30 f = lambda x: (x**2).sum() + sqrt(x**3).sum() 31 df = lambda x: 2*x + 1.5*x**0.5 33 32 34 f = lambda x: (x**2).sum() + sqrt(x**3-arange(n)**3).sum() 35 df = lambda x: 2*x + 0.5*3*x**2/sqrt(x**3-arange(n)**3) 36 37 lb = arange(n) 33 lb = zeros(n) 38 34 solvers = ['ralg'] 35 #solvers = ['ipopt'] 39 36 for solver in solvers: 40 p = NLP(f, x0, df=df, lb=lb, xtol = 1e- 7, iprint = 100, maxIter = 10000, maxFunEvals = 1e8)37 p = NLP(f, x0, df=df, lb=lb, xtol = 1e-6, iprint = 50, maxIter = 10000, maxFunEvals = 1e8) 41 38 #p.checkdf() 42 39 r = p.solve(solver) 43 # expected r.xf = [0, 1, 2, ..., n-1]40 # expected r.xf = small values near zero
