lbfgs_jpl_version/lsopt.2/liblsopt_ecco.a

!<arch>
/               1361899272  0     0     0       1106      `
5√√√√√√√√СРСРСРСРСРСРСРііііііі]]]]]]]|.├J├J├J├J├J├J├J∙Ў∙Ў∙Ў╡З╡З╡Зечечьё(ґ≤ґ≤ґ≤ґ≤ґ≤сНmpi_parm_ph_opti_lsopt_top_mpi_sgi_private_status_mpi_sgi_private_unweighted_mpi_sgi_private_mpi_sgi_private_inplace_mpi_sgi_private_char_mpi_parm_lsupdxx_mpi_sgi_private_status_mpi_sgi_private_unweighted_mpi_sgi_private_mpi_sgi_private_inplace_mpi_sgi_private_char_mpi_parm_lsline_mpi_sgi_private_status_mpi_sgi_private_unweighted_mpi_sgi_private_mpi_sgi_private_inplace_mpi_sgi_private_char_mpi_parm_hessupd_mpi_sgi_private_status_mpi_sgi_private_unweighted_mpi_sgi_private_mpi_sgi_private_inplace_mpi_sgi_private_char_cubic_dgscale_mpi_sgi_private_status_mpi_sgi_private_unweighted_mpi_sgi_private_mpi_sgi_private_inplace_mpi_sgi_private_char_mpi_parm_mpi_parm_xxstore_instore_mpi_parm_xxstore_dostore_xxstore_outstore_lswri_distwork_pddot_mpi_sgi_private_status_mpi_sgi_private_mpi_sgi_private_inplace_mpi_sgi_private_char_lsupdxx_1_lsopt_top.o/    1361899271  10100 26113 100644  61216     `
ELF> К@@ЯЪ
#(L|▄0²00ў`0©░%Г╦%ЬЮ% %0%+X%<─%M╗%^п%o~@▐,═цтЬ#"░'h,8I,ZDkt |■&█╪&Є`мxС▄Д,)P"9пI╟Y■+iю.yП.┴ .≥d'╘═люэпЛЛЭЭ(,D<TLl\└x▄)┬╦&·Єў(ЎднЮ$ч°У,,FX%V─%f╗%v<├\√|іЬ!І° е*т╗ ЦH%РЄ юЮ  p'0Л @ЬP `≤*p  ─д)░, ═П.╟8 ю &пD ЮH!ПP \|  ┬0п%?TPha|r@'┐эуЭЕє У╟)<&%П5тP`ahvp▀x═─╣┬й░ъ≤Т═         ╗ ╟3 ╦H ю] хr п┤ ь° Ю╠ Хф Пш ЬП 

/
D
 Y
(n
0┐
8≤
@ґ
Hб
Pв
XЛ
`hp+x@─U┬j░≤■═╘╗Ў╟с╦ХюЩхл#п4тEьZЮoХ└П≥ЬўцьМ
 
(,
0A
8V
@k
H─
P∙
X╙
`©
hт
pИ
xЧ
─┬(░=≤R═g╗|╟▒╦ію╩хппЕьЖэЮХ-ЛBТWЭl│√╚ю$у,Й4Ъ<D)L>TS\hd}l▓tї|╪└я▄Ф■Ш°є%╛:ЄO╪dдyл▌тёэ╦ДмЛБТВЭ!6K`$u,┼4÷<ЄDиLчTС\dl2tG|\└q▄├■⌡°╟єе╛зЄО╪дл.тCэXДmЛ┌Т≈Ь╛ажК (*,?4T<iD~L⌠T╗\Ґ`рhГpЭtx&─;┬P░e≤z═▐╗є╟╧╦нюЦхЬп
ь"Ю7ХLПaЬv▀═╣й ъ(Т0 8@3HHP]Xr`┤h°p╠xф─ш┬П░≤═/╗D╟Y╦╚╩   е
лв
ДС.РЪ8РЪA0Ojп╠ЕЭдщД╧l≤Н'<nРЪ├РЪ╒РЪЁРЪлРЪБТ+DUk┐⌠6H76H76H76H$76H       76H 76H 76H$ 76H$ 76H$ 76H$ 76H$ 76H$76H$76H$76H 76H 76H 76H76H76H76H76H$76H$76H$76H$76H$76H76$
$
8
976H76$
$
8
976H$ 76H$ 76H 76H76H76H76H76$H76$H76$H76$H76$H76$H76$H76$H76$H76$H76$H76$H76$H76H$     7UH┴ЕH│Лp
I╩L┴}хI╨L▀}hH╦ _═BL┴e╟I┴лH┴╣hЖЪЪHЎH▀M@H┴]пH┴ШH┴∙ЖЪЪ3рA┴РAD▀9L┴uюH┴E┬H┴E≤H╦П?L┴m╦L┴█═ЖЪЪL┴┘╗ЖЪЪD▀3гEьЪЪЪЪ┴UэHгE░H┴E═A┴R┴∙pЖЪЪР┘╟ЖЪЪ┴UЮH┴E╗┴UД┴UХгEЛ┴UПD┴}ТгEЬ┐>┘╒
H▀E(D▀(E┘М▌пI╧I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘XЧЪЪI╨Hг┘ЩЪЪHг┘XЧЪЪ/L┴∙`ЧЪЪAЪсI╩I╨L┴∙pЧЪЪI╨ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юM█KHг┘ЩЪЪL█┘hЧЪЪHг┘hЧЪЪ/AЪрI╩I╨L┴∙─ЧЪЪI╨ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юM█K0Hг┘ЩЪЪL█┘xЧЪЪHг┘xЧЪЪ/AЪрI╩I╨L┴∙░ЧЪЪI╨ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юM█KHHг┘ЩЪЪL█┘┬ЧЪЪHг┘┬ЧЪЪ$AЪрH╧H█ҐЩЪЪHЎH█∙@ЪЪЪ3юг┘@ЪЪЪXскQЪяI╩I╨L┴∙═ЧЪЪI╨ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юM█K`Hг┘ЩЪЪL█┘≤ЧЪЪHг┘≤ЧЪЪ$AЪрH▀MI╧HЎH█ҐЩЪЪ3юH█∙HЪЪЪL▀L┴┘HЪЪЪAЪяI╩I╨L┴∙╟ЧЪЪI╨ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юM█KxHг┘ЩЪЪL█┘╗ЧЪЪHг┘╗ЧЪЪ$AЪрH▀M I╧HЎH█ҐЩЪЪ3юH█∙PЪЪЪL▀L┴┘PЪЪЪAЪяI╩I╨L┴∙юЧЪЪI╨ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юM█▀░Hг┘ЩЪЪL█┘╦ЧЪЪHг┘╦ЧЪЪ$AЪрH▀MI╧HЎH█ҐЩЪЪ3юH█∙XЪЪЪL▀L┴┘XЪЪЪAЪяI╩I╨L┴∙пЧЪЪI╨ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юM█▀╗Hг┘ЩЪЪL█┘хЧЪЪHг┘хЧЪЪ$AЪрH▀M0I╧HЎH█ҐЩЪЪ3юH█∙`ЪЪЪD▀D┴┘`ЪЪЪAЪяI╩I╨L┴∙ЮЧЪЪI╨ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юM█▀юHг┘ЩЪЪL█┘ьЧЪЪHг┘ьЧЪЪ$AЪрH▀M8I╧HЎH█ҐЩЪЪ3юH█∙hЪЪЪD▀D┴┘hЪЪЪAЪяI╩I╨L┴∙ПЧЪЪI╨ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юM█▀ьHг┘ЩЪЪL█┘ХЧЪЪHг┘ХЧЪЪ$AЪрH╧H█ҐЩЪЪHЎH█∙pЪЪЪ3юD┴ґpЪЪЪЪяI╩I╨L┴∙ЪЪЪI╨ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юM█▀ПHг┘ЩЪЪL█┘ЬЧЪЪHг┘ЬЧЪЪ$AЪрH╧H█ҐЩЪЪHЎH█∙xЪЪЪ3юD┴ҐxЪЪЪЪяI╩I©ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘ЪЪЪI╨M█▀Hг┘ЩЪЪHг┘ЪЪЪ$L┴∙ЪЪЪAЪвI╦H█ҐЩЪЪHЎH█U─3юC█v┴M─AЪпE┘Ж▌в>H▀E0D▀8E┘Ъ▄>H▀E8D▀8E┘Ъ▌D=H▀EfОюf/┌аE┘М▌ёI╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘hЩЪЪH╩M█┼hHг┘ЩЪЪHг┘hЩЪЪH┴²pЩЪЪAЪсH▀MI╦HЎH█ҐЩЪЪ3юH█∙ЬЩЪЪH▀H┴²ЬЩЪЪAЪпH▀EhA╪гИ0H▀E f/┌аE┘М▌ёI╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘хШЪЪH╩M█┼─Hг┘ЩЪЪHг┘хШЪЪH┴²пШЪЪAЪсH▀M I╦HЎH█ҐЩЪЪ3юH█∙`ЩЪЪH▀H┴²`ЩЪЪAЪпH▀EhA╪гИE/H▀E H╨Рf/├аE┘М▌ёI╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘xЗЪЪH╩M█┼≤Hг┘ЩЪЪHг┘xЗЪЪH┴²─ЗЪЪAЪсH▀M I╦HЎH█ҐЩЪЪ3юH█∙юШЪЪH▀H┴²юШЪЪAЪпH▀EhA╪гИh.I╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ШЪЪ3юI╧Hг┘ЩЪЪHг┘ШЪЪL┴█ШЪЪAЪрI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ЧЪЪ3юI╧Hг┘ЩЪЪHг┘ЧЪЪ"L┴█ ЧЪЪAЪрH╧H█ҐЩЪЪHЎH█∙0ЪЪЪ3юг┘0ЪЪЪЪЪЪЪЪяH╧H█ҐЩЪЪHЎH█∙8ЪЪЪ3юг┘8ЪЪЪЪяH┴ъL█U═ARI╨L█]ЮASL█mХAUL█}ДAWH█EьP3юH█UПRH█u┬H█U╗H█MЬL█EТL█MЛAЪрH┐д0H╦D▀8E┘Ъ┘I╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ЧЪЪ3юI╧Hг┘ЩЪЪHг┘ЧЪЪ!L┴█ЧЪЪAЪр▀MьI╦HЎH█ҐЩЪЪ3юH█∙ ЪЪЪ┴█ ЪЪЪAЪп▀MДI╦HЎH█ҐЩЪЪ3юH█∙(ЪЪЪ┴█(ЪЪЪAЪпH╨H▀MhH▀E┬H┴B┐9└▐H▀E(D▀(E┘М~nI╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘PЩЪЪH╩M█┼╟Hг┘ЩЪЪHг┘PЩЪЪH┴²XЩЪЪAЪсH▀EhA╪гИ│+┐}Ь└▐H▀E(D▀(E┘М~nI╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘╟ШЪЪH╩M█┼╪Hг┘ЩЪЪHг┘╟ШЪЪ+H┴²╦ШЪЪAЪсH▀EhA╪гИХ*┐}Т▌g9КH╨H▀E┬H┴BЖEь└▄E┘Ъ┘бI╩IҐЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘ьШЪЪI╨M█▀тHг┘ЩЪЪHг┘ьШЪЪL┴∙ЮШЪЪAЪуI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ХШЪЪ3юI╧Hг┘ЩЪЪHг┘ХШЪЪL┴█ПШЪЪAЪрH╧H█ҐЩЪЪHЎH█∙┬ЩЪЪ3юг┘┬ЩЪЪXскQЪяI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ЭЪЪ3юI╧Hг┘ЩЪЪHг┘ЭЪЪL┴█ЭЪЪAЪрH▀█hЖЪЪI╧HЎH█ҐЩЪЪ3юH█∙≤ЩЪЪL▀L┴┘≤ЩЪЪAЪяI▀$I╦HЎH█ҐЩЪЪ3юH█∙╗ЩЪЪH┴█╗ЩЪЪAЪпг┘ЪЪЪH█}эH▀█hЖЪЪH█╣ЪЪЪL▀┘ЖЪЪH┴зM┴АЪ∙╗ЖЪЪH╦▀┘р┴∙░УЪЪ┘UI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ЬШЪЪ3юI╧Hг┘ЩЪЪHг┘ЬШЪЪL┴█ЭЪЪAЪрH╧H█ҐЩЪЪHЎH█∙░ЩЪЪ3юг┘░ЩЪЪXскQЪяH▀█ЖЪЪI╧HЎH█ҐЩЪЪ3юH█∙═ЩЪЪL▀L┴┘═ЩЪЪAЪяI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ЭЪЪ3юI╧Hг┘ЩЪЪHг┘ЭЪЪ%L┴█ ЭЪЪAЪрH▀█hЖЪЪI╧HЎH█ҐЩЪЪ3юH█∙╟ЩЪЪL▀L┴┘╟ЩЪЪAЪяH▀█hЖЪЪI╧HЎH█ҐЩЪЪ3юH█∙╦ЩЪЪL▀AL┴┘╦ЩЪЪAЪяI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘(ЭЪЪ3юI╧Hг┘ЩЪЪHг┘(ЭЪЪ%L┴█0ЭЪЪAЪрI▀$I╦HЎH█ҐЩЪЪ3юH█∙юЩЪЪH┴█юЩЪЪAЪпI▀L$I╦HЎH█ҐЩЪЪ3юH█∙хЩЪЪH┴█хЩЪЪAЪп▀3┘Ж▌ґHcфH┐Ь▄X5H▀UXH┐Б┘рtЖб┘C5jZ┴яH█yH;г▄15Взж┐Б+Р3рHcЖH┘иvH▀}XI╦П?L┴вHЪбH;яrТH╨(H▀UX)й)Dй)Dй )Dй0H┐аH;нrДH;ПsH▀UXH╧П?H┴РHЪфH;ПrТ┐Ґ░УЪЪu]I╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘8ЭЪЪ3юI╧Hг┘ЩЪЪHг┘8ЭЪЪL┴█@ЭЪЪAЪрH╧I╦H┴ъH╨3юH▀╣hЖЪЪAЪпI╦H┴ъL┴ФH╨H╧3юAЪпI╦H┴ъH╨H╧3юH▀uXAЪпH╦┐8u]I╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘HЭЪЪ3юI╧Hг┘ЩЪЪHг┘HЭЪЪL┴█PЭЪЪAЪрH╧H┴ъL┴ФH╨3юЪяР┘(ЧЪЪH┐дПРYюH╨I╩H█ҐЧЪЪH╧H█u╗I╦I┴я3юL█∙ЪЪЪР┘ЧЪЪL┴$AЪсH┐дРQE╗H╦РE╗▀┘р┴∙░УЪЪ┘▐I╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘XЭЪЪ3юI╧Hг┘ЩЪЪHг┘XЭЪЪL┴█`ЭЪЪAЪрH▀M╗I╦HЎH█ҐЩЪЪ3юH█∙пЩЪЪH┴█пЩЪЪAЪпH╦Рf/E╗vH▀EhгИ?H▀┘ЖЪЪH╧┐Ґ░УЪЪH▀H┴Q┘uI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘┬ЗЪЪ3юI╧Hг┘ЩЪЪHг┘┬ЗЪЪ+L┴█░ЗЪЪAЪрH╧I╧HЎH█ҐЩЪЪ3юH█∙hЭЪЪL▀AL┴┘hЭЪЪAЪяИвH▀U@▀EТ┴┘≤УЪЪ;└©H▀E(D▀(E┘М▌═I╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘`ЫЪЪH╩M█┼ЬHг┘ЩЪЪHг┘`ЫЪЪ!H┴²hЫЪЪAЪс▀█≤УЪЪH╩HЎH█ҐЩЪЪ3юH█∙`ЗЪЪ┴█`ЗЪЪЪсH▀EhгИG!▀EЛ┘ю┴┘░УЪЪ▌щ6▀┘≤УЪЪ;┘░УЪЪ█©H▀E(D▀(E┘М▌═I╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘┬ЬЪЪH╩M█┼(Hг┘ЩЪЪHг┘┬ЬЪЪH┴²░ЬЪЪAЪс▀█░УЪЪH╩HЎH█ҐЩЪЪ3юH█∙XЫЪЪ┴█XЫЪЪЪсH▀EhгИe ▀EП┴┘░УЪЪ;┘≤УЪЪ▌©H▀E(D▀(E┘М▌═I╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘hЖЪЪH╩M█┼@Hг┘ЩЪЪHг┘hЖЪЪH┴²pЖЪЪAЪс▀█░УЪЪH╩HЎH█ҐЩЪЪ3юH█∙─ЬЪЪ┴█─ЬЪЪЪсH▀EhгИ▒E┘Ъ┘ I╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘╟УЪЪ3юI╧Hг┘ЩЪЪHг┘╟УЪЪ*L┴█╦УЪЪAЪрI╧H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧E3ю3юHг┘ЩЪЪAЪяH╧I╦H┴ъH╨3юH▀╣hЖЪЪAЪпI╦H┴ъL┴ФH╨H╧3юAЪпHЎH╧H▀┘ЖЪЪH▀U┬H┴H┴Q┐>┘%I╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘юУЪЪ3юI╧Hг┘ЩЪЪHг┘юУЪЪ-L┴█хУЪЪAЪрH▀█ЖЪЪI╧HЎH█ҐЩЪЪ3юH█∙xЖЪЪL▀L┴┘xЖЪЪAЪяI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘пУЪЪ3юI╧Hг┘ЩЪЪHг┘пУЪЪ-L┴█ьУЪЪAЪрH▀█hЖЪЪI╧HЎH█ҐЩЪЪ3юH█∙─ЖЪЪL▀L┴┘─ЖЪЪAЪяH▀█hЖЪЪI╧HЎH█ҐЩЪЪ3юH█∙┬ЖЪЪL▀AL┴┘┬ЖЪЪAЪяI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ЮУЪЪ3юI╧Hг┘ЩЪЪHг┘ЮУЪЪ-L┴█ХУЪЪAЪрI▀$I╦HЎH█ҐЩЪЪ3юH█∙░ЖЪЪH┴█░ЖЪЪAЪпI▀L$I╦HЎH█ҐЩЪЪ3юH█∙≤ЖЪЪH┴█≤ЖЪЪAЪпЖEь┘ш.H▀E0▀г┘pЖЪЪ┘ю▌sD┴╣ЖЪЪ┴┘ЖЪЪD▀╣pЖЪЪL▀ґЖЪЪH╦┐8┘ГI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘юЭЪЪ3юI╧Hг┘ЩЪЪHг┘юЭЪЪ&L┴█хЭЪЪAЪрI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘пЭЪЪ3юI╧Hг┘ЩЪЪHг┘пЭЪЪL┴█ьЭЪЪAЪрH╧H█ҐЩЪЪHЎH█∙ПЩЪЪ3юD┴╣ПЩЪЪЪяHcH┘р~!H┐З▌3*HаБH╦L┴ФH▀}PЪпH╦H┴аI▀UH┴∙PЧЪЪD┴0▀Eь╗H┴Q┘╙,╗┘БЪuXЪ╣hЖЪЪATЪuHЪuH▀Ґ╗ЖЪЪH█E═PH▀UhL█UХARL▀M8L█]ДASH┴чL█ҐЮЭЪЪAWH╧L█EЮAUЪ∙═ЖЪЪH┐дPH╦┐8┘xI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘юЫЪЪ3юI╧Hг┘ЩЪЪHг┘юЫЪЪL┴█хЫЪЪAЪрI╦H█ҐЩЪЪHЎH█∙пЫЪЪ3юH╧Hг┘пЫЪЪH┴█ьЫЪЪAЪпH▀MhI╧HЎH█ҐЩЪЪ3юH█∙ ШЪЪD▀D┴┘ ШЪЪAЪяI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ПЫЪЪ3юI╧Hг┘ЩЪЪHг┘ПЫЪЪL┴█ЬЫЪЪAЪрI╦H█ҐЩЪЪHЎH█∙ЗЪЪ3юH╧Hг┘ЗЪЪH┴█ЗЪЪAЪпD▀}ДH╧HЎH█ҐЩЪЪ3юH█∙0ШЪЪD┴Ґ0ШЪЪЪяI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘0ЗЪЪ3юI╧Hг┘ЩЪЪHг┘0ЗЪЪL┴█8ЗЪЪAЪрI╦H█ҐЩЪЪHЎH█∙@ЗЪЪ3юH╧Hг┘@ЗЪЪH┴█HЗЪЪAЪп▀MХI╦HЎH█ҐЩЪЪ3юH█∙HШЪЪ┴█HШЪЪAЪпКD▀}ДAВг┘Ы┐}Ю┘ОH╦┐8┘]I╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ЮЫЪЪ3юI╧Hг┘ЩЪЪHг┘ЮЫЪЪL┴█ХЫЪЪAЪрI╦H█ҐЩЪЪHЎH█∙ЗЪЪ3юH╧Hг┘ЗЪЪH┴█ЗЪЪAЪпH╧H█ҐЩЪЪHЎH█∙(ШЪЪ3юD┴╣(ШЪЪЪяI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ ЗЪЪ3юI╧Hг┘ЩЪЪHг┘ ЗЪЪL┴█(ЗЪЪAЪрH▀█hЖЪЪI╧HЎH█ҐЩЪЪ3юH█∙8ШЪЪL▀L┴┘8ШЪЪAЪяH▀█hЖЪЪI╧HЎH█ҐЩЪЪ3юH█∙@ШЪЪL▀AL┴┘@ШЪЪAЪяI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘PЗЪЪ3юI╧Hг┘ЩЪЪHг┘PЗЪЪL┴█XЗЪЪAЪрI▀$I╦HЎH█ҐЩЪЪ3юH█∙PШЪЪH┴█PШЪЪAЪпI▀L$I╦HЎH█ҐЩЪЪ3юH█∙XШЪЪH┴█XШЪЪAЪпI╦H┴ъH╨H╧3юH▀╣hЖЪЪAЪпI╦H┴ъL┴ФH╨H╧3юAЪпH╧H┴ъL┴ФH╨3юЪяР┘(ЧЪЪH┐дПРYюH╨I╩H█ҐХЭЪЪH╧H█╣ПЭЪЪI╦I┴я3юL█∙ЪЪЪР┘ХЭЪЪL┴$AЪсH┐дРQ┘ПЭЪЪH╧H┴ъH╨3юH▀╣hЖЪЪР┘ПЭЪЪЪяР┘(ЧЪЪL█∙ЪЪЪРYюH╨ARI╩H█ҐЬЭЪЪH╧H█╣0ЧЪЪI╦I╧3юRР┘ЬЭЪЪAЪсH┐дH╧H┴ъH╨3юH▀uHЪяР┘(ЧЪЪL█∙ЪЪЪРYюH╨ARI╩H█ҐЩЪЪH╧H█╣8ЧЪЪI╦I╧3юRР┘ЩЪЪAЪсH┐дH╦┐8┘QH▀E(┐8▌╠I╩I©ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘░ЫЪЪI╨M█▀Hг┘ЩЪЪHг┘░ЫЪЪ(L┴∙≤ЫЪЪAЪвI╦H█ҐЩЪЪHЎH█∙═ЫЪЪ3юH╧Hг┘═ЫЪЪ%H┴█╗ЫЪЪAЪпI╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘`ШЪЪHг┘ЩЪЪM█┼D┴╣`ШЪЪAЪсH▀M═I╦HЎH█ҐЩЪЪ3юH█∙hШЪЪH┴█hШЪЪAЪп▀MЮI╦HЎH█ҐЩЪЪ3юH█∙pШЪЪ┴█pШЪЪAЪпI▀MI╦HЎH█ҐЩЪЪ3юH█∙xШЪЪH┴█xШЪЪAЪпH▀█ПЭЪЪI╦HЎH█ҐЩЪЪ3юH█∙─ШЪЪH┴█─ШЪЪAЪпРQ┘0ЧЪЪH╧H█ҐЩЪЪHЎH█∙┬ШЪЪ3юР┘┬ШЪЪЪяРQ┘8ЧЪЪРYE═H█ҐЩЪЪH╧HЎ3юH█∙░ШЪЪР┘░ШЪЪЪяH▀Eh▀┐Ь        uD┴╣pЖЪЪИ3     ┐ЬuD┴╣pЖЪЪИ"        ┐Ьu
H▀Ehг▀3┘Ж▌фHcнH┐Ы▄НH▀EXH┐Ю┘юt╗┘зjX┴бH█zH;о▄хВь3Ъф┐ЮВьфHcюH┘рv$РE═L▀EXL▀MH(хРAYЫРAЬHЪгH;ЗrИH▀}HРE═fюL█вIВюtzH▀}XL▀EH░H█t&H█©РAпРATпРA\п РAdп0fALпfATпfA\п(fAdп8fYхfYпfYьfYЮ)в)Tв)\в )dв0H┐бH;пr²КNH▀}XL▀EHH┴ЖA(пA(TпA(\п A(dп0fYхfYпfYьfYЮ)в)Tв)\в )dв0H┐бH;пrҐH;аs"РE═H▀UXH▀}H(хРYгРбHЪюH;аrКH┐Ы▄┐H▀UPH┐Б┘рtЖб┘njZ┴пH█xH;о▄\Взж┐Б+Р3рHcЖH┘юvH▀}PРAтР\вРвHЪбH;пrХI█дHВбtuH▀UPH█t&РAдРALдРATд РA\д0fADдfALдfATд(fA\д8f\бf\Lбf\Tб f\\б0)б)Lб)Tб )\б0H┐юH;фr√КNH▀UPA(дA(LдA(Tд A(\д0f\бf\Lбf\Tб f\\б0)б)Lб)Tб )\б0H┐юH;фrІH;ЯsH▀EPРAТР\ПРПHЪфH;ЯrХH╨I╧H┴ъI┴п3юH▀uPH▀MXAЪяH┐дПH╨I╩H█Ґ(ЧЪЪH╧H█╣≤ЗЪЪI╦I┴я3юL█∙ЪЪЪР┘(ЧЪЪL┴$AЪсH┐дfОюf/┘≤ЗЪЪ┐гH▀E@▀┘р┴∙ЖЪЪ▌D▀UП┘р█B┴┘ЬУЪЪ┘▌H╨┴EП┐:┘I╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘8ЖЪЪ3юI╧Hг┘ЩЪЪHг┘8ЖЪЪL┴█@ЖЪЪAЪрI╦H█ҐЩЪЪHЎH█∙HЖЪЪ3юH╧Hг┘HЖЪЪ'H┴█PЖЪЪAЪпI╦H█ҐЩЪЪHЎH█∙XЖЪЪ3юH╧Hг┘XЖЪЪH┴█`ЖЪЪAЪпH╧H█ҐЩЪЪHЎH█∙XЬЪЪ3юD┴╣XЬЪЪЪя▀MЛI╦HЎH█ҐЩЪЪ3юH█∙`ЬЪЪ┴█`ЬЪЪAЪп▀█ЬУЪЪI╦HЎH█ҐЩЪЪ3юH█∙hЬЪЪ┴█hЬЪЪAЪпИ⌡;┘ЖЪЪ~+┘ЖЪЪ┴┘ЬУЪЪ┴EПК┴EП▀EЛ┴┘ПУЪЪ;┘ЬУЪЪ┘jH╦┐8┘3I╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘░УЪЪ3юI╧Hг┘ЩЪЪHг┘░УЪЪ#L┴█≤УЪЪAЪрI╦H█ҐЩЪЪHЎH█∙═УЪЪ3юH╧Hг┘═УЪЪH┴█╗УЪЪAЪпH╧H█ҐЩЪЪHЎH█∙ ЖЪЪ3юD┴╣ ЖЪЪЪя▀█ПУЪЪI╦HЎH█ҐЩЪЪ3юH█∙(ЖЪЪ┴█(ЖЪЪAЪп▀█ЬУЪЪI╦HЎH█ҐЩЪЪ3юH█∙0ЖЪЪ┴█0ЖЪЪAЪп▀┘ПУЪЪЪю┴┘ПУЪЪ;┘ЖЪЪ~+┘ЖЪЪ┴EЛК┴EЛI╦I╧H┴ъH█╣╦ЫЪЪH╧3юH▀UXРAР^┘≤ЗЪЪРQюР┘╦ЫЪЪAЪяI╦H┴ъH╧H█╣╦ЫЪЪ3юH▀UPAЪпH╦┐8┘МI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘pЬЪЪ3юI╧Hг┘ЩЪЪHг┘pЬЪЪ+L┴█xЬЪЪAЪрH╧H█ҐЩЪЪHЎH█∙@ЫЪЪ3юD┴╣@ЫЪЪЪя▀MЛI╦HЎH█ҐЩЪЪ3юH█∙HЫЪЪ┴█HЫЪЪAЪпD▀}ПH╧HЎH█ҐЩЪЪ3юH█∙PЫЪЪD┴ҐPЫЪЪЪяКD▀}ПI╧H┴ъH╨H██ШЪЪ3юG█D?H▀uPD┴┘ШЪЪAЪяD▀EПI╨H┴ъH╨3юH██ШЪЪH▀uXG█LD┴█ШЪЪAЪрI╧H┴ъI╦3юH▀uPH▀UXH▀MHAЪяР█ПЭЪЪР^M╗H▀U H▀┘ПЭЪЪР█╟ЖЪЪРH┴E╗f/а┤'AЪфD;╣ЖЪЪ▌╞КЪЪD┴╣pЖЪЪH▀EhЪ█pЖЪЪгH▀E Р┘╟ЖЪЪЖEДР┘╨┐}Ю┘╟H┐д░L█ґPЧЪЪL▀∙ЖЪЪL█u╗L▀]L█ҐЮЭЪЪH▀UHH█}░L┴$H█E≤L┴\$L█E═L┴l$I╨L┴t$L┴|$ L┴d$(H┴T$0H╨H▀█hЖЪЪH▀uXH┴L$8H█MЛH┴t$@H┴|$HH┴ъL▀MH┴D$P3юL┴D$XL█EПL┴L$`H▀uhL▀M@AЪрH┐дpH╦┐8┘хI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘xЭЪЪ3юI╧Hг┘ЩЪЪHг┘xЭЪЪL┴█─ЭЪЪAЪрH▀MXI╧HЎH█ҐЩЪЪ3юH█∙ьЩЪЪL▀L┴┘ьЩЪЪAЪяH▀MXI╧HЎH█ҐЩЪЪ3юH█∙ЮЩЪЪL▀AL┴┘ЮЩЪЪAЪяH╨H┴ъ3юH▀uXЪр┐}ЮuВUДЖEьtгEДI©L█U═ARH┴ъL█]ЮAS3юL█mХAUH▀╣ЖЪЪL█uДAVH▀M@H█U╗L█EЛL█MПAЪвH┐д H▀E(┐8▄H╦┐8┘I╧H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧E3ю3юHг┘ЩЪЪAЪяIҐI╨L┴ґ╗ЗЪЪIҐЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘═ЗЪЪHг┘ЩЪЪM██╟Hг┘═ЗЪЪ$AЪр▀█pЖЪЪI╦HЎH█ҐЩЪЪ3юH█∙pЭЪЪ┴█pЭЪЪAЪпI╩M██хЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘╟ЗЪЪI╨Hг┘ЩЪЪHг┘╟ЗЪЪ$L┴∙╦ЗЪЪAЪс▀MЮI╦HЎH█ҐЩЪЪ3юH█∙┬ЭЪЪ┴█┬ЭЪЪAЪпI╩M██ЮЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘юЗЪЪI╨Hг┘ЩЪЪHг┘юЗЪЪ$L┴∙хЗЪЪAЪсH▀M I╧HЎH█ҐЩЪЪ3юH█∙░ЭЪЪL▀L┴┘░ЭЪЪAЪяH╧H┴ъH╨3юH▀╣hЖЪЪЪяР┘(ЧЪЪL█∙ЪЪЪРYюH╨ARI╩H█Ґ@ЧЪЪH╧H█╣0ЧЪЪI╦I╧3юRР┘@ЧЪЪAЪсH┐дH╧H┴ъL┴ФH╨3юЪяР┘(ЧЪЪH█²ЪЪЪРYюH╨SI╨H█ҐHЧЪЪH╧H█╣8ЧЪЪI╦I╧3юRР┘HЧЪЪAЪрH┐дH▀E(D▀(A┐Щ
▄*H╦┐8┘│H╩H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧E3ю3юHг┘ЩЪЪЪсH╩I╨H┴²ьЗЪЪH╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘пЗЪЪHг┘ЩЪЪL█▀ЬHг┘пЗЪЪAЪрH▀█ЖЪЪI╧HЎH█ҐЩЪЪ3юH█∙≤ЭЪЪL▀L┴┘≤ЭЪЪAЪяI╩L█▀ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘ЮЗЪЪI╨Hг┘ЩЪЪHг┘ЮЗЪЪL┴∙ХЗЪЪAЪсРQ┘0ЧЪЪH╧H█ҐЩЪЪHЎH█∙═ЭЪЪ3юР┘═ЭЪЪЪяЎЪЪЪЪL█▀(H╩H█ҐЩЪЪH╨┘┐ H╧3юL█┘ПЗЪЪI╨Hг┘ЩЪЪHг┘ПЗЪЪL┴∙ЬЗЪЪЪсРQ┘8ЧЪЪH╧H█ҐЩЪЪHЎH█∙╗ЭЪЪ3юР┘╗ЭЪЪЪяКH╦D▀8E┘Ъ┘TH▀EhD▀ E┘Д└DA┐Щ▄H╩H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧E3ю3юHг┘ЩЪЪЪсH╩I╨H┴²─ЩЪЪH╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘xЩЪЪHг┘ЩЪЪL█▀@Hг┘xЩЪЪ AЪрE┘Д▄CA┐Э┤ 
McДH╦BЪ$ЮЎЪЪЪЪL█▀▄H╩H█ҐЩЪЪH╨┘┐ H╧3юL█┘hВЪЪI╨Hг┘ЩЪЪHг┘hВЪЪL┴∙pВЪЪЪсH╧H█ҐЩЪЪHЎH█∙ПЬЪЪ3юD┴╔ПЬЪЪЪяH╩H█ҐЩЪЪHЎH█∙(ЬЪЪ3юH╧Hг┘(ЬЪЪ)H┴█0ЬЪЪЪсИ8 ЎЪЪЪЪL█▀\H╩H█ҐЩЪЪH╨┘┐ H╧3юL█┘XВЪЪI╨Hг┘ЩЪЪHг┘XВЪЪL┴∙`ВЪЪЪсH╧H█ҐЩЪЪHЎH█∙ХЬЪЪ3юD┴╔ХЬЪЪЪяH╩H█ҐЩЪЪHЎH█∙ЬЪЪ3юH╧Hг┘ЬЪЪ$H┴█ ЬЪЪЪсИaЎЪЪЪЪL█▀,H╩H█ҐЩЪЪH╨┘┐ H╧3юL█┘HВЪЪI╨Hг┘ЩЪЪHг┘HВЪЪL┴∙PВЪЪЪсH╧H█ҐЩЪЪHЎH█∙ЮЬЪЪ3юD┴╔ЮЬЪЪЪяH╩H█ҐЩЪЪHЎH█∙ЬЪЪ3юH╧Hг┘ЬЪЪH┴█ЬЪЪЪсИ┼ЎЪЪЪЪL█▀ЭH╩H█ҐЩЪЪH╨┘┐ H╧3юL█┘8ВЪЪI╨Hг┘ЩЪЪHг┘8ВЪЪL┴∙@ВЪЪЪсH╧H█ҐЩЪЪHЎH█∙ьЬЪЪ3юD┴╔ьЬЪЪЪяH╩H█ҐЩЪЪHЎH█∙ЬВЪЪ3юH╧Hг┘ЬВЪЪ&H┴█ЬЪЪЪсИЁЎЪЪЪЪL█▀лH╩H█ҐЩЪЪH╨┘┐ H╧3юL█┘(ВЪЪI╨Hг┘ЩЪЪHг┘(ВЪЪL┴∙0ВЪЪЪсH╧H█ҐЩЪЪHЎH█∙пЬЪЪ3юD┴╔пЬЪЪЪяH╩H█ҐЩЪЪHЎH█∙ХВЪЪ3юH╧Hг┘ХВЪЪH┴█ПВЪЪЪсИэЎЪЪЪЪL█▀°H╩H█ҐЩЪЪH╨┘┐ H╧3юL█┘ВЪЪI╨Hг┘ЩЪЪHг┘ВЪЪL┴∙ ВЪЪЪсH╧H█ҐЩЪЪHЎH█∙хЬЪЪ3юD┴╔хЬЪЪЪяH╩H█ҐЩЪЪHЎH█∙ьВЪЪ3юH╧Hг┘ьВЪЪ)H┴█ЮВЪЪЪсИЎЪЪЪЪL█▀lH╩H█ҐЩЪЪH╨┘┐ H╧3юL█┘ВЪЪI╨Hг┘ЩЪЪHг┘ВЪЪL┴∙ВЪЪЪсH╧H█ҐЩЪЪHЎH█∙юЬЪЪ3юD┴╔юЬЪЪЪяH╩H█ҐЩЪЪHЎH█∙хВЪЪ3юH╧Hг┘хВЪЪ(H┴█пВЪЪЪсИ.ЎЪЪЪЪL█▀<H╩H█ҐЩЪЪH╨┘┐ H╧3юL█┘ЬЖЪЪI╨Hг┘ЩЪЪHг┘ЬЖЪЪL┴∙ВЪЪЪсH╧H█ҐЩЪЪHЎH█∙╦ЬЪЪ3юD┴╔╦ЬЪЪЪяH╩H█ҐЩЪЪHЎH█∙╦ВЪЪ3юH╧Hг┘╦ВЪЪ-H┴█юВЪЪЪсИWЎЪЪЪЪL█▀H╩H█ҐЩЪЪH╨┘┐ H╧3юL█┘ХЖЪЪI╨Hг┘ЩЪЪHг┘ХЖЪЪL┴∙ПЖЪЪЪсH╧H█ҐЩЪЪHЎH█∙╟ЬЪЪ3юD┴╔╟ЬЪЪЪяH╩H█ҐЩЪЪHЎH█∙╗ВЪЪ3юH╧Hг┘╗ВЪЪ%H┴█╟ВЪЪЪсИ─ЎЪЪЪЪL█▀эH╩H█ҐЩЪЪH╨┘┐ H╧3юL█┘ьЖЪЪI╨Hг┘ЩЪЪHг┘ьЖЪЪL┴∙ЮЖЪЪЪсH╧H█ҐЩЪЪHЎH█∙╗ЬЪЪ3юD┴╔╗ЬЪЪЪяH╩H█ҐЩЪЪHЎH█∙≤ВЪЪ3юH╧Hг┘≤ВЪЪ H┴█═ВЪЪЪсИ╘ЎЪЪЪЪL█▀╛H╩H█ҐЩЪЪH╨┘┐ H╧3юL█┘хЖЪЪI╨Hг┘ЩЪЪHг┘хЖЪЪL┴∙пЖЪЪЪсH╧H█ҐЩЪЪHЎH█∙═ЬЪЪ3юD┴╔═ЬЪЪЪяH╩H█ҐЩЪЪHЎH█∙┬ВЪЪ3юH╧Hг┘┬ВЪЪH┴█░ВЪЪЪсИрЎЪЪЪЪL█▀|H╩H█ҐЩЪЪH╨┘┐ H╧3юL█┘╦ЖЪЪI╨Hг┘ЩЪЪHг┘╦ЖЪЪL┴∙юЖЪЪЪсH╧H█ҐЩЪЪHЎH█∙≤ЬЪЪ3юD┴╔≤ЬЪЪЪяH╩H█ҐЩЪЪHЎH█∙xВЪЪ3юH╧Hг┘xВЪЪH┴█─ВЪЪЪсH╩H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧E3ю3юHг┘ЩЪЪЪсИE┘М▌H╩H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧E3ю3юHг┘ЩЪЪЪсI╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘╦ЖЪЪH╩M█┼╪Hг┘ЩЪЪHг┘╦ЖЪЪ$H┴²юЖЪЪAЪсH╧H█ҐЩЪЪHЎH█∙хЖЪЪ3юD┴╔хЖЪЪЪяH╩H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧E3ю3юHг┘ЩЪЪЪсH▀]пL▀e╟L▀m╦L▀uюL▀}хH┴Л]цA┐Щ
▄⌡I╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘hЗЪЪH╩M█┼PHг┘ЩЪЪHг┘hЗЪЪH┴²pЗЪЪAЪсH╧H█ҐЩЪЪHЎH█∙≤ШЪЪ3юD┴Ґ≤ШЪЪЪяH▀EhA╪гИИРЪЪE┘М▌⌡I╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘═ШЪЪH╩M█┼8Hг┘ЩЪЪHг┘═ШЪЪH┴²╗ШЪЪAЪсH╧H█ҐЩЪЪHЎH█∙ЩЪЪ3юD┴ҐЩЪЪЪяH▀EhA╪гИ0РЪЪE┘М▌⌡I╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘╟ЭЪЪH╩M█┼ Hг┘ЩЪЪHг┘╟ЭЪЪH┴²╦ЭЪЪAЪсH╧H█ҐЩЪЪHЎH█∙ХЩЪЪ3юD┴╣ХЩЪЪЪяH▀EhA╪гИwЯЪЪH▀E(D▀(E┘М~nI╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘ЬЬЪЪH╩M█┼хHг┘ЩЪЪHг┘ЬЬЪЪ&H┴²ЫЪЪAЪсH▀EhA╪гИХПЪЪ3ЖИ"кЪЪH▀uP3юI▀дH┴фHЪюH;бrПИиуЪЪ3юИLАЪЪ3ЖИ°БЪЪH▀EhD┴╣pЖЪЪгИЙХЪЪH╨H▀EhD┴╣pЖЪЪD▀╣ЖЪЪг┐:┘©ХЪЪI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ПУЪЪ3юI╧Hг┘ЩЪЪHг┘ПУЪЪL┴█ЬУЪЪAЪрH╦McЖM┘ЖР╦~+H▀MPH╨(пРаHЪюI;фT
(аР_бrэM┘ЖH╧I╦HЎH█ҐЩЪЪH█∙═ЖЪЪР     (а3юР┘═ЖЪЪAЪпI╨H█ҐЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ЖЪЪ3юI╧Hг┘ЩЪЪHг┘ЖЪЪL┴█ЖЪЪAЪрH╧H█ҐЩЪЪHЎH█∙╗ЖЪЪ3юг┘╗ЖЪЪЪяИHГЪЪЎЪЪЪЪL█▀LH╩H█ҐЩЪЪH╨┘┐ H╧3юL█┘pЫЪЪI╨Hг┘ЩЪЪHг┘pЫЪЪL┴∙xЫЪЪЪсH╧H█ҐЩЪЪHЎH█∙ШЪЪ3юD┴╔ШЪЪЪяH╩H█ҐЩЪЪHЎH█∙─ЫЪЪ3юH╧Hг┘─ЫЪЪH┴█┬ЫЪЪЪсИЫЪЪH┐д░L█²PЧЪЪL▀UL█}╗H▀MHH█∙ЮЭЪЪH▀}XH█E░L┴,$L█E≤L┴T$L█M═L┴\$I╩L┴|$H┴T$ H╨L┴d$(H┴L$0H█MЛH▀╣hЖЪЪH┴t$8H┴|$@H┴ъH┴D$H3юL┴D$PL█EПL▀UL┴L$XH▀uhL▀M@L┴T$`AЪсH┐дp▀EьИ╔рЪЪH╧H┴ъH╨3юH▀╣hЖЪЪЪяР┘(ЧЪЪL█∙ЪЪЪРYюH╨ARI╩H█Ґ╟ЫЪЪH╧H█╣0ЧЪЪI╦I╧3юRР┘╟ЫЪЪAЪсH┐дH╦┐8┘MH▀E(┐8▌╠I╩IҐЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘8ЬЪЪI╨M█▀XHг┘ЩЪЪHг┘8ЬЪЪ(L┴∙@ЬЪЪAЪуI╦H█ҐЩЪЪHЎH█∙HЬЪЪ3юH╧Hг┘HЬЪЪ%H┴█PЬЪЪAЪпI╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘ЫЪЪHг┘ЩЪЪM█┼dг┘ЫЪЪAЪсH▀M═I╦HЎH█ҐЩЪЪ3юH█∙ЫЪЪH┴█ЫЪЪAЪп▀MЮI╦HЎH█ҐЩЪЪ3юH█∙ЫЪЪ┴█ЫЪЪAЪпH▀█ЖЪЪI╧HЎH█ҐЩЪЪ3юH█∙ ЫЪЪL▀L┴┘ ЫЪЪAЪяH▀M╗I╦HЎH█ҐЩЪЪ3юH█∙(ЫЪЪH┴█(ЫЪЪAЪпРQ┘0ЧЪЪH╧H█ҐЩЪЪHЎH█∙0ЫЪЪ3юР┘0ЫЪЪЪяH╧H█ҐЩЪЪHЎH█∙8ЫЪЪ3юHг┘8ЫЪЪЪяH▀E0▀┘ю┘4нЪЪH▀Ehг
И╨БЪЪH▀E(D▀(E┘М▌═I╨I╩ЎЪЪЪЪH█ҐЩЪЪH╨┘┐ H╧3юL█┘ЖЪЪH╩M█┼Hг┘ЩЪЪHг┘ЖЪЪH┴²ЖЪЪAЪс▀█░УЪЪH╩HЎH█ҐЩЪЪ3юH█∙╟ЖЪЪ┴█╟ЖЪЪЪсH▀EhгИ ИЪЪ░H█t&H█©#B▓║°г;8888    8080808 8 8 8 8 88080808 8 8 888888808008008 8 8 8 88 8008 08008008ЪЪЪЪ8808088     0 000088 88 888888 80080088 0 000088   888   8   808 8008     8 808080808 88 88 88 88 88 88 88 88 88 88 88 88 88 8П?П?ЪЪЪЪЪЪЪ _═B#B▓║°г;ПЪЪЪЪЪЪЪОЪП?==================================================         O P T I M I Z A T I O N         ==================================================  number of control variables.......  precision on x....................  precision on g....................  expected optimal function value...  maximal number of iterations......  maximal number of simulations.....  information level.................  number of updates.................  size of used memory...............bnc: cold start before instore is bnc: cold start after instore is  ERROR : m is set too small in instore ERROR : uncompatible floating point formatpathei-lsopt: warm start; read via dostorepathei-lsopt: warm; first dostore read: ff =    bnc-lsopt:                           xx =    bnc-lsopt:                           gg =   ERROR : inconsistent nupdate = pathei: nach simul:   nn, ff = pathei: nach simul:   xx(1), xx(2) = pathei: nach simul:   gg(1), gg(2) = pathei-lsopt: cold; first call simul: ff =   Itn      Step  Nfun      Objective    Norm G    Norm X   Norm (X(k-1)-X(k))pathei-lsopt: ++++++++++++++++++++++++  Itn      Step  Nfun      Objective    Norm G    Norm X   Norm (X(k-1)-X(k))pathei-lsopt: pointers updated for first pointer update after cold start; pathei-lsopt: dostore at iiter, jmin, jmax   number of iterations..............  number of simulations.............  relative precision on g...........  after optimization ifail..........  optimization stopped because :    convergence achieved within precision    cold start, no optimization done    the function could not be improved    maximal number of simulations reached    maximal number of iterations reached    the search direction is not a descent one    the initial gradient is too small    warm start file is corruptedpathei: vor instore   ERROR : epsg  =   ERROR : epsg  =   ERROR : epsx =   ERROR : nfunc  =   ERROR : itmax =   ERROR : n     =  ERROR : reading restart file  ERROR : inconsistent jmax =   ERROR : inconsistent jmin =   ERROR : inconsistent jmin = pathei-lsopt: cold startpathei-lsopt vor simulpathei-lsopt xx(1), gg(1) pathei: vor dostore pathei: vor gnorm0 pathei: gnorm0 = pathei-lsopt: entering iiter =pathei-lsopt: back from optline; ifail =    bnc-lsopt:                    sflag =    bnc-lsopt:                    tflag = pathei-lsopt: dostore 1,2 at iiter bnc-lsopt: dostore xxbnc-lsopt: dostore gg iiter, jmin, jmax = pathei-lsopt: iiter, jmin, jmax = bnc: max value of abs(gold) = pathei-lsopt: ys < 0; ifail =    bnc-lsopt: write control   cost function...............  norm of x...................  norm of g...................ifail = ifail = ifail =     optline parameters wrongifail = ifail =     the linesearch failedifail = ifail = ifail = ifail = ifail = ifail =     an input argument is wrongifail =     optimal solution foundifail =     user request during simulx░└0O├▐      ▄┐▌
█°Dцлмноф┐├▄█▌
▐       
╒0Olsopt_top.f
lsopt_top_йи|ц└иyц┤иsм█▓   nЪаtб
дЧбvЪцsц┴т┴фrб┤цце┤┌ц┌┌ф┌цфх▐fjj ÷÷╒╒╒  ≤хооя┬╒тм┬╒тш┬╒т\╪KрицгЬх▀mти▀mткяХи%хm▄фиєыдма^┐гЪёшЪ┐гЪх\Л╗╘н\ЕY└идЧо▌т▌пффЩе²т▐÷нпя▐÷нт▐÷нх\╪Л╗ииЩмЩцх■маи
Рр\┴я√и┌йЩбдцПегVр     тпязрлмаЛ╗Е\гъгШи^Ф^р▄╟!┤МЪ┤РЪГнQ⌡кЪ┬АЪYуЪ╧мЪ≥гЪG▐щЪ`нЪ@яЪ▐гЪФEгЪ┼ясд┬еЪибЪ┬wьЧЙ▓2м▓ЧВ     ёьЪ┌ірЛЖІГлцфЪгЪццЪ┴жцеЪ┘кццЫо╞р
гтх▄Dъ╪═⌠≤И^Е]Ц╩╔c┤МЪb┤ОЪу  о┴╩m┬з"жЩжЩжЩжЩжЩжЩжЩжЩжЩжЩжЩя$@┬╩ Шь
≈Щ     тхУ тхУ ткQmтф√┤лЪ┐фуMиҐвЪ┬Д▒╒жАЧ╗гИ╘^р▄╟'┤Ъ^м▌о╪÷шlsopt_top.f.strtab.symtab.data.bss.textmpi_parm_ph_opti_lsopt_top_.rodata.L_2il0floatpacket.23for_write_seq_fmt.rodata.str1.32__STRLITPACK_150__STRLITPACK_148__STRLITPACK_146__STRLITPACK_144for_write_seq_fmt_xmit__STRLITPACK_142__STRLITPACK_140__STRLITPACK_138__STRLITPACK_136__STRLITPACK_134__STRLITPACK_132__STRLITPACK_130__STRLITPACK_128.rodata.str1.4__STRLITPACK_120__STRLITPACK_118__STRLITPACK_116for_write_seq_lis__STRLITPACK_115__STRLITPACK_114for_write_seq_lis_xmitinstore___STRLITPACK_113__STRLITPACK_111__STRLITPACK_109__STRLITPACK_105__STRLITPACK_104__STRLITPACK_103__STRLITPACK_102__STRLITPACK_101__STRLITPACK_100.L_2il0floatpacket.19__STRLITPACK_99dostore___STRLITPACK_98dnrm2_mpi_allreduce___STRLITPACK_97.L_2il0floatpacket.20__STRLITPACK_96__STRLITPACK_92__STRLITPACK_88__STRLITPACK_86__STRLITPACK_85__STRLITPACK_84__STRLITPACK_83__STRLITPACK_82__STRLITPACK_78__STRLITPACK_77_intel_fast_memcpy__STRLITPACK_76__STRLITPACK_75__STRLITPACK_74__STRLITPACK_73__STRLITPACK_72__STRLITPACK_71__STRLITPACK_70__STRLITPACK_69__STRLITPACK_68__STRLITPACK_67mpi_reduce___STRLITPACK_65__STRLITPACK_64ddot___STRLITPACK_61__STRLITPACK_60__STRLITPACK_59__STRLITPACK_58__STRLITPACK_57dscal___STRLITPACK_56dgscale_lsupdxx___STRLITPACK_55optim_write_control_outstore___STRLITPACK_53__STRLITPACK_51__STRLITPACK_49__STRLITPACK_47__STRLITPACK_45__STRLITPACK_43__STRLITPACK_41__STRLITPACK_3__STRLITPACK_2__STRLITPACK_6__STRLITPACK_5__STRLITPACK_9__STRLITPACK_8__STRLITPACK_12__STRLITPACK_11__STRLITPACK_15__STRLITPACK_14__STRLITPACK_18__STRLITPACK_17__STRLITPACK_21__STRLITPACK_20__STRLITPACK_24__STRLITPACK_23__STRLITPACK_27__STRLITPACK_26__STRLITPACK_30__STRLITPACK_29__STRLITPACK_33__STRLITPACK_32__STRLITPACK_36__STRLITPACK_35__STRLITPACK_0__STRLITPACK_122__STRLITPACK_124__STRLITPACK_126__STRLITPACK_107__STRLITPACK_63.L_2il0floatpacket.21.L_2il0floatpacket.24.L_2il0floatpacket.22__STRLITPACK_62__STRLITPACK_39__STRLITPACK_38__STRLITPACK_80__STRLITPACK_79__STRLITPACK_90lsopt_top_$format_pack.0.1__NLITPACK_6.0.1__STRLITPACK_156.0.1__STRLITPACK_157.0.1__STRLITPACK_158.0.1__STRLITPACK_159.0.1__STRLITPACK_160.0.1__STRLITPACK_161.0.1__STRLITPACK_162.0.1__STRLITPACK_163.0.1__STRLITPACK_164.0.1__STRLITPACK_165.0.1__STRLITPACK_166.0.1__STRLITPACK_167.0.1__STRLITPACK_168.0.1__STRLITPACK_169.0.1__STRLITPACK_170.0.1__STRLITPACK_171.0.1__STRLITPACK_172.0.1__STRLITPACK_173.0.1__STRLITPACK_174.0.1__STRLITPACK_175.0.1__STRLITPACK_176.0.1__STRLITPACK_189.0.1__STRLITPACK_187.0.1__STRLITPACK_188.0.1__STRLITPACK_185.0.1__STRLITPACK_186.0.1__STRLITPACK_183.0.1__STRLITPACK_184.0.1__STRLITPACK_181.0.1__STRLITPACK_182.0.1__STRLITPACK_179.0.1__STRLITPACK_180.0.1__STRLITPACK_177.0.1__STRLITPACK_178.0.1__STRLITPACK_190.0.1__STRLITPACK_191.0.1__STRLITPACK_192.0.1__STRLITPACK_193.0.1__STRLITPACK_194.0.1__STRLITPACK_195.0.1__STRLITPACK_198.0.1__STRLITPACK_197.0.1__STRLITPACK_196.0.1__STRLITPACK_230.0.1__STRLITPACK_231.0.1__NLITPACK_1.0.1__NLITPACK_5.0.1__NLITPACK_2.0.1__STRLITPACK_232.0.1__STRLITPACK_233.0.1__STRLITPACK_234.0.1__STRLITPACK_235.0.1__STRLITPACK_236.0.1__STRLITPACK_237.0.1__STRLITPACK_238.0.1__STRLITPACK_239.0.1__STRLITPACK_228.0.1__STRLITPACK_229.0.1__STRLITPACK_226.0.1__STRLITPACK_227.0.1__STRLITPACK_224.0.1__STRLITPACK_225.0.1__STRLITPACK_222.0.1__STRLITPACK_223.0.1__STRLITPACK_199.0.1__STRLITPACK_200.0.1__STRLITPACK_201.0.1__STRLITPACK_202.0.1__STRLITPACK_203.0.1__STRLITPACK_204.0.1__STRLITPACK_205.0.1__STRLITPACK_206.0.1__STRLITPACK_207.0.1__STRLITPACK_208.0.1__STRLITPACK_209.0.1__STRLITPACK_210.0.1__STRLITPACK_211.0.1__STRLITPACK_212.0.1__STRLITPACK_213.0.1__STRLITPACK_214.0.1__NLITPACK_0.0.1__NLITPACK_3.0.1__STRLITPACK_215.0.1__NLITPACK_4.0.1__STRLITPACK_216.0.1__STRLITPACK_217.0.1__STRLITPACK_218.0.1__STRLITPACK_219.0.1__STRLITPACK_240.0.1__STRLITPACK_241.0.1__STRLITPACK_242.0.1__STRLITPACK_243.0.1__STRLITPACK_244.0.1__STRLITPACK_245.0.1__STRLITPACK_246.0.1__STRLITPACK_247.0.1__STRLITPACK_248.0.1__STRLITPACK_249.0.1__STRLITPACK_250.0.1__STRLITPACK_251.0.1__STRLITPACK_252.0.1__STRLITPACK_253.0.1__STRLITPACK_254.0.1__STRLITPACK_255.0.1__STRLITPACK_256.0.1__STRLITPACK_257.0.1__STRLITPACK_258.0.1__STRLITPACK_259.0.1__STRLITPACK_260.0.1__STRLITPACK_261.0.1__STRLITPACK_262.0.1__STRLITPACK_263.0.1__STRLITPACK_264.0.1__STRLITPACK_265.0.1__STRLITPACK_266.0.1__STRLITPACK_267.0.1__STRLITPACK_268.0.1__STRLITPACK_269.0.1__STRLITPACK_270.0.1__STRLITPACK_271.0.1__STRLITPACK_272.0.1__STRLITPACK_273.0.1__STRLITPACK_274.0.1__STRLITPACK_275.0.1__STRLITPACK_276.0.1__STRLITPACK_277.0.1__STRLITPACK_278.0.1__STRLITPACK_289.0.1__STRLITPACK_290.0.1__STRLITPACK_291.0.1__STRLITPACK_292.0.1__STRLITPACK_293.0.1__STRLITPACK_283.0.1__STRLITPACK_284.0.1__STRLITPACK_285.0.1__STRLITPACK_286.0.1__STRLITPACK_287.0.1__STRLITPACK_288.0.1__STRLITPACK_294.0.1__STRLITPACK_295.0.1__STRLITPACK_296.0.1__STRLITPACK_297.0.1__STRLITPACK_279.0.1__STRLITPACK_280.0.1__STRLITPACK_281.0.1__STRLITPACK_282.0.1__STRLITPACK_298.0.1__STRLITPACK_299.0.1__STRLITPACK_300.0.1__STRLITPACK_301.0.1__STRLITPACK_302.0.1__STRLITPACK_303.0.1__STRLITPACK_304.0.1__STRLITPACK_305.0.1__STRLITPACK_306.0.1__STRLITPACK_307.0.1__STRLITPACK_308.0.1__STRLITPACK_309.0.1__STRLITPACK_310.0.1__STRLITPACK_311.0.1__STRLITPACK_312.0.1__STRLITPACK_313.0.1__STRLITPACK_314.0.1__STRLITPACK_357.0.1__STRLITPACK_358.0.1__STRLITPACK_359.0.1__STRLITPACK_360.0.1__STRLITPACK_315.0.1__STRLITPACK_316.0.1__STRLITPACK_353.0.1__STRLITPACK_354.0.1__STRLITPACK_355.0.1__STRLITPACK_350.0.1__STRLITPACK_351.0.1__STRLITPACK_352.0.1__STRLITPACK_347.0.1__STRLITPACK_348.0.1__STRLITPACK_349.0.1__STRLITPACK_344.0.1__STRLITPACK_345.0.1__STRLITPACK_346.0.1__STRLITPACK_341.0.1__STRLITPACK_342.0.1__STRLITPACK_343.0.1__STRLITPACK_338.0.1__STRLITPACK_339.0.1__STRLITPACK_340.0.1__STRLITPACK_335.0.1__STRLITPACK_336.0.1__STRLITPACK_337.0.1__STRLITPACK_332.0.1__STRLITPACK_333.0.1__STRLITPACK_334.0.1__STRLITPACK_329.0.1__STRLITPACK_330.0.1__STRLITPACK_331.0.1__STRLITPACK_326.0.1__STRLITPACK_327.0.1__STRLITPACK_328.0.1__STRLITPACK_323.0.1__STRLITPACK_324.0.1__STRLITPACK_325.0.1__STRLITPACK_320.0.1__STRLITPACK_321.0.1__STRLITPACK_322.0.1__STRLITPACK_317.0.1__STRLITPACK_318.0.1__STRLITPACK_319.0.1__STRLITPACK_356.0.1mpi_sgi_private_status_mpi_sgi_private_unweighted_mpi_sgi_private_mpi_sgi_private_inplace_mpi_sgi_private_char_mpi_null_copy_fn_mpi_null_delete_fn_mpi_dup_fn_mpi_comm_null_copy_fn_mpi_comm_null_delete_fn_mpi_comm_dup_fn_mpi_win_null_copy_fn_mpi_win_null_delete_fn_mpi_win_dup_fn_mpi_conversion_fn_null_.note.GNU-stack.eh_frame.trace.rela.text.rela.rodata.rela.eh_frame.rela.trace
JAIЬL"h5     _i
zL pйтЕLx5?PLp─═M╠┬пз
КL░?MI≤pz▀L╚═ъMИ╗+LK╟┌M▄╦ЁҐнLНю%M/хV`qL▒пхMрьЫL4ЮgMxХ■·╞LоПMЬ/9LYl²MўL68I~M┬@шЕL    ( M MW 0≈ ╦ б LБ У *
M4
 p
N░
ё
м
NМ
x*O;─ZOk┬⌠PжIЛN░LOV≤|O├═╘JтчLЧ╦
m
w
L≈
╟╙
Э
J!+LKX^▐N╞`бЛOЩhN<pO ─O┼x╣O©─I,NL┬_!┴O ░юOй≤ЯN═$"UO_╗█O≈╟©Nъ╦Р# O*юSO]хЛюANaпt$·л╗Q╣ькQшьЕтТQьэI-NMЮ`%┼R л╦лбSсХАэI.NNЛa&▐O≥ТҐХХJN&Э9'cJmOwисLСH(=MGP╚╣Lу(Х)M)0┴L╘╪*СMЩ 7NWюj+■NЄхялшQХпЧQпт'I1JVNvь┴,╨OдЮКNХ-OOYП┤O▒Ь╧NыЛ.O$MOW╦IкNКTЧ/(NH\[0┘O√dхTшJJьdIwN≈l╙1тOЕtЬ2O%|LNl└3╘O╨▄м4ПOЗ■ N= °P 5z O▀ є· 6ю Oй ╛
!I!N=!ЄP!7z!O▀!╪·!8Ґ!Oн!дЙ!N
"л"9N"OX"т├"O░"э╦"Nь"ДК":#O##ЛL#OV#Тz#Q┤#ь▒#лї#QЇ#ьа#тп#RЮ#лЧ#л$S$Х*$э[$Rh$л≤$лє$U╣$Хф$эп$пЛ$RЫ$л%л*%U;%ХL%эV%пr%I▓%°%L╪%Эо%;&M&$&<C&M&Lm&═&M╙&я&Mш&'M'$7'MA',m'M~'4╞'M╧'<щ*лГ*V+л+S#+Х4+э≈+Iґ+Nм+lЮ+=
,O,t.,>M,O^,|q,?░,O║,└ю,Oй,▄С,OЩ,■V-Ii-N┴-D°-@ф-Oв-LЙ-A .O.T<.OF.\o.Oy.dю.й.Wч.л
/W/л1/ID/Nd/°w/B║/O╡/єя/Oш/╛0O0Є50QB0ьl0Qy0ь÷0X╛0`┐1Y║1ьТ1I2N'2э:2Ch2Or2Д²2Oї2Ло2ZЗ2[L3I_3N3Т°3Dі3LЇ3в3Ь4M4=4Ld4w4Eє4Mў4я4LЬ45F95MC5 j5Rw5л÷5л╚5U╪5Хм5эв5пС5R6л$6л/6U@6ХQ6э[6п┬6I⌡6N╩6(в6GА6LР67,L7MV74}7Lє7<Ї7HИ7MЗ7D#8L>8LQ8I┌8M⌠8T╡8IБ8N9t9J(9L99Y9x╒9╪9Lв9─Й9K:M$:┬@:MQ:░d:L⌠:Lў:≤а:MЙ:MШ:═;M(;╗;;Nj;L┘;╟≤;Oа;Mр;╦Н;MЪ;ю<PA<L\<хo<Q≤<M╘<пе<Mж<ьИ<R=L3=ЮF=So=M─=Х°=Mґ=Пю=TО=L
>Ь>UF>MW>s>M└>≈>Vф>LА>Т>W?M.?J?M[? n?X²?L╦?(к?YТ?M@0!@M2@8E@Zt@L▐@@╒@[к@Mэ@HЬ@M APA\KALfAXyA]╒AMЁA`оAMЮAhСA^"BL=BpPB_yBM┼BxіBMЇB─йB`ЫBLC┬'CaPCMaC░}CM▌C≤║Cb©CNъC╦ DN)D\EDODLoD`┌DcЁDMдDhЮDNEp?EIELiEH|EdґEMЎEPЬEFL"FX5FefFMwF`╠F╩FLшFhНFfGM0GpmGwGL≈G╗╙Gg6HIaHN│H╪■HhЎHПщHпIЬIOIдFINfIлyIiёIOЄIтДILЪI═Jj;JMLJ╗hJMyJ╟▄JkКJYЪJь`KRmKл∙Kл║KU╡KХцKэмKпИKI LLL3LFLlwLM┬L⌡Lm╨LдLLДLMM$M$KMMUM,MM┴M4ЄMMЎM<ЙMMШMDNM)NLvN─NL═N8ЁNnЙNMТN@КBB=Af@ ▐?(╦>0А=8
=@3<H\;P┘:Xў9 K,K00000@"IH"Ю #()(()0OLXx  |2xlo2Д├п╚Є▀╩Є▀═еT▄╗
Э▐ЪлШ╙╦>вЁИ ДсЙ     СКЙ0
lsupdxx.o/      1361899271  10100 26113 100644  10872     `
ELF>╦&@@ЯЪ!&?╠а,(ГЖ*2HHW(~H┼`H∙юH═ЄгшЛЩ $",64J<^D╞©  и
пш
Й,РЪ6÷п#AfrРЪ┼РЪіРЪЇРЪпРЪФЬ/HYo┤≈UH┴ЕH│Л─L┴ґ ЪЪЪIҐL┴ХM┴КH┴∙╗ЩЪЪH╨L┴╔ЪЪЪH┴ҐХЧЪЪL┴ОL▀`(H┴пL┴Ґ0ЪЪЪM┴ОM▀S0I┴сH┴█░ЩЪЪH┴╣хЧЪЪL┴█≤ЩЪЪI┴яL┴┘═ЩЪЪI┴пL┴∙@ЪЪЪI┴рL┴╔PЪЪЪI┴тH▀wH▀J@H╨H▀xH╦H┴²8ЪЪЪI▀_8H┴²HЪЪЪI▀_H┴²xЪЪЪH┴█XЪЪЪH╧I▀Z I╨H▀R@H▀@L┴╣(ЪЪЪM▀u@H┴UпH┴∙0ЧЪЪH╨H┴E═H┴┘ЬЩЪЪH╦L┴u┬M▀u H┴╣hЪЪЪI▀t$I╪M▀h8I╦L┴╣`ЪЪЪM▀q0I╧H┴uюH▀I8H┴╣@ЧЪЪH┴жL┴ґpЪЪЪM▀RL┴ґpЧЪЪIҐH┴M╦H┴█(ЧЪЪM▀{(I╩L┴U╗L┴∙ЧЪЪH▀
L┴ЙL▀PL┴∙ЧЪЪL▀∙═ЩЪЪH┴]░H┴}≤H┴ҐHЧЪЪH┴²XЧЪЪM▀@0M▀I(H▀~H▀²XЪЪЪL┴uхL┴}╟L┴E─L┴MьL┴█ЧЪЪL┴┘ ЧЪЪH┴█8ЧЪЪH┴ҐPЧЪЪL┴Ґ`ЧЪЪL┴╣hЧЪЪH┴²xЧЪЪM▀[ M▀$$M▀uL▀ҐhЪЪЪH▀┘xЪЪЪH▀²`ЪЪЪH▀╣PЪЪЪH▀Ґ@ЪЪЪL▀┘HЪЪЪL▀M┬H▀JA▀┘рL┴]ЮL┴╔ПЩЪЪL┴²ЧЪЪL┴╣─ЧЪЪL┴Ґ┬ЧЪЪH┴┘░ЧЪЪH┴█≤ЧЪЪH┴²═ЧЪЪH┴╣╗ЧЪЪH┴Ґ╟ЧЪЪL┴┘╦ЧЪЪL┴█юЧЪЪu3H╩H©H▀uH▀EH▀M РH▀РYH┴TК&H▀┘≤ЩЪЪ┐8┘ыH╦H▀UРРYH▀E(РXюРРYH▀∙ХЧЪЪР^а▀┘ю▌вHcхH┐Ы▄цH▀]@H┐Ц┘шtЖц┘ўj[┴зH█rH;н▄°Вшь┐Ц+ц3шHcюH┘рv,H▀}@L▀E8РAьHЎРYхWРъHЪцH;зrэH▀]8H█4сHВф└╒L▀M@(хL▀U8РAр(ИРA\рH╩РAdр H┴чРAtр0H┴ъfATрI┴ьfA\рfAdр(fAtр8fИfYуfYщfYЕfYУWWW'AW0A)яA)\яA)dя A)tя0H┐бH;п┌nЪЪЪИ}L▀M@(хL▀U8H█І(ИH╩fИH┴ч(у(щ(ЕH┴ъfAYрfAY\рfAYdр fAYlр0I┴ьWWW'AW(A)яA)\яA)dя A)lя0H┐бH;пr∙H;аs+H▀]@H▀u8РфH╨РYхW
РцHЪюH;аrщH▀┘╗ЩЪЪЖ└_H╦┐8┘LI╧H█Ґ░ЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘пЩЪЪ3юH╩Hг┘░ЩЪЪHг┘пЩЪЪ'H┴²ьЩЪЪAЪяH╩H█Ґ░ЩЪЪHЎH█∙ЮЩЪЪ3юH╧Hг┘ЮЩЪЪH┴█ХЩЪЪЪсИ╗H▀┘ХЧЪЪ▀┘ю▌сH≤H┴┘─ЩЪЪH┐Ь▄
H▀M@H┐А┘иtЖа┘
jY┴иH┴█┬ЩЪЪH█YH;²─ЩЪЪ▐Ц    Вых┐А+а3иH▀²┬ЩЪЪH┘шHcюv6H▀╣┬ЩЪЪH▀}@L▀E8H╩РAхWРоHЪаH;нrЮH┴╣┬ЩЪЪH▀M8H▀²┬ЩЪЪH█4ыHВф└■L▀┘┬ЩЪЪL▀M@L▀U8░H█t&H█©H╧H┴кH┴нH┴оРCбРCLбРCTб РC\б0fCDбfCLбfCTб(fC\б8WWWWC)аC)LаC)Tа C)\а0I┐юL;юr┼КWH╧H▀]@H▀u8(H▀█┬ЩЪЪ(х(п(ь(ЮWнWTнW\н Wdн0)к)Tк)\к )dк0H┐аH;хrеH;┘─ЩЪЪs.H▀²─ЩЪЪH▀u@H▀}8H╧РгWРфHЪюH;цrА┘р~HH┐дПI╨H▀²╗ЩЪЪ3юH▀U@L▀MPH┴$H▀ҐХЧЪЪH▀╣≤ЩЪЪH▀█░ЩЪЪL▀┘═ЩЪЪAЪрH┐дКxH▀┘╗ЩЪЪЖtlH╦┐8u]I╧H█Ґ░ЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘─ЩЪЪ3юH╩Hг┘░ЩЪЪHг┘─ЩЪЪ)H┴²┬ЩЪЪAЪяH╨H╩I┴п3юH▀ҐХЧЪЪH▀u@H▀M8ЪсH┐дПH╨I╨H█}ХH╧I╦I┴я3юH▀u0H█]ПРEХH┴$AЪрH┐дH▀E0fОиРf/а┌┌H╨H▀┘хЧЪЪг┐:┘fI╧H█Ґ░ЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ХЧЪЪ3юH╩Hг┘░ЩЪЪHг┘ХЧЪЪH┴²ПЧЪЪAЪяИH▀∙ХЧЪЪH▀EXHcH┘рHг╦~/H▀]XH▀u@H╧РDфЬHЪюРH;бTР_хР~ыH▀uXH┐дПH╨I╩H╧H█}ХH▀I┴пI┴я3юH┴]ХL█UПL┴$AЪсH┐дH▀EpH╧П?H▀UXH▀]hРР^H▀u`H┴H╧ _═BРH┴f/rH▀UhfH~юf/H┴vH▀EXH┴H▀EhРf/vH▀UhH▀█хЧЪЪH▀H┴гH╧H╨3юH▀ҐХЧЪЪH▀u@ЪяРEХH┐дПРYюH╨I╨H█ҐЪЪЪH╧H█╣ЪЪЪI╦I┴я3юH█]ПР┘ЪЪЪH┴$AЪрH┐дРQ┘ЪЪЪH╦Р┘ЪЪЪ┐8┘#I╧H█Ґ░ЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘хЧЪЪ3юH╩Hг┘░ЩЪЪHг┘хЧЪЪH┴²пЧЪЪAЪяH▀MhI╦HЎH█Ґ░ЩЪЪ3юH█∙ЬЧЪЪH▀H┴²ЬЧЪЪAЪпI╧H█Ґ░ЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ьЧЪЪ3юH╩Hг┘░ЩЪЪHг┘ьЧЪЪH┴²ЮЧЪЪAЪяH▀█ЪЪЪH╩HЎH█Ґ░ЩЪЪ3юH█∙ЪЪЪH┴█ЪЪЪЪсH▀┘ХЧЪЪ▀┘ю▌HcьH▀MHH┐Ш▄▐H▀uPH┐Ф┘ЖtВф┘wj^┴РH█zH;ъ▄eВчП┐Ф+ф3ЖHcюH┘рv+H▀}hL▀E@РH▀}P(хРAYПРXЯРВHЪфH;РrЕH▀u@H█<жHВг└дH▀uhH▀}@РH▀uPРвD(юР\вРlв Р|в0fLвf\вflв(f|в8fDюfAYхfAYьfAYХfDYгРяРdяРtя РDLя0fTяfdяftя(fDLя8fXяfXЦfXУfEXх)ж)dж)tж D)Lж0H┐бH;п┌QЪЪЪИ H▀uhH▀}@РH▀uP(ЬfЬ(о(ъ(ОfYвfY\вfYlв fY|в0РяРdяРtя РDDя0fTяfdяftя(fDDя8fXяfXЦfXУfDXг)ж)dж)tж D)Dж0H┐бH;п┌vЪЪЪH;цs*H▀UhH▀u@РH▀UP(хРYфРXаРбHЪюH;цrФH▀²≤ЧЪЪHВц┘ЬH▀²PЧЪЪHВц┘≤H▀┘ЧЪЪH╘┘QH╨H┴сI┴уI┴вH┴жI┴пI┴рH▀E═L▀]─L▀e╦L▀uпH┴BH╨H▀M╗H▀}ЮL▀MьL┴[0H┴сM┴e8I┴уM┴w@I┴вH┴NH┴жI┴x I┴пM┴J(I┴рH▀EюL▀]хL▀╔pЪЪЪL▀╣XЪЪЪH┴BH╨H▀M≤H▀}░L▀M╟L┴[0H┴сM┴e8I┴уM┴w@I┴вH┴NH┴жI┴x I┴пM┴J(I┴рL▀²@ЪЪЪL▀╔HЪЪЪL▀u┬H▀┘hЪЪЪH▀█xЪЪЪH▀Ґ`ЪЪЪL▀█PЪЪЪL┴[0M┴e8M┴w@H▀²8ЪЪЪL▀╔ЪЪЪL▀ґ ЪЪЪL▀╣(ЪЪЪL▀Ґ0ЪЪЪH┴BH┴NI┴x M┴J(H┴Л]ц3юИпТЪЪ3юИmВЪЪ3юИEЧЪЪH┴фH╨H┐Ф┐ЮHяНюH▀ҐПЩЪЪаФП3ю│нЪрИ{ЧЪЪH┴ч┴зH┐Ф┐БHяНрH╧3юаФР│нH▀Ґ8ЧЪЪЪяH┐ЦЧHг┘8ЧЪЪH┴²PЧЪЪИЧЪЪH┴ч┴зH┐Ф┐БHяНрH╧3юаФР│нH▀Ґ─ЧЪЪЪяH┐ЦЧHг┘─ЧЪЪH┴²≤ЧЪЪИ╪ЩЪЪH█vH█©88880808──ЪЪЪЪЪЪЪ─{ўGАz└? _═BП?pathei-lsopt: no hessupd for first optim.pathei-lsopt: first estimate of dd via .01*ffbnc: initial stepsize tact = bnc: norm of dd            = bnc: ifail = 4, dotdg > 0.x░└├█(▄▐g┐D▌Жцлмноф┐├▄█▌▐
≤       lsupdxx.flsupdxx_ф:иЪцгп:сфи:Йфи9йгф9иЪцхф:яфф:бЪбЪйфг╩м9ггй9огм:БЪцгц:УЧд
ц╩м:ЧЧа
и┴цЪгЁцMбШцбШдоВуцЁкP├цЪX⌡эЪйЪебЧCцЧСяЪтЪицЧPуЪ гЪБё
~≈ЪFlйЪMЮЪ√гЪцGз\И     ╪дЪ┼уилЪjф ц──фЪ▌▐┌цЪ┌йЪ│┐╔бЪ┤цЪ▌|Ъ▄     цЫ{Ъ┘{Ъ▄┴цЬб┘┬┐бЪ│├м┐├еЪ┘ГY┤игЧх▒░сZЪіFоОрэЪ╠яЪ▄яЪоОягЪvф÷ффFжlsupdxx.f.strtab.symtab.data.bss.textmpi_parm_lsupdxx_.rodata.L_2il0floatpacket.28.L_2il0floatpacket.31.L_2il0floatpacket.32.L_2il0floatpacket.29for_write_seq_lis.rodata.str1.32__STRLITPACK_5for_write_seq_lis_xmit.rodata.str1.4__STRLITPACK_4hessupd___STRLITPACK_3ddot_mpi_allreduce___STRLITPACK_2dnrm2___STRLITPACK_1__STRLITPACK_0for_dealloc_allocatablevar$101.0.1var$91.0.1var$81.0.1__STRLITPACK_11.0.1__STRLITPACK_9.0.1__STRLITPACK_10.0.1__NLITPACK_0.0.1__NLITPACK_2.0.1__NLITPACK_1.0.1__STRLITPACK_13.0.1__STRLITPACK_14.0.1__STRLITPACK_15.0.1__STRLITPACK_16.0.1__STRLITPACK_12.0.1mpi_sgi_private_status_mpi_sgi_private_unweighted_mpi_sgi_private_mpi_sgi_private_inplace_mpi_sgi_private_char_mpi_null_copy_fn_mpi_null_delete_fn_mpi_dup_fn_mpi_comm_null_copy_fn_mpi_comm_null_delete_fn_mpi_comm_dup_fn_mpi_win_null_copy_fn_mpi_win_null_delete_fn_mpi_win_dup_fn_mpi_conversion_fn_null_.note.GNU-stack.eh_frame.trace.rela.text.rela.eh_frame.rela.trace+`°ю╙юьюФю`#юAюOюdю▓╚юь─Б`─╚pЫP∙Pp=$P&p┐     ґ'ЎяspрPJP╧pЮ(0$?&_r°і)йт*Б Л)    $I &i D| 
р    `

*
 ╩
+е
К
У
* E$`&─$⌠а'к,Р&4%V'`<Йю`t>,─,л, %,%00@═$Ю  !&?═╠2╟TГ2d╞h©h═и≤═Жп√ьшn&    Й├&0
lsline.o/       1361899272  10100 26113 100644  15920     `
ELF>p:@@ЯЪ
 %OWfФЖ(HAPf▀,ТH`HюH*;L]n┌ √(╙0Ў8р@ФH7G   Q
Xc
r+РЪ5╟!=u┬▐·эЗРЪРЪ.РЪ?РЪXРЪn─■═ЇпАВUH┴ЕH│Л╟H╦H┴²8ЪЪЪI┴бH┴цL┴┘xЩЪЪI┴юH┴∙hЩЪЪH╨L┴Ґ0ЪЪЪL┴╣(ЪЪЪI┴фL▀{0H┴сL┴╔ЪЪЪH┴█`ЩЪЪL▀X@H▀H M▀b8L┴█pЩЪЪH┴ҐPЩЪЪL┴]юI┴цL┴╔PЪЪЪI┴тL┴}╟I┴вH┴MХH┴пM▀HH┴яH▀z@I┴пH╨L┴ґ ЪЪЪM▀n(L▀s0L┴M╦I╧L┴uхL┴╣xЧЪЪIЎH┴}─H┴вH▀Z(H┴]пH┴²(ЧЪЪH╩H▀@ H┴╣─ЩЪЪM▀I@Hc6L┴ґHЪЪЪM▀l$8I┴ЭH┴┘pЪЪЪH┴┘hЧЪЪI▀L┴█hЪЪЪH▀W L┴█@ЧЪЪM┴ЯL▀╣PЪЪЪH┴╣XЩЪЪHЎH┴U┬H┴∙ ЧЪЪL┴╣хЧЪЪH▀L▀╣XЩЪЪM┘ЖM▀SI╩H┴∙HЧЪЪ╨L┴U░I╨I▀|$IOжM▀d$H▀IH┴█XЪЪЪL┴e≤L┴╔ЧЪЪI╪H┴█XЧЪЪH┴ыM▀(H┴┘░ЧЪЪH▀EюL┴Ґ`ЪЪЪL┴ҐpЧЪЪM▀@M▀[8M▀R0L▀}─H┴┘пЧЪЪH█уL┴m╗L┴E═L┴]ЮL┴∙xЪЪЪH┴}ьH┴ҐЧЪЪL┴∙0ЧЪЪL┴²8ЧЪЪL┴┘PЧЪЪL┴ґ─ЧЪЪL┴Ґ┬ЧЪЪH▀6M▀d$H▀YH▀}░L▀E╦L▀UХL▀²HЪЪЪL▀m╟M▀yH┴┘┬ЩЪЪH┐ЮПH┴╣ЧЪЪL┴╔ЧЪЪH┴²`ЧЪЪH┴Ґ≤ЧЪЪL┴┘═ЧЪЪL┴Ґ╗ЧЪЪL┴∙╟ЧЪЪL┴²╦ЧЪЪL┴ґюЧЪЪH┴╔@ЪЪЪH+ЮH┐ДПH┴ЮI┴еH▀┘┬ЩЪЪH┐ЮПH+ЮH┐ДПH┴ЮH┴┘XЩЪЪH▀┘┬ЩЪЪH┐ЮПH+ЮH┐ДПH┴ЮH┴┘┬ЩЪЪM┘Ж▌0H▀┘`ЩЪЪЖtlH╦┐8u]I╧H█Ґ░ЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ПЩЪЪ3юH╩Hг┘░ЩЪЪHг┘ПЩЪЪ
H┴²ЬЩЪЪAЪяI┐Ч▌▓H╦J█УL┴ОH▀uHЪпH╧H╨3юH▀Ґ─ЩЪЪH▀uHЪяР┘ЪЪЪH┐дПРYюH╨I╨H█ҐЪЪЪH╧H█╣ЪЪЪI╦I┴я3юH█]ПР┘ЪЪЪH┴$AЪрH┐дРQ┘ЪЪЪH▀E Р┘ЪЪЪЖ└г┘ТЧЪЪЪЪЪЪH©H▀∙─ЩЪЪH█╣ТЧЪЪH▀MXL█┘ЬЧЪЪL▀MHЪ∙PЩЪЪI┐Ч▌QH╦J█УH▀Ґ┬ЩЪЪH▀uXЪп3рM┘ЖH┴ЦIMжH█уH┐ЮПH+ЮH┐ДПH┴ЮM┘Ж▌шI┐Ч▄ХH╧ЬЪЪЪ3рH▀}XI#нL▀EPРвР\вРlв Р|в0РAпРATпРAdп РAtп0fADпfATпfAdп(fAtп8fLвf\вflв(f|в8f\хf\зf\Лf\Ч)п)\п)lп )|п0H┐бH;я┌jЪЪЪI;нsH▀UXH▀}PРйР\оРхHЪаI;нrИH╨I╧H┴фL┴ИI┴п3юH▀Ґ─ЩЪЪAЪяH┴ьР┘ЪЪЪH┴дH┐дПH╨I╨H█ҐЪЪЪH╧H█╣ХЧЪЪI╦I┴я3юH█]ПH┴$AЪрH┐дH▀MH╩РґЬЧЪЪHЎР█ХЧЪЪfОюH▀UH▀┘ХЧЪЪР\)H┴Р\ґХЧЪЪРXМР^мWР.Р_хР█PЩЪЪf/И┌ H▀M0H╦П?H▀UM┘ЖР∙ЪЪЪH┴Р"Р▌вI┐Ч▄ЖH▀UXH┐БtHВб┘ъ╨H█BL;П▄мL┴П3иH+бH┐ЮHВьIфH┘рvMH▀uXH▀}PH█v(мH╩РYл(бРYбРYРAYLмРYкР^хР4оР\ЯР4нHЪаH;йr©I█LуHВа└┘L▀EX(кL▀MP(УH╧РYТ(аРY1D(бРDYбРA|уD(уfA|уH┴кfЖD(УfYВРDYтРDYТРDYfаH┴нfYПРDY6fEюH┴оfA^ПD(еD(БРDYдРDYБРDYРEяfELяfD\нРE\у(РfE\уfEрE)пfEYсРYРfDYпРE|у РELу0fE|у(fELу8fEЖfEюfEYВfEYаfDYПfDYю(бРYбfEДfE^тfЖfюfD^ЖfD^юРElяfElяРA|я fE\ЙРETя0fA|я(fETя8fA\ЧfE\пE)lпA)|п E)Tп0H┐бH;п┌▀ЧЪЪИPL▀EX(кL▀MP(УH╧РYТD(мРDYлРY1D(ЕD(ЩРDYДH┴кРDYЭРDYH┴нH┴оfЖ(аfAYtуРDY&РDY?fEи(ЗfEДD(рfEЪD(ЙfEYLуfEYdу fEY|у0РYЗРDYрРDYЙfаfYПfDYхfDYЮfDYЬ(бРYбfЪf^ВfEрfEМfюfE^йfE^ЕfD^ЬРEяfEDяРE\яРEtя fD\фРAtя0fE\яfEtя(fAtя8fE\ыfE\ТfA\ВE)пE)\пE)tп A)tп0H┐бH;п┌╩ЧЪЪI;фsTH▀MXH▀]PH█vH█©(мH╨РYл(бРYбРY
РAYLеРYкР^хР4цР\ЯР4аHЪюI;фr©H▀∙`ЩЪЪH▀E(Жг└┌H╦┐8┘жI╧H█Ґ░ЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘pЩЪЪ3юH╩Hг┘░ЩЪЪHг┘pЩЪЪ'H┴²xЩЪЪAЪяР┘PЩЪЪH█Ґ░ЩЪЪH╧HЎ3юH█∙пЩЪЪР┘пЩЪЪЪяH╩H█Ґ░ЩЪЪHЎH█∙─ЩЪЪ3юH╧Hг┘─ЩЪЪH┴█┬ЩЪЪЪсI╦A┐8└СИ│M┘Ж▌╗I┐Ч▄H▀EXH┐ЮtH╘┘Л╦H█PL;Р▄зL┴Р3иH+пH┐БHВзIжH┘ю├eH╩Р∙PЩЪЪH▀uPРH▀]X(кРкР\йРYбРYнРXхРкHЪаH;хrзH▀MPH█аHВц└ВР█PЩЪЪР∙PЩЪЪH▀MXH▀]P(Ц(ШD(сD(К(аР\БР\ЗРD\рРD\ЙР,цРDDцРD\ц РDtц0fаflц(ПfDDцD(хfD\ц(D(ЮfDtц8fДfЪfEрfEМfYЕfY4аfDYLаfDYdа fYDа0fAYЬfEYсfEYНfXФfAXЫfEXтfDXХ)$а)|аD)Tа D)lа0H┐юH;б┌&ЪЪЪИбР┘PЩЪЪРD∙PЩЪЪH▀MXH▀]P(к(Ц(СD(цD(хРA\йРA\БРA\РРE\бfDхfиA(яfДA(ИfЖA(ЫfEюfYцfYdцfYtц fDYDц0fYаfYlаfY|а fDYLа0fXйfXЕfXВfEXа)а)dа)tа D)Dа0H┐юH;б┌WЪЪЪI;жsDH╦Р²PЩЪЪH▀MPРH▀EX(пРпР\сРYкРYяРXяРпHЪбI;жrзH▀█`ЩЪЪH▀E0H▀U(Р┘PЩЪЪРгЪЪЪЪЖ└іH╦┐8┘⌠I╧H█Ґ░ЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘ьЩЪЪ3юH╩Hг┘░ЩЪЪHг┘ьЩЪЪH┴²ЮЩЪЪAЪяР┘PЩЪЪH█Ґ░ЩЪЪH╧HЎ3юH█∙ХЩЪЪР┘ХЩЪЪЪяL▀╔ЧЪЪH▀²`ЧЪЪL▀Ґ╗ЧЪЪИг┘ТЧЪЪH©H▀∙─ЩЪЪH█╣ТЧЪЪH▀MXL█┘ЬЧЪЪL▀MHЪ∙PЩЪЪH▀E(Ж└H▀┘xЩЪЪH█ҐПЧЪЪH▀∙─ЩЪЪH█╣ТЧЪЪH▀█┬ЩЪЪL█┘ьЧЪЪ▀Ъц┴²ПЧЪЪL▀█XЩЪЪЪ∙PЩЪЪI┐Ч▄H▀E@H┐ЮtH╘┘Ч╦H█PL;Р▄ЛL┴Р3иH+пH┐БHВзIжH┘юv&H▀]@H▀uPH▀Ґ┬ЩЪЪРоР\нРкHЪаH;хrИH▀MPH█аHВц└╒H▀M@H▀]PH▀╣┬ЩЪЪРфР\фРlф Р|ф0РцРTцРdц Рtц0fDцfTцfdц(ftц8fLфf\фflф(f|ф8f\хf\зf\Лf\Ч)а)\а)lа )|а0H┐юH;б┌rЪЪЪИ─H▀M@H▀]PH▀╣┬ЩЪЪH█Є&H█©РфРLфРTф Р\ф0fDфfLфfTф(f\ф8f\цf\Lцf\Tц f\\ц0)а)Lа)Tа )\а0H┐юH;бr·I;жs&H▀E@H▀MPH▀²┬ЩЪЪРсР\яРпHЪбI;жrИH▀]H╨Р┘ЬЧЪЪI╧L┴ИI┴п3юР\H▀Ґ─ЩЪЪH▀u@Р┘PЩЪЪAЪяH┐дПH╨I╨H█ҐЪЪЪH╧I╦I┴я3юH▀uH█]ПР┘ЪЪЪH┴$AЪрH┐дH▀E0H╨H▀MР┘PЩЪЪРРY
РY     f/а├┘H▀∙`ЩЪЪH▀┘hЩЪЪЖгtlH╦┐8u]I╧H█Ґ░ЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘PЩЪЪ3юH╩Hг┘░ЩЪЪHг┘PЩЪЪH┴²XЩЪЪAЪяH╨H╩I┴п3юH▀Ґ─ЩЪЪH▀u@H▀MHЪсH┐дПH╨I╨H█ҐЪЪЪH╧H█╣ЮЧЪЪI╦I┴я3юH█]ПР┘ЪЪЪH┴$AЪрH┐дH╦H▀UРРYf/┘ЮЧЪЪ├ПH▀∙`ЩЪЪH▀┘hЩЪЪЖгtlH╦┐8u]I╧H█Ґ░ЩЪЪЎЪЪЪЪH╨┘┐ H╧L█┘`ЩЪЪ3юH╩Hг┘░ЩЪЪHг┘`ЩЪЪH┴²hЩЪЪAЪяH▀█xЩЪЪH▀┘pЩЪЪ▀ ;█НH▀E0H╨H┴сРР]Рf/бvH╦РРYпM┘Ж▌ХI┐Ч▄│H▀EXH┐ЮtH╘┘k╦H█PL;Р▄YL┴Р3шH+пH┐БHВзIжH┘ю├║HЎH▀}PL▀┘┬ЩЪЪРH▀uX(ыРAьР\зРYбРYъРXьРчHЪцH;ьrыH▀²┬ЩЪЪH█4цHВф└H▀]XH▀uPH▀Ґ┬ЩЪЪ(ы(ЯD(иD(АР$ф(бР|фР\зР\РРD\йРD\БРDTф РDlф0Р,гРDDгРD\г РDtг0flгfDDгfD\г(fDtг8fшfdфf|фfDTф(fDlф8fюfЖfEиfEДfYэfYХfYВfDYюfEYйfDYьfEYЕfAYфfXщfAXПfEXкfDXЮ)ц)tцD)Lц D)dц0H┐юH;б┌ЪЪЪИХH▀]XH▀uPH▀Ґ┬ЩЪЪ(а(ИD(аD(ыD(ЙР\бР\ЙРD\бРD\зРфРtфРDLф РDdф0fEМf\фA(ЕftфA(ЩfDLф(E(уfDdф8fюfМfEюfEшfYцfY$гfY|гfDYTг fDYlг0fYНfEYаfEYэfXдfXОfEXбfEXщ)ц)lцD)Dц D)\ц0H┐юH;б┌'ЪЪЪI;жsCH╦H▀]PH▀╣┬ЩЪЪРH▀EX(ьРжР\зРYйРYсРXыРпHЪбI;жrзH▀E0ЪаH▀∙xЩЪЪL▀╔ЧЪЪH▀²`ЧЪЪL▀Ґ╗ЧЪЪР┴
ИнH▀┘hЩЪЪ┐8└┘ЭЪЪИ├H╦H▀UР┘ЬЧЪЪРРY
f/а├^H▀∙xЩЪЪH▀┘pЩЪЪ▀;█FH▀M0H╦H▀]I┐ЧР²ЪЪЪР)РY(Р)Р#▄ТH▀EXH┐ЮtH╘┘ч╦H█HL;Я▄лL┴Я3шH+хH┐АHВыIнH┘ю├ЗHЎH▀}XL▀EPРH┴Ж(йHЎРYл(цРYцРYРAYLщРYмР^хРA4ьР\ЯР4ъHЪцH;ьrЎI█\еHВц└┘L▀MX(мL▀UP(РH╩РYТ(аРY3D(цРDYцРA|еD(рfA|еH┴чfЖD(РfYВРDYтРDYТРDYfаH┴ъfYПРDY7fEюI┴ьfA^ПD(бD(ЦРDYдРDYЦРEYРEбfELбfD\нРE\е(СfE\еfEрE)аfEYсРYСfDYпРE|е РELе0fE|е(fELе8fEЖfEюfEYВfEYаfDYПfDYю(цРYцfEДfE^тfЖfюfD^ЖfD^юРElбfElбРA|б fE\ЙРETб0fA|б(fETб8fA\ЧfE\пE)lаA)|а E)Tа0H┐юH;а┌▀ЧЪЪИPL▀MX(мL▀UP(РH╩РYТD(йРDYлРY3D(БD(ЗРDYДH┴чРDYЭРDYH┴ъI┴ьfЖ(аfAYtеРDY'РEY8fEи(ШfEДD(сfEЪD(КfEYLеfEYdе fEY|е0РYШРDYсРDYКfаfYПfDYхfDYЮfDYЬ(цРYцfЪf^ВfEрfEМfюfE^йfE^ЕfD^ЬРEбfEDбРE\бРEtб fD\фРAtб0fE\бfEtб(fAtб8fE\ыfE\ТfA\ВE)аE)\аE)tа A)tа0H┐юH;а┌╩ЧЪЪI;нscH╦H▀]XH▀uPР░H█t&░H█|'(пH╦РYт(кРYкРYРAYTмРYуР^яР4нР\РР4кHЪаI;нr©H▀┘xЩЪЪЪбL▀╔ЧЪЪH▀²`ЧЪЪL▀Ґ╗ЧЪЪ┴И3M┘Ж▌бI┐Ч▄жH▀E@H┐ЮtH╘┘ю╦H█PL;Р▄ўL┴Р3иH+пH┐БHВзIжH┘юv#H▀]@H▀uXH▀}PРнР\оРкHЪаH;хrИH▀MPH█аHВц└°H▀M@H▀]XH▀uPРцР\цРlц Р|ц0РфРTфРdф Рtф0fDфfTфfdф(ftф8fLцf\цflц(f|ц8f\хf\зf\Лf\Ч)а)\а)lа )|а0H┐юH;б┌rЪЪЪКvH▀M@H▀]XH▀uPH█Є&РцРLцРTц Р\ц0fDцfLцfTц(f\ц8f\фf\Lфf\Tф f\\ф0)а)Lа)Tа )\а0H┐юH;бr·I;ж┐~H▀E@H▀MXH▀]PРяР\сРпHЪбI;жrИM┘ЖH▀UH▀┘ЬЧЪЪH┴~%I┐Ч▌IаФH╦L┴РH▀}PH▀uXЪпH▀┘xЩЪЪH╨П?H▀M0L▀╔ЧЪЪH▀²`ЧЪЪL▀Ґ╗ЧЪЪгH┴IВг┘hHВц┘IВд┘╣H▀┘@ЪЪЪH┴дH╨H┴сI┴уI┴вH┴жI┴пI┴рH▀E≤L▀²xЪЪЪL▀eЮL▀╣hЪЪЪH┴BH╨H▀MьH▀}┬L▀MпL┴[0H┴сM┴e8I┴уM┴w@I┴вH┴NH┴жI┴x I┴пM┴J(I┴рH▀E═L▀]хL▀e╗L▀u─H┴BH╨H▀█XЪЪЪH▀ҐpЪЪЪL▀█`ЪЪЪL┴[0H┴сM┴e8I┴уM┴w@I┴вH┴NH┴жI┴x I┴пM┴J(I┴рL▀]╟L▀╔PЪЪЪL▀uюH▀E░H▀M╦H▀}ХL▀█HЪЪЪL┴[0M┴e8M┴w@H▀²8ЪЪЪL▀╔ЪЪЪL▀ґ ЪЪЪL▀╣(ЪЪЪL▀Ґ0ЪЪЪH┴BH┴NI┴x M┴J(H┴Л]ц3иИЎЕЪЪ3юИwЙЪЪ3рИLНЪЪH▀MX3юH▀²┬ЩЪЪH▀аH┴цHЪюI;фrПИ╛ДЪЪ3рИ│ЯЪЪ3рИВЪЪ3иИ┼ШЪЪ3рИ╛ЩЪЪH▀MX3юH▀]PH▀аH┴цHЪюI;фrПИЯЩЪЪL┴ФD┴БH┐Ф┐БHяНрH╧3юаФР│нH▀ҐЧЪЪЪяI┐ДЧHг┘ЧЪЪL┴╔ЧЪЪИЧЩЪЪH┴ч┴зH┐Ф┐БHяНрH╧3юаФР│нH▀ҐHЧЪЪЪяH┐ЦЧHг┘HЧЪЪH┴²`ЧЪЪИ╔ЩЪЪL┴ЧD┴ЗH┐Ф┐БHяНрH╧3юаФР│нH▀Ґ░ЧЪЪЪяI┐ГЧHг┘░ЧЪЪL┴Ґ╗ЧЪЪИKЩЪЪH▀MH3юH▀аI┴TеHЪюI;фrОИmБЪЪH▀┘hЩЪЪг      ИЩЪЪM┘ЖИ═ЭЪЪH╩РИNВЪЪH╩РИ░СЪЪH╧РИлЙЪЪ░H█t&88880880─+┤ынВО?Ю?{ўGАz└?Э╘ЯрMbP?@млллллЛ?@П?lsopt: simul.lsopt: fails wolfe test 1lsopt: fails wolfe test 2< .01lsopt: trial stepsize =lsopt: using steepest descent, t value x░└╟!├┐(▐▌▄_█Вцлмноф┐├▄█▌▐
ч╟!lsline.flsline_тфXиЪн╘фYиЪГ╗фXбЪй╘фYиЧф╘иYцїфYфЧ┌в╗фXО╗пYиЪц╗йWф╘фXк╗цXфЪч╘тYтЪ┤цЪзЪй╘мYфЪИЪ╘Kх_з\Д╙Yгц┤гyхи╛Х<Mц┴гЪигЩ┐йЬцбЗг┼цкЧи█▌бЪДфцЧбеЧьМўв=мцб=Кцц=ц²ь=лцц=сцд=цгфцЭ┐кЪ┘хлне2▓ы]и╛╛л@пцигЧяцЧйгЧб©гЪ┼ОффЪ┌еЪ░\ИBгЪ┼ЦффЪ┌еЪ░\в █┌гЪ⌡Uэа?:эаЕ?┐┬т"цщ├тЫГв█▐кЪLь╘W-А╘W@┬тфКгЪР
і├█т    еЖ┌    iф├ффАФф1ф(фф      чЕwЪяМг·▒К
lsline.f.strtab.symtab.data.bss.textmpi_parm_lsline_for_write_seq_lis.rodata.rodata.str1.4__STRLITPACK_5_intel_fast_memcpydnrm2_mpi_allreduce_ddot_.L_2il0floatpacket.33.L_2il0floatpacket.27.L_2il0floatpacket.29.rodata.str1.32__STRLITPACK_4for_write_seq_lis_xmit__STRLITPACK_3.L_2il0floatpacket.32__STRLITPACK_2.L_2il0floatpacket.28__STRLITPACK_1.L_2il0floatpacket.30__STRLITPACK_0.L_2il0floatpacket.31.L_2il0floatpacket.26.L_2il0floatpacket.25for_dealloc_allocatablevar$116.0.1var$106.0.1var$96.0.1__STRLITPACK_9.0.1__NLITPACK_0.0.1__NLITPACK_2.0.1__NLITPACK_1.0.1__NLITPACK_3.0.1__STRLITPACK_15.0.1__STRLITPACK_16.0.1__STRLITPACK_10.0.1__STRLITPACK_11.0.1__STRLITPACK_12.0.1__STRLITPACK_13.0.1__STRLITPACK_14.0.1mpi_sgi_private_status_mpi_sgi_private_unweighted_mpi_sgi_private_mpi_sgi_private_inplace_mpi_sgi_private_char_mpi_null_copy_fn_mpi_null_delete_fn_mpi_dup_fn_mpi_comm_null_copy_fn_mpi_comm_null_delete_fn_mpi_comm_dup_fn_mpi_win_null_copy_fn_mpi_win_null_delete_fn_mpi_win_dup_fn_mpi_conversion_fn_null_.note.GNU-stack.eh_frame.trace.rela.text.rela.eh_frame.rela.trace
5`╛юиючЧ`]ю┼юєюрю5%D'dw     ╚(ф)пЫ*%r╘(нь+     *$5WPipE─ё─(        ─┘
─щ
%П
'(#\,f0┐,■8їе%4≤K≤╧%л'Л@Ъ
8,BH┐сЕ+*0:hxґ%╪'эО#+GQ*bs≥┬ж%Е' ^░}hЕ≤=≤б`hyp∙─Т─y─юpЕ─1(╛юБ`6i -╣ -!-s!p├!≤≥!≤ &,&00@╦%Ь   %╟!Oх(═W2h)hФ2п)(7Ь)GЬ)═Q≤*Ю
x,~XЖ00        c&:      r>:0
hessupd.o/      1361899272  10100 26113 100644  7904      `
ELF> @@ЯЪ!&QYh▌²ґ ╪0И Э"5I ](q0┘8≥@ґHаPрTЦXТ\BR  \
cn
},РЪ6═
?wктзРЪРЪ9РЪJРЪcРЪy▀÷╚бшЛ*UH┴ЕH│Л0H▀EH┴²hЪЪЪL┴}┬I┴оL┴u─I┴жL┴ґpЪЪЪI┴УL┴╔xЪЪЪI┴ЭL┴M░A▀Ж└XH╦┐8┘EI╨H█ҐПЧЪЪЎЪЪЪЪH╨┘┐ H╧L█┘пЧЪЪ3юI╧Hг┘ПЧЪЪHг┘пЧЪЪL┴█ьЧЪЪAЪрI╦H█ҐПЧЪЪHЎH█∙ЮЧЪЪ3юH╧Hг┘ЮЧЪЪH┴█ХЧЪЪAЪпA▀I╦HЎH█ҐПЧЪЪ3юH█∙0ЪЪЪ┴█0ЪЪЪ┴MХAЪпI╦H█ҐПЧЪЪHЎH█∙8ЪЪЪ3ю┴²8ЪЪЪAЪпE▀MI╦HЎH█ҐПЧЪЪ3юH█∙@ЪЪЪD┴█@ЪЪЪD┴█HЪЪЪAЪпD▀█HЪЪЪ▀MХКA▀E▀M;ы}Aы┴ь;ы▄>L┴РH┐БA┴пA┐ЮL┴}╟D┴E╦┴Uю┴MХ┴]ьL┴mЮA┴еL▀}░H▀UЖ└/H╨┐:┘I╧H█ҐПЧЪЪЎЪЪЪЪH╨┘┐ H╧L█┘HЪЪЪ3юH╩Hг┘ПЧЪЪHг┘HЪЪЪH┴²PЪЪЪAЪяH╩H█ҐПЧЪЪHЎH█∙XЪЪЪ3юH╧Hг┘XЪЪЪH┴█`ЪЪЪЪсH╧H█ҐПЧЪЪHЎH█U≤3юD┴m≤Ъя▀MьH╩HЎH█ҐПЧЪЪ3юH█U═┴M═ЪсH▀M╟I╦HЎH█ҐПЧЪЪ3юH█U╗▀┴]╗AЪпH▀UЮD┴К▀D;Й~D┴К+зI╧L┴ГL┴ЧH╨3юH█MПD█DD┴EПAЪяH╨I╧L┴ГL┴ЖL┴ЫI┴п3юAЪяH┐дПH╨I╩H█}хH╧H█uпI╦I┴я3юL█UТРEхL┴$AЪсH┐дI╧L┴ГL┴ЧH╨3юH█MЬD█DD┴EЬAЪяHcшH╧H▀UпIc4$H┘ЖH┴TыЬ▌╪H┐Ч▄b▀Uю┘рt┐}╦┘Q╨┴сH█KH;Я▄=┴ЯВзя┐Б+й3рHcиH┘шv&РEп(хРAYвРAжР\яРAжHЪбH;сrъРEпI█ъfюHВб└·РAъРA\ъРAlъ РA|ъ0fALъfA\ъfAlъ(fA|ъ8fYхfYьfYХfYЬA(чA(dчA(tч E(Dч0f\яf\Цf\УfD\гA)чA)dчA)tч E)Dч0H┐цH;ы┌mЪЪЪКx░░H█©A(ъA(\ъA(lъ A(|ъ0fYхfYьfYХfYЬA(чA(dчA(tч E(Dч0f\яf\Цf\УfD\гA)чA)dчA)tч E)Dч0H┐цH;ыr▒H;нs&РEп(хРAYоРAнР\яРAнHЪаH;нrъAЪмD;mХ█ШШЪЪL▀mЮ▀]ьL▀}╟E▀ME┘и└▒I╦L┴ГH╨H╧3юH▀u░AЪпA▀$┘ю▌[HcхH┐Ы▄OL┴ЖH┐Ф┘ЖtВф┘8j^┴РH█zH;о▄&ВчП┐Ф+ф3ЖHcюH┘рvH▀}░РAЖРYВРAЖHЪфH;РrГH▀u░H█<жHВгtpH▀u░РжРLжРTж Р\ж0fDжfLжfTж(f\ж8fAYжfAYLжfAYTж fAY\ж0A)жA)LжA)Tж A)\ж0H┐бH;пr√КRH▀u░A(жA(LжA(Tж A(\ж0fYжfYLжfYTж fY\ж0A)жA)LжA)Tж A)\ж0H┐бH;пr╡H;аsH▀U░РAфРYбРAфHЪюH;аrГE▀?A;ъ▄лL┴ПH┐Ю┴ф┐ФL┴mЮ┴u╟┴E╦┴]ьL▀m░H▀EЮD┴Ш▀D;Ь~D┴Ш+ьI╧L┴ГL┴НH╨3юH█MюD█DD┴EюAЪяH╨I╧L┴ГL┴ЖL┴ИI┴п3юAЪяH┐дПH╨I╩H█}хH╧H█uпI╦I┴я3юL█UТРEхL┴$AЪсH┐дHcшI╦I╨L┴ГL┴НH╨3юH█MдРADьЬD█LD┴MдР\EпРEпAЪрIc$H┘ш▌ H┐Ш▄м▀E╦┘юt┐}╟┘╪╦┴аH█QH;з▄╗┴зВьб┐Ю+п3юHcрH┘иv#РEп(хРAYLеРAXфРAфHЪюH;аrБРEпI█DмfюH╘└▌РALмРATмРA\м РAdм0fALмfATмfA\м(fAdм8fYхfYпfYьfYЮfAXнfAXTнfAX\н fAXdн0A)нA)TнA)\н A)dн0H┐аH;й┌yЪЪЪКhH█t&A(LмA(TмA(\м A(dм0fYхfYпfYьfYЮfAXнfAXTнfAX\н fAXdн0A)нA)TнA)\н A)dн0H┐аH;йr²H;сs#РEп(хРAYLуРAXжРAжHЪбH;сrБAЪгD;}ь▌QЩЪЪH▀²hЪЪЪL▀╔xЪЪЪL▀ґpЪЪЪL▀u─L▀}┬H┴Л]ц3иИШЪЪ3юИжЭЪЪ3рК≈H█v88    88   pathei-lsopt: in hessupd; jmin, jmax, mupd:j,jfin,jmin = pathei-lsopt: in hessupd; loop x░└═
├┐▐▌
█
▄:
цлмноф┐├▄█▌▐
Е═
    hessupd.fhessupd_Шб9шDмиИъЯцбЪ▀╚Д╪дЪ┼К░├дЪI²гЪ└цЪ┐лЪ┘вЪeмЪ²гЪвл╘X░эЪ`нЪDяЪ░гЪкэцЖбЪ▀╚Д╪дЪ┼лъфЩхиЩбPдЦугЪ└дЪ┐кЪyсЪYмЪдЦугЪлАфшф
гhessupd.f.strtab.symtab.data.bss.textmpi_parm_hessupd_for_write_seq_lis.rodata.rodata.str1.4__STRLITPACK_3for_write_seq_lis_xmit__STRLITPACK_2.rodata.str1.32__STRLITPACK_1__STRLITPACK_0dostore_ddot_mpi_allreduce_hessupd_$ALPHA.0.1__STRLITPACK_7.0.1__STRLITPACK_8.0.1__STRLITPACK_9.0.1__STRLITPACK_10.0.1__STRLITPACK_11.0.1__STRLITPACK_12.0.1__STRLITPACK_13.0.1__STRLITPACK_14.0.1__STRLITPACK_15.0.1__STRLITPACK_16.0.1__NLITPACK_0.0.1__NLITPACK_1.0.1__NLITPACK_3.0.1__NLITPACK_2.0.1mpi_sgi_private_status_mpi_sgi_private_unweighted_mpi_sgi_private_mpi_sgi_private_inplace_mpi_sgi_private_char_mpi_null_copy_fn_mpi_null_delete_fn_mpi_dup_fn_mpi_comm_null_copy_fn_mpi_comm_null_delete_fn_mpi_comm_dup_fn_mpi_win_null_copy_fn_mpi_win_null_delete_fn_mpi_win_dup_fn_mpi_conversion_fn_null_.note.GNU-stack.eh_frame.trace.rela.text.rela.eh_frame.rela.traceJ#]%}░      ╨&кч
&
0&Aa&k Я#%$(7a&r0┘
ё&Є8м&в@В&H3'CP_Ti(┬T▓)═\ўXя'АP'PXн'чPЗT(#T-);\IXoy'┴P $,$0000000@X#≤!≤  &≤═
Q8`Y2≤@²2ь BЬRЬ═\≤Х─┴c  хnя     }И0
cubic.o/        1361899272  10100 26113 100644  2528      `
ELF>═@@
       ЯЪ          $1▐÷╘╟
╩й* UH┴ЕI╩РAI╦Р(бР\Р\РAYР^ьР2HЎ(н(фРA L▀UH▀EРAX     Р\к(ИРXаT.f/ЕrB(А(чРYАРAYfEОюР\ЦfD/дvfD/vH▀H┴fA/П┐cИdРQДКm(ч(ИР^ыРAYfEОюfA/хР\КrfA/ХrРQАРQМРYЕК9fD/а┌fD/е┌ЗHЎ(АI┴ПW&AW(РQДРQМРYЕР?HЎ(ъ(Т(ОI╦Р\Р\)fD/цW6w(Т(ГР\!AT Р^ЛРYХfA/ХvР\РXфРYР^пРXРКZH╨РXхР\фР\9РAX РYЬР H┴я(ыРA\"TРYЦT:f/ГvР\РYбР^аРXРКH▀H┴РРA_Р]РКfD/vH▀H┴fA/ПrI▀H┴H┴Л]ц░H█t&░H█|'─ЪЪЪЪЪЪЪ@П?x░, ├ф
▓ cubic.fcubic_цнишЫибЫ┌дгГхбцЗх▌└цЪ├▄▐еббцЧ└ицЪ┬█∙╒ц┴бЪ┌лцЭцлЩн┼≥ицЧ┐┐дЧБ∙ех┴лД▄пcubic.f.strtab.symtab.data.bss.textcubic_.rodata.L_2il0floatpacket.8.L_2il0floatpacket.10.L_2il0floatpacket.11.L_2il0floatpacket.9.note.GNU-stack.eh_frame.trace.rela.text.rela.eh_frame.rela.trace (;Щ$7┼ ,00@h ╗╗$╗ 1х0▐Ь÷ЬH╘@≤ ьж╟ў╗╩Vйn0dgscale.o/      1361899272  10100 26113 100644  3896      `
ELF>Ь@@
       ЯЪ !&>q┌⌠є╣
'.
9H,─5FjфРЪчРЪЗРЪРЪ$РЪ:РЪDVjv█іЇмЕУUH┴ЕH┐ЛPI╧L┴}хI┴оL┴mюI┴УL┴e╦I┴тL┴ЧH╨H╧3юH┴]пL┴uьM┴фH┴ШAЪяHгEХ3юHcH┘р~#РADеРYюРAYгHЪюH;бРXEХРEХrщH┐дПH╨L▀UХH█}╟L┴U╟I╨H╧H█uХI╦I┴я3юL█]ЮL┴$AЪрH┐дI╦I╧H┴ъH█uХL┴ЗH╧3юРAР^EХРEХAЪяHc3рHгEПL┴ЬH┘и~fОюf/rI▀6H┴0HЪбH┐юH;яrХH┘и╦~"РAдРYюРA^гHЪюРXEПРEПH;аrчH┐дПH╨L▀UПH█}╟L┴U╟I╨H╧H█uПI╦I┴я3юL█]ЮL┴$AЪрH┐дHc╦H┘р~jH╧fОюР  РA\дЬ(АРA^\гЬРA^dгЬРYш(ИР^]ПРATеЬРYрРXБР\ЦР^ЛРAlгЬf/еrI▀I┴LгЬHЪюH;б~╗I╦H┴ъL┴ЧH╨H╧3юAЪпH▀]пL▀e╦L▀mюL▀uьL▀}хH┴Л]цH█Є&ЪЪЪЪП?x░l─├▐       █
▄ ┐▌цлмноф
▒─       dgscale.fdgscale_ги&тзь&мзб&г┴░еЪ┴аЪмцЪ┐цЪ╣иънЧбд йЧдиеуЭобЫидЫмцЪ┐цЪ╣ляГшбЪ≈йЪе┤гХЮdgscale.f.strtab.symtab.data.bss.textdgscale_dostore_.rodatampi_allreduce_.L_2il0floatpacket.5dscal___NLITPACK_1.0.1__NLITPACK_0.0.1__NLITPACK_2.0.1__NLITPACK_3.0.1__NLITPACK_4.0.1mpi_sgi_private_status_mpi_sgi_private_unweighted_mpi_sgi_private_mpi_sgi_private_inplace_mpi_sgi_private_char_mpi_parm_mpi_null_copy_fn_mpi_null_delete_fn_mpi_dup_fn_mpi_comm_null_copy_fn_mpi_comm_null_delete_fn_mpi_comm_dup_fn_mpi_win_null_copy_fn_mpi_win_null_delete_fn_mpi_win_dup_fn_mpi_conversion_fn_null_.note.GNU-stack.eh_frame.trace.rela.text.rela.eh_frame.rela.trace
,6▄╒╛╨ьБЖu▀∙ён8HR ,0000@╗     Х!Х&Х─>h 
┬┬┬'≤╗T.Э        ╟9╛Hд0instore.o/      1361899272  10100 26113 100644  7424      `
ELF>@@@ЯЪ!&Zbq│▒$
═0╦
%@<O@u┴² ╠(е0ы8М@HP)X=`Qhepyx█─║└(╨╛нЄБ╪Жд
лп2ьFЮZХn~   ┬
▐ 
╘,РЪ6РЪ?юH╞гьНЬ)^UH┴ЕH│Л`H╦D▀H╧I┴бL┴mПI┴уH▀UH┴]ьH┴ШL┴}ЮM┴оL┴uХM┴фL┴eЬI┴Тг[AгB\D┴XгЪЪЪЪ┐9uTI╨H█Ґ@ЪЪЪЎЪЪЪЪH╨┘┐ H╧L█E╦3юI╧Hг┘@ЪЪЪHгE╦L┴MюAЪрI╧H╦H╨H█Ґ@ЪЪЪL┴█ЪЪЪI╧H┴┘ ЪЪЪЎ[H┴∙0ЪЪЪH╨┘┐ H╧L█┘ЬЧЪЪ3юA╨I╩Hг┘@ЪЪЪL┴∙ЬЧЪЪL┴∙ЪЪЪL┴²ЪЪЪHг┘ЪЪЪ     Hг┘(ЪЪЪ
AЪяI╧I╨H╨┘┐ H█Ґ@ЪЪЪH╧L█Eп3юA▀1Hг┘@ЪЪЪH┴]пAЪр┘ю┘зH╧H█Ґ@ЪЪЪHЎH█Uх3юL┴eхЪя┘ю┘╚H╧H█Ґ@ЪЪЪHЎH█U╟3юL┴m╟Ъя┘ю┘|H╧H█Ґ@ЪЪЪHЎH█∙ПЧЪЪ3юL┴╣ПЧЪЪЪя┘ю┘GH╧H█Ґ@ЪЪЪHЎH█∙хЧЪЪ3юL┴ҐхЧЪЪЪя┘ю┘H▀MI╦HЎH█Ґ@ЪЪЪ3юH█∙юЧЪЪH┴█юЧЪЪAЪп┘ю┘ьH▀M I╦HЎH█Ґ@ЪЪЪ3юH█∙╦ЧЪЪH┴█╦ЧЪЪAЪп┘ю┘·H▀M(I╦HЎH█Ґ@ЪЪЪ3юH█∙╟ЧЪЪH┴█╟ЧЪЪAЪп┘юuhH▀M0I╦HЎH█Ґ@ЪЪЪ3юH█∙╗ЧЪЪH┴█╗ЧЪЪAЪп┘юu2H▀M8I╦HЎH█Ґ@ЪЪЪ3юH█∙═ЧЪЪH┴█═ЧЪЪAЪп┐ЬЪtk┘ю~]H╩I╧H╨┘┐ H█Ґ@ЪЪЪH╧E3ю3ю▀3Hг┘@ЪЪЪAЪяH▀]ьL▀eЬL▀mПL▀uХL▀}ЮH┴Л]цH▀EгI╧I╨H╨┘┐ H█Ґ@ЪЪЪH╧E3ю3юA▀1Hг┘@ЪЪЪAЪрH╦┐8u┴I╧H█Ґ@ЪЪЪЎЪЪЪЪH╨┘┐ H╧E3ю3юHг┘@ЪЪЪAЪяI╩H█Ґ@ЪЪЪЎЪЪЪЪH╨┘┐ H╧L█┘пЧЪЪI╧3юI╨Hг┘@ЪЪЪHг┘пЧЪЪ$L┴∙ьЧЪЪAЪсLc#H╧HЎH█Ґ@ЪЪЪ3юH█∙8ЪЪЪD┴╔8ЪЪЪЪяH╩H█Ґ@ЪЪЪHЎH█∙ЮЧЪЪ3юH╧Hг┘ЮЧЪЪH┴█ХЧЪЪЪсH╩H╧HЎH█Ґ@ЪЪЪ3юH█U─Hc[┴]─ЪяI╧H█Ґ@ЪЪЪЎЪЪЪЪH╨┘┐ H╧E3ю3юHг┘@ЪЪЪAЪяI╨H█Ґ@ЪЪЪЎЪЪЪЪH╨┘┐ H╧L█E═3юI╧Hг┘@ЪЪЪHгE═L┴M╗AЪрH┴ыI╦I╞лHЎH█Ґ@ЪЪЪ3юH█U┬H┴M┬AЪпH╧H█Ґ@ЪЪЪHЎH█U░3юD┴e░ЪяH╧H█Ґ@ЪЪЪHЎH█U≤3ю┴]≤ЪяИ-ЩЪЪH█Є&88
888   00    06H$
H$78 8       8pathei: in instore OPWARMIunknownformattedsequential * pathei rectmp   opening direct access with recl = x░└ю├%█┐▐▌▄═цлмноф┐├▄█▌▐
▓ю       instore.finstore_йибибЭФМе┤┐┘XўцчУ"цчУ"цчЯ"цчЯ"цчЖ"%CьищDн╪ЭSбЩицЩwinstore.f.strtab.symtab.data.bss.textmpi_parm_xxstore_instore_for_write_seq_lis.rodata.rodata.str1.4__STRLITPACK_11__STRLITPACK_10__STRLITPACK_8__STRLITPACK_7for_open__STRLITPACK_9for_read_seq_lisfor_read_seq_lis_xmitfor_closefor_write_seq_fmt.rodata.str1.32__STRLITPACK_2for_write_seq_fmt_xmit__STRLITPACK_1__STRLITPACK_0for_write_seq_lis_xmit__STRLITPACK_12.0.1__STRLITPACK_13.0.1__STRLITPACK_14.0.1__STRLITPACK_15.0.1__STRLITPACK_16.0.1__STRLITPACK_17.0.1__STRLITPACK_18.0.1__STRLITPACK_19.0.1__STRLITPACK_20.0.1__STRLITPACK_21.0.1__STRLITPACK_22.0.1__STRLITPACK_23.0.1__STRLITPACK_35.0.1__STRLITPACK_24.0.1__STRLITPACK_25.0.1instore_$format_pack.0.1__STRLITPACK_26.0.1__STRLITPACK_27.0.1__STRLITPACK_28.0.1__STRLITPACK_29.0.1__STRLITPACK_30.0.1__STRLITPACK_32.0.1__STRLITPACK_33.0.1__STRLITPACK_34.0.1__STRLITPACK_31.0.1.note.GNU-stack.eh_frame.trace.rela.text.rela.eh_frame.rela.trace
21k4▀⌡     ©
исК5+
n2x6⌠ ю7я(О707/8S7d@▄7√Hф7пP7
X67@`l7vhё2ґ8хp
28/xO1^4~─⌡9╩╛л└ь:Є2:C╪Vt2~:┬д╗4хлЕ4Х<;Jпh;yь▐;═Ю 3,3000@═1Ю!Ю&ЮюZ═Пb2░
P
2Ю
(n~═┬╗≤@╣▐У У      ╘
0
dostore.o/      1361899272  10100 26113 100644  4776      `
ELF>(@@
ЯЪ
!&HW$fw,├4∙@╪кзИ-$A,THgP{▀∙ °їІ
,РЪ6РЪ? nє╡ЬUH┴ЕH│Л E3юL┴eПI┴ЭL┴uьH┴]ЮL┴}пI┴оMc4$M┘ЖL┴mХI┴уMOфH┴╣ЮЧЪЪH┴eЬJ█┘H┐ЮПH+ЮH┐ДПH┴ЮH┴цAЖE└1M┘Ж▌⌡I┐Ч▄▓D┴Р3ю┐БЬHcрH▀█ЮЧЪЪ░H█t&РLаРаРlа0Рdа fLаfDаflа8fdа(fZьfZяfZЭfZУзЧ)┐)|┐H┐юH;бr╘I;жsH▀┘ЮЧЪЪРпРZюС⌠HЪбI;жrЙH©I╧H┴ЗH╧L┴█ХЧЪЪI╧A╨L█┘ЮЧЪЪHcGH█ҐpЪЪЪI╞фH┴┘ ЪЪЪ3ю▀rH╨┘┐ I╩I╪IҐHг┘pЪЪЪL┴∙ЮЧЪЪL┴∙ПЧЪЪL┴²ЬЧЪЪHг┘ЪЪЪL┴╔ЪЪЪHг┘ЪЪЪL┴ґЪЪЪAЪяE▀I©AЪиI╩H┴]╦L┴ШE3рM┘ЖI╪H█ҐpЪЪЪMLРH╨┘┐ E╞OH╧DL█E╟Mcи3юIЪаIаФA▀sL┴∙pЪЪЪL┴u╟AЪтH╩I╧H╨┘┐ H█ҐpЪЪЪH╧E3ю3ю▀sHг┘pЪЪЪAЪяИTH╦I╧H┴аH©L┴█0ЪЪЪI╧H┴Ґ`ЪЪЪH█ҐpЪЪЪHcPL█┘(ЪЪЪI╞жH┴∙hЪЪЪH╨┘┐ ▀qH╧3юA╨I╩IҐHг┘pЪЪЪL┴∙(ЪЪЪL┴∙8ЪЪЪL┴²@ЪЪЪHг┘HЪЪЪL┴ґPЪЪЪHг┘XЪЪЪAЪяE▀I©AЪиM┴ШE3рM┘ЖIҐH█ҐpЪЪЪMLРH╨┘┐ E╞OH╧EL█EюMcи3юIаФIЪаL┴uюIЎA▀uL┴∙pЪЪЪH┴]хAЪжI╧I╨H╨┘┐ H█ҐpЪЪЪH╧E3ю3юA▀qHг┘pЪЪЪAЪрA▀$┘р▌©HcбH┐Ь▄ыH▀█ЮЧЪЪH┐А┘иtЖа┘аjY┴нH█~H;г▄╞Вый┐А+я3иHcрH┘ЖvH▀ҐЮЧЪЪfОюСZ▀РоHЪаH;нrЙH▀█ЮЧЪЪZЁ)ЯZLЁ)LЯZTЁ)TЯ Z\Ё)\Я0H┐фH;РrяH;пsH▀█ЮЧЪЪfОюСZ⌠РяHЪбH;пrЙH▀EЬH┴дH▀]ЮL▀eПL▀mХL▀uьL▀}пH┴Л]ц3рИжШЪЪ3рКЁH█vH█©OPWARMDunknownunformatteddirectOPWARMDunknownunformatteddirect8
888A8
888Ax░└ ├▄▌┐▐█фцлмноф┐├▄█▌▐
w        dostore.fdostore_Yй╙MтЪ█гЪхyIхyE\█нЪ╔тЪ█гЪъфЦнdostore.f.strtab.symtab.data.bss.textmpi_parm_xxstore_dostore_.rodata.str1.4__STRLITPACK_7.rodatafor_open__STRLITPACK_6__STRLITPACK_5__STRLITPACK_4for_write_dirfor_close__STRLITPACK_3__STRLITPACK_0__STRLITPACK_2__STRLITPACK_1for_read_dir__STRLITPACK_11.0.1__STRLITPACK_12.0.1__STRLITPACK_13.0.1__STRLITPACK_8.0.1__STRLITPACK_9.0.1__STRLITPACK_10.0.1.note.GNU-stack.eh_frame.trace.rela.text.rela.eh_frame.rela.trace"     ,3 o
y┐вДШ! HNX"s P≤╒
╞ю Ш 
dz·     ©#шЕ" $ ,000000@`
═!═&═ H2юHf        X{` ▀` ═∙
xx
б°:═їзІР0      outstore.o/     1361899272  10100 26113 100644  4584      `
ELF>h@@
ЯЪ
"'@O^
v~кчЯ (0*8=@PHdPxX▄`═hЄpхьБ   ИТ
-РЪ6─m█≤╙аUH┴ЕH│ЛПH╦L┴}ЬI╨I┴ЪH©H┴]ПH┴СL┴uЮM┴нL┴mьM┴еL┴eХI┴тH┴M╗I╧H┴ҐxЪЪЪH█ҐЪЪЪ▀0H╨┘┐ H╧L█┘PЪЪЪ3юI╩Hг┘ЪЪЪHг┘PЪЪЪL┴∙XЪЪЪHг┘`ЪЪЪ
L┴²hЪЪЪHг┘pЪЪЪ       HгE─AЪяH╧I╦H╨┘┐ H█ҐЪЪЪ3юHг┘ЪЪЪ▀1AЪпI╧I╨E▀?H█ҐЪЪЪD┴}┬H╨┘┐ H╧L█E┬E▀9D┴Ч3юHг┘ЪЪЪAЪрH▀H╩HЎH█ҐЪЪЪ3юH█U░H┴M░ЪсI▀$H╩HЎH█ҐЪЪЪ3юH█U≤H┴M≤ЪсH▀M╗I╦HЎH█ҐЪЪЪ3юH█U═▀┴]═AЪпA▀MH╩HЎH█ҐЪЪЪ3юH█U╗┴M╗ЪсA▀H╩HЎH█ҐЪЪЪ3юH█U╟┴M╟ЪсH▀MI╦HЎH█ҐЪЪЪ3юH█U╦▀┴]╦AЪпH▀MI╦HЎH█ҐЪЪЪ3юH█Uю▀┴]юAЪпH▀M I╦HЎH█ҐЪЪЪ3юH█Uх▀┴]хAЪпH▀M(I╦HЎH█ҐЪЪЪ3юH█UпH▀H┴]пAЪпH╩H█ҐЪЪЪD┴ЧH╨┘┐ H╧E3ю3юHг┘ЪЪЪЪсH╩I╧H╨┘┐ H█ҐЪЪЪH╧E3ю3ю▀sHг┘ЪЪЪAЪяH▀]ПL▀eХL▀mьL▀uЮL▀}ЬH┴Л]цOPWARMIsequentialformatted8
88            00    0x░l─├▐┐▌█▄"цлмноф
n─
outstore.f outstore_йи
цСи
бСи
ъС
┘Ж$цМХцМХцМХцМЙЫDьoutstore.f.strtab.symtab.data.bss.textxxstore_outstore_.rodata.str1.4__STRLITPACK_2__STRLITPACK_0for_open.rodata__STRLITPACK_1for_rewindfor_write_seq_lisfor_write_seq_lis_xmitfor_close__STRLITPACK_3.0.1__STRLITPACK_4.0.1__STRLITPACK_5.0.1__STRLITPACK_6.0.1__STRLITPACK_7.0.1__STRLITPACK_8.0.1__STRLITPACK_9.0.1__STRLITPACK_10.0.1__STRLITPACK_11.0.1__STRLITPACK_12.0.1__STRLITPACK_13.0.1__STRLITPACK_14.0.1__STRLITPACK_15.0.1.note.GNU-stack.eh_frame.trace.rela.text.rela.eh_frame.rela.trace
(     R!v
┴ьБ"#;
b$l
 █$≈
(╦$б
0Е$О
8$
@8$B
He$o
P▓$°
X©$и
`Й%
h$.%I
p  , 0@░
п"п'п─@2P vpxхХьХ┬БppЮИО
0Т70      lswri.o/        1361899272  10100 26113 100644  1416      `
ELF>х@@ЯЪ        $1AKR   a
*UH┴ЕH┴Л]цH█Іx░,├ф
>lswri.flswri_цк5lswri.f.strtab.symtab.data.bss.textlswri_.note.GNU-stack.eh_frame.trace.rela.eh_frame.rela.trace ,000@8xx$x1┬A┬HKп@ mR}a∙0lsopt_top.f/    1361899270  10100 26113 100644  24763     `


      subroutine lsopt_top( nn, xx, ff, gg, simul, optline
     $                  , epsx, fmin, epsg
     $                  , iprint, itmax, nfunc, nupdate
     $                  , dd, gold, xdiff
     $                  , loffline
     $                  , ifail )

c     ==================================================================
c     SUBROUTINE lsopt_top
c     ==================================================================
c
c     o uses a set of control variables, their adjoint variables,
c       and a cost function value
c       to compute an improved set of controls with respect to the
c       cost function via a
c       variable-storage Quasi-Newton method
c
c     o Reference: J.C. Gilbert & C. Lemarechal
c                  Some numerical experiments with variable-storage
c                  quasi-Newton algorithms
c                  Mathematical Programming 45 (1989), pp. 407-435
c
c     o started: ??? not reproducible
c
c     o changed: Patrick Heimbach, MIT/EAPS
c                
c     o Version: 2.0, 24-Feb-2000: 
c                    (Version 1.0 is considered as version
c                     starting from which changes were made).
c                   - severe changes in structure including various
c                     bug-fixes to enable multi-optimization runs;
c                   - routine lsoptw incorporated into lsoptv
c                   - optimization iteration loop restructured
c                   - complete restructuring of handling
c                     cold start cases
c                   - number of 3 control flags for error handling
c                     (indic, moderl, ifail) reduced to 1 (ifail)
c                     and homogenized with routine lsline
c
c     o Version: 2.1, 29-Feb-2000:
c                   - handling of case ifail = 6 changed;
c                     leads to stop of optimization
c                     (similar to case ifail = 8)
c                   - logical lphprint included
c
c     ==================================================================
c     SUBROUTINE lsopt_top
c     ==================================================================

      implicit none

ccc#include <blas1.h>
      include 'mpif.h'         ! cheung
c-----------------------------------------
c declare arguments
c-----------------------------------------
      integer nmax
      parameter( nmax = 1372312408 )

      integer nn, iprint, itmax, nfunc, nupdate, ifail

      double precision  xx(nn), gg(nn)
      double precision  dd(nn), gold(nn), xdiff(nn)

      double precision  ff, epsx, fmin, epsg
cph(
      integer phniter0, phniterold
      double precision phff
      COMMON /PH_OPTI/ phniter0, phniterold, phff
cph)

      external simul, optline

c-----------------------------------------
C declare local variables
c-----------------------------------------
      logical cold, lphprint, loffline,sflag,tflag,isforward
      parameter (lphprint = .true.)

      integer mm, mupd, jmin, jmax, indic, isize, REAL_BYTE
      integer i, iiter, ifunc

      double precision   fupd
      double precision   r1, tmin, tmax, tact, gnorm, gnorm0, eps1
      double precision   fold, ys
      double precision   dotdg
      double precision   temp,temp1,temp2
!cheung
      integer mylen, mystart,myend, pidlen, myindx(2)
      integer ierr,myid,nprocs
      integer status(MPI_STATUS_SIZE), requests(2)
      integer time_array_0(8), time_array_1(8)
      real    start_time, finish_time
      common /mpi_parm/ myid,nprocs,mystart,myend
!cheung

      external DDOT, DNRM2, DSCAL
      double precision     DDOT, DNRM2

c-----------------------------------------
c parameters
c-----------------------------------------

      double precision rmin
      parameter( rmin = 1.e-20 )

      character*(*) iform
      parameter(iform='(i5,2x,1pe8.1,1x,i5,4x,1pe11.4,3(2x,1pe8.1))' )

c     ==================================================================
c
c-----------------------------------------
c initialization
c-----------------------------------------
      cold  = .true.
      fupd  = 1.0e+10
      indic = 0
      tmin  = 0.
      tmax  = 1.0e+10
      tact  = 1.0
cph(
      phniterold = 0
cph)
      iiter = 0
      eps1  = 1.0
      ifunc = 0
      ifail = 0
      gnorm0 = 1.
      sflag  = .false.
      tflag  = .false.
c-----------------------------------------
c initialization for dd and dds
c-----------------------------------------

      jmin  = 1
      jmax  = 0

      mm    = nmax
      mupd  = nupdate

      REAL_BYTE = 4
      isize = REAL_BYTE

c-----------------------------------------
c print information
c-----------------------------------------
      if (myid .eq. 0) then
      if (iprint .ge. 1) then
 
         print '(2x,a)',
     $         '==============================================='
         print '(2x,a)',
     $         '===         O P T I M I Z A T I O N         ==='
         print '(2x,a)',
     $         '==============================================='
         print '(a,i11)'
     $        , '  number of control variables.......', nmax
         print '(a,e9.2)'
     $        , '  precision on x....................', epsx
         print '(a,e9.2)'
     $        , '  precision on g....................', epsg
         print '(a,e9.2)'
     $        , '  expected optimal function value...', fmin
         print '(a,i9)'
     $        , '  maximal number of iterations......', itmax
         print '(a,i9)'
     $        , '  maximal number of simulations.....', nfunc
         print '(a,i9)'
     $        , '  information level.................', iprint
         print '(a,i9)'
     $        , '  number of updates.................', nupdate
         print '(a,i9)'
     $        , '  size of used memory...............', 3*nn
      endif

c-----------------------------------------
c check arguments
c-----------------------------------------

      if (nn .le. 0) then
         if (iprint.ge.1) then
            print '(a,i6)'  , '  ERROR : n     = ', nn
         endif
         ifail = 1
         goto 999
      endif

      if (itmax .lt. 0) then
         if (iprint.ge.1) then
            print '(a,i6)'  , '  ERROR : itmax = ', itmax
         endif
         ifail = 1
         goto 999
      endif

      if (nfunc .le. 0) then
         if (iprint.ge.10) then
            print '(a,i6)'  , '  ERROR : nfunc  = ', nfunc
         endif
         ifail = 1
         goto 999
      endif

      if (epsx .le. 0.) then
         if (iprint.ge.1) then
            print '(a,e9.2)', '  ERROR : epsx = ', epsx
         endif
         ifail = 1
         goto 999
      endif

      if (epsg .le. 0.) then
         if (iprint.ge.1) then
            print '(a,e9.2)', '  ERROR : epsg  = ', epsg
         endif
         ifail = 1
         goto 999
      endif

      if (epsg .gt. 1.) then
         if (iprint.ge.1) then
            print '(a,e9.2)', '  ERROR : epsg  = ', epsg
         endif
         ifail = 1
         goto 999
      endif

cph(
      print *, 'pathei: vor instore '
cph)
      endif

      if (myid .eq. 0 ) print *,'bnc: cold start before instore is ',cold,sflag
      call instore( nn, fupd, gnorm0, isize, mupd, jmin, jmax, cold,
     &              sflag, tflag, ifunc, tact )
      if (myid .eq. 0 ) print *,'bnc: cold start after instore is ',cold,sflag


cph(
      phff = fupd
cph)

c-----------------------------------------
c check warm start parameters
c-----------------------------------------
      if (myid .eq. 0 ) then

      if (ifail .ne. 0) then
         if (iprint.ge.1) then
            print '(a)', ' ERROR : reading restart file'
         end if
         ifail = 2
         goto 999
      end if

      if (isize .ne. REAL_BYTE) then
         if (iprint.ge.1) then
            print '(a)', ' ERROR : uncompatible floating point format'
         end if
         ifail = 2
         goto 999
      end if

      if (mupd .lt. 1) then
         if (iprint .ge. 1) then
            print '(a)', ' ERROR : m is set too small in instore'
         endif
         ifail = 2
         goto 999
      endif

      endif
c-----------------------------------------
c cold start or warm restart ?
c-----------------------------------------

      if (cold) then
c--- start if cold start ---
         if (myid .eq. 0) then

         if (lphprint) then
            print '(a)', 'pathei-lsopt: cold start'
         end if

         print *, 'pathei-lsopt vor simul', nmax
         print *, 'pathei-lsopt xx(1), gg(1) ',xx(1), gg(1)
         end if 

         isforward = .false.
         call simul( indic, isforward, nn, xx, ff, gg )

         if (myid .eq. 0) then
         print *, 'pathei: nach simul:   nn, ff = ', nmax, ff
         print *, 'pathei: nach simul:   xx(1), xx(2) = ', xx(1), xx(2)
         print *, 'pathei: nach simul:   gg(1), gg(2) = ', gg(1), gg(2)
         endif

         do i = 1, nn
            xdiff(i) = 1.
         end do

cph(
         if (myid .eq. 0) then
      print *, 'pathei: vor dostore '
cph)
         endif

         call dostore( nn, xx, .true., 1 )
         call dostore( nn, gg, .true., 2 )
         call dostore( nn, xdiff, .true., 3 )
               
         if (myid .eq. 0) then
cph(
      print *, 'pathei: vor gnorm0 '
cph)
         endif

         temp = DNRM2(nn,gg,1)
         call MPI_ALLREDUCE(temp*temp,gnorm0,1,MPI_DOUBLE_PRECISION,MPI_SUM,MPI_COMM_WORLD,ierr)
         gnorm0 = dsqrt(gnorm0)
 
       if (myid .eq. 0) then
cph
      print *, 'pathei: gnorm0 = ', gnorm0
cph)
       endif

         if (gnorm0 .lt. rmin) then
            ifail = 3
           goto 1000
         endif

cph(
         phff = ff
cph)

         if (lphprint .and. myid .eq. 0)
     &        print *, 'pathei-lsopt: cold; first call simul: ff = ', 
     &        phff

c--- end if cold start ---
      else
c--- start if warm start ---
         if (mm .ne. nmax) then
            if (iprint.ge.1) then
               print '(a,i6)'
     $              , '  ERROR : inconsistent nn = ', mm
            endif
            ifail = 2
            goto 999
         endif
         if (mupd .ne. nupdate) then
            if (iprint.ge.1) then
               print '(a,i6)'
     $              , '  ERROR : inconsistent nupdate = ', mupd
            endif
            ifail = 2
            goto 999
         endif
         if (jmin .lt. 1) then
            if (iprint.ge.1) then
               print '(a,i6)'
     $              , '  ERROR : inconsistent jmin = ', jmin
            endif
            ifail = 2
            goto 999
         endif
         if (jmin .gt. mupd) then
            if (iprint.ge.1) then
               print '(a,i6)'
     $              , '  ERROR : inconsistent jmin = ', jmin
            endif
            ifail = 2
            goto 999
         endif
         if (jmax .gt. mupd) then
            if (iprint.ge.1) then
               print '(a,i6)'
     $              , '  ERROR : inconsistent jmax = ', jmax
            endif
            ifail = 2
            goto 999
         endif
         
         if (lphprint .and. myid .eq. 0) then
            print *, 'pathei-lsopt: warm start; read via dostore'
            print *
         endif

         call dostore( nn, xx, .false., 1 )
         call dostore( nn, gg, .false., 2 )
         ff = fupd

cph(
         phff = ff
cph)

         if (lphprint .and. myid .eq. 0) then
             print *, 'pathei-lsopt: warm; first dostore read: ff = ',
     &        ff
             print *, '   bnc-lsopt:                           xx = ',
     &        xx(1),xx(2) 
             print *, '   bnc-lsopt:                           gg = ',
     &        gg(1),gg(2)
         endif

c--- end if warm start ---
      endif

c-----------------------------------------
c print information line
c-----------------------------------------
         if (cold) then

            temp = DNRM2(nn,xx,1)
            call MPI_REDUCE(temp*temp,temp1,1,MPI_DOUBLE_PRECISION,
     & MPI_SUM,0,MPI_COMM_WORLD,ierr)
           
            if(myid .eq. 0)  then
      if (iprint .ge. 1) then
         print '(2a)', '  Itn      Step  Nfun      Objective    '
     $        , 'Norm G    Norm X   Norm (X(k-1)-X(k))'
      end if
         print iform, iiter, tact, ifunc, ff, gnorm0
     $           , dsqrt(temp1), 0.
            endif

            if ( itmax .EQ. 0 ) then
               ifail = 10
               goto 1000
            end if
         end if

c=======================================================================
c begin of loop
c compute x(k+1) out of x(k) + t*d, where  t > 0.
c=======================================================================

      do iiter = 1, itmax

         if (lphprint .and. myid .eq. 0) then
            print *, 'pathei-lsopt: ++++++++++++++++++++++++'
            print *, 'pathei-lsopt: entering iiter =', iiter
         endif
c-----------------------------------------
c store old values
c-----------------------------------------
      do i = 1, nn
         gold(i) = gg(i)
      end do

      fold = ff
cph(
      phniter0 = iiter
      phff = ff
cph)

c-----------------------------------------
c compute new dd and xx
c-----------------------------------------
c

!cheung

!      if (myid.eq.0) then
!      CALL date_and_time(values=time_array_0)
!      start_time = time_array_0 (5) * 3600 + time_array_0 (6) * 60
!     &           + time_array_0 (7) + 0.001 * time_array_0 (8)
!      endif          ! myid = 0

      if (cold) then
!get steepest descent controls
      call lsupdxx(
     &     nn, ifail, lphprint
     &     , jmin, jmax, nupdate
     &     , ff, fmin, fold, gnorm0, dotdg
     &     , gg, dd, xx, xdiff
     &     , tmin, tmax, tact, epsx
     &     )

      endif

!      if(myid.eq.0) then
!        CALL date_and_time(values=time_array_1)
!        finish_time = time_array_1 (5) * 3600 + time_array_1 (6) * 60
!     &              + time_array_1 (7) + 0.001 * time_array_1 (8)

!        write (*, '(8x, 1a, 1f16.6)') 'Clock (sec.) for lsupdxx :',
!     &         finish_time - start_time
!      endif

c-----------------------------------------
c check whether new direction is a descent one
c-----------------------------------------
!      if ( ifail .eq. 4) goto 1000

c-----------------------------------------
c optline returns new values of x, f, and g
c-----------------------------------------

c
      if (.not. cold ) then
      call optline( 
     &     simul
     &     , nn, ifail, lphprint
     &     , ifunc, nfunc
     &     , ff, dotdg
     &     , sflag, tflag, tact, epsx
     &     , dd, gg, xx, xdiff
     &     )

      if (lphprint .and. myid .eq. 0 ) then
      print *, 'pathei-lsopt: ',
     &     'back from optline; ifail = ', ifail
      print *, '   bnc-lsopt: ',
     &     '                   sflag = ',sflag
      print *, '   bnc-lsopt: ',
     &     '                   tflag = ',tflag           
      endif
      endif

      if (.not. sflag .and. ifunc .eq. 0 ) then

      if (lphprint .and. myid .eq. 0 ) then
       print *, 'pathei-lsopt: ',
     &     'dostore 1,2 at iiter ', iiter
       print *,'bnc-lsopt: dostore xx',xx(1),xx(2)
       print *,'bnc-lsopt: dostore gg',gg(1),gg(2)
      endif

      call dostore( nn, xx, .true., 1 )
      call dostore( nn, gg, .true., 2 )

      temp = DNRM2( nn, gg, 1 )
      call MPI_ALLREDUCE(temp*temp,gnorm,1,MPI_DOUBLE_PRECISION,MPI_SUM,
     & MPI_COMM_WORLD,ierr)
      gnorm = dsqrt(gnorm)

c-----------------------------------------
c print information line
c-----------------------------------------
      temp = DNRM2(nn,xx,1)
      call MPI_REDUCE(temp*temp,temp1,1,MPI_DOUBLE_PRECISION,MPI_SUM,0,
     & MPI_COMM_WORLD,ierr)
      temp = DNRM2(nn,dd,1)
      call MPI_REDUCE(temp*temp,temp2,1,MPI_DOUBLE_PRECISION,MPI_SUM,0,
     & MPI_COMM_WORLD,ierr)

      if (myid .eq. 0) then
      if (iprint .ge. 1) then
         print '(2a)', '  Itn      Step  Nfun      Objective    '
     $        , 'Norm G    Norm X   Norm (X(k-1)-X(k))'
      end if
      print iform, iiter, tact, ifunc, ff, gnorm
     $     , dsqrt(temp1), tact*dsqrt(temp2)
      endif

c-----------------------------------------
c check output mode of ifail
c-----------------------------------------
           if ( ifail .eq. 8
     &     .or. ifail .eq. 9) goto 1000

c-----------------------------------------
c maximal number of simulation reached
c no goto in order to update Hessian
c-----------------------------------------
         if (ifail .eq. 6) ifail = 0

c-----------------------------------------
c NOTE: stopping tests are now done
c after having updated the matrix,
c so that update information can be stored
c in case of a later warm restart
c-----------------------------------------
c-----------------------------------------
c  compute new s, y
c-----------------------------------------
         do i = 1, nn
            xdiff(i) = tact*dd(i)
         end do

c-----------------------------------------
c compute ys
c-----------------------------------------
         do i = 1, nn
            gold(i) = gg(i)-gold(i)
         end do

         temp = DDOT( nn, gold, 1, xdiff, 1 )
         call MPI_ALLREDUCE(temp,ys,1,MPI_DOUBLE_PRECISION,MPI_SUM,
     & MPI_COMM_WORLD,ierr)

         if (ys .le. 0.) then
            ifail = 4
            if (myid .eq. 0 ) then
            print *, 'bnc: max value of abs(gold) = ',maxval(abs(gold)) 
            print *, 'pathei-lsopt: ys < 0; ifail = ', ifail
            endif
            goto 1000
         endif

cph(
cph-----------------------------------------
cph at this point it is clear that xdiff 
cph provides a true optimized solution;
cph i.e. take new gradient gg to compute new dd
cph-----------------------------------------
cph)

c-----------------------------------------
c update pointers for hessupd
c-----------------------------------------
      if (nupdate .gt. 0) then

         if (jmax .eq. 0) then
            jmax = jmax+1
            if (lphprint .and. myid .eq. 0) 
     &           print *, 'pathei-lsopt: ',
     &           'first pointer update after cold start; ',
     &           'iiter, jmin, jmax = ', iiter, jmin, jmax
         else
            jmax = jmax+1
            if (jmax .gt. nupdate) jmax = jmax-nupdate

            if (jmin .eq. jmax) then
               if (lphprint .and. myid .eq. 0) 
     &              print *, 'pathei-lsopt: pointers updated for ',
     &              ' iiter, jmin, jmax = ', iiter, jmin, jmax
               jmin = jmin+1
               if (jmin .gt. nupdate) jmin = jmin-nupdate
            end if
         end if

c-----------------------------------------
c compute sbar, ybar store rec = min 4,5
c-----------------------------------------
         r1 = sqrt(1./ys)
         call DSCAL( nn, r1, xdiff, 1 )
         call DSCAL( nn, r1, gold, 1 )

      if (lphprint .and. myid .eq. 0) 
     &        print *, 'pathei-lsopt: dostore at iiter, jmin, jmax ', 
     &        iiter, jmin, jmax

         call dostore( nn, gold, .true., 2*jmax+2 )
         call dostore( nn, xdiff, .true., 2*jmax+3 )

c-----------------------------------------
c compute the diagonal preconditioner
c use dd as temporary array
c-----------------------------------------
         call dgscale( nn, gold, xdiff, dd, rmin )

      endif

c-----------------------------------------
c test convergence and stopping criteria
c-----------------------------------------
      eps1 = gnorm / gnorm0
      gnorm0 = gnorm
      if (eps1 .lt. epsg) then
         ifail = 11
         goto 1000
      endif

c=======================================================================
c return of loop
c=======================================================================
     
      endif  !end sflag check
      
      end do

      iiter = iiter - 1
      ifail = 5

c-----------------------------------------
c loop exit
c-----------------------------------------
 1000 continue
!      nfunc = ifunc
      epsg = eps1

c-----------------------------------------------------------------------
c compute dd(i+1), xx(i+1) based on latest available gg(i), xx(i)
c for offline version
c-----------------------------------------------------------------------

      if (.not. sflag .and. ifunc .eq. 0 ) then
         call lsupdxx(
     &        nn, ifail, lphprint
     &        , jmin, jmax, nupdate
     &        , ff, fmin, fold, gnorm0, dotdg
     &        , gg, dd, xx, xdiff
     &        , tmin, tmax, tact, epsx
     &        )

       endif
c        Save xx(i+1) to file for offline version.
         if (myid .eq. 0) then
             print *, '   bnc-lsopt: write control ',
     &       xdiff(1),xdiff(2)
         endif

         call optim_write_control( nn, xdiff )
         
c-----------------------------------------------------------------------
c save data for warm start
c-----------------------------------------------------------------------

      if ( ifunc .eq. 0 ) sflag = .not. sflag
      if (cold) sflag = .false.
      call outstore( nn, ff, gnorm0, nupdate, jmin, jmax, sflag, tflag, ifunc, tact )


c-----------------------------------------
c print final information
c-----------------------------------------
      if (iprint .ge. 5 .and. myid .eq. 0 ) then
         print *
         print '(a,i9)'
     $        , '  number of iterations..............', iiter
         print '(a,i9)'
     $        , '  number of simulations.............', ifunc
         print '(a,e9.2)'
     $        , '  relative precision on g...........', epsg
      end if


      temp = DNRM2(nn,xx,1)
      call MPI_REDUCE(temp*temp,temp1,1,MPI_DOUBLE_PRECISION,MPI_SUM,0,
     & MPI_COMM_WORLD,ierr)
      temp = DNRM2(nn,gg,1)
      call MPI_REDUCE(temp*temp,temp2,1,MPI_DOUBLE_PRECISION,MPI_SUM,0,
     & MPI_COMM_WORLD,ierr)

      if (iprint.ge.10 .and. myid .eq. 0 ) then
         print *
         print '(a,e15.8)'
     $        , '  cost function...............', ff
         print '(a,e15.8)'
     $        , '  norm of x...................', dsqrt(temp1)
         print '(a,e15.8)'
     $        , '  norm of g...................', dsqrt(temp2)
      end if

c-----------------------------------------
c print error message
c-----------------------------------------
 999  continue

      if (myid .eq. 0 ) then

      if (ifail .ne. 0) then
         if (iprint .ge. 5) then
            print *
            print '(a)', '  optimization stopped because :'
            if (ifail .lt. 0) then
               print '(2x,a8,I3,a)', 'ifail = ', ifail
     $              , '    user request during simul'
            else if (ifail .eq. 0) then
               print '(2x,a8,I3,a)', 'ifail = ', ifail
     $              , '    optimal solution found'
            else if (ifail .eq. 1) then
               print '(2x,a8,I3,a)', 'ifail = ', ifail
     $              , '    an input argument is wrong'
            else if (ifail .eq. 2) then
               print '(2x,a8,I3,a)', 'ifail = ', ifail
     $              , '    warm start file is corrupted'
            else if (ifail .eq. 3) then
               print '(2x,a8,I3,a)', 'ifail = ', ifail
     $              , '    the initial gradient is too small'
            else if (ifail .eq. 4) then
               print '(2x,a8,I3,a)', 'ifail = ', ifail
     $              , '    the search direction is not a descent one'
            else if (ifail .eq. 5) then
               print '(2x,a8,I3,a)', 'ifail = ', ifail
     $              , '    maximal number of iterations reached'
            else if (ifail .eq. 6) then
               print '(2x,a8,I3,a)', 'ifail = ', ifail
     $              , '    maximal number of simulations reached'
            else if (ifail .eq. 7) then
               print '(2x,a8,I3,a)', 'ifail = ', ifail
     $              , '    the linesearch failed'
            else if (ifail .eq. 8) then
               print '(2x,a8,I3,a)', 'ifail = ', ifail
     $              , '    the function could not be improved'
            else if (ifail .eq. 9) then
               print '(2x,a8,I3,a)', 'ifail = ', ifail
     $              , '    optline parameters wrong'
            else if (ifail .eq. 10) then
               print '(2x,a8,I3,a)', 'ifail = ', ifail
     $              , '    cold start, no optimization done'
            else if (ifail .eq. 11) then
               print '(2x,a8,I3,a)', 'ifail = ', ifail
     $              , '    convergence achieved within precision'
            end if
            print *
         else if (iprint .ge. 1) then
            print *
            print '(a,i9)'
     $           , '  after optimization ifail..........', ifail
            print *
         end if
      end if

      end if
c-----------------------------------------
c end of subroutine
c-----------------------------------------
      return

      end

lsupdxx.f/      1361899256  10100 26113 100644  5032      `


      subroutine lsupdxx(
     &     nn, ifail, lphprint
     &     , jmin, jmax, nupdate
     &     , ff, fmin, fold, gnorm0, dotdg
     &     , gg, dd, xx, xdiff
     &     , tmin, tmax, tact, epsx
     &     )

c     ==================================================================
c     SUBROUTINE lsupdxx
c     ==================================================================
c
c     o conceived for variable online/offline version
c       computes - new descent direction dd based on latest 
c                  available gradient
c                - new tact based on new dd
c                - new control vector xx needed for offline run
c
c     o started: Patrick Heimbach, MIT/EAPS
c                29-Feb-2000: 
c
c     o Version 2.1.0, 02-Mar-2000: Patrick Heimbach, MIT/EAPS
c
c     ==================================================================
c     SUBROUTINE lsupdxx
c     ==================================================================
c


      implicit none
      include 'mpif.h'

c-----------------------------------------
c declare arguments
c-----------------------------------------
      integer nmax
      parameter( nmax = 1372312408 )
      integer nn, jmin, jmax, nupdate, ifail
      double precision    ff, fmin, fold, gnorm0, dotdg,dnorm
      double precision    gg(nn), dd(nn), xx(nn), xdiff(nn)
      double precision    tmin, tmax, tact, epsx,temp
      logical lphprint

      integer pidlen,myindx(2)
      integer status(MPI_STATUS_SIZE),ierr
      integer myid, nprocs,mystart,myend
      common /mpi_parm/ myid,nprocs,mystart,myend

c-----------------------------------------
C declare local variables
c-----------------------------------------
      integer i
      double precision    fdiff, preco

      real*8,allocatable:: vv(:)
      integer, allocatable:: displs(:)
      integer, allocatable:: counts(:) 
      double precision     DDOT,DNRM2
      external DDOT,DNRM2

c     ==================================================================

c-----------------------------------------
c use Fletchers scaling
c and initialize diagional to 1.
c-----------------------------------------
c
      if ( ( jmax .eq. 0 ) .or. (nupdate .eq. 0 ) ) then

         if (jmax .eq. 0) then
            fold = fmin
!            if (lphprint .and. myid .eq. 0) 
!     &           print *, 'pathei-lsopt: using fold = fmin = ', fmin
         end if
!         fdiff = fold - ff
         fdiff = .01D0*ff
         if (jmax .eq. 0) fdiff = ABS(fdiff)
         
         preco = 2. * fdiff / (gnorm0*gnorm0)
         do i = 1, nn
            dd(i)    = -gg(i)*preco
         end do

         if (lphprint .and. myid .eq. 0) 
     &        print *, 'pathei-lsopt: first estimate of dd via ',
     &        '.01*ff'

c-----------------------------------------
c use the matrix stored in [diag]
c and the (y,s) pairs
c-----------------------------------------

      else

            do i = 1, nn
               dd(i) = -gg(i)
            end do

            if (jmax .gt. 0) then
               call hessupd( nn, nupdate, dd, jmin, jmax, xdiff, 
     &              lphprint )
            else
               if (lphprint .and. myid .eq. 0) 
     &              print *, 'pathei-lsopt: no hessupd for first optim.'
            end if

      endif

c-----------------------------------------
c check whether new direction is a descent one
c-----------------------------------------

         temp  = DDOT( nn, dd, 1, gg, 1 )
         call MPI_ALLREDUCE(temp,dotdg,1,MPI_DOUBLE_PRECISION,MPI_SUM,MPI_COMM_WORLD,ierr)
         
         if (dotdg .ge. 0.0) then
            ifail = 4
            if(myid .eq. 0) print *,'bnc: ifail = 4, dotdg > 0.'
            goto 999
         end if

c----------------------------------
c declare arguments
c----------------------------------

      tmin = 0.
      do i = 1, nn
         tmin = max( tmin, abs(dd(i)) )
      end do
      temp = tmin
      call MPI_ALLREDUCE(temp,tmin,1,MPI_DOUBLE_PRECISION,MPI_MAX,MPI_COMM_WORLD,ierr)
      tmin = epsx/tmin

c----------------------------------
c make sure that t is between
c tmin and tmax
c----------------------------------

      tact  = 1.0
      tmax = 1.0e+10
      if (tact.le.tmin) then
         tact = tmin
         if (tact.gt.tmax) then
            tmin = tmax
         endif
      endif

      if (tact .gt. tmax) then
          tact = tmax
          ifail = 7
      endif

c----------------------------------
c compute new x
c----------------------------------
         temp = DNRM2(nn,dd,1)
         call MPI_ALLREDUCE(temp*temp,dnorm,1,MPI_DOUBLE_PRECISION,MPI_SUM,MPI_COMM_WORLD,ierr)
         dnorm = dsqrt(dnorm)

      if (myid .eq. 0) then
      print *,'bnc: initial stepsize tact = ',tact
      print *,'bnc: norm of dd            = ',dnorm
      endif

      do i = 1, nn
         xdiff(i) = xx(i) + tact*dd(i)
      end do

c----------------------------------
c save new x to file for offline version
c----------------------------------

 999  continue

      return

      end
lsline.f/       1361899257  10100 26113 100644  7069      `


      subroutine lsline( 
     &     simul
     &     , nn, ifail, lphprint
     &     , ifunc, nfunc
     &     , ff, dotdg
     &     , sflag, tflag, tact, epsx
     &     , dd, gg, xx, xdiff
     &     )

c     ==================================================================
c     SUBROUTINE lsline
c     ==================================================================
c
c     o line search algorithm for determining control vector update;
c       After computing updated control vector from given gradient,
c       a forward and adjoint model run are performed (simul.F)
c       using the updated control vector.
c       Tests are then applied to see whether solution has improved.
c
c     o Reference: J.C. Gilbert & C. Lemarechal
c                  Some numerical experiments with variable-storage
c                  quasi-Newton algorithms
c                  Mathematical Programming 45 (1989), pp. 407-435
c
c     o started: ??? not reproducible
c
c     o changed: Benny Cheng, NASA/JPL
c
c     ==================================================================
c     SUBROUTINE lsline
c     ==================================================================


      implicit none
      include 'mpif.h'
c----------------------------------
c declare arguments
c----------------------------------
      integer nmax
      parameter( nmax = 1372312408 )
      integer nn, ifail, ifunc, nfunc
      double precision ff, dotdg, tact, epsx
      double precision xx(nn), gg(nn), xdiff(nn),dd(nn)
      logical tflag,sflag,isforward
      logical lphprint

      external simul

      integer pidlen,myindx(2)
      integer status(MPI_STATUS_SIZE),ierr
      common /mpi_parm/myid,nprocs,mystart,myend
      integer myid,nprocs,mystart,myend
c----------------------------------
c declare local variables
c----------------------------------

      double precision xpara1, xpara2
      parameter( xpara1 = 0.001D0, xpara2=0.9D0 )

      double precision    factor
      parameter( factor = 0.2 )

      double precision    barmax
      parameter( barmax = 0.3 )
      double precision    barmul
      parameter( barmul = 5.0 )
      double precision    barmin
      parameter( barmin = 0.01 )
      double precision    tmax
      parameter( tmax   = 5. )

      integer i 

      double precision    tg, fg, td, ta
      double precision    fa, fpa, fp
      double precision    fnew, fdiff,ftemp
      double precision    z, test, barr
      double precision    left, right, told
      double precision    t,a,b
      double precision    xtemp(nn),gtemp(nn),gold(nn)
      external DDOT,DNRM2
      double precision     DDOT, DNRM2

      double precision    temp,gnorm0
      real*8,allocatable:: vv(:)
      integer, allocatable:: displs(:)
      integer, allocatable:: counts(:)

c----------------------------------
c check parameters
c----------------------------------

      if (  (nn.le.0)
     &     .or. (xpara1.le.0.0) .or. (xpara1.ge.0.5)
     &     .or. (xpara2.le.xpara1)  .or. (xpara2.ge.1.0) ) then
         ifail = 9
         go to 999
      endif

c----------------------------------
c initialization
c----------------------------------

c=======================================================================
c begin of simulation iter.
c=======================================================================

         if (lphprint .and. myid .eq. 0) 
     &        print *, 'lsopt: simul.'

c------------------------------------
c compute cost function and gradient
c------------------------------------

         gold = gg
         temp = DNRM2(nn,gg,1)
         call MPI_ALLREDUCE(temp*temp,gnorm0,1,MPI_DOUBLE_PRECISION,
     & MPI_SUM,MPI_COMM_WORLD,ierr)
         gnorm0 = dsqrt(gnorm0)

         if (sflag ) then

!get estimated cost from hessian update control, this is the trial step
         isforward = .true.
         call simul( 0, isforward,nn, xdiff, fnew, gg )
        
         xtemp = xdiff
! new direction vector dd

         temp = DDOT( nn, xdiff - xx, 1, gold, 1 )
         call MPI_ALLREDUCE(temp,b,1,MPI_DOUBLE_PRECISION,MPI_SUM,
     & MPI_COMM_WORLD,ierr)
         dotdg = b
     
c-----------------------------------------
c quadratic line search
c-----------------------------------------

         a = fnew - ff - b   
         t = -b/(2.*a)
         t = dmax1(t,0.D0) 
         
         if (t <= 1.d-2 ) then
         tact = 1.
         xdiff = xx -  2.*(.01*ff)*gold*tact/(gnorm0*gnorm0)
         tflag  = .false.

         if (lphprint .and. myid .eq. 0)
     &   print *, 'lsopt: using steepest descent, t value ',t,'< .01'

         else

         xdiff = (1. - t)*xx + t*xdiff
         tact = t
         tflag  =  .true.
         endif

          if (lphprint .and. myid .eq. 0)
     &        print *, 'lsopt: trial stepsize =',t

! no hessian update , get new adjoint run
         goto 999

         else

!no line search
!read new ctrl vector and new cost gradient and new cost
         isforward = .false.
         call simul( 0, isforward, nn, xdiff, fnew, gg )

      if (tflag) then

!read in trial step adjoint output
         call simul( ifunc+1, isforward, nn, xtemp, ftemp, gtemp)

         dd = xtemp - xx
c-----------------------------------------
c apply 1st Wolfe test
c-----------------------------------------
         fdiff = fnew - ff

         temp = DDOT( nn, dd, 1, gold, 1 )
         call MPI_ALLREDUCE(temp,dotdg,1,MPI_DOUBLE_PRECISION,MPI_SUM,
     & MPI_COMM_WORLD,ierr)

         if (fdiff .gt. tact*xpara1*dotdg) then
            ifail  = 7
         if (lphprint .and. myid .eq. 0)
     &      print *, 'lsopt: fails wolfe test 1'

         endif

c-----------------------------------------
c 1st Wolfe test 1 ok, apply 2nd Wolf test
c-----------------------------------------
         
         temp = DDOT( nn, dd, 1, gg, 1 )
         call MPI_ALLREDUCE(temp,fp,1,MPI_DOUBLE_PRECISION,MPI_SUM,
     & MPI_COMM_WORLD,ierr)

         if (fp .lt. xpara2*dotdg) then
            ifail = 7
         if (lphprint .and. myid .eq. 0)
     &      print *, 'lsopt: fails wolfe test 2'
         endif

c-----------------------------------------
c 2nd Wolfe test 2 ok, donc pas serieux, on sort
c-----------------------------------------

        if (ifail .eq. 7 ) then

        if (ifunc .lt. nfunc) then
        t = dmin1(tact,tmax)
        if (t < 5.D0 ) t = .5*t
        xdiff = (1. - t)*xx + t*xtemp         
        tact = t
        ifunc = ifunc + 1 
        goto 999
        endif

        endif !end ifail = 7

      else

        if ( fnew > .999D0*ff ) then

          if (ifunc .lt. nfunc) then
          tact = tact*.5
          xdiff = xx -  2.*(.01*ff)*gold*tact/(gnorm0*gnorm0)
          ifunc = ifunc + 1
          goto 999
          endif 
     
        endif

      endif  ! end of tflag

         dd = xdiff - xx
         ff = fnew
         xx =xdiff
         ifunc = 0
         tact = 1.
      endif  ! end of sflag 
c-----------------------------------------
c end of routine
c-----------------------------------------
  999 continue

      return

      end

hessupd.f/      1361899259  10100 26113 100644  3641      `


      subroutine hessupd( nn, mupd, dd, jmin, jmax, xdiff, lphprint )

c     ==================================================================
c     SUBROUTINE hessupd
c     ==================================================================
c
c     o controls update of descent vector using available
c       approximation of Hessian Matrix based on gradients of
c       previous iterations
c
c     o Reference: J.C. Gilbert & C. Lemarechal
c                  Some numerical experiments with variable-storage
c                  quasi-Newton algorithms
c                  Mathematical Programming 45 (1989), pp. 407-435
c
c     o started: ??? not reproducible
c
c     o changed: Patrick Heimbach, MIT/EAPS
c                24-Feb-2000: 
c                   - changed some variable names to be consistent
c                     with routine lsoptv, lsline;
c
c     o Version: 2.1.0, 02-Mar-2000: Patrick Heimbach, MIT/EAPS
c
c     ==================================================================
c     SUBROUTINE hessupd
c     ==================================================================

      implicit none


      include 'mpif.h'
c------------------------------------
c declare arguments
c------------------------------------
      integer nn, mupd, jmin, jmax
      double precision dd(nn), alpha(100), xdiff(nn)
      logical lphprint
      integer status(MPI_STATUS_SIZE),ierr
      integer myid, nprocs,mystart,myend
      common /mpi_parm/ myid,nprocs,mystart,myend

c------------------------------------
c declare local variables
c------------------------------------
      external DDOT
      double precision     DDOT

      integer jfin, i, j, jp
      double precision    r
      double precision temp

c------------------------------------
c initialization
c------------------------------------
      jfin = jmax

      if (lphprint .and. myid .eq. 0) 
     &     print *, 'pathei-lsopt: in hessupd; ', 
     &     'jmin, jmax, mupd:', jmin, jmax, mupd


      if (jfin.lt.jmin) jfin = jmax+mupd

c------------------------------------
c compute right hand side
c------------------------------------
      do j = jfin,jmin,-1

         if (lphprint .and. myid .eq. 0) 
     &        print *, 'pathei-lsopt: in hessupd; loop ',
     &        'j,jfin,jmin = ', j,jfin,jmin

         jp = j
         if (jp.gt.mupd) jp = jp-mupd
         call dostore( nn, xdiff, .false., 2*jp+3 )

         temp = DDOT( nn, dd, 1, xdiff, 1 )
         call MPI_ALLREDUCE(temp,r,1,MPI_DOUBLE_PRECISION,MPI_SUM,
     & MPI_COMM_WORLD,ierr) 
!         r = DDOT( nn, dd, 1, xdiff,1 )

         call dostore( nn, xdiff, .false., 2*jp+2 )
         alpha(jp) = r
         do i = 1, nn
            dd(i) = dd(i) - r*xdiff(i)
         end do
      end do

c------------------------------------
c multiply precondition matrix
c------------------------------------
      if (mupd .ne. 0) then
         call dostore( nn, xdiff, .false., 3 )
         do i = 1, nn
            dd(i) = dd(i)*xdiff(i)
         end do
      end if

c------------------------------------
c compute left hand side
c------------------------------------
      do j = jmin,jfin
         jp = j
         if (jp .gt. mupd) jp = jp-mupd
         call dostore( nn, xdiff, .false., 2*jp+2 )

         temp = DDOT( nn, dd, 1, xdiff, 1 )
         call MPI_ALLREDUCE(temp,r,1,MPI_DOUBLE_PRECISION,MPI_SUM,
     & MPI_COMM_WORLD,ierr)
         r = alpha(jp) - r
!        r = alpha(jp) - DDOT( nn, dd,1 , xdiff, 1 )

         call dostore( nn, xdiff, .false., 2*jp+3 )
         do i = 1, nn
            dd(i) = dd(i) + r*xdiff(i)
         end do
      end do

      return

      end

cubic.f/        1361899270  10100 26113 100644  2199      `


      subroutine cubic( t, f, fp, ta, fa, fpa, tlower, tupper )

c-----------------------------------------
c arguments
c-----------------------------------------
      double precision    t, f, fp, ta, fa, fpa, tlower, tupper

c-----------------------------------------
c local variables
c-----------------------------------------
      double precision sign, den, anum
      double precision z1, b, discri

c-----------------------------------------
c Using f and fp at t and ta,
c computes new t by cubic formula
c safeguarded inside [tlower,tupper].
c-----------------------------------------
      z1 = dble(fp) + dble(fpa) - 3.d0*dble(fa-f)/dble(ta-t)
      b  = z1 + dble(fp)

c-----------------------------------------
c first compute the discriminant
c (without overflow)
c-----------------------------------------
      if (abs(z1).le.1.) then
         discri = z1*z1-dble(fp)*dble(fpa)
         if (discri .lt. 0.d0) then
            if (fp.lt.0.) t = tupper
            if (fp.ge.0.) t = tlower
            return
         else
            discri = dsqrt(discri)
         end if
      else
         discri = dble(fp)/z1
         discri = discri*dble(fpa)
         discri = z1-discri
         if (z1.ge.0.d0 .and. discri.ge.0.d0) then
            discri = dsqrt(z1)*dsqrt(discri)
         else if (z1.le.0.d0 .and. discri.le.0.d0) then
            discri = dsqrt(-z1)*dsqrt(-discri)
         else
            if (fp.lt.0.) t = tupper
            if (fp.ge.0.) t = tlower
            return
         end if
      end if

c-----------------------------------------
c discriminant nonnegative,
c compute solution (without overflow)
c-----------------------------------------
      if (t-ta .lt. 0.0) then
         discri = -discri
      end if

      sign = (t-ta)/abs(t-ta)
      if (sngl(b)*sign .gt. 0.0) then
         t    = t + fp*(ta-t)/sngl(b+discri)
      else
         den  = sngl(z1+b+dble(fpa))
         anum = sngl(b-discri)
         if (abs((t-ta)*anum).lt.(tupper-tlower)*abs(den)) then
            t = t + anum*(ta-t)/den
         else
            t = tupper
         end if
      end if

      t = max( t, tlower )
      t = min( t, tupper )
      
      return
      end

dgscale.f/      1361899261  10100 26113 100644  2197      `


      subroutine dgscale( nn, gold, xdiff, diag, rmin )


c     ==================================================================
c     SUBROUTINE dgscale
c     ==================================================================
c
c     o computes new preconditioner and writes it to OPWARMD
c
c     o started: ??? not reproducible
c
c     o Version: 2.1.0, 02-Mar-2000: Patrick Heimbach, MIT/EAPS
c
c     ==================================================================
c     SUBROUTINE dgscale
c     ==================================================================

      implicit none


      include 'mpif.h'

      integer nn
      double precision    gold(nn), xdiff(nn), diag(nn)

      integer i
      double precision    r1, rmin, den

      double precision temp
      integer status(MPI_STATUS_SIZE),ierr
      integer myid, nprocs,mystart,myend
      common /mpi_parm/ myid,nprocs,mystart,myend

c-----------------------------------------
c read diagonal
c-----------------------------------------
      call dostore( nn, diag, .false., 3 )

      r1 = 0.
      do i = 1, nn
        r1 = r1 + gold(i)*gold(i)*diag(i)
      end do

      temp = r1
      call  MPI_ALLREDUCE(temp,r1,1,MPI_DOUBLE_PRECISION,MPI_SUM,
     & MPI_COMM_WORLD,ierr)
      r1 = 1.0 / r1

      call DSCAL( nn, r1, diag, 1 )

c-----------------------------------------
c update the diagonal
c (gg is used as an auxiliary vector)
c-----------------------------------------

      den = 0.0

      do i = 1, nn
cph(
         if (diag(i).LE.0) then
cph            print *, 'pathei-lsopt: in dgscale; diag = 0 for i=', i
            diag(i) = rmin
         end if
cph)
         den = den + xdiff(i)*xdiff(i) / diag(i)
      end do

      temp=den
      call  MPI_ALLREDUCE(temp,den,1,MPI_DOUBLE_PRECISION,MPI_SUM,
     & MPI_COMM_WORLD,ierr)

      do i = 1, nn
         diag(i) = 1./
     $        (1./diag(i)+gold(i)**2-(xdiff(i)/diag(i))**2/den)
         if (diag(i).le.0.) then
            diag(i) = rmin
         endif
      end do

c-----------------------------------------
c write diagonal
c-----------------------------------------
      call dostore( nn, diag, .true., 3 )

      return
      end

instore.f/      1361899270  10100 26113 100644  1780      `


      subroutine instore( n, fc, gnorm0, is, m, jmin, jmax,
     &                    cold, sflag, tflag, ifunc, tact)
c
c         arguments
c
      integer n, is, m, jmin, jmax, ifail, ifunc
      double precision    fc, gnorm0,tact
      logical cold
      logical sflag, tflag
      integer*8 rectmp, tmp1, tmp2
c----
c
      common /xxstore/ itape, ntape, isize
      integer itape, ntape, isize
      common /mpi_parm/ myid,nprocs,mystart,myend
      integer myid, nprocs,mystart,myend

      ifail = 0
      itape = 91
      ntape = 92
      isize = is
      cold  = .true.

      if (myid .eq. 0) print *, 'pathei: in instore '

cc one record = two arrays of real*(isize)

      open(  itape
     $    , file   = 'OPWARMI'
     $    , status = 'unknown'
     $    , form   = 'formatted'
     $    , access = 'sequential'
     $    )

      read( itape, *, end=800, err=900 ) n, fc, gnorm0,
     $                                   m, jmin, jmax,
     &                                   sflag, tflag, ifunc, tact 
!Note: bnc read itape defaults to cold if OPWARMI has 0 size and 
!      skips the next line, due to end condition

      cold  = .false.

 800  continue
      close(itape)

      if( myid .eq. 0) then
      print*
      write(*,'(a,i10,a,i2)')
     $        '  opening direct access with recl = ',n,' * ',isize
      print*
      endif

      tmp1 = n
      tmp2 = isize
      rectmp = tmp1*tmp2

cph(
      if(myid .eq. 0)  print *, 'pathei rectmp ', rectmp, n, isize
cph)

!      open( ntape
!     $    , file   = 'OPWARMD'
!     $    , status = 'unknown'
!     $    , form   = 'unformatted'
!     $    , access = 'direct'
!     $    , recl   = rectmp
!     $    )

      return

 900  continue
      close(itape)
      ifail = 1

      return
      end
dostore.f/      1361899270  10100 26113 100644  1137      `


      subroutine dostore( n, x, store, j )
c
c         arguments
c
      integer n, i, j
      double precision    x(n)
      real*4 tmpx(n) 
      logical store

      common /xxstore/ itape, ntape, isize
      integer itape, ntape, isize
      common /mpi_parm/myid,nprocs,mystart,myend
      integer myid,nprocs,mystart,myend
 
cph(
cph      print *, 'pathei in dostore ', 
cph     &   store, n, ntape, j
cph)

   
      if (store) then
        do i = 1, n
          tmpx(i) = x(i)
        enddo
!bnc may need to fix this
      open( ntape
     $    , file   = 'OPWARMD'
     $    , status = 'unknown'
     $    , form   = 'unformatted'
     $    , access = 'direct'
     $    , recl   = isize*n
     $    )
      write( ntape, rec=(j-1)*nprocs+myid+1) tmpx
      close(ntape) 
      else
      open( ntape
     $    , file   = 'OPWARMD'
     $    , status = 'unknown'
     $    , form   = 'unformatted'
     $    , access = 'direct'
     $    , recl   = isize*n
     $    )
      read(  ntape, rec=(j-1)*nprocs+myid+1) tmpx
      close(ntape)
        do i = 1, n
          x(i) = tmpx(i)
        enddo
      endif

      return
      end

outstore.f/     1361899270  10100 26113 100644  636       `


      subroutine outstore( n, fc, gnorm0, m, jmin, jmax, sflag, tflag, ifunc, tact )
c
c         arguments
c
      integer n, m, jmin, jmax, ifunc
      double precision    fc, gnorm0, tact
      logical sflag, tflag
c----
c
      common /xxstore/ itape, ntape, isize
      integer itape, ntape, isize

      if (myid .eq. 0 ) then 
      open(  itape
     $     , file   = 'OPWARMI'
     $     , access = 'sequential'
     $     , form   = 'formatted'  )
      rewind (itape)

      write( itape, * ) n, fc, gnorm0, m, jmin, jmax, sflag, tflag, ifunc, tact

      close(itape)
      endif

      close(ntape)

      return
      end
lswri.f/        1361899270  10100 26113 100644  1470      `


      subroutine lswri( isize, phniter, NN, xx, gg, lphprint )

c     ==================================================================
c     SUBROUTINE lswri
c     ==================================================================
c
c     o writes formatted output of updated controls and 
c       adjoint variables
c
c     o started: Patrick Heimbach, MIT/EAPS
c                01-Apr-1999
c
c     ==================================================================
c     SUBROUTINE lswri
c     ==================================================================

      implicit none

C     ========= Routine arguments ======================================

      INTEGER   isize, phniter, NN
      real*4 xx(NN), gg(NN)
cph      DOUBLE PRECISION      xx(NN), gg(NN)
      LOGICAL   lphprint

cph      CHARACTER yitnum*3
c

C     ========= Local variables ========================================


C=======================================================================

c      if (lphprint) 
c     &   print *, 'pathei: formatted output of xx, gg, in lswri: ',
c     &   phniter
c
c      open(95,
c     &     FILE='PH_X.dat',
c     &     FORM='UNFORMATTED',
c     &     ACCESS='DIRECT',
c     &     RECL=NN*isize)
c      write(95,REC=phniter+1) xx
c      close(95)
c
c      open(95,
c     &     FILE='PH_G.dat',
c     &     FORM='UNFORMATTED',
c     &     ACCESS='DIRECT',
c     &     RECL=NN*isize)
c      write(95,REC=phniter+1) gg
c      close(95)

      end
distwork.o/     1361899272  10100 26113 100644  1552      `
ELF>P@@ЯЪ"'7GQX  g
-pUH┴ЕD▀D┴п≥▀6ВЧ┘рt,D▀B█+б≥ВЧA┴A┴аE╞кЪнD;чA█y┴9B█u.E┴К,D┴п≥ВЧA┴A┴аD▀ЪнE╞кD;чA█y┴9B█uE┴КA┴H┴Л]цx░,p├kф
_p
distwork.f distwork_цйцбкЪ├деЧ┐├х┘ацЧбеЧ┐├бЗдdistwork.f.strtab.symtab.data.bss.textdistwork_.note.GNU-stack.eh_frame.trace.rela.eh_frame.rela.trace ,00000@8x"x'xp7ХGХHQ0`░sXg0distwork.f/     1361899270  10100 26113 100644  1000      `


      subroutine distwork(myid,numprocs,nmax,mystart,myend,my_veclen)
!     There are numprocs of MPI processes calling this routine
!     Distribute different sections to the processes
      implicit none

!     Input variables
      integer :: myid,nmax,numprocs
!     Output variables
      integer :: mystart,myend
!     Local variables
      integer :: rest, my_veclen

      rest = mod(nmax,numprocs)
!      if(rest.gt.numprocs/2) then
      if(rest .ne. 0) then 
        ! If each process takes my_veclen, total will be bigger than nmax
        my_veclen = (nmax+(numprocs-rest))/numprocs
        mystart = 1 + myid*my_veclen
        myend = mystart + my_veclen - 1
        if(myid.eq.numprocs-1) myend = nmax
      else
        ! If each process takes my_veclen, total will be less than nmax
        my_veclen = nmax/numprocs
        mystart = 1 + myid*my_veclen
        myend = mystart + my_veclen - 1
        if(myid.eq.numprocs-1) myend = nmax
      endif
      end subroutine distwork

pddot.o/        1269896888  9435  26131 100644  8016      `
ELF2░@@pddot.f.strtab.symtab.textpddot_pddot_??unwSTRLITPACK_0$0STRLITPACK_5$0$0for_write_seq_lisSTRLITPACK_6$0$0for_write_seq_lis_xmitpddot_$STATUSES$0$0pddot_$REQUESTS$0$0STRLITPACK_7$0$0STRLITPACK_9$0$0LITPACK_4$0$0LITPACK_3$0$0LITPACK_2$0$0STRLITPACK_10$0$0LITPACK_0$0$0LITPACK_1$0$0STRLITPACK_8$0$0STRLITPACK_11$0$0LITPACK_5$0$0STRLITPACK_12$0$0for_emit_diagnosticddot_mpi_waitall_mpi_reduce_STRLITPACK_13$0$0STRLITPACK_14$0$0STRLITPACK_15$0$0STRLITPACK_16$0$0LITPACK_7$0$0LITPACK_6$0$0LITPACK_8$0$0LITPACK_10$0$0LITPACK_13$0$0LITPACK_12$0$0LITPACK_9$0$0LITPACK_11$0$0.IA_64.unwind_info__udt_pddot_??unw.IA_64.unwind.bss.rodata.sbss.sdata.rodata.str1.4mpi_sgi_private_status_mpi_sgi_private_mpi_sgi_private_inplace_mpi_sgi_private_char_mpi_win_dup_fn_mpi_win_null_delete_fn_mpi_win_null_copy_fn_mpi_comm_dup_fn_mpi_comm_null_delete_fn_mpi_comm_null_copy_fn_mpi_dup_fn_mpi_null_delete_fn_mpi_null_copy_fn_.trace.rela.text.rela.IA_64.unwind.rela.traceЯЪ  &`s┘⌠▄0≤Єbеpe══і   A d .  ? _     $Р (Д ,ж 0Q 4 8є <І @х Dз HЗ LЛ PC T XQ \4 `% dґ
2
       елв
Й       Ruq┘▀≤╪РЪтРЪЕРЪЧРЪ$<Rc|⌠÷Ё@иT─═ХЧK`░`─>#пAH XNЮ≤   д  )И╟ЭСЪO a└H@░ ═` b≤╟B`b└`,Апpб1юHX `2! ю3B─g└XАl╟3@@pA!═e`░p8!─H@     p└     ░a!`C0B═     ь└
═4t
░X8ю(А─юX@PА═Pr% ≤╦0# L≤Хp0#
p@8
`H▓@Б$A
@Рa@б=A
@рp<A
@bh0A
╟04H.АPP\!ю0B X└
X$H2!#PV!p$cX$@х#HZp$═H≤пA(0 BП HЯOЭ▄Х
$═(`ю░!0,H`@X░!Ю7`h$ЮH0╙   п░!ПH─@└$PH░ ╦$╟H─░х$ pN═■  ЯGЧ▄ h А ≤  `8А╟б1
═0tП(АпЯ(@HtА
─,tЮ А╟  pBpl,▌
14`W─P┬!P@`├!пТB@ Ь└     pp!ПэB`     П└─▌!ЮBP▐0╛BPB`├!пB@ (└     p┌!ПB`      └─▌!BP@▌1╛BPB`├!пТB@ Ь└     px!ПЛB`     П└─▌!ПBP▐0╛бPB`├!пB@ (└     pt!ПДB`      └─▌!BPV!PL!ю┬BюsI─     pH!─ B─sI─XБ=?#пдs~FP l≤юлB═     ░└     Pj!╛B`     ═└PxZ!Ю╦B 
░└     h^!ююB
`└V!╟дB@ ╟└P8b╛B 8H══ЩЪHN═Hб1     ПLь03 ХJю─   HА─Пб18tА
`$t
Пt
Юtь<А═qб!░l,hdАпX``А
╟h8
═`8
░,I
─(I
p@H
`HH
П04
═<А`ДPa─═XB░$П░b ░─$Пх`BX"@9` h!Ю@B─Y░┬ !ПHB
 └PV!H$4▓`░$ H`BX@)`     h!ЮB─Y░┬▄!П$B
 └PV!ь$п░b@░х$─х`BX4@9`     h0!ЮdB─Y░┬2!ПlB
 └PV!┬$4▓`░x$ЮH`BX @)`     h!Ю<B─Y░┬▄!ПDB
 └PV!
Ьp.─p≤B@@I─ ь$░H`@
I─пА?#─ц~FPЬхT─@ь03 X└
Ь|.─pH@░└ю$░лBPюb!╛Bю░└╦l! п$пH─░х$╟HPЬхT─╟1B░
pb╛B═
╙@U─pH 6ью1B─└`пЮ
Ф╠(Б╧╦Д╟╘Й╡╙aФ│ю      fort: (%d): Subscript #%d of the array %s has value %d which is less than the lower bound of %d
fort: (%d): Subscript #%d of the array %s has value %d which is greater than the upper bound of %d
8   XXYYXXYYnwork = 
Й     pddot.fpddot_юю─юЛааЪюМаюЧюОаюЩаОб┌юЪюМюююЛаююЧююЩаЬйЪ├ьЪ┼аа     цЩююЩюЗкааЭюаЭю─аЭюбЭьЭШюдЗбюЩц┘юЪ└─бЪ──ца─абЗцбЩбаЩюбЩцдЩббЩббЩхю┌юЪ_е─юЪ┼аюЩюЧдаЩб2*▓*╒I0ю*рI1П2
Я*22
`*a*q*─*│*┌*░*▒*▓*═*┌I2бI2I2BI2┌I3╡I4РI5p*r2бI2Ю*Б2
2I2P* R2╒I2ю*!б2
I2@*#A*"RI3q*─*┌I4╟*%╠*(╡*'ю*$а*&бI5__     _'/0'/00        @Ж6H/~
@ @`Ў@┘p┌Ч⌠0≤Ю═ЖіЖpґ2fеvПлf╟вHЙ^0pddot.f/        1269896881  9435  26131 100644  1677      `


      FUNCTION PDDOT(nmax,X,xn,Y,yn,myid,my_veclen,numprocs)
      implicit none
      include 'mpif.h'

      integer :: xn, yn
      integer :: i,nmax,my_veclen,nwork,mystart,myend
      integer :: ierr, myid, numprocs
      integer :: requests(2)
      integer :: statuses(MPI_STATUS_SIZE,2)

      real(kind=8) :: answer, total, totend, DD
      real(kind=8) :: x(nmax), y(nmax)
      DOUBLE PRECISION PDDOT

      DOUBLE PRECISION DDOT
      EXTERNAL DDOT
      total = 0.0
      totend = 0.0

      nwork = nmax/(my_veclen*numprocs)
      write(*,*) "nwork = ", nwork
      
      do i=1,nwork
        !
        mystart = 1 + (i-1)*(my_veclen*numprocs) + myid*my_veclen
        myend = mystart + my_veclen - 1
!        write(*,*)"myid mystart myend = ", myid, mystart, myend

        answer = DDOT( my_veclen, X(mystart:myend), xn, Y(mystart:myend), yn )
!        write(*,*) "myid answer = ", myid, answer

        call MPI_WAITALL(2,requests,statuses,ierr)
        call MPI_REDUCE(answer,totend,1,MPI_DOUBLE_PRECISION,MPI_SUM,0,MPI_COMM_WORLD,ierr)
        total = total + totend
      enddo


      DD = (my_veclen*numprocs)*nwork
!     Finish the rest
      answer = 0.0
      mystart = 1 + DD + myid*my_veclen
      myend = mystart + my_veclen - 1
      if(myend .gt. nmax) myend=nmax
      if (mystart.le.nmax) then
        answer = DDOT( my_veclen, X(mystart:myend), 1, Y(mystart:myend), 1 )
!        write(*,*) "myid answer = ", myid, answer
      endif
      call MPI_WAITALL(2,requests,statuses,ierr)
      call MPI_REDUCE(answer,totend,1,MPI_DOUBLE_PRECISION,MPI_SUM,0,MPI_COMM_WORLD,ierr)

      PDDOT = total + totend

      return
      end function PDDOT

lsupdxx_1.o/    1270584817  9435  26131 100644  9952      `
ELF2`#@@lsupdxx_1.f.strtab.symtab.textlsupdxx_1_lsupdxx_1_??unwSTRLITPACK_3$0STRLITPACK_4$0$0for_write_seq_lisSTRLITPACK_5$0$0for_write_seq_lis_xmitSTRLITPACK_8$0$0STRLITPACK_6$0$0STRLITPACK_7$0$0STRLITPACK_9$0$0STRLITPACK_10$0$0STRLITPACK_11$0$0for_emit_diagnosticSTRLITPACK_2$0STRLITPACK_12$0$0STRLITPACK_1$0STRLITPACK_13$0$0STRLITPACK_14$0$0STRLITPACK_15$0$0STRLITPACK_16$0$0STRLITPACK_17$0$0LITPACK_1$0$0LITPACK_0$0$0ddot_STRLITPACK_20$0$0STRLITPACK_19$0$0STRLITPACK_21$0$0STRLITPACK_22$0$0STRLITPACK_23$0$0STRLITPACK_24$0$0STRLITPACK_26$0$0STRLITPACK_25$0$0STRLITPACK_0$0STRLITPACK_18$0$0hessupd_.IA_64.unwind_info__udt_lsupdxx_1_??unw.IA_64.unwind.rodata.sdata.rodata.str1.4.rodata.str1.32.trace.rela.text.rela.IA_64.unwind.rela.traceЯЪ
.z█ё╠╘b≤pe╧Mp_#D Д(Ж0V8h<z@▄D╛H·L╨PкTэXН\р`юdhl>p,tю5
о    (> 0#P `*ъ
ФЯ
#`^│╨qHаP─ю0|FHю─  л       X`!п3B`d└\6ьп,0 @д  :╟1B─)4Б┬QN`l0 `й P`!Ю─2Bd└ьP80 Ю│c└ ░@!0│1B┌d└P@0 ─┌b└ 8!0─B └8(Ю°B─└
ПH0 &─ьBПD@     ПP(~≤ХBх !@H░     ┬A!╟H─     `└░0░@ ═`     P≤≤д0#     ┬└H2≤░ЭСЪOPH\═└1B `└@b!░H$■≤PP~╟BЮ├░░  H@
ХJЛq─хBP~@LpHю≤Б 0Z@BЮ<ф└!юH0B(X
 H
─1▄ 
FH
p  ┴°─
\Da
─С└.A─Ц─@A
─сx<@A
─ц╦\@A─Ё8A
─ёh4A
─⌠`0A
─┐P(I
─sH$A
─c┴i─sAх░CX ▌!═@BЮЮX @ ▌!0@Bю░П ▌!юA        B@Z  b░Ь?~F@░bЮ80 ═HX┬$Ю%`b@0 п    ╙h    $ю6`bH@6░   8$`H`V░<b░┐BЮ&└
($pд2 Йx0Пbи└ю0!@Brb Пb╗Ц
ю
хn!pЮB░@rb░эBЮю└ @rb░пB╗Цю Pf!╟хB     ╟└px!ПдB─ ┬└     хn!pЮB═     ─└Йx0 ПbPX!PBPf!╟ЄB x└     `b!п°B      p└px!ПХBPX!PBPf!╟хB ╟└     `b!п≤B      ═└px!ПдBPX!─@BPf!╟ЄB x└     `b!п■B      p└px!ПХBP@ HюПB `└
` .h@F≤═`ЧЪHЮD╟A1B@░ х─!═┐0B─r░≤$─p(@РГЪ÷░B┬4░@     ═`     P2≤Х0#     п└D6≤░HB` ≤└`2!PPB!Пю0B `└@t!░HH`░$■≤░H
@@<0#PX!П@╦▌!pcю└Ь$─═Cх ▌!─AB`ЮX ≤ ▌! AB@░─ ▌!Ю@        BюЮX  2b@Ыo~F ░0bph0 Ю═Xx$═%`&b D0 П╙0       $`6`$bЮH6░   п$░H V░ bп┐B─&└
ю$8д2 З~
0Эb─И└0║@Bzb ЭbХШ
ю
Хh!@ПB░@zbппB─Ю└ @zb░хBХШю Pb!╟юB     ≤└p~!П╪B─ x└     Хh!@ПB═     p└З~0 ЭbPX!PBPb!╟≤B h└     `^!пХB      8└p~!ПЬBPX!PBPb!╟юB ≤└     `^!пДB      ░└p~!П╪BPX!─@BPb!╟≤B h└     `^!пЮB      8└p~!ПЬBPx nПЭB `└
Юx<"@8l≤═`ЧЪHH`─0C
DP──C@├!юH@     ░H└!╟BP ▄≤░ ■ `C░
═$п  `└
└║╒! 1lB░@Cь!HU╟
╙╦6ью0B═
ПA──└4┬1░C░ └!AB`пX      ─ └!П@B─┌└P$░@B`Z (bЮЬO~F`░$bP,0 pXh$ю$`"b(0 ─╙x$ 6` b0H@&└         $ЮHЮ░b$д2 dbpпb`8ёю4j≤─@BхP^║╠╦B       ┬└и`Z║я°B# ─└яph║ЯЄBRX!@BhP^А╟░б 0└i`ZАп▄б! (└qphАПлбR     ` D╟т00─═└
X!═0 пX(*═PЪЪHп├┬ 
 пф _═B `eЫ
5▄└┼≤─▒▄@A
─║─a─qH$A
─│(A
─a│A─А8A
─Я A
─а0A
─╠p8I
─║`0A
─▒5i
─Z$H
Ю─· Л─j≤≥#│ ─ц│є C╘─j≤у03░┴(
░┴$s░H┬v≤YМ03 *`B░HЛ03─═ЭЪJП └!пAB`пX ь └!░H@░≤ └!═AB`Z <bюЫ~F┐└:b░0 ─ЮX╟$@'`6bd0 @    ╙╦    $@6`0b`H6░   а@     !pHW░&b0┌     B═░ p$@B═░4b ┌BЮ░(─!@H h1°я`в─(!p@Ь░8─L!ьB`╟└ Ь`!ЮЮB@═└Хb!юхBP
ЬюP─°я` `└─pBЬ░8─N!ьB`≤└ Ьn!ЮЮB@═└Хr!юХBP
ЬюP─░пb `└─pBЬ░8─Z!ьB`╟└ Ь`!ЮЮB@═└Хb!юхBP
ЬюP─°я` `└pBЬ░8─\!ьB`≤└ Ьn!ЮЮB@═└Хr!юХBP
ЬюP─юпb `└─BЬ░8─^!ьB`╟└ Ь`!ЮЮB@═└Хb!юхBP
ЬюP─п°я` `└─ BHv`@Ь░8─j!ьB`≤└ Ьn!ЮЮB@═└Хr!юХBPЬюP─░Л00 `└@ Jp@░08 │└ Fb0┬B@└!@пB
`@$H@@≤═`ЩЪHЫЪH`б!═ю3B ░ @є$╟H     a└П░@ ═`     ░≤@01# ЫГЪ÷x≤PX!@ЬЪHP└!╟▄B     └     `H!п─B      (└pD!PX!ЬЪH`*ЮФ╠)Б╧╦Х╠╗Д╟╙Й╡╚a╟│юaХ║ю`fort: (%d): Subscript #%d of the array %s has value %d which is less than the lower bound of %d
fort: (%d): Subscript #%d of the array %s has value %d which is greater than the upper bound of %d
80888XDIFFXDIFFDDDDGGGGDDDDGGGGDDDDXXXXDDDDfold - ffpathei-lsopt: first estimate of dd via pathei-lsopt: using fold = fmin = pathei-lsopt: no hessupd for first optim.
j`lsupdxx_1.f
lsupdxx_1_бю;иею;юобю;гею;баЪ─≈юаЧююЩюЪеЫцююЧю┘бЭббЧГ┌─аЪ─аЪ─юЪ│юЪ└юЪ└юЪ└бЪ─дЪ│аЪ│юЪ│дЪ─фЪ┤аЪ╗юЪ┘юЪю┘юЪ┌юЪ²г
─аЪ─аЪ─юЪ│юЪ└юЪ└юЪ└бЪ─дЪ│аЪ│юЪ│дЪ─фЪ┤аЪ╗юЪ┘юЪееф      ──ю     юЪюЖбаВюЪх1ды─бЪ─юЪ─юЪ│юЪ└юЪ└юЪ─юЪ┌цЪ⌡юЪ└юЪбюХюаеПаЬцЦЬа   юВб
─юЧюВбцююЧ└─цЧб─юЪ│аЪ─юЪ─аЪ└юЪ└юЪ│юЪ│аЪ│юЪ─аЪ─бЪsгЪ┐цЪяркЭ2(*
BI/a*rI0"2B2
▒*═*║*ю*▓I1рI1I1RI1▓2'╟*I/"2%1*RI0╒2б2
* *!*@*I1RI1▓I1рI1Q    *R *b I2
2
"
2q
*┌
*Б
I1"I11
2
2
2║
*╠
*б
* п
*!р
*#Б
*"RI1╡I1I1rI1рI1BI1б2)я*I/BI3__`_'.0'.
@PЮ.0─@z╟@ёp┌П╠Ю╧Х─ю2hо2x░ъpФxpЯХ"H0#0
lsupdxx_1.f/    1270584807  9435  26131 100644  4007      `


      subroutine lsupdxx_1(
     &     nn, ifail, lphprint
     &     , jmin, jmax, nupdate
     &     , ff, fmin, fold, gnorm0, dotdg
     &     , gg, dd, xx, xdiff
     &     , tmin, tmax, tact, epsx
     &     )

c     ==================================================================
c     SUBROUTINE lsupdxx
c     ==================================================================
c
c     o conceived for variable online/offline version
c       computes - new descent direction dd based on latest 
c                  available gradient
c                - new tact based on new dd
c                - new control vector xx needed for offline run
c
c     o started: Patrick Heimbach, MIT/EAPS
c                29-Feb-2000: 
c
c     o Version 2.1.0, 02-Mar-2000: Patrick Heimbach, MIT/EAPS
c
c     ==================================================================
c     SUBROUTINE lsupdxx
c     ==================================================================
c


      implicit none

c-----------------------------------------
c declare arguments
c-----------------------------------------
      integer nn, jmin, jmax, nupdate, ifail
      double precision    ff, fmin, fold, gnorm0, dotdg
      double precision    gg(nn), dd(nn), xx(nn), xdiff(nn)
      double precision    tmin, tmax, tact, epsx
      logical lphprint

c-----------------------------------------
C declare local variables
c-----------------------------------------
      integer i
      double precision    fdiff, preco

      double precision     DDOT
      external DDOT

c     ==================================================================

c-----------------------------------------
c use Fletchers scaling
c and initialize diagional to 1.
c-----------------------------------------
c
      if ( ( jmax .eq. 0 ) .or. (nupdate .eq. 0 ) ) then

         if (jmax .eq. 0) then
            fold = fmin
            if (lphprint) 
     &           print *, 'pathei-lsopt: using fold = fmin = ', fmin
         end if
         fdiff = fold - ff
         if (jmax .eq. 0) fdiff = ABS(fdiff)
         
         preco = 2. * fdiff / (gnorm0*gnorm0)
         do i = 1, nn
            dd(i)    = -gg(i)*preco
         end do

         if (lphprint) 
     &        print *, 'pathei-lsopt: first estimate of dd via ',
     &        'fold - ff'

c-----------------------------------------
c use the matrix stored in [diag]
c and the (y,s) pairs
c-----------------------------------------

         else

            do i = 1, nn
               dd(i) = -gg(i)
            end do

            if (jmax .gt. 0) then
               call hessupd( nn, nupdate, dd, jmin, jmax, xdiff, 
     &              lphprint )
            else
               if (lphprint) 
     &              print *, 'pathei-lsopt: no hessupd for first optim.'
            end if

         endif

c-----------------------------------------
c check whether new direction is a descent one
c-----------------------------------------
         dotdg = DDOT( nn, dd, 1, gg, 1 )
         if (dotdg .ge. 0.0) then
            ifail = 4
            goto 999
         end if

c----------------------------------
c declare arguments
c----------------------------------

      tmin = 0.
      do i = 1, nn
         tmin = max( tmin, abs(dd(i)) )
      end do
      tmin = epsx/tmin

c----------------------------------
c make sure that t is between
c tmin and tmax
c----------------------------------

      tact  = 1.0
      tmax = 1.0e+10
      if (tact.le.tmin) then
         tact = tmin
         if (tact.gt.tmax) then
            tmin = tmax
         endif
      endif

      if (tact .gt. tmax) then
          tact = tmax
          ifail = 7
      endif

c----------------------------------
c compute new x
c----------------------------------
      do i = 1, nn
         xdiff(i) = xx(i) + tact*dd(i)
      end do

c----------------------------------
c save new x to file for offline version
c----------------------------------

 999  continue

      return

      end