- Timestamp:
- 07/02/08 03:19:46 (4 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
trunk/openopt/scikits/openopt/solvers/BrasilOpt/algencan_oo.py
r1069 r1090 72 72 else: nh = 0 73 73 nb, nbeq = p.b.size, p.beq.size 74 p.algencan.nc, p.algencan.nh, p.algencan.nb, p.algencan.nbeq = nc, nh, nb, nbeq74 #p.algencan.nc, p.algencan.nh, p.algencan.nb, p.algencan.nbeq = nc, nh, nb, nbeq 75 75 76 76 m = nc + nh + nb + nbeq … … 320 320 i = ind - 1 # Python enumeration starts from 0, not 1 321 321 322 if i < p. algencan.nc:322 if i < p.nc: 323 323 c = p.c(x, ind=i) 324 elif p. algencan.nc <= i < p.algencan.nc + p.algencan.nh:325 c = p.h(x, ind=i-p. algencan.nc)326 elif p. algencan.nc + p.algencan.nh <= i < p.algencan.nc + p.algencan.nh + p.algencan.nb:327 j = i - p. algencan.nc - p.algencan.nh328 assert 0<= j <p. algencan.nb324 elif p.nc <= i < p.nc + p.nh: 325 c = p.h(x, ind=i-p.nc) 326 elif p.nc + p.nh <= i < p.nc + p.nh + p.nb: 327 j = i - p.nc - p.nh 328 assert 0<= j <p.nb 329 329 c = p.dotmult(p.A[j], x).sum() - p.b[j] 330 elif i < p. algencan.nc + p.algencan.nh + p.algencan.nb + p.algencan.nbeq:331 j = i - p. algencan.nc - p.algencan.nh - p.algencan.nb332 assert 0<= j <p. algencan.nbeq330 elif i < p.nc + p.nh + p.nb + p.nbeq: 331 j = i - p.nc - p.nh - p.nb 332 assert 0<= j <p.nbeq 333 333 c = p.dotmult(p.Aeq[j], x).sum() - p.beq[j] 334 334 else: … … 389 389 i = ind - 1 # Python enumeration starts from 0, not 1 390 390 391 if 0<=i < p. algencan.nc:391 if 0<=i < p.nc: 392 392 dc = p.dc(x, ind=i) 393 elif p. algencan.nc <= i < p.algencan.nc + p.algencan.nh:394 dc = p.dh(x, ind=i-p. algencan.nc)395 elif p. algencan.nc + p.algencan.nh <= i < p.algencan.nc + p.algencan.nh + p.algencan.nb:396 j = i - p. algencan.nc - p.algencan.nh393 elif p.nc <= i < p.nc + p.nh: 394 dc = p.dh(x, ind=i-p.nc) 395 elif p.nc + p.nh <= i < p.nc + p.nh + p.nb: 396 j = i - p.nc - p.nh 397 397 dc = p.A[j] 398 elif i < p. algencan.nc + p.algencan.nh + p.algencan.nb + p.algencan.nbeq:399 j = i - p. algencan.nc - p.algencan.nh - p.algencan.nb398 elif i < p.nc + p.nh + p.nb + p.nbeq: 399 j = i - p.nc - p.nh - p.nb 400 400 dc = p.Aeq[j] 401 401 else:
