Changeset 1171

Show
Ignore:
Timestamp:
08/01/08 12:42:48 (4 months ago)
Author:
dmitrey.kroshko
Message:

minor change in example

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved

  • trunk/openopt/scikits/openopt/examples/nlp_bench_1.py

    r1129 r1171  
    2828 
    2929############################################################### 
    30 colors = ['b', 'k', 'y', 'g', 'r'] 
    31 solvers = ['ralg', 'scipy_cobyla', 'lincher', 'scipy_slsqp', 'algencan'] 
    32 #colors = ['b',  'r'] 
    33 #solvers = ['ralg',  'ALGENCAN'] 
     30colors = ['b', 'k', 'y', 'g', 'r', 'm'] 
     31solvers = ['ralg', 'scipy_cobyla', 'lincher', 'scipy_slsqp', 'ipopt','algencan'] 
    3432############################################################### 
    3533 
    …… 
    3836    solver = solvers[j] 
    3937    color = colors[j] 
    40     p = NLP(ff, x0, c=c, h=h,  lb = lb, ub = ub, gradtol=gradtol, ftol = ftol, maxFuevals = 1e7, maxIter = 1e4, plot = 1, color = color, iprint = 0, legend = [solvers[j]], show=False,  contol = contol,  maxTime = maxTime,  maxFunEvals = maxFunEvals) 
    41    
     38    p = NLP(ff, x0, c=c, h=h,  lb = lb, ub = ub, gradtol=gradtol, diffInt = diffInt, ftol = ftol, maxFuevals = 1e7, maxIter = 1e4, plot = 1, color = color, iprint = 0, legend = [solvers[j]], show=False,  contol = contol,  maxTime = maxTime,  maxFunEvals = maxFunEvals) 
    4239 
    43     if solver == 'ALGENCAN':  
     40 
     41    if solver =='algencan': 
     42        p.gradtol = 1e-2 
     43    elif solver == 'ralg': 
     44        p.debug = 1 
    4445        pass 
    45         #p.debug = 1 
    46     elif solver == 'ralg': 
    47         #p.debug = 1 
    48         #p.iprint = 1 
    49         p.gradtol = 1e-8 
    50         p.ftol = 1e-9 
    51         p.xtol = 1e-7 
    52          
    53    
     46 
    5447    r = p.solve(solver) 
    55     for fn in ('h','c'):  
     48    for fn in ('h','c'): 
    5649        if not r.evals.has_key(fn): r.evals[fn]=0 # if no c or h are used in problem 
    5750    results[solver] = (r.ff, p.getMaxResidual(r.xf), r.elapsed['solver_time'], r.elapsed['solver_cputime'], r.evals['f'], r.evals['c'], r.evals['h']) 
    …… 
    6861xl = ['Solver                              f_opt     MaxConstr   Time   CPUTime  fEvals  cEvals  hEvals'] 
    6962 
    70 for i in range(len(results)):  
     63for i in range(len(results)): 
    7164    s=(ljust(lower(solvers[i]), 40-len(solvers[i]))+'%0.3f'% (results[solvers[i]][0]) + '        %0.1e' % (results[solvers[i]][1]) + ('      %0.2f'% (results[solvers[i]][2])) + '     %0.2f      '% (results[solvers[i]][3]) + str(results[solvers[i]][4]) + '   ' + rjust(str(results[solvers[i]][5]), 5) + ' '*8 +str(results[solvers[i]][6])) 
    7265