U Pc9@s>dgfddZddZddZddZdd d d gfd d ZdS)Nc Cs dd|D}|}d}d} tt|D]} t|| || || d|| d|| <||f|}| d7} ||krz|}q&||  || <t|| || || d|| d|| <||f|}| d7} ||kr|}q&|| || <q&tt|D]} || || <q|| fS)a ## given a point, look for a better one nearby one coord at a time f is a function that takes a list of floats (of the same length as point) as an input args is a dict of any additional arguments to pass to f delta, and point are same-length lists of floats prevbest is a float point and delta are both modified by the function cSsg|]}|qSr.0xrrJ/home/sensei/ce-5319-webroot/MLBE4CE/chapters/optimization/hookeNjeeves.py sz&_hooke_best_nearby..)rangelen_value_in_bounds) fdeltapointZprevbestboundsargszZminfZftmpZfevirrr_hooke_best_nearbys( ,,rcCs`tt|D]N}||||dkr6||d||<q ||||dkr ||d||<q dS)z? shifts the point so it is within the given bounds rr Nr r )rrrrrr_point_in_boundss rcCsLd}tt|D]6}||||dkr.d}q||||dkrd}q|S)z< true if the point is in the bounds, else false TrFr r)rroutrrrr_is_point_in_boundssrcCs ||kr |S||krtS|SdSN)Zhight)vallowhighrrrr s r ?gư>ic st}d|d<d|d<dgt|}dgt|} |dkrJdd|D}ndd|D}d d|D}d} d } d } |D]X} | d dkrtd | d <nt| d | d <| d dkrtd | d <qvt| d | d <qvzzPt||dd|D}dd|D}tt|D]0}t||||<||dkr ||<q }||f|}| d 7} |}|} | |kr||kr| d 7} dd|D}t||||||\}}| |7} d }||kr|d kr|} tt|D]}||||krt|| ||<nt||||<||}||||<t|||||||d ||d ||<q|}t||||||\}}| |7} ||krqd }tt|D]<}d }t||||dt||krqnd }qq||kr^||kr^|}fdd|D}q^tt|D]}||| |<qWnTtk rz}z4t \}}}d|d<t |dt |j |d<W5d}~XYnXW5| |d<| |d<| |d<| |d<X|S)a In this version of the Hooke and Jeeves algorithm, we coerce the function into staying within the given bounds. basically, any time the function tries to pick a point outside the bounds we shift the point to the boundary on whatever dimension it is out of bounds in. Implementing bounds this way may be questionable from a theory standpoint, but that's how COPASI does it, that's how I'll do it too. TsuccessmessagerNcSsg|] }ddgqSrrrrrrrszhooke..cSsg|]}|d|dgqS)rr rrrrrrscSsg|]}|qSrrrrrrrsrz-infr infnitZfevalsfunrcSsg|]}|qSrrrrrrr scSsg|]}|qSrrrrrrr scSsg|]}|qSrrrrrrr"srcsg|] }|qSrrrrhorrrHsFz. line number: ) dictr floatrr absrr Exceptionsysexc_infostr tb_lineno)r Zstartptrr$epsilonZitermaxrresultrZendptZfminZnfevZitersboundZxbeforeZnewxrZ steplengthZfbeforeZnewfZevalsZkeeptmpreexc_typeZexc_objexc_tbrr#rhookes      4 & 0 r4)rrrr r4rrrrs+