lxTc"@ sddlmZdddddddd d d d d ddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~ddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddd d d d d ddddddddddddddddddd d!d"d#g"Zdd$lZdd%lmZdd&lmZmZd'e fd(YZ d)e fd*YZ e d+krd,GHej ed- Zej eeZx:eeeD]&Zeeeeeee8ӵ?g ߵ?g,6?gMTO[d(?gśf>?g)b?g8:#s?g$8?gSn?gL.?ge?gҥz?g1?g ؞K?g ?g*?gQ:?gH^;?gS&~TO?gϗv?gk}?g^3?g#"◷?ga?gRι?gnT?gim]з?gOr?g?ȷ?gU3"?g$?ggW.?g.,B1?gHA?g?/ğH?gOv?go,?g♰?gWƊ?gB?gҿBة?ggjT{?g!\?g̏o?grCo9ظ?g0x.Sݸ?gc?g~*0O?g ?gZw?gJ+M?g߷+?g{?g-e ?g?gK+ ?ggiiRttarget_evaluatorR(RRtnRRRRRRtltutinttnbdtnumber_of_lbfgs_iterationstnumber_of_function_evaluationsRtrunt minimizert run_lbfgsbR (R RRtminimizer_type((sT/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/lbfgsb/boost_python/tst_fit_cdf.pyt__init__\s    icC s|jj}tjd|ddd|jd|jd|jdtdd d d d | }|j\}}xtr|j |j||r|j\}}qo|j r|j d |j d|j }|jd|qo|jr|jddPqo|jddqoW|S(NR#tmiR$R%R'tenable_stp_inittfactrgcAtpgtolgh㈵>tiprinttstptcsdtvaluetsuffixs FINALR+(RtsizeRR+R$R%R'tTruet compute_functional_and_gradientstprocesstrequests_stp_initt adjust_stpt relative_step_length_line_searchtcurrent_search_directiont$set_relative_step_length_line_searcht is_terminatedtcallback_after_step_no_countingtcallback_after_stepR (tOR3R#R+Rtgtnew_stp((sT/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/lbfgsb/boost_python/tst_fit_cdf.pyR,us2          cC s||dkst|j|jjko8dknsCtdt||d}t||d}t||S(Ngig?i(RR8Rtabstmin(R R4R5t large_shifttmax_stp((sT/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/lbfgsb/boost_python/tst_fit_cdf.pyR=s 1tcC s#|jd7_|jd|dS(NiR7(R(RB(RDR+R7((sT/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/lbfgsb/boost_python/tst_fit_cdf.pyRCscC sdS(N((RDR7((sT/public2/ftp/pub/ccp4/6.5/ccp4-6.5/lib/py2/scitbx/lbfgsb/boost_python/tst_fit_cdf.pyRBscC s:|jd7_|j|jtj|j}x4t|jD]#}|jd|j|||s.o Hd $