Changeset 1090

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

some changes

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved

  • trunk/openopt/scikits/openopt/Kernel/runProbSolver.py

    r1089 r1090  
    146146        else: 
    147147            setattr(p, 'n'+s, 0) 
     148    for s in ['b','beq']: 
     149        if hasattr(p, s): setattr(p, 'n'+s, len(getattr(p, s))) 
     150 
    148151    if p.probType not in ['LP', 'QP', 'MILP', 'LLSP']: p.objFunc(p.x0) 
    149152 

  • trunk/openopt/scikits/openopt/solvers/BrasilOpt/ALGENCAN_oo.py

    r978 r1090  
    7171            else: nh = 0 
    7272            nb, nbeq = p.b.size, p.beq.size 
    73             p.ALGENCAN.nc, p.ALGENCAN.nh, p.ALGENCAN.nb, p.ALGENCAN.nbeq = nc, nh, nb, nbeq 
     73            #p.nc, p.nh, p.nb, p.nbeq = nc, nh, nb, nbeq 
    7474 
    7575            m = nc + nh + nb + nbeq 
    …… 
    304304            i = ind - 1 # Python enumeration starts from 0, not 1 
    305305 
    306             if i < p.ALGENCAN.nc: 
     306            if i < p.nc: 
    307307                c = p.c(x, ind=i) 
    308             elif p.ALGENCAN.nc <= i < p.ALGENCAN.nc + p.ALGENCAN.nh: 
    309                 c = p.h(x, ind=i-p.ALGENCAN.nc) 
    310             elif p.ALGENCAN.nc + p.ALGENCAN.nh <= i < p.ALGENCAN.nc + p.ALGENCAN.nh + p.ALGENCAN.nb: 
    311                 j = i - p.ALGENCAN.nc - p.ALGENCAN.nh 
    312                 assert 0<= j <p.ALGENCAN.nb 
     308            elif p.nc <= i < p.nc + p.nh: 
     309                c = p.h(x, ind=i-p.nc) 
     310            elif p.nc + p.nh <= i < p.nc + p.nh + p.nb: 
     311                j = i - p.nc - p.nh 
     312                assert 0<= j <p.nb 
    313313                c = p.dotmult(p.A[j], x).sum() - p.b[j] 
    314             elif i < p.ALGENCAN.nc + p.ALGENCAN.nh + p.ALGENCAN.nb + p.ALGENCAN.nbeq: 
    315                 j = i - p.ALGENCAN.nc - p.ALGENCAN.nh - p.ALGENCAN.nb 
    316                 assert 0<= j <p.ALGENCAN.nbeq 
     314            elif i < p.nc + p.nh + p.nb + p.nbeq: 
     315                j = i - p.nc - p.nh - p.nb 
     316                assert 0<= j <p.nbeq 
    317317                c = p.dotmult(p.Aeq[j], x).sum() - p.beq[j] 
    318318            else: 
    …… 
    373373            i = ind - 1 # Python enumeration starts from 0, not 1 
    374374 
    375             if 0<=i < p.ALGENCAN.nc: 
     375            if 0<=i < p.nc: 
    376376                dc = p.dc(x, ind=i) 
    377             elif p.ALGENCAN.nc <= i < p.ALGENCAN.nc + p.ALGENCAN.nh: 
    378                 dc = p.dh(x, ind=i-p.ALGENCAN.nc) 
    379             elif p.ALGENCAN.nc + p.ALGENCAN.nh <= i < p.ALGENCAN.nc + p.ALGENCAN.nh + p.ALGENCAN.nb: 
    380                 j = i - p.ALGENCAN.nc - p.ALGENCAN.nh 
     377            elif p.nc <= i < p.nc + p.nh: 
     378                dc = p.dh(x, ind=i-p.nc) 
     379            elif p.nc + p.nh <= i < p.nc + p.nh + p.nb: 
     380                j = i - p.nc - p.nh 
    381381                dc = p.A[j] 
    382             elif i < p.ALGENCAN.nc + p.ALGENCAN.nh + p.ALGENCAN.nb + p.ALGENCAN.nbeq: 
    383                 j = i - p.ALGENCAN.nc - p.ALGENCAN.nh - p.ALGENCAN.nb 
     382            elif i < p.nc + p.nh + p.nb + p.nbeq: 
     383                j = i - p.nc - p.nh - p.nb 
    384384                dc = p.Aeq[j] 
    385385            else: 

  • 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: