|
Revision 1148, 0.5 kB
(checked in by dmitrey.kroshko, 5 months ago)
|
minor change in llavp example
|
| Line | |
|---|
| 1 |
__docformat__ = "restructuredtext en" |
|---|
| 2 |
|
|---|
| 3 |
from numpy import empty, sin, cos, arange, ones |
|---|
| 4 |
from scikits.openopt import LLAVP |
|---|
| 5 |
|
|---|
| 6 |
M, N = 150, 15 |
|---|
| 7 |
C = empty((M,N)) |
|---|
| 8 |
d = empty(M) |
|---|
| 9 |
|
|---|
| 10 |
for j in xrange(M): |
|---|
| 11 |
d[j] = 1.5*N+80*sin(j) |
|---|
| 12 |
C[j] = 8*sin(4.0+arange(N)) + 15*cos(j) |
|---|
| 13 |
|
|---|
| 14 |
lb = sin(arange(N)) |
|---|
| 15 |
ub = lb + 1 |
|---|
| 16 |
p = LLAVP(C, d, lb=lb, ub=ub, dump = 10, X = ones(N), maxIter = 1e4, maxFunEvals = 1e100) |
|---|
| 17 |
|
|---|
| 18 |
|
|---|
| 19 |
r = p.solve('nsp:ralg', iprint = 100, maxIter = 1000) |
|---|
| 20 |
|
|---|
| 21 |
|
|---|
| 22 |
|
|---|
| 23 |
print 'f_opt:', r.ff |
|---|
| 24 |
|
|---|
