Changeset 428

Show
Ignore:
Timestamp:
02/06/10 14:17:49 (1 month ago)
Author:
dmitrey
Message:

minor llsp changes

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved
  • PythonPackages/OpenOpt/openopt/kernel/LLSP.py

    r427 r428  
    82 82 #    return dot(r, r) 
    83 83 ff = lambda x, LLSPprob: LLSPprob.objFunc(x) 
    84   def dff(x, LLSPprob): 
    85       r = dot(LLSPprob.C.T, dot(LLSPprob.C,x)  - LLSPprob.d) 
    86       if not LLSPprob.damp is None: r += LLSPprob.damp*(x - LLSPprob.X) 
    87       if LLSPprob.f is not None and all(isfinite(LLSPprob.f)) : r += LLSPprob.f 
      84 def dff(x, p): 
      85     r = p.matMultVec(p.C.T, p.matMultVec(p.C,x)  - p.d) 
      86     if p.damp is not None and p.damp != 0: r += p.damp*(x - LLSPprob.X) 
      87     if p.f is not None and all(isfinite(p.f)) : r += p.f 
    88 88     return r 
    89 89  
    90   def d2ff(x, LLSPprob): 
    91       r = dot(LLSPprob.C.T, LLSPprob.C) 
    92       if not LLSPprob.damp is None: r += LLSPprob.damp*eye(x.size) 
      90 def d2ff(x, p): 
      91     # TODO: handle sparse! 
      92     r = dot(p.C.T, p.C) 
      93     if not p.damp is None: r += p.damp*eye(x.size) 
    93 94     return r