Changeset 1533
- Timestamp:
- 10/10/08 09:14:54 (2 months ago)
- Files:
Legend:
- Unmodified
- Added
- Removed
- Modified
- Copied
- Moved
trunk/openopt/scikits/openopt/examples/oofun/speedup.py
r1531 r1533 15 15 r += z ** (1 / k**1.5) 16 16 return r 17 f = lambda z: ((z-aN)**2).sum() 17 18 18 N = 150 19 K = 500 19 N, K = 150, 500 20 20 solver = 'scipy_ncg'# try also scipy_cg, scipy_ncg, ralg, algencan etc 21 ftol = 1e-6 22 xtol = 1e-6 23 gtol = 1e-6 21 ftol, xtol, gtol = 1e-6, 1e-6, 1e-6 24 22 iprint = 5 25 23 … … 30 28 """ 1: using oovar & oofun """ 31 29 v = oovar('v', size = N) # start value will be zeros(N) 32 f = oofun(lambda z: ((z-aN)**2).sum(), input = v) 33 g = oofun(CostlyFunction, input = f) 34 30 ff = oofun(f, input = v) 31 g = oofun(CostlyFunction, input = ff) 35 32 p = NLP(g, maxIter=1e4, iprint=iprint, ftol=ftol, xtol=xtol, gtol=gtol) 36 33 r = p.solve(solver) 37 34 38 35 """ 2: classic """ 39 f = lambda z: ((z-aN)**2).sum()40 36 g = CostlyFunction 41 37 p = NLP(lambda x: g(f(x)), x0=zeros(N), maxIter=1e4, ftol=ftol, xtol=xtol, gtol=gtol, iprint=iprint)
