&]\c@`sdZddlmZmZmZddlZddlmZddl m Z m Z m Z m Z mZmZedddd Zedddd d eZd efd YZdS( s} Unit tests for Krylov space trust-region subproblem solver. To run it in its simplest form:: nosetests test_optimize.py i(tdivisiontprint_functiontabsolute_importN(tget_trlib_quadratic_subproblem(tassert_tassert_array_equaltassert_almost_equalt assert_equaltassert_array_almost_equaltassert_array_lesst tol_rel_ig:0yE>t tol_rel_bgư>tdisptTestKrylovQuadraticSubproblemcB`s5eZdZdZdZdZdZRS(c `stjdddgdddgdddggtjdddgd}tddd d d fd d ddfd}|j|\}}t|tjdddgt|tttjj j ||j |dttjj ||d}|j|\}}t|tjdddgt|tttjj j ||j |dttjj ||dS(Ng?gg@g@g@g@txitfuncS`sdS(Ni((R((sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyt$stjacc`sS(N((R(tg(sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyR%sthesscS`sdS(N(tNone(R((sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyR&sthesspc`s j|S(N(tdot(Rty(tH(sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyR'sgg?gb5m1ݿgp돬ӳȿ( tnptarraytKrylovQPtsolveRRtTrueRtlinalgtnormRtlam(tselft trust_radiustsubprobtpt hits_boundary((RRsKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyttest_for_the_easy_cases4     ' 'c `stjdddgdddgdddggtjdddgd}tdddd d fd d d dfd}|j|\}}t|tjdddgttjjj||j |dttjj||d}|j|\}}t|tjdddgttjjj||j |dttjj||dS(Ng?gg@g@g@RiRcS`sdS(Ni((R((sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyRNsRc`sS(N((R(R(sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyROsRcS`sdS(N(R(R((sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyRPsRc`s j|S(N(R(RR(R(sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyRQsgg?g( RRRRRRRRRR (R!R"R#R$R%((RRsKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyttest_for_the_hard_case@s.    ''c `s?tjdddddgddddd gddd d d gddd d dgdd d ddggtjdddddgd}tdddddfddddfd}|j|\}}ttjjj||j|d t |d!d"d#d$d%gt |t dS(&Ngno?gu?g,wn)?g/Kg뚘G?g켮wO@g٢yn*?gr<8ٿg&Î?g6 @g ĽCԿg8K?gc?g >ڿg-Q?g*g?RRcS`sdS(Ni((R((sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyRsRc`sS(N((R(R(sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyRsRcS`sdS(N(R(R((sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyRsRc`s j|S(N(R(RR(R(sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyRsggiG7㰱?gl$WgZg-O|οgm/( RRRRRRRRR RR(R!R"R#R$R%((RRsKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyttest_for_very_close_to_zeros*   '  c `stjd tjdddddgd}tdddddfd d d d fd }|j|\}}|j\}}t|jdt|dS(Niigư>g?RRcS`sdS(Ni((R((sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyRsRc`sS(N((R(R(sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyRsRcS`sdS(N(R(R((sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyRsRc`s j|S(N(R(RR(R(sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyRss TR Solving trust region problem( RteyeRt KrylovQP_dispRt readouterrRt startswithtrepr(R!tcapsysR"R#R$R%toutterr((RRsKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyt test_disps   (t__name__t __module__R&R'R)R*R3(((sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyR s  * '  (t__doc__t __future__RRRtnumpyRtscipy.optimize._trlibRt numpy.testingRRRRRR RRR,tobjectR (((sKlib/python2.7/site-packages/scipy/optimize/tests/test_trustregion_krylov.pyts .