a
    \:b*                    @   s  d dl Zd dlmZ e Ze e_dej_dej_dej_e	ddgej_
dej_d	ej_d
ej_e	ddgej_e	g dej_dej_dej_dej_dej_dej_dej_dej_e	g djddddej_e	g de_e	g de_eje_de_de_de_e	g de_e	g de_ de_!e	g de_"ej"ej e_#de_$e	g de_%e	g de_&e	g de_'e	g djd d dde_(e	d!d"ge_d#e_)e	g d$jd d dde_*e	g d%jd d dde_+e	g d&jd d dde_,d'e_-d(e_.d e_d)e_/d*e_0e	g d+e_1e	g d,e_2e	g d-e_3d.e_4d/e_5d0e_6e	g d1e_7e	g d2e_8e	g d3e_9e	g d4e_:d.e_;de_e	d5d6ge_<e	d7d8ge_=de_>d9e_?d:e_@d;e_Ae	d<d=ge_Bde_Ce	d>d?ge_De	d@dAge_Ee	ddge_Fe	g dBjd d dde_GejGej> e_He	g dCjd#d#dde_Ie	g dDjdddde_Je	g dEjdd#dde_Kd*e_0d'e_Ld e_d e_e	g de_MdS )F    N)BunchZweight            ZNAz	s(weight))  r   r   r   ga8ㄞ@gP8@i  gqǟ@gBq@g	@ga8㐧@gX/8@  r	   r	   r	   
      gZ>      )ygDs-}?')_+Ϭ?gM, ?Yj?Mƪ?aGBs?r?xn?甇llM{?Ԏr2v?r   gj,Yl*?sx\uJ?g/r?
dEcB?!\Ͼp?1WF)_pa?i6sN?hl%?r   r   g(w?O$~sǿgЙzSVg?}k?	e?g^[7}g}lE?g#9?TԐ6Cg[, ?r   r   g7z6?DTҿ=ܬ?Nӑܡ?W!?Yˍv?>Lv?p=eo?r   g0r?gЙzSVg?r#   g<5?sK4ǿgK=]pg,p mw`gofyFD)cYWr   r   r    r$   r*   g>9N?2k5?V/?INm?@ɷ`3c?r   r   r!   r%   gK=]pr-   gY~?OE_;%,߂?o?f/4X?r   r   gS[7}r&   gCp mr.   r2   gޞN,!i?[d4CQYN^g?r   r   glE?r'   r+   r/   r3   r6   gx_A6o?
Hy磷Yd?r   r   g#9?r(   gofyFr0   r4   r7   r9   g{B@?݆r   r   r"   r)   r,   r1   r5   r8   r:   r;   gu+^C?F)order)(cA=@U_/	?S4.?g>0D4@ggC)gjQI+g;-gFgg
gd[g|5j""gsX!g%gmg^IO_@g6+*?g<;"g 	gg`Y1ug'Q;Ug	1 dg{û*)竵z^5@rB   gڥm2@gy7@gc3@gz{6@gCo5@g|'	5@gLuX3@gzT1@*B;7@rC   gtaÔ4@gY4@gD4@ggX0@gj.@gZka4-@g<H@VA@us:@@ a3tA@rF   皱̽9@)1^?@ڂ&p?@rI   gD]h9@'8@gֹ؊3@gk;$E@gPB@gtv&B@g`e?@g3ؠ>@gL˴n=@g0鼮=@g}(:@g߬$:.\:@g!c-:@g[q+b9@giQY#QM7@gH5&ȵ8@g\j:@rD   rE   g*"7@X"~-@rK   gҖ~,@g܎i,A@g`Q@@gT@@g1W<@@g'dg(*@@Ol07@rL   g!Ǻ&7@gmJ4@eK9@~c9@rM   rN   g=$:A@g#)U9@gg4@g @?@g'ю7@g4,bn6@gnQ7@g=k6@grM:-@gpi-@g'c6,@g5K+@gi-2@`Sը@@g߰?@@g-.?@gvӧ9@g.ʚ6@g s9@gO%F3@g֝3@g-_3@g39@g=8@\6@rP   g3qA@g*tC@gDV
@@gx6B@@g\L$>@g@ؗ@@gSχ>@g(	ئ?@g8>@gEu=>@gi8@g#Z8@gDχ2@g^[&1@gj_3@gQ;2@g:CX0@g5a1@g9ma6@gX# 7;@g?Av4@gEV:9@g, 5@uoS:@gZK4@g8I9@g1VuHX5@rQ   g}}2@rO   rG   rH   rI   gsh=@rJ   gP;
?3@gɳ?&83@緰ɀ1@rR   g,n1@gP|6/@gZc5@gA5@gQC5@g*Y8z5@gd
Y2@gda2@gM7u>@gH'1na=@gPM9;@g]$5+;@gDT9;@g -8@g\8@g"xd5@g'(:@gtF	8@gg,8@gǕ<hn4@giH<@@g?sGF?@g\>8?@g\ >@g#Jr<@g(:W9@gRt>@gMe=@籄뎞B@rS   gZ6=@gxێ=@geF==@g*;>)=@gٝV<@g`aT7@gVyg@7@gQw+5@gG'ި5@gbvD5@gJ(bX4@g
U"G4@g딪4@g&jj:)9@g3<$@@J_8@rT   g@D8@g}
2@gey2@g!1@g
v1@g+]ʸC@gKu'6;@g6;B@gg@1;@g~:@gdA@gaN!8@gl^::@g՛_:@gVH5@g	Z@@gTÐ8@g=4@g~gM4@gV}4@gVf2hw4@g־5+Q2@ga"1@g0/4@g\2@gYw3@g[:@ge*y4@gPy[9@g_@)r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r      )   rV            rW   rW   rW            r\   rV   rV      r[   r[   r   /   &   r_   %      rX   ra   ra   ra   rX   rX   rW   1   ra   r_      rc   rc   rc      rd   rd   rd   rX      rX   r_   r_   rX   r   r      rc   ra   ra   ra   ra   rZ   rZ   rZ   r[      rg   rg   rg   $   rg   rW   ra      ri   ri   ri   r[   r[      rj   rW   r`   ra   ra   rX   r\   re   rW   rW   rW   re   re   r\   r\   ra   -   ra   ra   ra   ra   ra   ra   ra   ra   rd   rd   rZ   rZ   rW   rW   rZ   rW   rW      rW   re   rW   rl   rW   re   rW   rl   rY   r`   rX   ra   ra   ra   rX   rW   rW   rZ   rZ   rZ   rZ   rV   rV   rV   rV   rW   rW   ra   rg   rg       rl   rg   rX   rX   rl   re   r\   r\   #   ra   ra   ra   rd   rd   rc   rc   "   r_   r_   rl   rl      rp   rg   rg   rX   rX   rX   rX   rX   rX   rp   rc   rd   rd   rd   r]   rW   r]   rW   r`   rd   r`   rd   rd   r`   rg   rX   rX   rW   !   re   r\   r\   rX   rX   rZ   rZ   r\   rW   rY   rg   rW   gnaUp@)gi?     ?g     ?gZ     ?g'?gwhA?gm?g,%?g;c?gE?g釐:?g3O?gHX?g{}#?g"C?gkТo?gPeru>?gpj(?gWfBf?gvn?gu?gIQ?gP?g{-˒t?)g?g    ?gh    ?rr   g
?gvl?g1Q"`?g6?gkG[?gd,?gIâ}?g^{?g@?gC?gJB-	?g(*
?g\Op?gYk?gg=/Փ?gk?g3k6?g?gx56?gҖZ?)H?rs   g{ӥ?g&?g!@.?gK!S)6?gM'5۲?gk?gry?g	7J`#?/{?rt   gjp?gNeG?g)}?gsn  ظ?g8?gq?g"!7?綆`?R1r?"sZnH?rw   AʡK?QB? pL?rz   g7ӛ)?$?g4"?g>?g/l> ?g꽸ø?gQ?gubI?g]`!á?g3Ŗ?gxLz$?ğ5?g^aNUR?gzƼ`8?gd
Z?g1'b۩?g?ru   rv   ga_?{:&?r|   g/\sF?gT2ȥ?gGf?g7ßy?g?g	K,˰?ȕ"e?r}   gsZM&?gۗb?Tϼۧ? G>Ѩ?r~   r   g?gъ5ܖ?g!?ga7{=?gn6;SԼ?g/R?gUN;?gND`?gD?g(W&;?g m
k?g7+?gdZ1o?@O6L?gCEvJ?g1?gO1?g?gc ?g@^4B?gS':3?gPO?gɓY?gϦ?犥v'?r   g8*?gX?gk˱^?g
7?g_?g.A?g3-?g9?gIxĐ?g!_7?gɓ?gW'?g&((?g'f?go='?g}K$?gagN?gx$é?gF[\?gı=Ϲ?grӿdt?g޹|?g>iU?.?gd?g7|?g:7?r   grFEd?r   rx   ry   rz   g8?r{   g1?g$5?礳?r   g{:?g"?g(2?g/f?g閰?gDpt?g?g%?g+T	
?gkH1?gԾܾ$}?gQ"j?gJ{<߉?gLcpM?gY&?g&?gl,?g`-l?g]kR?gr?g\p?g{rk?g'n?gt?gOUA1?g"x?g雵F?g5^?S7?r   g6\({?gQ z?g`?gǆR?gXۖ#?g8~yKu?gW?gtĽϴ?gjr?gTPմ?gP,X?g%Bj?g.^ڻ?g?gUQ[?qG?r   gT?g!D|ή?g#?gTl?g0R+?g\Y>%?g^1CI?g7џ/?g̅c?g䨂OӨ?g+o?g3D?gj?gXYiy?g\6?g|!P?gH5?g޴?g"A?gFNu?gÈz?gd䱤?g3NsZKt?gY5?gg̋}?g:ׄ!?gX7?g0%?(@  gHm~,r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   gFH+)g脾?gn.,r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   gI5_ gz@r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   g=,Vg^#gYԙ?g9Ϟs?g^{ԿgҥٿgJ7kr   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   g     `<gN+?gE
Ε9;g?$?gCe©g?g:?gyf>ӞgrǿgΟHQ?g>?gsgU=Ho?gegP@Rÿgdsx?gCVdc7?gD`6Rֿg0*3?gaPgcʑ闿gRu[g>WL)?gY<Pg0&!)ܿgg?g(˿gMG?gb[ڬ՜?gT?g1g|g|v>­D?gLo=?g
pixܿge?gbS势ٿgQ[Agm7)Tl?ga,ܜ?gѩֿgKIt?r   r   r   r   r   r   g     gP*}p?gzɿg6r?g%?g~ϐ斍?g] j9ܦgyṖ?gm1?gr# ?ggJUg>ċ$?gI5ݿgǱR߿g7.Dhgm	пg?r   r   r   r   r   r   r   gg\{^?g)-n?g_?=gdeÿg.l?g("?g1V@r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   g<g\$?g:R?g5
r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   g<gQS(ѿg(TYT?gUV&ghѻڿgeʑgV?gjbUg.C?gF'z߿gqݔ+?gϻ?gqo~q?r   r   r   r   r   r   r   r   r   r   r   ggQcngE$пgd<'?g94ؿg`g)£p?gTFNӿgLs?g[ڿg;q@r   r   r   r   r   r   r   r   r   r   r   r   r   g<g"] g&`?Ŀg4#p?gےkVпg[gZez/gڂ'g@M6r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   gg~5J{g;Q?D?gpz@?g㸻DgTsjsjĿg/ɿg\ ϿgN2~g{!ٿgJg[ZP|g2b4#ǿgϽߑ9g0jp g;3tbgTkܿgE8uƿgKߘ?gg#ɤo?r   r   r   r   g    g30Їݿg|a?g~ire?g_(?gpɿg׷U1ÿgy-:`?g"RWZ?g1ѿg!`ٿgwZ᤿gG}pӿgI9gpdik?g/jޯпgQ*пg߷g.'
?ga3(Q?gwDr   r   r   g      gSgl ٿgCF?g7"
~?g\2?gjʴ=gVRڞȿg?gun?gs^gKj"6̿g<ÿgh"RĿgM?g+Ї?r   r   r   r   r   r   r   r   r   g     gBE?g]`?g4hbx*?g,8p2?g?gFuKY?g
!?g>BVg>o3?g$<2?ggO~oEgtoձ?g4mdԿg *gP@M?g%*?gk"?gћDD$?gHJ?giҕgUW?r   r   g<g/y?g.Qu
Fi?g"J\?g!F͆?gǊ˯?gvn?gfg
ڿg3e)?g?gػw'g9?g[tSӣ?g[r   r   r   r   r   r   r   r   r   r   g     <gTŤ?gĿg`1y4?gdv?Y?gNJ r?g.?gPÿg?gHР?gڔ0m2ͤgZ*S?r   r   r   r   r   r   r   r   r   r   r   r   g     g=?gɠ4g<þg=OToԿgǦf?gξϿgE?gi"g{%r   r   r   r   r   r   r   r   r   r   r   r   r   r   gԼg.g;g"PVYg4r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   ggvHUqgFJbۿgKzA?g[TD?gw$r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   ggqg<aXܿgBo-?gkv4?g
'g\*kӿgV+>}?gA,gvI඿gFt ?g@cnF&ݿg}g8.!?grzl?g	yHr   r   r   r   r   r   r   r   gg&ֿgRiĿg|.?go?gH]m?gM kƿg5&_)g揾
geĔ7>?g:W?gܴSȿg0{yʙgug?g3.y?g-yxpgco<?g?gT!7V?g}gN8gt6g?g(aT?r   gg-Q7ÿgOg	SъF?gd?gŪ[8;w?g\sa[g(Og#>`?g-]l?gr?gfȜRfg etgvI~A?g:?g}CF?gl<g+gS?r   r   r   r   r   rX   g??g@{LǕz@   (@  g>(?gG[J;w Pag5a\0gV='?gCӿ?g	cI?gPſgn%ȿKX9_(ο-n7jտg 	Aӿgpbz)gF(7޿Qb,|̿jKտ߸tJ쿿g'r۹-?Lμ?g_2?g8?T?gaz?gYO+j7?gG[J;gqV?gIlϼ?gnG?b.uѿ-}	CgT-a;g08?I~`?s>&?^	!?g1Ec?%^?gK8?gF?J4K2?gxmX?g9ZwU촿gCG#ʿg]zѿgɅL&vԿ^gԿgn1ֿr   gIlϼ?gNn?g5?9Dl$ÿgǭgjE:gZlm?g!B?4A22?g/<?guhs?g,՘%|?gp1?gN4W?gUIS?g1F?gAlꚑghg;%^gīMglgN_Egg5a\0gG?g5?g) )?gv!gmR	ܹg_*}?gjr?M,?g.R6h?g_M?gĕl-d0?gE_\?gz<?g(=q?gE1?g#?g?e`ƒg-
gg{v;g(GĿpϥ¿ ÿg/ÿg	W='?r   r   g!gvS 	 @gOS/00}W?goۿ犗}eԿgMtWjԿgϜۿ6g-M8gܿgcҿg[JO`ԿgUdﹿg/3F?gN?:?gp6"?gi!?_?;P?g?gӿ?r   gg^R	ܹgOS/00gil]@g݀G˿̈M?g9y?g&*	Jc?;'S?g^B?5w?gCzgY䖳X#uk0¿g5hg	0 g0P)c?SsJ?gt
I?gPa_?@W?g>tg?g&cI?g-a;gE:g]]*}?r   gG݀G˿gX?g+̈?gӀ9}K?g-s??gױCp?gZT?g1?g!.?ghE?g|ΠqEMgu!_WgO㪿gPgPcEJga+
=?g]WXx?gj8?gPſg,18?gKZlm?gr?gܙoۿr   g=̈?gʚ3?zbz?xd?c>֠?M2Ԑ?gfڠ?ge?XN?*Y6C?g#'b?gGpgme\I̱gneC1޼Vp!ga05LgGW߄g&gn%ȿr   g!B?r   r   g'y?gӀ9}K?r   g,t?.ѽZ?q . ?gy;:?gn?ޑ?g.Wk?ſF?:1sN?gRJJ1ng03 lgFCͿᰠпiYBiv¿gar   r   r   gCR6h?g;tWjԿg*	Jc?g-s?r   r   g7;x?gg5e?gD*"T?gUm?gZN?r?gI?g>[g_q͐gA>ɿ$+reϿgO Tg 2/眗Wh騿r   r   g羇<?g_M?gۿr   r   r   r   g/g?gsU?g/l(*?gLض?g6_˕?B ?矞u?g0 g#cֹgrѿgOd8^ٿj_ۿ`AۿJFۿg 	AӿgEc?gths?g|l-d0?r   g^B?gױCp?r   gy;:?g,e?gsU?g^0 @g/a&տg?g(Zb??瑝jv?gIT訿绢9[n_XҿgOۿ[߿.>pſgh8=;fgpbz)r   gԘ%|?g_\?g-M8r   gmT?gfڠ?ge?g *"T?g/l(*?g/a&տg8(@48q' g<ܞ,1?2l1m?%g1?gih	¿g٦/ҿg%>cmֿgmUq׿	jngLƱgF(7޿gK8?gYp1?gz<?geܿggzg1?g[e?r   gUm?g˙ض?g?r   g֗l@g6grԳ?j`;e?gggŻ1ggrJӿgbW߿g33zտg6aտr   g4F?g7N4W?g՜=q?g$cҿgv䖳gL.?r   gRWk?gZN?gH_˕?g(Zb?gEܞ,1?g6ggQcz@U%9w?/lo¹?yO8f?gngP׿Pz;t˿罽7[ϹgJhr   r   gIS?gi1?gmJO`Կr   ghE?r   r   r   r   r   r   r   r   g$Y 6?E?gGkjĿgo^P]ֿ`$cܿg+bؿQtԿGLԿr   gmX?g1F?g#?gyUdﹿg5hgΠqEMg'b?r   gI?r   r   r   r   r   r   g1?th?gad\gp#dɱ篬~ĿɨĿrt}IΩg9(r۹-?g9ZwgBlꚑg?e`ƒgL3F?g	0 g!_WgGpghJJ1ng>[g:0 gnIT訿ghggr   gVkgad\gB&z?g A<$?gc7?g.Oa?gxg(EgA?gNk ?r   r   gLhg
gr?:?gP)c?gO㪿gfe\I̱g3 lg_q͐g#cֹr   r   gŻ1gnr   gpg@<$?gKcb%/?zY?g[?g*?OgS&3?g75?gQ_2?ggCG#ʿg%^gש{v;g^6"?r   gBQg&eC1޼r   g>ɿgwrѿr   g٦/ҿgUrJӿg׿g]^P]ֿr   g?7?r   g}Q?g|Ԯ?g{x?g.W,?g-E?g8?g]zѿgīMg)GĿgi!?gt
I?g%cEJr   r   r   gad8^ٿg*Oۿg %>cmֿgbW߿r   r   r   g Oa?g[?g|Ԯ?gZB?g?gLB?gPk+?r   gۅL&vԿgglr   r   gta_?g+
=?g05Lr   g~O Tr   r   gUq׿g3r   g+bؿr   g}g(Eg+?g{x?g?gO?gk?gRҳ?gaz?r   g^Er   r   r   gq]WXx?gcW߄r   g~ 2/r   r   r   r   r   r   r   gA?r   g*.W,?gLB?gY?g??gGO+j7?gn1ֿgtg//ÿg?g?>tg?g8?g&g*ar   r   g29=;fgUƱg6aտr   r   r   g۠k ?g75?gQE?g>k+?gRҳ?r   gF 	KG	@(@  g o?Nn*1..5Q٧3瘛<T?97?lSς?祜nʿ3οI;>пFPf׿*mԿΧ!y1޿狖myfο_?ٿi~'ٿǶ?;X?y*+d?b?'v?-mȦ?*v@?r   ga?1)?狲п?A}ӿڣv~ҡz?_ VvC?L?6?sV?yj?7-K?i?73 ??s_ֵ?_Q8Z<!ϿmBӿ_vԭ׿8?ӿV"y׿r   r   gW]8?n8?D
ÿה(8.6O'v?(?=)Jl?xݮ?7}/?p/MM?w}\d?獆Hk?'G?m&?x|瞣 ˁYtٿ@0:ĿqQV	躿r   r   r  g?Zmf. " [?[\=e?YY?StX?6Yܼ?>lgW?.IU?^2Q@?.<o6~?@a>?z?4ɕiV׬7Yſ.8¿G"ȿGS(ÿr   r   r  r%  gJl''@(/Z$H}?&#@'6߿1Z?
̿G wԲȿSbDʿ1\JؿlvܓؿZ@Ϳ;HԿTg*Sֶ?q;?[7G?'?~9.(?`?c,[R&?r   r   r  r&  r9  g./A[@j\$pQ??d*?[Q%8?P}?Z/;?9`o?79z翭蚷_5yu\-5翿qHƴIZ~$?nC='?(K?,V?`U?o?r   r   r  r'  r:  rL  g!/(?s?߮? CfE?xs6?\:=?_?ҕҊ?&O?Ƴ.
P?=nHaj>U<˷ 	8(B6qfqK';@pr   r   r  r(  r;  rM  r^  g a?{?ȹ0f?$??H?Z^yzD??猊#9?S2*?M@?Y{^]iOÿ4 P盀
`CRIXWUr   r   r  r)  r<  rN  r_  ro  g;<Ƈ?$VjiD?d?S̾y?M?ߞu?0 i?<?D,+?TlƤƿ3s=SͿB)rpտ.b7wԿ翯I/	̿bIƿ;59ſr   r  r  r*  r=  rO  r`  rp  r  gDRU(?̳̚:?p>α?Xp0S?X%fj#α?0
/?(H>?`鍸~N疐^ZοޜGӿv>^=[Űr   r  r  r+  r>  rP  ra  rq  r  r  g?vh?j,{s?Q?qSА?q-$?|!gb?QÿӲտ=4aݿtwݿaK(t3ۿr   r  r  r,  r?  rQ  rb  rr  r  r  r  g@{տQe?"} ?c	E?@75?_4$F6ĿL}Rֿgl20`t玻W:NLӴr   r  r  r-  r@  rR  rc  rs  r  r  r  r  gUU3Vk@p{e' b?lv.?]-,?j	hIƿ,ѿbڿc_@ֿy$sm̞߿g>r   r  r  r.  rA  rS  rd  rt  r  r  r  r  r  ga@+֤oh?׏?)=V擿gWrR|ѿ@(W>߿cw}ֿ|8Կr   r  r  r/  rB  rT  re  ru  r  r  r  r  r  r  g@PK8\?yϱӸ? wЮAШ3/Z+ÿmɿQyuYl+|ðpXr   r  r  r0  rC  rU  rf  rv  r  r  r  r  r  r  r  g:*OTH?d?9)\Cƿ|qJϿQ_ۿTqW4ܿZqտտr   r  r  r1  rD  rV  rg  rw  r  r  r  r  r  r  r  r  ge?0	k?-K
jk3&Q+?¿"ma}r   r	  r  r2  rE  rW  rh  rx  r  r  r  r  r  r  r  r  r  g_z?XK?N8f?Cci?]&/B?(r,?8c?r   r
  r  r3  rF  rX  ri  ry  r  r  r  r  r  r  r  r  r  r  gPl? [??l@?)Ŵ?3'-??r   r  r   r4  rG  rY  rj  rz  r  r  r  r  r  r  r  r  r  r  r  gNU?QN=?)ۙ$?va?XP[?r   r  r!  r5  rH  rZ  rk  r{  r  r  r  r  r  r  r  r  r  r  r  r  gjE+?fs?u?=&?r   r  r"  r6  rI  r[  rl  r|  r  r  r  r  r  r  r  r  r  r  r  r  r  g_Պ?3f6?r   r  r#  r7  rJ  r\  rm  r}  r  r  r  r  r  r  r  r  r  r  r  r  r  r  g~|*@3h?r   r  r$  r8  rK  r]  rn  r~  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  g$z@(@  gFwg ՘g^Z蒿gSg;m@(Jg<Z_0g<O4g(־g>jZWg$O?güZ)?g3AN1?gl=~4?g@'%?gl⛮Q?g.PfpJ/?g9(gh^D/gʷV>gn""?gXly,?g7C$?g۵pk#>gP>gyyD8?g*y*?gϽYF>gF@_I?gtgҼ8g\*;g 	?gļ'(?ghN?g3$?g}4?gGKف!?gBJ
?gZ+i?g}l?gFRFg!R[?g$!D?gwȓ?gmhXW?g-yºz?g9RC^?gM/'ZF?gBt]GDg'f g~ gb?gihg!I
-?g5J;T?g"ڪQ\?gt/QgJ_@g ~4gZ7g(gGUgJ4yHg L/g"g@7jg}Hg<i4N?gSMg%ZNtgdxv[gL*qgXRSgׯє<g5	߀?gex,}t?gN2Byg/*ƙ?g`Ӧ॰cg+kg&]D~7g.r|gw?Vpgo]S8g?gUP蔡?gYn|?gr.?g)(	Mp?gF1k?g5ʃgpG
g(QJg@ZFg!#[x?g]sV\A?g\7?g5DGg?g$R?g%ŮPgC gyXn t?gd	 ~wgJgB0Jog拰mg?^k?gd|nu?gb2?g']Rg9/VgA3g突Xg9,gRnd-g9//gퟨhE?g	
a?gYګgg^T?gm#agu`v"Og;;gaL?g[c<6e?gK}qg%-`a0?gbVj?g0}>?g??gԴ{̆g
&+gn9Aga"X?g52?g7Hq?g-j&Jg?g5/k?gG?g
a?g*Ln/?gC}gOhag*Ƥ?g#?g>t?g>Xh[\?g,tm<Eg&?gl/?gYgUYC`g.ϝRg{ZcgLM砿gJםg9XU?gޜ I?g>U?gzga[-|g	QRvzgP΢sgsըg92g_Wwg,~)gd^5_?g"<g?ği_giētg_fįa?gHwC㿇?gMRߏgq0acٓge'Lg:㳅g`(gw6;Z?g/"8?g(e?g?gPtBбgP(~W>g&xg9mh9g|`qgEQ?grw(-'?gRig=xgMM?gM0g:aGtgZ'B`gIӸgJ8U'?gs@3goE3hW9g9v"y?g.>BP?gB"R?gagS |ZgQ(]?gO"	?g5=ҁ?g/8*a稿g?i֥]g}a￵}g1tgE\D~?gylu?g[U4?gsc8?gXPUkg;Ϙ?g^Ʉb?g:D;`?gؚ	n?g4;4pObgPU8g2?g	8b,?g1[3?g"?g)uJgS,Xb,gtwp?gpj?g͔<uſgW@ȴg,ǔ?g$H̔?gL=?gm&*4|g=Ƀ8g&ݠgEҝg2Jgչ׋?gcȱg502igt 2xkYgz-@gN=
?gT?32g|`;}g74pcg,㯔gy#D?gy?gل6pgª㼿gehg3CմgT̗ݣgjH)t?g0=k?gI?g
+OgVigB=g2pg>#{?g 2hcgM`?gz&ñj_?gB=9j?gN~2{g%$?g w
$?g
>~Ԡg?@ؽg8?g=?gbSRƿgq?g{?gSi?̆?g
dg傶Hg;Scg앂gEE{?g身|?gO}?gc&o5?g aLg-PcgKEg?gd5J,g8gg%6ltg3qIg?gE겦g'&w?gM.g'{ gժ%Q?g=a FS?gL;,?gݾ=\@?g!af-g_{-g09?g;?gJT =?gh?g) JgGg<"5Jg¿gg0?gIeؾ?gKL?g%}JM?gc֊>_?gػQe⺿gCEեgߥ>k?gܷXƝ?g6ÿg"pgDo(f?g[gM`paġ?g[-C?g6Bg<2G¿g~x¿gVZn<?g
ٔ?g{gDgGDsÿg[T8ƿgRQwg%hg1v:>gÿgy ee렿g7<?gX,3?g0g?g*u<?g~+壿gPg,`?gCltkѿgR?g)0?g~ǿg:ٳgEK
ՋgMNglg?gf:7?g,|?gZI.?g*ǆOgJf?gұ. gg"ԩ?gY?gXVຠgZ{_gv}Z~<?gQ?gSҿg.S?gdp|<跿g_g% ץG?g˥ݭgĉ4?g[?g?f?g>`ڛgaT`?gs?g!V9?g_/?gcN]U?g2ff
¿ge_Ҏg븘J4?g5|пgL%?gWK?gs|5?g0~:?g;FJ?gfǿg>o	?gmſg*?g^5iHGͿg0pI?gjk?grȱ+gl׋?ﺿg/?gϻb a?gS-Bgflfg:EMg1·gkA?gg7?g_z!?gzΏhg Pڒ?gHl}?gE _,ug!?gxH)YgjKʲg/LgmZ?g{gl ?g71ҿgr;!`gt	4k:?gbz?gXv?gO?gL?gk@n?gyPhAegΩtſg-?gdtgyyƿg~o?gM^ʋ?ge5_n?g[CNg"1fпg,e2g^[_ҐgZd4Ggp̖gU?gRp/9g0l?g)l^ԥ?gFb?g^ӵ?gЙg!נ?gnOg>$bgtMgajga(
ſg!Gog~i*A'?gԜ?g4HHgek lgFGwԥgt;?gOۮl2?g7r?gRi?gT!?g >G?gzz?g>
?g=Rg6g.8?gBdZKg5g0+ugE~{?gADH?gpT?g@Ic ?gW>gn1?gž)kg?gIz?gC%?g~zCg[#˿gh	пg:b̿gs˗uпgCjȿgGTǸg(Fӿgߥ#Eg?gTf??goQjG_?gpOwȃ?g3I?g Oz?g?g6gh}?g.GQ?gp-Xؿg}Hѿg/F:?gi<Z?gyy	"ں?gW
[?gp!gr?g< k?g6(?g{|:(?gzx9?gŤJK?g-Da?gn?gT*#?g	(S}g?gͯ?gӨ}g jRg㲄ȿgx0пg41:b̿g8-)ƺտg5IpֿgXd{dҴ?g0{}g!pCCgIF:g.RH!g'xʕ?gvDxg #gg϶bH?g޺iϿgvW?gF7"g}OV8?gC;V׿gQWg=ՃgM:?g
|?gxZc?gwW[֘?gc&투?g6j?g/?gğB?g(E늮gvşgozg`X.%?g<gY?gj7gF(g-X
ڋNggIXg[	geg{-	
g2C= )1geT
g";c?g=$?gՙr?gv?gǵf?gc?gN>8?gY@gQ`l@g9@ZGCV)r   g
?gy*ڙ?gqL?r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )r   r   r   r   )r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   g^r3#g@   magic)r>   r?   r@   rA   )4ķ?PZd??r(x?gNßT@gva?gZ7M?g%T?g?g":C?gJ1y5?gJ;||g?g?eϱ @g(@g9?g]>Y?g|.K?g-w?gb	u?gyP?g?g $?g51 @g1Z?)E)7:@	ER !@?:|̏?)A[۳;J{2JA<fUɞ*?{8T?g `@g`NЧVQ@(vU<w*@W=g]rK ?   gXJ?9BӁ$@Yw+@0A<#f"?oYUR@guHp@(@  g}ˬ?CtvJ¿¿?qx?{_?-֕)Ħ,F#7̻w+p㳿 ұ 'vj¿b;A⽗B֩
tnSu^Y<Kc^1\?4Z$?4.?痩-?~OZ?V\{?qh?r  gE?N4d"|?k+
?YȐJ/#V?E<(p5_?iLݠ?	??+?楜Wߤ?m	l?aEʵ?|iU?Yܡ'?`?2vB0Y,_pIͪds|{^Rk'b\	@Vr  r$  gr\z?v?
hƠl_Gs;i{Xܣr?xƓ]V?|?狪j?z*֤? QR~?㵍?pJA?0T߻?ﴗҬ?v*s燯hxU:pD>!^S-0砃|Cr[ޖr  r%  r:  gLx?]I_vq\4_?!Va?yf?翧'?2ߕ?s.77?8/G?&?.<ynY?/c[?ot?C7r4ۦ	t77&O	痉j8lx񑠿r  r&  r;  rO  g*jW@ۤ^LkD?຿<i%?B vOPæT8爭T	32XȒ^穼]g=KqSFh?_pK?MA?0Aӄ?uc؎?Iq4;?2jqٮ?r  r'  r<  rP  rc  g=<12?2LK?y#?׆8+?s\퉍?W׸?U`f?4u\:\g& >t=il޽犳p?ͯƘ??ڪ?ARd?ٜÛ? [/?r  r(  r=  rQ  rd  rv  g۞Y?mZň?m`3\?N+6?2lzT??@3|?
(?Aru?ah]?9(e^wjZvkK척<{D᧖r5n:N
stb?Kr  r)  r>  rR  re  rw  r  gwvR?hUo?{A 
?t ^m?'p?ڈ߸?&1v?u^xo?8J?T[ZZ?l]buܓ]^N>̶/ҡtdJ!8n/N<0zr  r*  r?  rS  rf  rx  r  r  gɜ ?6?P%?c@?_ΒB?eQc?S!1?ysG?J?ҡ>HLҫU]8Z[~1diӧx'9Z9\r  r+  r@  rT  rg  ry  r  r  r  gD&b=F?;Fv?g3a?ھYRo?C?璜vY?4?TXޔ*dkE̩cb繚&m{A	͹Sr  r,  rA  rU  rh  rz  r  r  r  r  g\J?ƌ6JPG?>z:?c5?(\?1c7? ?ŴO[┿ڠq qESV痝_inЫr  r-  rB  rV  ri  r{  r  r  r  r  r  gx/Y?!H;+?繀AЮ?x-9ɶ?#S5?JĘC>֡|?;}u6Wث,_OQr  r.  rC  rW  rj  r|  r  r  r  r  r  r  g@?/޿l&?xF"?ҭl>?uR[FGYVBܶ=鲿2\R	Ch jɿr  r/  rD  rX  rk  r}  r  r  r  r  r  r  r  gB`?fc̿(Tޫ?a>2V?ӺOpT oy+zT[oLPY\N;a<'r  r0  rE  rY  rl  r~  r  r  r  r  r  r  r  r  gw7?ݓR2?<~?홳d"JJC~bж}Ȝv¿}6X.r  r1  rF  rZ  rm  r  r  r  r  r  r  r  r  r  r  gv4A#?1?b?1`|H矟}>罬)?ӷNwAOƶr  r2  rG  r[  rn  r  r  r  r  r  r  r  r  r  r  r  grh2?爁tZ?9%*H4,f&{/Eɘ웿|؝j0r  r3  rH  r\  ro  r  r  r  r  r  r  r  r  r  r  r  r  g#"?<<t?	?;gʝ?U?砸u>Vx?iGm|?r  r4  rI  r]  rp  r  r  r  r  r  r  r  r  r  r  r  r  r  gƣ%?c7C?A6Q?6<?#o?mǪ%?r  r5  rJ  r^  rq  r  r  r  r  r  r  r  r  r  r  r   r  r  r  gs^DD?	{	~?D?>9y?&~E?r   r6  rK  r_  rr  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  gv>?6]?ʴ}?"-r|?r!  r7  rL  r`  rs  r  r  r  r  r  r  r  r  r  r  r  r	  r  r  r  r  g~ml?1LG~M?r"  r8  rM  ra  rt  r  r  r  r  r  r  r  r  r  r  r  r
  r  r  r  r  r  g5ӹ͌z?ce:jO?r#  r9  rN  rb  ru  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r   gQ:+w?)r>   r?   r@   rA   r  r  r  r  r  r  r  r   r  r  r  r  )r   r   r  guHp?r
  r  )r  r	  g&%@g0&@ghdn(@g&m(/@r  r  )NZnumpyZnpZstatsmodels.tools.toolsr   Zmpg_bsZsmooth0ZtermZbs_dimZdimZarrayZp_orderZbyZlabelZspmZknotsZrankZnull_space_dimZdfZS_scaleZvnZ
first_paraZ	last_paraZreshapeSZcoefficientsZfitted_valuesZlinear_predictorsZdevianceZnull_devianceiterZweightsZprior_weightsZdf_nullyZ	residualsZsig2Zedf_allZedf1ZhatRZnsdfZVeZVpZrVZgcv_ubreZaicZgcv_ubre_devmethodZcmXZassignoffsetZdf_residualZmin_edfZ	optimizerZp_coeffZseZp_tZp_pvZresidual_dfZchi_sqZs_pvZscaleZr_sqnZdev_explZedfZ
dispersionZ	pTerms_pvZpTerms_chi_sqZ	pTerms_dfZcov_unscaledZ
cov_scaledZp_tableZpTerms_tableZs_tableZsp_criterionparams r*  r*  Klib/python3.9/site-packages/statsmodels/gam/tests/results/results_mpg_bs.py<module>   s   

*
,
5

Fpr B  D B  D B  D	 B  D