Changeset 1024

Show
Ignore:
Timestamp:
06/15/08 13:35:45 (6 months ago)
Author:
dmitrey.kroshko
Message:

code cleanup

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved

  • trunk/openopt/scikits/openopt/solvers/UkrOpt/ralg_oo.py

    r1021 r1024  
    143143                        newPoint = p.point(x) 
    144144 
    145                         if newPoint.mr() > PrevPoint.mr() or newPoint.f() > PrevPoint.f() or abs(newPoint.f() - iterPoint.f()) < 15 * p.ftol or p.norm(newPoint.x - iterPoint.x) < 15 * p.xtol: 
     145                        if PrevPoint.betterThan(newPoint) or newPoint.f() > PrevPoint.f() or abs(newPoint.f() - iterPoint.f()) < 15 * p.ftol or p.norm(newPoint.x - iterPoint.x) < 15 * p.xtol: 
    146146                            iterPoint, hs = PrevPoint, hs_prev 
    147147                            break 
    …… 
    178178                    g1 = g2-constrGradPrevIter 
    179179                    if sum(p.dotmult(constrGradPrevIter, g2))>0: 
     180                        #TODO: add debugmsg here 
    180181                        doDilation = False 
    181182            elif prevIterPointIsFeasible: 
    …… 
    195196            if not self.needRej(p, b, g1, g): 
    196197                if all(isfinite(g)) and ng > 1e-50 and doDilation: 
    197                     g = (g / ng).reshape((-1,1)) 
    198                     vec1 = self.__economyMult__(b, g).reshape((-1,1)) 
     198                    g = (g / ng).reshape(-1,1) 
     199                    vec1 = self.__economyMult__(b, g).reshape(-1,1) 
    199200                    vec2 = p.dotmult(g.T, w) 
    200201                    b += p.matmult(vec1, vec2) 
    …… 
    205206                #hs = max(p.norm(xPrevIter - x), hsmin) 
    206207 
    207             #xk, fk, rk = iterPoint.x, iterPoint.f(), iterPoint.mr() 
    208             p.iterfcn(xk,  fk) 
     208            p.iterfcn(xk, fk, rk) 
    209209 
    210210            s2 = 0 
    …… 
    225225                            s2 = key 
    226226                            break 
    227                     #p.msg += ' /checkme/' 
    228  
    229227                p.istop = s2 
    230228 
    …… 
    372370            ind = where(v != 0)[0] 
    373371            r = dot(M[:,ind], v[ind]) 
    374 ##            if norm(r-dot(M,v))>1e-17: 
    375 ##                print 'diff', norm(r-dot(M,v)) 
    376372            return r 
    377373