Changeset 1217
- Timestamp:
- 08/10/08 06:24:47 (3 months ago)
- Files:
-
- trunk/openopt/scikits/openopt/tests/chain.py (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
- Modified
- Copied
- Moved
trunk/openopt/scikits/openopt/tests/chain.py
r1199 r1217 35 35 36 36 37 #MaxForces *= 1.001 38 #MaxForces -= array([ 504.84391737, 353.6779311 , 220.55447721, 119.6013147 , 39 # 37.84759491]) 37 40 38 MaxForces -= array([ 470.83299096, 319.6860156 , 186.5286786 , 85.57322313, 3.84866201])-150 41 39 lengths = 5*ones(n)+cos(arange(n))#array([4, 3, 1, 2]) … … 51 49 ub = [inf] * (n+2) 52 50 51 #lb = [-1e5] * 2 + [0] * n 52 #ub = [1e5] * (n+2) 53 53 54 54 ######################################################## … … 161 161 162 162 p = NLP(f, x0, c=c, goal = 'max', ftol = 1e-16, xtol = 1e-16, plot=0, maxIter = 1e4, contol = contol, iprint=1, maxFunEvals = 1e10) 163 p.diffInt = 2 * [1e- 4] + n * [1e-7]163 p.diffInt = 2 * [1e-5] + n * [1e-7] 164 164 p.Aeq = 2 * [0] + n * [1] 165 165 p.beq = 100 # sum of additional masses should be this number 166 #p.h = lambda x: x[2:].sum()-100 166 167 #solver = 'scipy_slsqp' 167 168 solver = 'ralg' 168 solver = 'ipopt'169 #solver = 'ipopt' 169 170 #solver = 'scipy_cobyla' 170 171 #r = p.solve(solver, debug=1, maxIter=1e4, showLS=1, show_nnan = 1) 171 r = p.solve(solver, debug=0, maxIter=1e4, show_nnan = 0) 172 173 172 #solver = 'algencan' 173 r = p.solve(solver, debug=0, iprint = 100, maxIter=1e4, show_nnan = 0) 174 174 175 175 176 176 #----------------------------------------------------- 177 #solver: ralg problem: unnamed goal: m inimum177 #solver: ralg problem: unnamed goal: max 178 178 # iter objFunVal log10(maxResidual) 179 # 0 nan 2.01180 # ...181 # 122 1.854e+02 -2.10179 # 0 1.034e+01 1.98 180 # 100 2.448e+01 -8.33 181 # 116 2.448e+01 -8.33 182 182 #istop: 4 (|| F[k] - F[k-1] || < ftol) 183 #Solver: Time Elapsed = 91.98 CPU Time Elapsed = 91.3184 #objFunValue: 185.4275 (feasible, max constraint = 0.008)183 #Solver: Time Elapsed = 1.7 CPU Time Elapsed = 1.65 184 #objFunValue: 24.480519 (feasible, max constraint = 4.68276e-09)
