Show
Ignore:
Timestamp:
07/02/08 03:19:46 (4 years ago)
Author:
dmitrey.kroshko
Message:

some changes

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • trunk/openopt/scikits/openopt/solvers/BrasilOpt/algencan_oo.py

    r1069 r1090  
    7272            else: nh = 0 
    7373            nb, nbeq = p.b.size, p.beq.size 
    74             p.algencan.nc, p.algencan.nh, p.algencan.nb, p.algencan.nbeq = nc, nh, nb, nbeq 
     74            #p.algencan.nc, p.algencan.nh, p.algencan.nb, p.algencan.nbeq = nc, nh, nb, nbeq 
    7575 
    7676            m = nc + nh + nb + nbeq 
     
    320320            i = ind - 1 # Python enumeration starts from 0, not 1 
    321321 
    322             if i < p.algencan.nc: 
     322            if i < p.nc: 
    323323                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.nh 
    328                 assert 0<= j <p.algencan.nb 
     324            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 
    329329                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.nb 
    332                 assert 0<= j <p.algencan.nbeq 
     330            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 
    333333                c = p.dotmult(p.Aeq[j], x).sum() - p.beq[j] 
    334334            else: 
     
    389389            i = ind - 1 # Python enumeration starts from 0, not 1 
    390390 
    391             if 0<=i < p.algencan.nc: 
     391            if 0<=i < p.nc: 
    392392                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.nh 
     393            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 
    397397                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.nb 
     398            elif i < p.nc + p.nh + p.nb + p.nbeq: 
     399                j = i - p.nc - p.nh - p.nb 
    400400                dc = p.Aeq[j] 
    401401            else: