From 8e2b7a923db673677087f69b1df0cd1a1a7c7a4f Mon Sep 17 00:00:00 2001
From: Benjamin Root <ben.v.root@gmail.com>
Date: Tue, 17 May 2016 09:58:04 -0400
Subject: [PATCH 1/9] Defer to 2D scatter() for handling of 'c'. Closes #5974.

---
 lib/mpl_toolkits/mplot3d/axes3d.py            |  22 +++++++++++-------
 .../test_mplot3d/scatter3d_color.png          | Bin 0 -> 41598 bytes
 lib/mpl_toolkits/tests/test_mplot3d.py        |  11 +++++++++
 3 files changed, 24 insertions(+), 9 deletions(-)
 create mode 100644 lib/mpl_toolkits/tests/baseline_images/test_mplot3d/scatter3d_color.png

diff --git a/lib/mpl_toolkits/mplot3d/axes3d.py b/lib/mpl_toolkits/mplot3d/axes3d.py
index b1032f45091e..e4ff16cf9410 100755
--- a/lib/mpl_toolkits/mplot3d/axes3d.py
+++ b/lib/mpl_toolkits/mplot3d/axes3d.py
@@ -2209,7 +2209,7 @@ def add_collection3d(self, col, zs=0, zdir='z'):
 
         Axes.add_collection(self, col)
 
-    def scatter(self, xs, ys, zs=0, zdir='z', s=20, c='b', depthshade=True,
+    def scatter(self, xs, ys, zs=0, zdir='z', s=20, c=None, depthshade=True,
                 *args, **kwargs):
         '''
         Create a scatter plot.
@@ -2233,7 +2233,9 @@ def scatter(self, xs, ys, zs=0, zdir='z', s=20, c='b', depthshade=True,
                       that *c* should not be a single numeric RGB or RGBA
                       sequence because that is indistinguishable from an array
                       of values to be colormapped.  *c* can be a 2-D array in
-                      which the rows are RGB or RGBA, however.
+                      which the rows are RGB or RGBA, however, including the
+                      case of a single row to specify the same color for
+                      all points.
 
         *depthshade*
                       Whether or not to shade the scatter markers to give
@@ -2262,13 +2264,15 @@ def scatter(self, xs, ys, zs=0, zdir='z', s=20, c='b', depthshade=True,
 
         s = np.ma.ravel(s)  # This doesn't have to match x, y in size.
 
-        cstr = cbook.is_string_like(c) or cbook.is_sequence_of_strings(c)
-        if not cstr:
-            c = np.asanyarray(c)
-            if c.size == xs.size:
-                c = np.ma.ravel(c)
-
-        xs, ys, zs, s, c = cbook.delete_masked_points(xs, ys, zs, s, c)
+        if c is not None:
+            cstr = cbook.is_string_like(c) or cbook.is_sequence_of_strings(c)
+            if not cstr:
+                c = np.asanyarray(c)
+                if c.size == xs.size:
+                    c = np.ma.ravel(c)
+            xs, ys, zs, s, c = cbook.delete_masked_points(xs, ys, zs, s, c)
+        else:
+            xs, ys, zs, s = cbook.delete_masked_points(xs, ys, zs, s)
 
         patches = Axes.scatter(self, xs, ys, s=s, c=c, *args, **kwargs)
         if not cbook.iterable(zs):
diff --git a/lib/mpl_toolkits/tests/baseline_images/test_mplot3d/scatter3d_color.png b/lib/mpl_toolkits/tests/baseline_images/test_mplot3d/scatter3d_color.png
new file mode 100644
index 0000000000000000000000000000000000000000..1ded1b33885072645a41bb261781f0bcebe28f0b
GIT binary patch
literal 41598
zcmeEuhdb48`2UfakWDt3#|R0@4B6YUS5_IvKK8E6kd-3qL^$^5AcU-tU1nDHCM*1&
z)93sB{S{x=<@)sLy3YGN@B6u**ZsO*_v<NITT_Xckd6=nfe@=GE9gQX*enpp4I=_P
z@Hgq+;cnnRIG&GH^a;REAc1WZ`2U-3%Eq1$$nB??U#voz0tX0$1)`!Lr|<h=Yc}AU
z$fP;*-tM<K`+Gyuem0MUqCGqC9FweiUy1W@<ZF(S_G(iZDkz3QzZ+0-pS>VdHKmG%
z#>9{!;U<cN9Iz4D$9lQp>2Get@FtPIj>~w-mItNb=!97>{urumx?laM{;G2K{mItW
zx2>%mv%&Xg!R-BiE~L|{tE=y?imjeJlDQj$BvAPG^FT0z2||YXLEvMPvSEHu#Q*n)
z{|mzZWy1fS0*%&nEK)KuQbnD=2@dIkL?72#!Eu3wV-b+)bbJwfxqD$$KDigf{K)}1
zzx#1=Gk0>cT|@(qLqTDTjz~5tijG{7EO0(@_hSD#HpEpl<McXczS~V)p><xVbw1J|
z-QDU-&daDMuMWt|s6u-kN>+j!KkDKJ>|d;V5q1b#rQ3f9=Q#bGz(PP8!74e_VP75n
zf`lxphnJ(pB0=sZ85vI;W+fU(?wHZiFM<lBIPQJ0x3(3m!SSa`I2?;k#2v~S^L-bL
z;JhQOA&(XGvB`XK<RHHcDp;;B{ac<;D!T8zsZC1!Yw^qA{H+nszT|Bk$}UUZ@|2!+
zQ!{>6HY8dI!peqvW0X%2Zsi)u+aXA*NM<(qYDwIkbivQ!@5!Rm3zn$Snvtk-iemOf
zAME;Y=AlFH;CizyfdW?{Ztw)?ox}!tGBO<ph0I9hdUA>Yp~eIxLGd_LGn+2%2`v$x
zdk*V{bZ09rnK=8c;aEIkkK6#-2oSLBc1&zgUoo*!u(rOZn-hAO!ZUyJblSW5Vmwe9
zgfPKSbr}{fNc41C^vTCzO)ED1+Ys*?5qcxV%V_e4zsO`_r*BY@=@`EYS5(ly(a3G4
z{%dqYlA*iZza_#Uoi^?X*uJ%s+m<nkq%5BeNT`>?y%Rs<T5%bP*CerCo10WN85Jdt
z_Ryw~k`43FhLVzHst5>bZQ2S<`q4&&_1cCHKYux3PeGH#{K!2`?^BYE3LeLSYIL)z
z2lw+HsCU{6a=)#`XVr#at5%+Z_<OMpD{JF0IUtxa2JqsSScv3bKe_w5dr7svQ)Al0
zaO`7L8#_O@3La+x9>*5!v4TRRAPB$2NbU_!tNKA3LE2x7i`LH~Q4-N!Jnxm~Rq;5(
zn;&~T=lJ+sn=OPtAeEF1K?8CZ8J-k}gSDmc?FCEdlK2m2(Y0g$hK7bVE5YZEJ@1$2
zSRPk>%8SXt`LPXS3pQ+>@3s=ODko<{l2~DKL^e_4cDMU4^p0kx(%iDK=6GNF^T~&2
zZ(op;i7*nXcAl`u&Q=FYS4<bP6NHN!scpiNQa^Wry;8Qs<h)0e#BHH<P3quSiQCh4
zWQt^b{a<_~T?TN@w6)o2n)>IH4Ysn553xEgl{$xLuHX9nvBrG2H|D!jkjj~nt9})5
zdO1=naRzHirfk@1qF9(gk)58wxk1$r^Rb|)-G`h(YFAmyy96DiWt(3^2H|1RQAb}u
znkehm40hDPeX|O@GL<q%>Tus;<B`3Kn!0@!0O=&A{aFk$=@PSzv+R?8ko7nm^T}76
zXh98%(RIp>g}vpQrDoV8PkK{$tYo;!x=l#bHQ_nj6=3%gr@*IlR;hRdEz2j}81Qu{
zpGq=yRF#Zym0FL?*!SkJ6UqOmW2D)d9F|RNQH_d9?*hw_;4K}z`>>KK4o5Mw+gA8#
zYYd_7G4$b>Y~@yt!j;J%C^%VoV7oRL5^pT~SxDv%6d$ggT)Z1v4WO_^86@Jx$dPsy
znaZA-&vrdz`)pvr^Ar=Q(LEPcg-q7noMp7Y*QNbAoCiem;_)wz{fS3Kgfa#MX^Eik
zt*DEz7COBREwH<wO=#?T-QyqN<<}&^1{G2)OWdA6sr>0I_+H&R976UEQdly$qz^AA
ze%Q@pK05@@>9GI^R$v1*q?D603fU$JFE$z;x{JxykG^&<k6=%*yV&GZ5-m6uhx}WZ
zB05;>aXFEL|2~j1sjRJj3P5iY4S?}7n7R6A$*{z2O@(4uUYIp0P7yI>EH$#=7)tgi
z=uSjC=9@4(&##x4m~(o#;$iHUsgzZBf;B7$)^5{;W|jyUInTvx0)h<0$}x`ZQt}@F
zK_LBAf=meyD?c(SHlh1#nPnS^DKT3-^%{VpgVvZ|lr&kd$esO^pY9AVpGGaXMGe?{
z8qgd<3q`1>QYXe4WgRtCstlvWX*vXFjMsb48LMmpQUw8KF=9>~oXUZH?Lyk1m@1qF
zoSe3<E=7hnYvB~Ci2U)F?4$MJERYJgW8+U}Ro8z3%cI9*QOL_g#^Eq_x7%0Cl}~nn
zk2KRKy(1Mh4@fk?PI{o?>|dg9RuzWieoTXSV*D@G+`%;9+$L^-1WPnf_{B=to19$8
z9rH_^1ia-_+IQ-l0H!`X*I80Pr)Ya)D~5vfQoarLQTdi4;Wx1iV&`c9<9vc#E14s*
z`zbkj%H`zy)s6guWaD8V4m=tz?6eV29i^L|v(w^Wfvo~5Gngk3Jw30R8UJtyWVlgH
zantGfa*T?Hsnz~7BGsruLd@e;z~cdaz^C^lUCRgM@rhm*M#=MXD5BBKqp~_zj2L`^
zVYW_{XG-}y!n~pE7RbL$&4ejk!}@3$E<}ocUhxIklLUY=lc_ynH-{|QR54Ykk6U-O
z(xQApoCL`YpdrG}_`@&Kg&MpVJ3Nl3wAQsIUF#;8ryFAq-O}=t@V4M#zR-rmoWes8
zjz%61_{V-n;v$aHFl!h%8uPKCFF;VZZ(-6y-6T`VgPd*B@8}}v+v;H^)2#w+4qUZV
zcWo*C_h#9n<5J~e)>YZWAm%#MV8isAkO-#uC94*T4iSxS_SFrE^kk{!%7!P+crlUE
z8JK0<0MA9Gnyq5xiIG#{{)RwH4}bMxkIStd^8FlV6rYiG2OmDU`~sHJ3i7Y;5oX!^
z3Ie_sVU266WD%@SK9YT1a<@Y;eJCN#<aELpm_J#f>i4AMj+w{g|276Ipd7m56PpjH
zV{t!?u3JD;cP{`4e%=L;ZzY<MX}+3!t3$9R_^~F=R+$x1WJd#mRC)2mLt;u?fFdJ$
zVmo^Z`j3RpNSo%PLY7q&Xd}R<GX^-z4rxguYhP;>fc;R<`=|%bN&A;ZIZW1{zMr9f
zS7?j@1pR|JOUiczD?UwSx4PY}Xz@6cZ-O;*VLkM^`8J9Km#9EgTm6$k<}}!L3~He;
zH3q6-aaV4C7jTtv$5)^N8&mW3Z|KY1-YBcndagyB6f3G-f-Zgp;yfM+@-S0_GAE|(
zZfSO;ay2NrpzMQ!^uv$PL=2`|p??QIQFX&Phufx{9O(}ZA2OK=&b}tD`S)m~icg(0
z`BQ|tep4ClJDA4T{^vPH?A%x}06<KMZ{uNu^+=F|^+=c`J_*;MY<Fj9D5IbB`%yQ$
zhEht3p^&d+`_g+$M+m&bx6KU!P^|zb1Ca)}1=f}s3f9KRQ9g(*+syb0J{^T5cx#8j
zcEQ)DL;5)z3Zd=z>#C_kqHktpE1w&lAhG1D3bR1A#bC<$tEeMv_fgG&4Q{%dcygft
zVfKBS>lbrJ*cEZDyW}NxsN!Zi*%k9V_8Q1CS5wM%jEO0`R)y|J0yxtV0-@(N1nGK!
zcLxW{c{*a%@4qq~b#j?j<fR95&}EKbEh{f4!@*J?Ro{+d@!ohE{IE!~!1KP7=`VTk
zsE#xi0)XplIG9JlsuEc$EHU}U&>fvUrB$E@@??fSADSFY)sy_eqCE1aGc~>+WU)7a
z+uyu=a{r`{g2LrCDzdIs27cyD1J5B~nfqM7Rkj0?9^MUd@=p(fsqJb&H}BTN_!rzU
z0BoNc8{`H$Ou;(|AOP8EXlNe!EqS1$-W>rBBhF1$%<f?NMwnYi&6M-E`CBs{vZjl;
zCwC#-W8+f%@CRp)z*^B5?%YK}^$t7g*0&D9PB-yRTj7VO#2jd{P>0|mWAg?(bL;SP
zocWVmtAPb86$NlPXn0bMJ>2<R8MZOiX&6}c68R4}2|y6nG5mLX+%IpIGH5<>NqkT8
z_P>hb>?$0noaP$pk@?za9~0HjpXbidHVHbyh2mo++dm4b6cmdN+;T~b%V9e(wh;I*
zUX0=k@_nh^efEG<_J9IV1z;c#Hv}68?-sX78n6CWST#Sp7YVBr6%&~6M0#xOKl<Yo
z>-9!60NOv}fn1_AZ~S(l&ZJr-`6<lWVz0x3Baa*)45NNcV$LRdU!(XrsB8kWKawuB
z`^19_thFB_*83iM>~zG|q>}$?DBDrIf%ot^DDk>Cu}P8fJR{25M~9>c!lnI&6cNVM
zy3eIa|APw8{>YQXFeO(dm|6fzh30yZWaE9#id_0IHGcUHF9z%%6ae!W8Jn09Ae&nA
z%h88cf>E24Q-v{B%9@7#n~_LV^WWp;9n*u8yp{6NG69BU0ZB8`DmpR)R`_AG9zFN(
z%Xs-mZv1KeEwXu&n@=0f*a&XS_^oEHe0Fnu0N#o}Z9-Yluli;OX&v0A`e)}QcSOh2
zbsGoA)gn~yG&y#gt#JEu5y!5)`@GyLYE)$LnSX1{(>Q<YCoI>swJqSqtm1JX;|D|r
zGA%fE{e}YeehvjwxBhWzH5iD^Hvm;)Ot*W?PdpyGL(t@by8!Zabgr`Dt?-;NuLr4_
z_s^bM_%$I`q*C=#^jSj$M(Mci5r}(920Z1b?B%Y-<Ai{8d@+0i?}WT8e5D|_?^yfA
zO?+k=p|X4?5&z<U&<Yi}SfWeZ)0iB&T!%h9`@05@?h>D>fO_>UE0H|@BB<}7tJGC=
z@bi-?Nss`5vlO47swRBYfSw&Uvp5oCg248+YsEY+qGMrDLq6Dg{*Ocgv}byQE+TB*
z3+?mTp;%d)90w~D{^IL^eMc(!NJQDH-2~z9T)Y^Fcim0u>>&W1O*`n~&or{7ZD6Um
zBO8*PjuF0400`g^Z%U?fzcGO0F2Ac}_rb5`y9F}gNJV+7!7Fnu4|E#HE9Ef;JE}3J
z>SOK)1-!el(8<aw2FIMWi{R#H6;hnjy-1nY*_S`FyMj(;Blmq@1F*gqydUeku4Xy)
z-ZOL%w4KT}9(U#3<tdZz?_XPocXRzI(0l4rQK#XXw7eY}c4;1baL&5+ecX7ZGBw|c
z#|s`(>Y;5-ix)Vbpog%$I!IgMclEm^efWOpSHjg=!rsxS)>6xE&^?!#hOmp>R_2rW
z5Sh8X2(e#?@&NPEdoL(Kty7jin216KvpC)0^tMNBL|+Q&n_!V0aB4W1|LUJ(RPVAu
zg0AMI4u>PJye!Cm1>fwG!y;9=FRsOYj3ET_`TC-jS>|tt%xjvKHAUaWXs*3Z3awxA
z%;%&G-2j!}gYxU}Tl=ouFM~|{P+wE?3BRqrjCXz08=|vaA6$9~RS(i@^0)W&PA-1b
z_BtAR?fbuhy}Uxwv=F!kF46CQ75gl7<*75Tm=b|SkjM$>$5Yy`VUdCIZ+z^@lS|vJ
za@6;7J4S;Wl|w*>BK*E4*wj4J?tz1vsHXZY)ny9;^AG9<xxC$#PV4r@liS=e6=nkN
zetwN7KByTVWeot;lxZ-lEN6D^0>WsCgyQ02Jf$p+{ET{9wW>QyMbjI)hcTdwjF%=-
z9V|HSz?Cq!SUp+s=#^&i5}7lPri9fbT3MjpIP#J}HP$+fI_uELJ2tb5=cpDXcq(;D
z-Z3|sTr|N^@*_7gEJrB?2fB;U@0i||XH6SWBefqTBr_jo>?ykKDvg$z<f+_{dc<V<
zQEduy-ugHkS*zQ)9KHNu>(p7r+bL9TiaIaH9l@T(gDzZ|6*(w>B&)pBPqO85b9(m#
zU}B95RgxiN`iDZGDoi4<oVjvQ4{W|Q{qzhsWb~UJZq8x<nZpl1DXqz#^)~naG2As|
zJ+*0GYCzEOhR4WDM9Q8YAc?_h;x{?1X>xBt+T}T%#v4_<Kf4$poXOF;?f6A2V>t>I
z91G9hZ0t1--7yk4*~PCmei^XyaF^kmL}Xe&$EVRE(cA6HjFL7{bo!wGL(HB=U;aZM
z=+Th|rD5@73#=kNN=48HSHK^YrLG_Q8wsSv$#mJiDa-~7G+q<%H1{e0-U6RGO;0)?
ztjah5U(|4YwFQgf$}-BLdTCpF%13?coBjR#c-(bBs&(W+2W8coXa)M3zEw<%0BbHS
zBarS~ev=#meZi1Iz{X8BAlvoK;wr0Vw%)6ANT+(lK>E-1Rr3iqFQl=dL40#nR}3gC
zsP_bX2Xr_45KJ%He=XDDj%SrkHr`+>C^`0*GBY!?5~S6?3W$Dbt1Zlwumscc!xwnS
zoftcfxUD)>rnaj#K^XNTyk4z4yssUg2dCnNUr!77V5Uehb?BJfTsm!!rm+*0af`CK
zv0lkQ?4*W2>{ZQnsGp{YQl!%ps*%CQ%}uAW(C;cBJ7DVS0+K_qu4cvvCL~Rl5Q0gf
z$bV$l&!U+(R)ELJoQT{ZjUgopp9GD^0ev0TmuC2Kkg;No5~MSf?Vwd^DE3x6kXDu#
zpZ8Z;85=kx+rzl(?Gkt}f)?@6Va;^9k$2&RH%<MObF9t0*lajCRr!6GP9;H$qaTma
z7S`i}@4F-#d13NExah?F4cdS3?p2H>j#$b+#kK%IcY6s<^5Di2^xRfekWHdaVW;Qv
zc37ODsy4ntY<-_!?Y-5|DPW|~4N8ap|EMgNlHnNrBY)hFk}*v+Io^}&vX}A8KRLo+
z;{o|iSn$<>BM^|t&+ey9h*y^EBNpx#E)w=gohSfUv9k-3t#YH`lb{0Cn4Hpq-FI;G
zT(8LM1y0``E#pToiyo)$?2R#L;Bc6!IM0dpa9;=Y+-i8+&yjWgcQx(mQ4vOI$dsEI
z13en}TfY9OO-+H_KK2dC4S$pJjpaWny{aZ$rweYr2ekSc?Yiqh$Ui-Vs!2tr2dk{5
zN;ijwOnv$4u)u0&MvnxQ<w2x1fo#zeaMmxz(KsmY5<Uqln5js{n)%m4mA0G}Q<Lfh
z+=-~HH<mXM^BFo^Y?Yx*sX6xK;XS;6-gS@8sPs?5PcOO!yL);|CT@jzFToKc1Z7N#
z9Kp2ckn51XuVH3ZnA{{huzL8!+fMhGCuDe^9lN2|R`<LPjbh9YF@HCfNH+rb*w<(_
zZWZs4G(Byed79SeziP`fAII!kzrfnV52?`QB*3cDtDyZ4cc}g;0@i@*Io8IEzT7B0
z<d*P;Q9h*Zjd%^ianoq!`Xc1u0c}a}Vt@_S>qX60+xBjUzF>ibHhkrWtNp6uE+z=7
z75?QH>hR<Z&DR%@J;_Spk-yKRchMbi-*tTu7ZPs3l}KmnVJZ}k=O+cd({$oL-{_=5
z;-kEya~@050$R9vOxLmCwwnIlj7NnUy@~o(dEzg7aMm|4p`+SE*ns2)dBX(+j1>8@
zL&ymI%@`3@$Ub;!$!)y{xV$`Y(^+A9VSAlYEy7wonn5A=(nPTZ(gCBD)T9CX>O(0B
zJPvbTnOJLP(tqhfsR$1Bbo-x+AO)S5|0)Kn-g?CQ4@X>^c3$%aX@#*PWkQQ#8daJC
zjTZ*O+=-)gUN=>zw!Vutr=_JKZbzGdz9f;|cr~%EUv0{};zc`Ta&l5t;|S}4UIk9$
zy_i6TO(xILU1N__)s>UiYi1tWZ0ta6G?d$ya%ggWkiP{Va22R{cGXHaVo|Q#ZcQeU
zi6P$Jo03VmzYivRj&COS;(Q$DzM`WV^#j>=5X3wrec`+`z2Q-%HP$SO!%^@~^DYtI
zboim|K(|%JWdN1!tPD^jR{ct_a+3uixtbJVkw*KFKexUmAK&~qOt7@?dv9E|3!<QH
zXXYVAe1W7anC^^07i>{Am*lX+)vV<5;4sFNl-e5yRcVT+FJz1yhVj{XFZ-%Ws_DMy
z`UpXC2PNt-8s~Gv?3}(%n5HRqUDr93QUFCO^1J03DlRofZ>74hy2ND%6sXZ~^2RcS
z{1a<>+qJVs3xPCMvu+hFBH2REe|o`Peyi5oU0q$ccrWC3h8_^zwNTViNltcdjC}}E
z=cI5JrO)StXe5s=E64@YMm$>Q$P>ehA5VNi?jBcmmwSEW0iP4@^7%O(;E{R&UI9^(
zxLcZNH3f>Tpaf9$5@NJQUvYk2@nbS59=DiE|FH9<zGTv=WODQLgZd5;j<ezN$VPIM
zeoo&*V}$m>$~}6wS*9l+FBxxhB@^5Wu1r4(^UPzs<<hRU3YdDSVa*^Q(QS2b1syM4
zW<nVD8#EVgFC6i#+IZdkTR-wk04Jv$NMR7?js?=syB#6};?JqWs+&Cu{2gy1)JC&E
z${&Q2`>Fl11w!y=_0dktzptX1|F+Y+nNv70i}l}-J|Js?$thb9_XJ9>A|f@Bze_?3
z*(2$0-T-~ri+=(a=)-un6<_r_r=!27no0F-SeI?MqcT4H#)klX^yBbzt_RvKW4|w?
zfi|E%`Yk9X!+GOtFP09vIi39Mki)V(E7%`0$@p<T;Y~@)AIUZZ@sBzz^G%5nTU-`1
ztXxl!GzR343|FY$^=oa!PWVxtwU{SQbOlN>ISK8<*=WfTLUr&CL88~)HsbNgO;en%
z$AN1+E*@V}W>ETFPu`D_NCf`A4ftE<YQKnL_sKRjdU}`4@>$zSlP9wgA2lc<m=j!@
z1o*^T7ud4OW#zK9b*c+Ur7EqT93}m{EC-Gmg|TF-OGP9a?>bl851UUnM-5a9P2e(L
z;aBSP5o864eF*016R!(&{XLn-lvwdgZSKgH?NEdi{(&wA(TZq8wRwBuaJu3I`S(a4
z2?GWvoIO_T?^$I1C3}K6E1><E+73!IduZvVhMU6l5_0#W<WutbGz{!CZqALx6}n6(
zM3G(BpI&coEZcMK;r}tbiy%-vIBx>+0Wo4NjgC%Tvm#;mr`-U(E586)h`yow$(Gbd
zmBusGKo&J8@eoeGh-^Zos+4h&h4toYj+SuoL_$~n8PRURQ75B^uUpu>hNkgY0d6eq
zdn*c;x?6iZ(4D00UGa$S@9+N<DVGlLZVEc0EZBHA`R{V0vj+s3RK@udf#VKe#4wui
z_N7!w*Nybf6-2s#i53`B;f$9Bz}~B~y?#PA;Y8sOoA$Q$l8%+XP<fg}oL8}GJ%@c4
zU1P-OHfioNF2(M66OSd2@xN)urDPJ%V>?V0&Q={XckfYScC?nXD;Ab2Eg2x*_~K8F
z47?ZY{&N2BfNR@{Qqg2PTac`RcA{C_H?>6leNWE11r7;ZY|@bk0!fAHx*DzuyF7wS
z)mQY7-y^xhH_>hk1!LAhXRBgC-}?SudfdD7>&UD4Q-X?aRTj>Vknnzg@{Jho7_crZ
zK3;`ZOU-k8ix}HcxtmV1g-{|{Wj~L*?Z^IcpcMm=X&EE)<27nL^fuSmY?rWCcbMK_
znoimO?mP8lj=1F3x`VNn+u?@F50T1C^ih_<rAh1<d<FN~4xhuA$B0OF4OAm#zusV`
zwjFlMtiKO_`ZW<_DOcHl>ZwCZgf7u@{I3^)_pa3Tj8}2ZA!n)I$G(FN6_N)U`p2Gn
zqT}-&o<$W%zn)|*x`syCMjK%)i_<UGtCq>D4il?<_fBPe6oogBpXhW_@<{;2Q#Y%j
z^8KefmPB34oQ+glV8XM-S$L?Rx)Y~+T<P55^3($j%Pb0($TaO1)iA>hD=nuP!<6A#
z(hk2gT<L6Tr>>jUDNp0(-3aknLnW~6#^rA6Rs0v27K?j3sVoBF0C>Q-D&na~U9AK4
zIA^q<)>%fna_oHRM;3w#lhVcCvB|x8<HZVgs=OYz!y)lYIrF>S2n#sWWAa#1=2hQr
z92~otN|52k$Klcw+ZNFZ6IO&BdmPE`YA#i8vq-Ovfoj8z#~2`-BCWj%srczkS9I9(
zabrNtsf)EyW6M{v%#N*e#9e|DmlvKIH9#}OH5eJ(%8bwIoL<Px$oacX^prMXMJBr9
zXU+BB-l>Hh?esj-h=sXlaSvF$NI)XuHV$G#zAg!UZSo=GD6H#pri0Z*!!E^d_9FiR
zE!0bMw)jht28HpmCZe#(QY-afggrU?{;>KnO~^jYKBlSte6iCJM=u6OCUo>N7}<$u
z$en>Nt>3q**@!Qt<+PYKG|GkwGo{r-^46;!%-w6`4iTiS`OIwzTcvEQcCZfe4LnXG
zQU%f@nB@9O#_Q!W2W$y7UtjJj(MB{J4nA_94?0<R)Oz@6DP*VRWbTP*IBQE`ViTYw
z&Q1~lP`*4L=R6CKv7tx~@3O|4t`lM)0ez*hX{i!5RlJ_XLOy)D0gZl4it}#3o=BSV
zE~hFu{y@zXes3j_`K79)!7^_5!p)opfNUKc42%GAnOw*mEa=}iD<>8wSqz!jTq+BL
zn_%%g-je>gc`V34vVSOZedv2Mth8UG6+~D+AN(kIfhHJh{ttXULW#)_N9JltCVuxc
zezz_A`Yc;a%6SALudizXX4`hg``5+dy)^pt8QgogbvBRcWBu;vx4Oq71BbEgB$;@H
z?5V!yNTg<P9$Y$2*jCMM1C7RFpxX<58K9cxy>tFK0SFfLE!*%%w;QD|*R_`B0(MIS
zqVib5b~e{!5kNH#VC*#Nu}flLOKwX;Y8{!DD_=hSVRSnMxmrWmtuE%AuH05P)b08=
z!qc0XG{>~vcH#jS{MKA-?5)6|lRz!T4kY#SVA023?5WsqK6P`}jUyPig2QD%Sc8IQ
zf4uWz!FSR)NzQf5)YD>f#Wg$nt9@(&<1`;hioYpb<+X5`(n}m_*S|tu{Zswy<ii+f
zlT!*t38S9U5_D7~&hzh5a)R}Jk%JZ1++d53+eEA8o9eX0!_x_@?jh7tx7&*Ro1b8K
z-As7eAQ)fPj4rcZzu^r3SQEeJ5ki>C`4>gHZh&OUDJqRpF*gpP0IdPsF})C39*i&*
z$+DNk+i9+z+ccRh``CM>bE7F|w2?OSjkmrfR*{o4m{5YLg!_jXf!`=R7LC_jF>GK!
z^H1G>@=<5zcD@QUd(^9+M{&p64vcXpCX-Vxr#_#<u&>vcjT=5|UN;{4y|g@|s{g~1
z=Tvp$%}?ivw&#WX;0q4CRNGOp&tG$)-)V7Dz;$KQ(?ff^oki%nE~d&So%ja}M77Ki
z-i1BpV@ZH#Ig4f(P&XPLpSd&qI-d%)Q>-YNonJuXm6cZ;Dxs2U*}fYX@ACiu#goXQ
zm@+;{jBOqh!i$R+36(Q<a|6GfUtSFg>6J=H9HsGj?4DeJ5$s#Ld_2xJKgAt$D`T=x
zqs!S(+{~FiLRa7d6-Io5*+ZL5<h2B#XoO2X-L{4#sWYZ)7Pwc}a(3qnKU5@c8}3$D
zG!M8vE!+Iey7q<Qj(`6Y70L+=`!L;`m=w_nE%<EwUj@tBUfZ8Gu}E%BM)K(mlehVo
z#1$K%-2Ai9?RQ005@6FSHxFWF6tx%qEF2R9#x*Fk@~3>)BqI%Ef$88Mrmd?&$Z_Nx
z;`Y(P?NS_yamrN49;tE%q}g(rs+0;%8NcxO&Lfx&AvG+_$A09oFz&7AyBb<M)&h0o
zQuRVvm=1VJ<&&ufCj~GIz%#lK@NoI{4VUyn6uh)Mn9QF-zsW$KXScS@fEP6LU}y$%
zfJAoiPR5$G%e41S9lECfB&fuoIN2LV@A${Yh;jOc;giu$RCrQV$Zy}~S;v#}DP>~u
z--iTC?nSLtZ^{A%-iIr~06{Q$b1xK5B%74Hw(o1^=o%wh7TQBh8_<$9F@OUL2q_GN
zX_{oXq#KPl2CQaIOEWtw6u#w-0q^|`@USTdJxEfCaiIr!>EgbGqDTJP4=>*LeXF{<
zhFNIchnf!Q@9K0%M5r=|P;O=$Z|ZM+0Tm5QmPX?HlD{GWVAwB=xfd;taa{qvTL(Y=
zTyBp2Rk+Z1P)|B_u3BD$qDxifU^AAwu|MQwOb^e@=h}(ymhar**h9{4?e!MtDUZfw
zLW0a!SOgBN+n)*!23TC3vsRzhny?ae%<x<H+tT8F5iB|K(dUU_1E!mQ<MUe0dNy4}
zW5u1`%~?IIJEsNLUN1K{kI-?J=lQjDa>{A8rJ*=-K&l-CdRWa}mvkr|i=P&L^7O`_
zef3mb6l~j<iZ}`-@XgTf0Y5xV49WRHcUfT2c>ZK~^Y|B9%FL3!q6D-`rRAOsOZIYp
zE#*ag<`=K+;asg|0!tAM3I&*`>%t%N=V?3#Pp!+n>%$1f&yKa>yBA<)hM|Mhzl&e5
z;uS4xCpHl%gK7e4s|z3q2a{cb<Y@Q9Tr!A496VdCpqEhAzS_88#Ih(5aFScSu{kI#
z0yih;#q}0o-k~|%_nqlk)~#Ji70arT=I&Ui7}+oMa*pz1p=qR#2m^!&xLX8?hId_M
zw5C}&n$I4{lD;G8*zS16JyI^CLSra962<S^<iSt~j84j?C`~%z4p~whK*RS(FFeuj
z(iF#+I23enh7PT!MM0eNZxg!a=%<NhJ1tJJ5spX7S6W-=%qJ}33%7Z5m0|c7^ts(=
zXI<|GmT<-Z9W8M!7=nx>p~Q*=iI!4tr^e*N9$EBKkv2*{NewO4RCKh~H|=7R=oGA=
z&_D8*ZFRL4p?fNRL=R8*B^m#!ms+wp@ZbJ}QtcpjOv&eGI?EmN9aOdKq0~m%V(;dc
zHD1N-R@$5H;85fZ+kNm-Gn}ig87f(#g&;exU0=LjaPwI$MW=ISedOKE2Np2L%t!*+
zs6imL)#Q!TRk6X9y@YOu=L~1qW|f$0Y;exC0xkFJrVX%J^hP(@!Rx1IcYJFv-e(Qn
zHJaEt*_;rmRGZ}aWTgWGNA#_}uVdzxwRHB-@a{>a5g#41X_BhgHx`1$Sr3JDleCdi
zlKmi=9Ur8wk|PqCb<zwkPw)ienOTv2N|gCzv3DUs<{ATLXQBC4iaI!xhy`)YW*RgQ
zB2|m;hc`>PC9mH)aL3M3V5%f+FA;wR28_V!LS*CW!9lZ+dPk!*`&N><-$lg8U{oR4
zB|f>9_#QUtzw`H`%Gh)^`RlTP##_7gs*)B4V=B%XRAN(WY-TA#w}Zb?A89hzB*JHq
zrxEMJJENa<#8&l2hT|zRBmJXV^bL*Tjn71GxHRzxY;JLW7!eVAk4C89yXDf+)I~)K
z#m4OFjMK;X47sY<Mi&Qp4aDEKqhrfEyPGA(n%L<?v|#j_!GEJ|Wrah1^yR9bJ67R4
zv(iZ?V1pZePN`KPIW)~P0T;fPRkgk+EqPfJ7IIl0*??AM-dwSJnSIfyVNZ&T0q6h|
z@Gw-dsHg?v`~P$#cK~z9N4t;~Y>niB!ZM4coG9{OA7S47GuC|6c{=`<TZf=R`DA}k
zLQ)r6@nG?N_+!PA$+q#-uzm3_wp9TnQ*~j*#arWK<90zV@%$dUi^tGnf=pf(F9Go6
zI9oz$?O7ayqXRH{*P&z_r_}0#<Y^w1S^X17P76gM>kfmMi^m;TPaDdb?ql^ONtPcR
zEFS7=pVrh%qs?z|CKFuHGbrv;qm8{J5|)PsgcNrcH^W(;UlQ8P`SsKN7d+imt)rEj
zi~r3{9dDFcLv*_dl_$!q{Pe{&5iOPnjwa=k#U(0s7KEyEqjON7+emM{1<$}THpLw_
z7P4a}58}5dP^A#uxUpCSnc&x<DCQ=@<3`X1j!xGR5{gaQc=IksBmrih7Y>f*YzCuW
z1SK6Gl=89Ra%?-e&h5y-KNyiwrSf_Bn2T97Zh1)cT35COO1t{g<~;<FAU=C+@%!a@
z%1se?em(`2F~1{<TUrgaM5_23;IAMpceP*G-x@L3lbkY*lINzotUmWzkDr~PJ*0aa
zOvC%Y+EQz8fFV5)O+iM|gM`sFrZDmg_))b%g|4_kju~(R1E0VLmwXh;F>0CZ(3UKC
ziZ-3(>xgds0jB{%B>3a`>nGEhWo2cr$NUre`w2ekoS&e)N{e+Tp6Pdeini-jZMWv~
zl1R(CkGIk6$u$O~maab)(ATA^tMNr8<7}x3!;PBeRJuA-ozU|A>4#F(Q|jZr`CE0<
zb<<l}Blvl0F)j+ofYHuor$g}WD=T(!A9F=9EwyqJE%T!=4VU78wcaBh{SRho`{#rF
zgN88U%*}oQVJ&Hx!Bg^#@l*ct?(q0^E^~xt=<utP;}5A>=nT%RSRNyNFkQ{#*!L=2
zW4ClJ-W1oIi3^ZqGV~bL7lL4EY8?V5kXjj#jOt#~M^L>l_H8CPHtJE6c{p~GM?7_P
zzUCz%szEvV@RQXrh7F1)(M%yJ95JRz(ZKogT2iCh=<Z6ao`)^kXk-bSCzYdYU(&dZ
zoTdzzCqfogsBg_AoP4_<y1LC=-OlnQj(+_5s_sKm<e1K;uEA)hr(MfLmo+dxBH;tf
zMK7L_I&6N=(+Ku(O}DoTq8a&`|LEcKVzyb;e5yBkoY&pd?|sZIHVjB{_9X9Jxt?tT
zB=~WE2(lZvN;)}t5h9k2?HaWC=`G7kl)K|-r*O>7p?~P!q4FvwkbUXmqAdW-lAVkb
zqrI!1d@9S);DSSiU6XuAoYFZ#5A|H+F*vQ6`bEHJL^yQcEr{4?qAs@4YI<5`zU;Mr
z2>-^HfZRS8O9W-(vsqr?tQjpY1oH}IdPQA;A}`?6s+6jml%ZppF8L;J3~U*Id8BSz
z*XijRwNS2VhW-lGdX-r3gf00~9}C)c#O6aLRA78l$rjs5iItjj`}E2XUFr*Pder7x
zxPNvv_kO1ArtT5rC%6HqR>Ko};rHCbeB}a)E_pDT_cS+q(7qi*ydwksk%pzM;e<9h
zhql{<`BiK!L6Spiv$X=7CSE~71B?~_R%prwac97fU*b-SEjNKA5yd2}8&v-!0!6xV
zzfgIHMsE@iPo2MyOSeA0H!|%-5;~!LW7XH#cW`6j$ERrI;PwmVzbrdVyV7Rn=0do!
zs+f|+I{#}q&=52L7^99ZCyixSMkT=y;2hQ)N6jS};X8!A=~MiTTFxzIyT@#yDW&7U
zI|}{r$ECFU&vo%+Dn#vhM<ZMR`Cnjj#NV>K+%VkNlR4INpPly1I%H_xhP!qYTO30J
zCw+owKChU5R(v0DZyggJrC=yA1c=&w`9wwMUX~(>sKe-BYFGkU?6+e>_G1OjS*X)w
z4kr;F*>BHZ=X3!Nm=DkX?5+nwEs9ZWgF(PjMi;j?lJ=QWOZtrc-=G<r`4l3(AFcK4
zBl9Kq+}WVdjURI_n1fHB{@p0Qa!l^v8~`dBJgPvL&F;B4!`BDQC>^*5|LTB}&dZ<9
z5kQR}U@s==WjmDEj9R86{meD4tR~ev{h+W{V63~r-KJ^WB{@j8Feqsf2`JEH^X~_z
zE%O&rl1Ja=Pgg9j>?9Ar-S_>w-@X6CCIR%;BH|_^z?e7^vwPuAw&&mS)xKKCW8{-A
zgpEcUd%6GyWw-~6Po(W0ABU}EeoMGR0o(juCg*2f2;0=CA{I`+<SX*vLm_R*{HXxo
z?KCWrK-z%Nuy;P$%2Xb5xk=ANwH^iH2!v5H3kyfjyU+bS_zO_cNb~uVR2dfpVK@Zn
zGmynn1k6as>3n7#C0}GJM{4`Y1yL(@b7=feZ&SBGnp)tYScvDWJdzqiFcY)Wi$Bzj
zERS9Sqak`{ZqN$&Xtg|H-;acnV!vSlxd57Bk6p;@RzIFQR=Qj!=dtmYz_x0x3I4I|
zQihqR?wCH0#GENGm1nzFfEtu`wS0_tqg(nuv!tlEcGb@aGvOr#g)h06U0M^BC;5Qi
z*miF7*!xBzXJMegXjut(tmr3?3Rl`EX62l^^=0K@C#@l&o12%1F475mzr_Mfa9PdM
z-_@JXA)YsyEiNu1zYc4gf{|%ykJ!UhE;vH#&&re6|I9vo)uAS-R3Fg^ItGc2%Izbs
zpu5boZ2rjv9<HwdvE+{Vqg!uKES=@|DYCb-&6g;ue$1T+F!3VZyHMTzbJKi*z&|AU
zR$u@jj7c0J$ip(rS21FP=((o<+-=|qG2dXCD|e-OTd2U0wEn)vP6HK=w0{f13pvUH
zf7~#vZimk#(cVSA*+nJ?{<Ml@KtdD2@M0)Fi|p|d`2bGrpJHt|(p>4y+bvt^_v$I!
z!Ua1j)D2pkz*;KX%A7MfElO}4#QfQ{?GQ{%U##nA&#^-j;^I3f`A`_n6r;&wB>Nya
z=;17CCz#5vc%wzFjEHJIU}V`ZMKbN-`^^Zrwzak34{7lzHmmCB@0(a4Hc!`fDH{q9
zZm#EMI^acL2?UE)>sA1*&F;JB#A6((8*hZ8Rce>SA$4)NKowY+&!5^yZ9EBf<~kJ6
z;4Zx5xB700_Dz6q$t0)06gHSVB=@=%h(#;UCUJE*$lZi}<n4aT%c74DD%Jy-k2f>`
zcln3PajsAi@>`|I&)yGfM$@I<cR{G#W@nu`1DO_n<C)?IB)@$6?q`NV<QTtsq1SAQ
z$%JFVCK^oV!Hht6Z%5rk#ex(;AUpV<k&SB{xXkt^veCRwx|E7K&6jQ_iGIWGg<eeS
zRnZW}w+M@Y)wYcJt;{j}Wo=jTN?o=dUiQU+{;!AC=2au*U{nEK4Erv!4md+hoSYzR
zb=x&?`nLUb1K8^#!qxbEe#PJ2&0^4&5U;E%CiZ8?oGPMflr1XDWBiR`y*z-FpsM_o
z7^~(+LEZc5Y;gFzuoTqQ34LPO!G#NHz^ZpI#17lI@+B+v`5_Ni?X$N2SmgDXiJ_Xc
z@(TX<C>H7d(i~6;*ex$FcXxMhKMel|$s@lhOGFDy^B*3!s7-YA%+6^s;^$1fw9K#k
zUN%nnrYghwQd%%`LMF@RLoyms*w(V0tXmf9toLsRG4FGq2>U6msMKqz?oVS;Ud4Kj
zv?$u-7$aLv4b#%l>wOtubU2a80CtH;Cc>24k8O>MOxb8v)MiyDlGpsDa-EIa-~PZV
z<_)8T#rn98r0?*ufVq3ENotl$s#59_%&VZnHqr0Sv7qtN^B7RFIs%!E1c;5Ot&HpE
zcqboJn9+eF6lzh79r1Jv=Z_9$KCrD!Vo!f9glo+@Zf>e3|1mMfer$(9oMGU>?V&TJ
zR0PUR+2$9K0h5tZuV%lx{GMhaiO(T$J7Y3ViZd-_OR3lNV;stZ7i5iX=0_aHz@|9g
z-?a^2*BGg#EGdD!G~Ltvr*JIi$d5zVH=l51{hqFi6C~nQtM8*vgLd<};^y5_{5?jD
z3Y_a>FK~&bO6s*Gr>W-(Nba#L_jDy^><)}+rsOLN)y6z$LO)%bi=&&fz?2vd8#pdR
z4n!<L<C<b_8yw}bu7!S`j8@EZ>BZggD>8LpoPlNiI`r%oSkuQa{@SWX@o6eJEWKXE
zpsp;q{rK@tKk>sUmk-vaZld3TQQ?vF@Mkz0mU0@B#Hcz;{>Y_wHPA%+@o~ReJG!~q
z*-ByT$!tTDUNO5~1q6XjI)CDhK>TN=1$N<&OZx6z%3C)KV2O?KofC!VPpBGLYb(%6
z8%6$v$N%JethBfyr8qP!-=YG9K@ZoYPkQzG&E5K}yR1DPjh*H`v#^slh4~Wxjd^w2
z4}zY|sDI>k6q$v_LAzC(eUR(cdWyfHnex`P-`p7>ZT>C4N3C9|7qrRW#+~|uRrIAX
zJELZfoejZTDyCBPB4E4;Z-igH9Vv`~;35$}9PiiyETHq-h-UEI_tUpyI4FB4M=Fy3
z_7nLef-GWTuu}OGs-Y}4xV81r)+qI)_Lk;N;nOQht^$`{aJoo#dgs`CbMm&9X7k5-
z<at?UaAW&08t7e1S{R|MLOEfg%gV)Pk5{SCJgZ*C-w?q}P5Rr?>5>B0<MMw`0;R#_
zDF7Ldw0!<uaRC<3Kxj$yCmAP1(K($DH-eWS{M#5mBjF4J^wt3I^1*Q_{k|M0O(b7p
z!0VfZ86Omw@~5eBG}VY9N|F?R1Af)NQZ-6}{qQ+1xLNRh_e5`qo=l$z5Buoiax8qV
zIey$+dJqB`lE1DGp>I)Q!s5Gi)p7f-?|@S}Q|g~|bXH!z3&*C)GK=|@<cIW;AFqbA
zCkT}(6+zH+F+qFaC$IbmMXGhs^-kXiuIy;7Zk+k8OI-ia69jM_0~ckBMYr50pb5E$
z##t;3{?Ko*qznVz*Vlvda@Pshw8H~8H}__Dzy9V5A>Z_A#$IddtlC~4dhr3BI)itb
z!qO5uSy)2JhsUB1eC>IRRJs)k_vMC*!5x}ijWycPEg#`*&4PV*-Sr=1erBkdecsI?
zw@E{li0eP|$dl702t=N+F9LwZm+Q~F@_S-ePxRzJ5=PIle`#z{nL0mvvto|lNPI6f
z8gP-)9mu~6FPNKZ!L1<`hObnUw<LOt<&9gk^ilIgF9LV+wohS-q@A#&(sq1l$X$+7
z{*9C2YuD>Ite*E?@82;`XW+DHq77)uV39H38$m?^8#JY2LA&Ye&+iQ<@A?8x72A@&
z9Ah_XNAtH`(MCg*G*=xQo7;OwfBgq63FSkJ>hgGhwW@=2cXPLpboRm~4|m<Tn7wS%
z+i^wvw;vNj_kw_tJ{^cf3yN!8ailbzc9I9hFj|@U4d#zB8pi>SR*(xoR{*|k(9z7G
zzVJLE;^4AE-Lg25plv<N5OmpT_TIb~Zek#V7L<c0+ituC)CQQ>Nh00dJk&SdYE~J@
z0@8MVDGd%R1sEq_8kb^6UH*b^l<&0*;%<GvzPZ49KoRUMJWtMRUk88oDW|BU2O2G@
zqZ5>Bdw5^w#g&qDcn@NUg?;0gG`PRM?IM@i(bwx-Z~Iap!#4<dO``qxey+{k;hIW7
zwDSE;60ZPid80JjRX>?#O5_>{L3l2VR;b>K31#V}K8)AX`%V&<+v}-(0!!MP{!IL*
zTRLLl(OlTo?zJ>qfB3b$%*$h$i)U@Nx<7Ye3o~UVX%;qBEK<Ao$s(TGyOQtEqh@Gq
zjd!`QNuij=+X9x{U8QTUkCscLqNL&xmihgxHtad~8d?C@=LG}=7Pn>lucV;g;i>G+
zCl(FYE3ISakBUr3q$$|!0B}UdWJuf2(i??p+F=CY`$S7y5h1<juY@@h;6b|VlbuFy
zU+wz_9yr1Kg_7tA<GXE)-j1hF8jKYcvqQ&k3_DB(m#uCSr7F!K)OIwIW?pGVYA`F%
zc)`BPpJN`#b04?3?;Mwk%V4!6l)=>2B89Mfb<9)Qh1=-2kVCAEfeO$^1LqAr4itat
zpZ?ugPLY-R+1}DOq&U$cmuE$m`pQRjWUeS#>EES7yGN8oOrex8r8l2GR?!opbZ_hB
z*Yp~a;((e3MxoM^IG`Ch9Z3_s`h`Dy@pqN1<tTq|{;&PZkd<-~p}>s_XgSO{%wO1u
zCg~VvSKtCVuMHp-*RX3%#Ai12YL#k?mhc%`c0rbnR;l09<bK(I$Qie*EiWe#fBNA+
zoXLbG&2`JfpCw$2TTRQ--e_UM17_N{*9$T%forUA0~IgYyLQDFbR75_9vcQ#V$mPa
zFauZeiBN5GkQRSp@K}1e&S0w>U>ZrusP|s+x8m22*i>|c8(z(tR9c@Q>Q`N@p--oy
zI4Rl$L%DA_(<f(N8vkAK2@$h0Rz0!OUaz)%rRQ|lYwXBL{SRm2t(jzEim=gu<mLb9
zF@Q0C4nx_{zi7NhR<pEh6PlZ{eSE_A9!ILAlMv%hfBicu^Xe@B9xEM5eLH@|WxZYQ
zNu2Ai#p>*?^YUv1htCYsTW6n+gQnqD?->7?CzjRf{lv+PmskYjz*gyS?;a$v{z+*Q
z>xUR#7C#hJc<5tar=Sy6>0m^7cgZp$+&4YLJ6?h8TZXCKt?nZZ9f#DHg?x&-^~*zw
zCL0D5$9A<Q3-OP33Re)?Vt~qT=W85&<c`T+7HfQfa(vZg+>I&%SmK;(ma0cjsQd8C
z?kR~v{XEH1Q32_(zdPaE1$2~v0IB3&+u$TKba@4;f|v^4(s$$xD6!G258aL_U-jEx
z?G({4v$V94&YqMCbQA*cl5eYI@@pbb!dN`W-F?^yow|XgEc~D}hRt)9-EaUkmG_I(
zoJ%tQ4oT(JS$5>{MEw`H55EMCUH{h$U>Y*F`SkGKHHVv(aLON%Q|$uBvLptjd?$LU
z<TNygfpRBwQAy@p=kK?nzeGE&*TF~rQh&QF|KehkLBW0fk*s&@QM1wM%76f6#E3Zv
z^_HP~{*uDt6KEr6v@7Dz(<IHIdjpFye+56do^ujijfRAut%k(!?KCa4tnd?9oS>6~
zEdUEK557Ik5`#35c}LcceoBk+%euJixb3!_EGh*rEA4TG9B}<>d@(cX(|pTwj+ubn
z7T;fC_<L>Gw|`xH-o4mMN^8zrPm`q+W~17<Qcwc6GC5Mn4?HG}xlA(J*?*C{35=-0
zjZ-fPxIcF##&F%oz*~9l;pg|WXE&`ITnoc;zV)DB)M1+#mnIe$0xrC=1xsWH61JL@
zf|?A(07pjC)Xo>!{Z%Q?7r96DmMN(ge!$3wm-0$8sH09?B$?{{^=|86+Z!8<$`|R+
zU7Sr#ex{b$r;rH?ZcR{P4x&tuKIgtV(z>SH4L)hww~_f{Lmc<cYvS*C)*8;@r+Xu%
zo!JAcp7|~>QIZSzE>!PtLA|yO7EZw+Q=F<`@H|bgP@+tf?BStZNs>sAi(%<@xz@z1
zA&ZHNJO4CGi?cq%#qsjs2XhB+>mfg+Nm((Lz>%a6kJ(*_J`O7v{leC^3D4ZJ#Bxr6
zjkYU|E=Po0pQP2vEaWY%Z1Z^)zRg=7l1>Ff2OEp2RWKEI7M(u3pK9<&MP@cD>`k4*
z9;<FtR?wfHpCVIG>*iS_=Y}93#YD66IB(R{M1+GoqdZ2aQEOpJoVDu-*WByrR^p}i
z&q$lk$M@#`+AG}i*BQOb+a*-Y4n%{}#-+EP(N|gg)UjXn?20fb@>W5!7F-Ff7<u!l
z$bAJNY<O%m&RU!cuCT;0Q#>>p_phQ*q`fQ-d36N}`%6}<^tGV#QOq~B`(0m<&3bB$
z-sLl&fcDx~hk@e~atl7~uC`i9GkjE=D`5Pu!)P?~IN^n2f}-o;dMKcUO5WJ94UR|8
z)l1M5zEvDQ*&;oR=K6kJeoX~#0?nxZjQT?Q_4*=&_VZz;(REYE<Q~B$YdASDxH2ig
zH`PY+oytgdm0)0aOVnZ2N0l0bMNk_NBi12-{Sle@t%<sXaa3D-HNpn97Ics-D0w^a
zL)p{xw&yReKyf@CfF%l$@iG)lB-K1kRWLM18`zhI@};SPyMP8e%TER{eU_>P3oKKA
zQmD*Uqa!!-gT91;hfQY{mtnUE<#K14!l&r5YTg(L8eaz8-|?mzx^@g+StXl$fP@1B
z;`}qrn>ymc?!AU-{i$}OLu=+ogew+N(q`Q~%?P15J-KH3Y8|G${Q5$=;%D1+KuEcs
zTu<e0flDzvT-@^ZQ+e%Q=jZ;e_<TrDAW;}R_pdMFgGnt6@e1D)=7hz=es+V|{*SsK
zLE_R4cYVc`%9{)QmL$8(SMo%2moD$Lq+St-4K+3a;#lem_xGK0=!f1cy61{(E1I!y
z_K%wKuP24WUU_Ek)E^r~i_{1brMhzak}v}KL+e^3+^~i4@^%RA*Kw-w-M^Z?w`GWC
zyiG^tN3*~J(VJ6*$`*vwI`2W^)GbHXfUOkHF6meN{ZE-uzpxKK*ngNhqwV(Z>iXlx
zh?GnlT7u6f0c3%@C!$hGn&OC51rpW4%gel05VGUr2Y`8G&Q;SOGrS07H4)B<eu+Xw
zz9sRg6WVU#3|U@geZJX)<8M!1r5jE{Rni701#pZ(=Q%;@f%_5#s;!Bf`ka3Mf{c%~
zDQKLa6p?xfTBJB&B$KsTd770{@Z~|#I*i=B$fMLbsSt}5N&&LaKp0m|u|Tt8h)4LI
z1ScbGTwD7U)tfv{1BKUI*`2>yb`^IMl|tJgSM+m`-2qqo?I*WVD^=NoHEFyS5cX1S
zroej`u9=VV9_E}2Jn2?vFP?jBkguO@!+8Z##JLzc|7gj{^2Kq&wa~Yr5T;>YA_cOD
zU4I(?kRv6x`Ae5j9$F>66;Fe)O7L$eI0eqTZbPG8b{86}V!m7Mh0}-#(mCG8V;Lia
zO1I<12RL%l=ouQ@mibq_^-`CY%Y48{`$y)8^_}JEjW0BxNl7Wb58$;71>#*M1mav}
z%#Ho{QFtfPKqYV<9|9~OR!ypGuGrrg>yMHUm&Sb5z|K1Y>N>nvOmrOkO|yHdxougp
ze0M=vL^vu-#sNiVv#Xc=cx<=39Uo#Pj)Np<xxQ?@6G|F-2_pxz$=IM-l?|U0;_aWb
z)}V^?Pj@vwuvqnI1k<q8cyJEGQ4~G$OOD9oMiO3@0#VAmQAHXW!Ei{xGe=6+bj=J}
zavn`H`S?<f1!%MdXw+c87&B|~I|vEhWSJ!M3A?hDD!;Nx{qy<3Pp~OAH=MJfO=a|2
zJDqu@s_pc1p;uheSbWse{t4XS`&-qacC5a8?W!SSwf(s!<e3I4`oY?IpFnpBe8JA(
zmkRhILTBB_Qw>rSyE2O2<%nRm+J@~04E%cs;+m427|%QnxHRC9o*zFQvc`JE6hMSY
ze&64BAhIG%9T&sQihGk+CJwCyPWs~Y>0e>wLF@^h|C&5ijK5V!2tr14gpT<Z8?W@r
zhT|*nQZjU7xmkd1di9h5^3S)@)t>gQ5M1rXzq@FZdcN!{*ml0{8hMQ^133h3-r6x~
z;EB)O=3VI=n_Wd?{nOX5?GF)!=#nI_1(Ot?vj55RnpE|6M}L}J$<Mn>8yke75}=u9
zC}tsmIFFX4J9GcA)#ltiz=a&1F|=6R&RPYYhn6GWc})FWj5e_-kNB97@4)7JTucVa
zLVmNTW4*A<Jycmzp<8B&JUGx7q1sfy>FQKS3;YhGyUFFM-8-MKnXT$`Ss>tAsiyZ-
zYCs0ytp9b=GQ`wa>4Ae$T-{%CEr@$V6<(cW7ZSh>RM1E7pAu4!0}n=Rgeq2)nJtO4
z>uZ#HMF!|!L@9wCF%I7j@&9=G%D5`OU~9S?ltvn*TRIO&AG#Z)K|n&f1(ZCbAl)q>
zT_Q+|NOuWHw{+co{_nl-2Y-Ii^E`V$vuD<<Su;CHPrR4Y41}!;8Zdn1F$y!)g6w1^
zI)(6L3%t8{&o}d01ZLv!9*aXpj_ZwdQ1>d4nd(r=9}+KeFsPYuqXLQy5_VQ<_(?vM
z&sRc#vGb;}1hJOg>c18sA@RgjE$4d<O0Phv1-i@Oa~s&n-V^4q046~t1M(ZNc%ihj
z*M<}P76@cwrx@=+c_D@?h&)|owd*AIIPl-Av96YY>h^qkBqFy(gJb0qO^e7<*$PGr
zGSCP-{?E~A9Fzd9dDjHDwMDWuCMi_WX$|E5TGo|)V0J0?5e*W!@0?<Co|4Vx$Ncz%
zKAZ9UHXyw?dc+19Xd^P10X=yqu0fJ_dD|7g)^+`od)g5)>Vx-1-~2fH^4ryr-WQ+S
zg^*0dh^tShmqXCF31)kh1(9d*CLKejOd50p<+8fCaK6V$aIcpq*B69{X@VE5p2}Xl
z(8E84yC`jOJUi3sU&^iM(alfGBQa1QN>d&pu`EiU)k`%Qv@WWt6WmD1gYLPXBWT^j
zrf7>)K|(Y=cUS#{ffw}{pvO<hL|JWt7YDi6{%LmMAfn=o=HKU;c*GN5bzt;}&Gs2e
zrfn4oPwZ+j>2&0lJ|T=D0?Ovej&-#mf~H=p#u)G;R5|5dzX=ibyH-|qoaswGAVh^>
zX@noPatdau-r0ioI?>mNvudT)H5LKa!%^91pAUSz&xGX9gnxG3DW~T=iiN&;yvQ0%
zxe}*GNHhyFWc?^6o#Y2&L?)tU&iB+>Yp;e?*MRH+x?rNlkLK)`VX2Go{Ru%i=-eOA
z5c8>eOP{GT>S`BLUV3T+Q<1Bv&&yT%Q0W_$f`5HaAg%#v&1lbHW>(Cn=@n2_k!K(q
zX8L|jRcwLq9&zH*;@>UVq8ob58ux?5wKujUed^N$P^U&nKA4d5sQ$XN*^7))+q(Ss
zOni?92pr0>K(J2{VZQ12;?)yCc;TlKsX4YwHT+_TqW?t&&7Y&*^2k{M4)W%MDoqvU
zC>mdk`^$E8eRhIhHs5UOs~$mbX087w+gMvm!;v}n?1Ca&%Tmo1^{FK2@X!?<*Hy$j
zXSCnbvT|U#_NT@{BU1|T*TtvFdii!hshdcE76MPSt2E9b7mF;n<xluheJX)FgU@h)
zehOq!O$gHxKG^@P6EF7Og4Tc7hyc{od~>c)Zp<D<@=fvAmb`oJw;rk1>`<Ue3q%nI
ze^&l@M`7ylauc48H=cf2ND^0sH88RU&}{RY_~O^=Dt?%RZMQ0YXAakX0%t_w8i)tT
z<0Q-NDVdJebUJ{x#rBqhnHEA+G~!QFyf8)OQYvWBRZEVSPw>bdw;sRUeJ3R_F{Srk
zMa1d1!U^dJHV4hbIpnn&bm)wytcs@dsTj{;h>E6ee$!f3C!U3Aja{}rOWi+mo@puP
zgMUCGSB4)8!Nt&+gS#tjH##m&Yf}$n)~ev-TEk?g=R%;?;uBCJbe=xhv75;A4pQUK
z6sp(Bd+r(iwQ9`nYm(sRZ6IiDt4YyQ1Z&*L#S*G&Hn~FvMKuw;UhU7hSxT2D6K)=+
zOmNKrE`|<BR2$?!5#L*V0L45#%48oW1WyO4I3WE6+8CRubsYEkWc@0or>Sk^_C}WZ
zjGeFU-m1Cz!C3?riO9MS3!fT`iaFb#ZMDmr=qlD#cLTQsgtdN~T{*2_P=dW`CwH2_
zm}KKX`!<<QhIP)_*Of+R6lwLLT8%-Ten!WSYxET~jm}_<n^6T#&~+Mjb|A&rSq{F~
zWs6QnyRqeXQmSA|YNMxJ3Pf#<2v))14kcRE&1HCesKm_R6j#)~`I%NVQlnbBD^T97
z1e-zMc+KT!u+tR}ii{0?EyVZbP4PNA?*p8c<Wg`^Jq-~PM-#ivzclUAIWPy`>A+RF
z+3UHF)!yB}MOakoM*V#o9<=wS`w{&Sp?U70{e``~5pd%RkHzOGiqo*DGi6Q}sPcl6
z-_ZoK-RaYB_FZPno)sL@Jf2}g6^wu!EeCOPT1WsvxCtz|*+b`Rp80~aASIyY5Zbl|
z6f0=nhdi_`M=XP0=IuC9f{s<xZ}p@3m)-dEyO)G_v+)AwyjsziaqZnwytInN&V$m5
zy7Q^+s9$Dz_SEN7({rhJ6=9!qH3a7Uo$#OJ>6QJ`fwj*JB7~R2>1Hn%EM0Oylk}I4
zYdz3c;CTv=Xun2fbfg;8i3dI0^>35XaRvnN{~TTDio4KAuCtkiTn|EBCaPW9oNqrK
zwi8Fu|9jYagdB%SVu)=CpSK?qlfMrE`*qf<tMUOHJ&*!vEqwjTp`lEI(rv05btgJL
z-q;(ob<TJq#swU>&m~<x{jNR}kn1($a{XsreLTeXb=YOY6nZA7%=bDq7(0w_rep(Q
zW*$?r+?^K<&ht}e*-dA+`NPU&P;u1tJ+oJZ?+{?g?40Odh=Fg?sslo)n)0F|=jfED
zYm?cGr^Hq#9qcZWSxPhY<@_y3nh?N~`dlJ;$&~!J`OihhxY`EVpHsWcs`A6v57GTm
z9?9N%?*&Q1bPeb*dKVZFYbPkns4K^BoAFm|K2@k1gzKYQ-;OxPkDyW}?7XKAg;jp7
zs$}W=b6-)=Y%(>-)-VcKJfE%S|JzOAIEZb#Y~5I$oL|XP&OfL|Q^HGum0htT@!}7r
zzfJ{B-S!VI`Kxcp+xLe@E-R5Au~&Do`uIj^-%Vr&Q@(WY8S#Y0kG>OIF56U4ULb;c
z%SpNKyBL9m7>FzO3c-FHxJ!LEw=^S-i~dzp3n*5MZIc-}FoUN=^4=cEAl7th{wuI^
za&cu+NHhR-uS3dL1V}_O{gjLOshdM<UM_3D{9<z2@TNfpO@wYN4)O!^f${fV_uXWP
zx!4&_{v1ws5B?kPZihtKeiqP}Sm+d;3rV<5j9O~!y-XfUf;`-x3c;*647V3N)J>RD
zXiBw<VrzJuQ+m(DKx-HPOJE4mW#?tb)w$WDC&IUbrj<CLSIg-M!<UGH%M6+?m;F+f
zZ!-Q_K%DQwHEhEzn)#$#{*Nv}Rfi_}5W;zsyIY;rR!&fPd~M)H{nJCFsuGVGjA7~n
zQ~?kUP<*P4M$oDB9jK6eUF~bllkHD*bu-h1y5qfUu#BxGMZUBwcg%tW+tBSFJO~0V
ze8b2D?p=@4uoqE%_RM^mxTGte8>V_rI!=e!(SAy?nMJd-l*m#uqEtYN@~$PW;lq*j
zKYW|h!y<XD*nWTTra#MjIa>0zUQ(9qJs6JgA0Y&miSOOG?<;ry2MfcVgAKLf7yTK%
zcE#I-w8c!9=g(^xC&k}$<J1~ft%AW<KFMM6ntI~A5H5eZ<L{)YXQ3pCl?Os43x1~a
z$v{_80J(Nds$g`Ng3Tl`aycspIuB+ux`SB<38nu&)hK1M{4FQVOJnlM&M1z4b3gs&
z$J?A+<2yj{@OR!!)qsIQYqk9}JtY$7FT@_M)|;!CJRK&fKtja|q96tmu#PEH>p(kW
zRuy>=un;Kle`df`hBHTOfNUvS<!_bmyKMA+w_k?mc6_yo(x6>X)Hp>>Et9C*(sW+e
zFbQ%ZzE$j8^rVSpq7fj;7j<ujTk^l3SK>MQbe&qG8lK|mgBmtSOwdnO!|<miBpg^k
ztQ5EoeX->-V_<EnJ@L}1C2!SLC?;q5Wz`{lfX|ey+AwkVt2`b2mHvHYv(TN^DAQBe
z0jqp}bgr7(PLdv8Ipc@k2_&Ri8Tt}uCO;UJ)?#*QV@iV(Lk6GClcsu>==qFJo0kXw
zbuh+kvA`!5+CvNdqDsc!TL4N;rtG6ll9l4Zz4HFWWD%v)e0k(&7bfP|M(6${#gop#
z3YsK<fRn&S>`UuJK|$osVR&~ejS}>Dge;)iti$l8MwAK<<bLa-8PS)n{1gX8&Vv<<
z^t1}qF!)Kh^An<ax;y$Un{hI{^4|<-I)k?lM3rK<`ZN_z17!%n!JND=agejOUU(y?
zvG98if>%MIVMU?!tPI{be@RFzd%`b<5qdS)K$fHMBSW6B^6@Um!XKc+&g;yOH3DGL
zW4A|-)MIDak;k}iO7{i8%&dGNAwHafX`OhhLg4$%Wl)$-{Wxb527YX|UNV-2I?T-T
zHu#X%7gi!zOXPU}B~$O1Y||p3mKNI~2|(@Da+sPlK^}+5cE8ne^WS)2es8gH=*r})
zI<V4YTC?sJVG(s+5N~-E!A0A?MMD;bsDjW(?Tzpe)%S5Pi0*D7h;n;iL<ByH4H@KF
zWH`yC4I=L?$bf45lP8rs6h?1vbjtei%7;D$%5UnV1qBW8^&T0Omztmu0r&g%s$!lg
zAGfvJd8P%b+cCY}#s2JN>P|tFmP37jwpomZ^l%T1%mVro@8%DBUj!-15JFAy6PP3a
z+M%57lo0yT#AGH-2ohJ&Ab>6<Ltv?fQVGUPir5&i9bZdFS0CHzjp+^_dn0gmPA|F!
zzjF8#6&`FE6SvkO!O}v1fT}he!-@O2uWY1NM{0}DLDibiD6q;&f-zpl%He&78FA;-
z_vqBjgq{2$?fJ@%^ODkBWa;5EcldSe8}oa;M^@cNT!yhsnTV39*NI_#wWp;wL-S>j
z1X}cCcPV<WCL+<H!EP|bGUUIpW+z_l;AMawfkyen*1m$#GDRLm?vLXhrwcegOP(NI
zcL9HNfY73x>f>zIpc@sFZ#GrMPpc%-@Fync1DX1-Rvr#8z>tx@&mQPp7o54-vHP^3
zh@%0VXp+as=_)0D&|$=V6vkGGs?;#P#M9FCJ)u{I^Qj7Hu-in8?5!Vdr<+C2)DuY%
zJP(LeQX}e@tL5?!3p&2<G>SMsuLC&>GrWi}NUHWEXCJu@muR>>@47P{Gz!vU#O?az
zkGrv2vA?eJrTVOYd3*dgC@FoE*@+2%8Jrops}Fj#cMnIeJXQk85L%47nR8NEK7Q8G
zwwd_RNQf9RdyFX6j-1}EL^lv(2O=22)?u^*>@G6`1l36n_;uyU7b2Sy{|o2y$9oNR
zFH}`rp3nZLUx&^IC^W2ld??owzxcU@0l0Yru{-WBrpRe&QbU*2^H{opmf7s?x=hZX
z@k1z%_~p6v#1E5G8&m`8l^?>-LA~TF2VY_}i$sskKgPtcwK`wGEf{Zqcl_Ck517yP
zRaCB0AO6WF|7}0M)X=+p=9eT?@%=34#_A`<zWaWuX}-F<<*!l)t#6^y_F!>-T2bF!
zP8uxT4KF8l>52EMll*cvJpD(3?S+(H&2O+O8bV>b0A9OcxDV+r`;^QJG-8>&i<Jt*
zng+sQpoBmLYAHXFLUYWm7`CdsRpBHV9WNZw*OplTctU?E8J(9Z=%T+NOkniWV$9ss
zk#<_6dT(*q+d3?mIWpJ*{eY71KCe4{YbrmVnt-DKFxEzpmLKCw*F;c-D+gnM3>71y
z8}>hA8k>Ui-CsE5#wwC3TS0TGCTm|ov-fx#G<LhV<;p|#z5Ui@HQSHXPqk#g*N~IB
z*!}Xx`sP<l8%^}tPfwHdVgL+D+!2gyj*=>+w5y*D`L6t6l{1AB<?7F;-id5-1P83F
zd@8;0?a2C8QibS{X=X&oE>pq^r*1XQZQOQ;V&ZU7D_}aVT$1-szzg&Wf1np?#69y4
zBCU@PzzNGV+LTX?IBuwsJNAC$%r(Dsp%o2~FFf2RI-=5gMPEXU8`+aIel!2;G5|<1
zOAd<8XI>2RoXZY+bhLlva68118*hZDC+{oL2R=A7SVCb>%>e_WI?~2DM(Cl?E#hS?
zM~Rg!3YcT9Wqnu$jeB=T`z|7F{dCXW-4@U7Y`^(0o2bK1N<stxdVap5k@mP;0*9<?
zH%%Xl$AAhiGGg2t&-sy3PQyP&O_60c)LLAH7#|VjPQ89PDb*|`s)MKpTNjLknZ!6~
zVI6-iihW_pFSrECR8Onhz?ot_QIoAI02VF8j4-&jY9crlr!~$~K{NQ-t7Y2iBBdtx
zs1@U=^L84)%MJNze#YGYqCS`zQ6WRtlhJ<Uv{&mCbn2{JDyy<M5ZCg2%ab`{rpjz&
zLT9oUSA@{ItGKmcy#E0zNwk{iKW232lET%z!K4eIQUlIyQlu4)*7GHawmnjU|LBHQ
z;-XFSc6nrbPs+}%g81pZ1g*HROBi?%;6_NuKStY=z5h(M%vgi-skrRNkUHFVv<%G$
zUVi>_VI2@oRal5BXF+Cj!Ki5DrXi;%fe(6s!3xVu1LT-kX?@m!Vq17)@_p!*`nk{(
zgY$$ie7ZN0728Z)A!B=V?o{St%Yk8Ah;d=tMqLdTGx{vT0W3JOfT9H~WOR<w!ecp@
z4=WT`a6hE(9)a}DR0^iDpoGRG<qHMG>w0eQS$o)z1A8NakGVosEYTPW8DfX5rj!*4
z=OS_=6F`za!jNALy3Ch+zLIQyl01Inz0z>>nwSX#??WzoekiTG-ap2%B=Xp#==Cw|
zsP5N`Fm6r;B}<YfGOOd#XdFWCKVDPY5ZgQZW^7Rm`nqUoX$)3cTJNGX-=Ls`CjYk9
z@}>n79h3X}&q+yd+n-O#P25YI1D18*?c%#HjYGlKj(O(Ai|#qr!V0#wEUj&ANjBm!
z0-J*uzIG*_L${mVj%!SDwY3o>p8cLia~pc`n_1<{<qr|YF>xg??gp#78p;7#QUsK+
z?$0)~S*vH)>c+-YO<eu#qBiUcu%L@Pbp@^?&07uLIYs)%rR+n-K_;%Q6~sNF%Em@P
z6ePstVRS)=0cg__>jy1B1VJWq`p?rx_S|4N+~3z+e|B(R_>$#(D2M(dK0JTUA?GpA
zAC&}s{}_fJtbj7VQz9$6$Np_y;8)#e)x&t6wpTXj6WaaXS{po=sdvTQ%f}WmbrC}h
z=AOE4jG8`(SC*Fke(X^3h7CO06pJdl&n+INj>?yOZE7m0s~bN<BohXy{ib+55&QY@
zD91*K$O1O1pRc$Y?6sD=nKaZ-oG8!$E$sI%Vf1<7t&Ky}78ORKv63e#DJezmlFUmy
zUp~b>dLnzCSc+T@RNMsy8PDkoQOSFId#e=q`uVldw3lWAJ$meGqn1U*d^F~{hx)S{
zKROMxLcJQ(?4bJ%0CqPpcP-{|(H+?m6;M|E&v1kHIY{;~S&p<AiKM|#JV10L4D?EN
z(m6;}g!eSQ#S~ciwQ6&Zk<F^~J+>l}auv&S3?>=iTc!50l8sEM5edlJvM*?6F8G};
zzuogz;QWG^^%>M(;JY2g{^B?PhyXTfmI%pEW#s%GT~4))_?|}aB+E?;6?Nv$m6<13
zgw~||TTWB{v{H}J@^}UfZU;eUwmCJTin=Pu$bl|#*QpPSBBw!sUSmT;Ey_oCP}K8@
z%S0mZ%%gSP7;XmJx>h0bJ^GdeFLG>b%wqmBe>$)qg~r6>TVHIH{~b@h#51M$_^$?M
zYmc{Fi9IDW>uJ=)#-O%$HMQMyzv?*1dy$r%_($J+43EQ53!;6aetLm|Py{1#q;|#v
z@c|<u0dazU$0e3MVcrxTL6g5;B79&I)YS?z)h*AQpTct~?UTvnpwK}DG)JO1GUC^-
zU$f&`hb~cl@=^Gb;Wf=<(p}FUR5NoA=8zB2B@y&_k-ueOpN-qQHd@`iVI8^huoim1
z7zDE?Fd3lt|3582=5tw#Y1w2umgrh1wk)P-(*;G#l8Xy3PP)R+KwOkvAfe&?hrp&m
zUHw9|IU)=Cd8CfUVJ^ZS%Yl9Ag+UWu(cUzjh{tm!-%E+pqlN~$FTMi31BnDZkqf63
zi$k;hs}~VIHn_jMVjf=?__*)VWH}v4Di%dPPm#RU9%{5+O>~$`^`&-@s(i@47SyMs
z{u!L0wRo#4mrTckoB+W_U4nf|bPvdxtc-M<Pph4gO|=}z*(ri7G)gC@$|NVR9!uiQ
zJ=F1%czPp?dLyAy1;5LSr<++1nzPHM5RhWy;)=&4di$f?X?S<k6jzBT#rbtz62vr&
zhZZq6ZLe&iaZyB{PD@vH`)9B^Cw#lrC8egq8cZ}zo^KZ;K!)TH@@>RJyK_y6k`GPe
z>&yv4U__ANq1+d3h!Dn5b?=w-&d%8jiit5)ALpm0P;*AN$|T-9#XY91SA+Q_%V>5j
ze~S?~#&QDcVB%J(g(aejlSLmlH#cttKzZSwVxwl{^6v$NoItHR-x$8#R^AA(T^Z@7
z0Ab;*^Dh~$s8Ttu5rhoE@e-k3^~4V!J}`8-J_hDOpGB?p&><lL!2$8eh?hF@S1_hK
zIVHtMWLpVta9bP4Ydy!;{?%N{Z$B*X+1ZFg29+*iq%jk&)rfq!p_7Z%HcA^o3<hu8
z7W7eRsi~=jZng`ko9S6e?zEXlYuOMKz4LN2mA}-{85!U=2<v}#$LhT}IKH>?Yr#ht
z8X~Bwk{df-d(O*WG1L={W4&+Z()f`_s=h#_FP2Ev_ydIA_c7tF!QCy35*%hZ;9FSE
zHnq?=GzLxU6dj|zVCgB&pU(G~99^2`7qjDY#3S*jEz!dNJ9D8~O!{?B(t#l*t?Iv>
z5u+w29AQw+`jqpimUt$nsO%sK{{s2$CP&iMOWDt^NWGXW2tze|qREEhFGeDnd{`TV
z(<d_j3|5jNh*Q#stynXdpEWo!YxSX~xXiwqB({m&?CZL<A?Ni}q$e<IEJ9<<Q-5wT
zr<)&b=DwSm=?&9=Ej5cQ(*k$oENO^QPSc?LWY0z7;VY7+Ay_TVAvvr&&7kIXk{GSW
z@T%C~WH2LhshttZl%xIjB)`GhXk`SisO%!~UGzNhdaB0A95epEK1Ihe&H;8-R+;~t
z=u3ZzoWS@>FBsWy<Hl&^P-6`4;jX}dQ*>4kbF^2Mhld9@vWzh3?icIcN9NB&i~<|6
zb(U9G$Gc8M)`Uv?X|V*dX}g<Nb+e(alyN?QSb?Tj!Q`X%vav#x4Jia+Xn0t-V-Ug`
z=}<3Qrt_9BA&9(AyEJGEtBR$MRo3Fd+wAoA{y7_)vVJ2|nBn<KZARnkJKMhC-byFs
z`yiGt!1-0VN<@^J58&T9q0;t6gnml>)_hKzbQpjYA+x=`jX+=&X!GeDBQit^fg-_o
znxVfkTXRSPiG^~ZrJi2E-NKYP_HS8C@JJ+%0k3f<*=S{pI8k(NMTK&Z9X+xUtsw2l
zg&NVI6H9bP7B?E>+>O|mp<oO;{ZKh!Fa}oa1tntCluUQFanjvUt8}ZGrNeYVb2HhK
zCr@~Yt^W9mU+jkjt;=6MJa|Kp5k2t9LkFS=pC3*1dlo@<wcIQ0m0Mh#;^%_CXFW-E
zujS?%!_c>qa&U6-a-2)`A05hq;8>58;a>A5sA@O9{+UpoL)!O{=ybry!?rMVt_B0n
z&qnHb3D6g?q^Fqt+sMta%QJ`8j5#AJu+AiTO<Rr9pggQRhF#q_hO@E}9vsTgRhwXV
zcuM<Uw(;O@xV*NtYV+0q=MWGmNahEo;9#=PU9CVwUo5IOGS$oCSXVJQ3_vIKvz?x)
zasm^m>DCuX-%*1>#A$&Kss=u2bO=0iSEgH~nFB#lZr&h1s9pXsl_TnKF04oPJ;I|M
z1{YSu>g@(GfHg=wb_rHBFSJ_hoL+K<vVFa|({8cw)v7z`>dPXrQ9gC&Ae{ph3h&%P
zqGJ|pg{c63sk@vYs5f{nV$_Nv{jK3fjjAn6^oeaoF~}G`ufL@B(PSnk@Qjzsj|dgt
zOhsTckV^KNUVp}~Dr>zUXt%j}{J4${W-%iF_;Gx&j7ueiuwj8<IXThi)x;J07Q4wN
z@@F=T@>XG=>apSi&&@KU!RbIsu_%6Ee^fFli*YbxVWe7mf`2%R@pQgQ;&y5@Am12B
zUb25D5mkDSV%h9bEPdOTu+MlO7l(KM!lB#%pDtYDd)?P|s-mo{<BCYR>%!3j?0@6V
zPHC43=huHXZWU9BTEs<@)t)H|_CHe8+p)LsmMoVTTaS(3pEV7^&I_H`G<MQKq(84W
z%KA0Cy!`8uaPzaD<VdtHZ3+V81XyPa^Yv}3M*eN0(zL9^2p2f0gZ3qP+Z~s*c+^NY
z|Agy)_;0D%^;049y!)cXHKuD5d(dYu>i#=~ox@%0SxD04rB-JuQISz&lH#*OcJkLx
zIe1tOwelr{>LqyBS00U%oy%r&Q=JA8DdAc$e$}2Y7Cpbaf`U5J{xdY&Qt4N92!$*F
zqqzqJw;>Vuzysb%U}45NLNe?<Aww>;7{*Fp>@f)K2*1&Af=erqL#5G8%g!cEa84yQ
z>uSZ#UROYX?>$dsG=kaQ`JbK!3^6zGkL=tgyUe8;q!zM($m56*=l>8XLk9u{_yYbH
zbGBD2kN0N&XCp5`tRRyN6f3W&(9qHf+ue1*A?GbLmxz3R>Hl4|gf5>(U9$bPqRvWZ
zfK_0uTVAgJ(3>=5bi8e?hQ-y@?y%$A`x|0`^T1LA{7>0UL%t-Wqy?3g>=IE6P57Pn
zxA9JVk{9~9l|LJ2{yufEa;-DL=|hg5>$6@jq^?>0_2G)iXIBDL?9FVTpB)B9(OWoM
zhEN-vv>vX9gMkuRpgxOnckEZW6B`K@9Ti^xKVL|J2EU--pO;<NxDE~uU!3N6ot9dk
zAaJ|=fm9l|jK>*CqXv@`oW@v_CpdAC{G1jLQLPeCOFieF#?G`6|49|V-#hkd>)}Up
zvYV4rUEgZ9VB;1Wr<vIIa%*|<g!0!$5;;(O2-XYb)EJ6?TUz!RhK(<nhF}!+);3@t
zB?#7!3NBJhA4kBDW2<=ar%>^85(M+3<qrK$f8qK0dH3`2G6Xo9(7b($VC|V_-pP02
zOY-CiaN(Md2J!VYBxn3xEKsuvXQ!sOR(^tCLPUhY-6Aj;5kLjvSI;>(P!VQg9zBWI
z`ME{;7F)~(A4JDwF)LKcaAYO}iO7?~W%APNsu~*`)0x8isovrQUiCbizZ)9dgL|z(
ze3NSwkxV$?H&L^Zis|60Bd3u1Nt)C?9aBex%9#c|8=+V@vzN~PU=^j6RLRegb<$-T
zq_2*FO_ua}kO=o>0%G*Mq_(kfW&B1;&VA;~H#Uu7?m%4Q^L}ovsH%!%5B&D`GKAjQ
zTTKIlXz<!!8ybF$<Usk?J{aIT?$(Kk;!AnCIE_B7fA9Z+<yH7F`g5D^O0D&vrNZB&
z+(Yc2lI@y9?(j~p{V~wO+$I{cgC2auPI`&cd0&5EsENReJY(U_XUC50xtXATX)Ug;
zhC&<dYX+j`L$3;r^`~=%oXphCPbp4?y!L}}Flxb!619MgoL_h(q5y4Rl5+R)Eua&{
z%8xBARi&9=(z$N`rt?4YTj|*zE8u~fh*1-saw>Z_5R*nKCNT?N6O|&kVm<Yde@-G}
zWRz*@{VZ1AB<VU~eqkZcp2GL?FzakZdBpeXJ`f~_&6s*7CMGu)+SIq;>Wp(PHXHFu
zaz8Q<b4|_z`I4N*HLI1pW&!#X&HhUnrdeU{Xam*_k#K_LyfI~O8A20CV!pKmE}%OM
z-R<%I((1>DCs*4A@&BUrlbBUyR8`>?&22x<Xw+F#dJOo$QK5Iov>=;3=gOCj#A~Tu
zYwzfw1_if(n}35lM^-<5QNpB$P4UmA%Hsl1;;7_5eEcY`Q=Q1Mo8~KV$xTMNunFU(
z6$K89HWppZQB8aV`ueUW^F3whr>`>XElQDYW-&Lx5R8z`;WVrHulCJXe?FdoMpq0r
zz3SfaBIR0excu+WXPS4*E4YIm{7X1<#K3(QKJoy{FcO8V9L?k$4w#}hB}hwGBN13l
z7v$&D7n86|IsZM)%~Mm-$jrbNrcEre2b<>RIcaLKPshC=moVh_jZNCHx`jhw*yC9Q
zwdT&pybV8fuq$`hXKt5E-s_iLkK&vi|HjjApND^&t}-ITK;ibka%C|^^7-7<DNvS6
z{<0nEa9#v&qY<4ktMZjd2aQ9^djqX{qZ0Pb&|12c3aDxQ*labhjfoz0yKsW-{e6!q
z{x2#Q8nZkPX-8#~j0<7ob92$3_-rk~C$MvH_&!@}wf;aUk%P}5AM@p`uI$_K8!lk-
zX{K~2w0)J4hF+;^(>>)gJ>_{CH1a?3=&$F0XDHA?U6h(?O**$uUn4hYP0FRX1VIft
zfcZU+=st@Uv54!zCej+d`njyN(Dde&tSmAbE=9Obh3;BuL7e-xQrT~~pMXJ++i8OT
z@$coO5*GPh)nX}ZNcX-Y8@c>pw>tr$35?_R!?(+Lul~F+-XFCyAn55h=LhXl^j&}4
zabEJrOvPuI%GH`Bf}526!Ps=iGhGi)^o-q5t6Tn3z(=U>|Db;b#})2Y2=s%(`{LJ^
zRgG=7410-q6=6}n7>Pl6*ZdafQ2{a|ggMcJu|Gqkq`0mJE-WUoC%b~w;+=jKYxg}4
z;aHXjFghR$C|6gK^FemQp05t%^enc#<?%i=7Q5a^YCissEq;<+aM#y0vO#Wx-MaWR
z*m_-75&*w~f*y}Hzv?=g9d&$q8>sxIw2o;GWlmg$me=VyW`oHj%VZyN!<5COgH_v1
zBqfO<ewYE5M%}6E^&BBq#Hom^HEBcS0DG>pqhpW2jB!}J<OA7oDo4*iA_H~1Vg4JR
z+lhFrj7;PN6<%7GA>YED-N90^sB&okRAP=ADkC5ufb}ivB_g%3k8YhYiK}m|_gy$y
z;P(*f7Wth9<k~1B_JTnc_W=pcbt^sPk=yI>lJekFsR?hARO;7;-ii&tN2u7gh$ASz
zs;i7#TmtkE)u3<T@X=*KVId<Ytwud#aL1-dU#q_ru0Q}A<?ty0D5Eo`^cURvTF=xF
zRxY>^UWKZwPW)b4iak7Zxv!TRNqu@6klfw1f2odaNQbiQyT@Vtc0{0gGu6m-ud3;(
z^*7~1coNc6iK1*$zJ<1fM(Z%mK!DNf<dtR+Z^w^C7hVB&6owd`T<Rr88cE>k9+S3v
z`nG$PG^t8lbVrQZ+i@p5QB1uAM#pL$?WH>~inQK?`OWWBoi?M<cW}vz8&KL?c88+A
zJNZF+a(@)mMSF<7`7DqERs!PqOB<;|?doAigLn@{wi%wiX9L!x(Qb%T_|en**^Acx
zBH4CBu!Tm?w_4u2umkDZT?UQ*$MR2g{uJ(s)z))g0Xi`6iy@zF^XAg6GK|QRjTCv;
z0Es@yOcfiB1yDMYz+jm^2aTehrf!rWr7}_vPa=BAllS1+YU<(SAbYYwL5m4Y)xcfZ
zsx3oxedB8{?4Xi6$!^z_pCXpheXaQ0R~DEm$}RAmt~*6(`HhWFYq}mpK&nohrkKnN
z=qCdgJgLfR;J1_0()y`BbY!0^ysF!bvbje-B|xJ6Ys#}{p&DS;!Xvv;+~>94VTqq8
zU^S>yW7@&{OP7rAwbY;Nx@>Eetb|fA)rY&yG_$_gm&;MYZ~v~|?A2Tv|BY-P6klP4
z7?|UToH8L;Oc&S^o#T@Y^E1X>s<1y>(Nhy>tWdrFS<Gb9-@$v4%!H+~B(_RB*MSw;
zbM=EOFhtKR;K`fz6SUpkU1QICNHpBj?+y+u7;7OXC@_YCETqFa%pg9SZ$XYehkhvM
zRe=VGvOBmVrd5=q;jK0FIV*Zz{P1U7``7H7vqOic5yGr1k+d4(d!3Q-09=3dDsr=9
zzTOs}vu&RjAf`&w4vH{ze6J0WXkK`7IqWQrPC!Wl_5uI=MgC$pElc6%ggZpEbY#L5
z9F~#33+o~tl3^ji>A_R}3JIe6%MjxApeb|&n@flD-KlgDPu5?OkN#qJ$2|bCjO?5f
zY=%OctK=g}(KnJ*R|iN<{rp64E)Mp}!@5I|(E*(rqreE$oR$ST2bB@u0cJg3mZKCZ
z^2}1YS^byQE=iY4436}>6FV}?dKO|)Lzf}*lB<`uZdb<^uQE^1<OIF?gt~>Wj00;B
zTJ~6!8oy{_gxF^Lpw-)s_Pq(ZHTs#(!@zaOlt^wPHBirMa+9Xadj{MG0S3zN5@=3y
zN$pSHf5yzvQ{n@99;3MJlnA>%jb1g$EH6Wm*n?B?C$aKp`FV3L*J6Y}9qwv$BsHMG
zpx2{DLg&GOU34Gjcn8kQ@h{OM^FO`Ob3S<xJb2!+*@(bMg;!8efQ;}D#&+RrZ3g&F
zJw3hnoDCyt$dp@7a6N;WpAxRaoS3uP1^&RP0=W}li|sQ0eu_PqV@~7tWhzyjhE67(
zbPe_q9?dskYn7<U0W&;)xIJ1=w~8BF!e*)iv#Cx=&SjT55TpF-m?m~lzp=7-@$4m1
z_sRLpN&HI}sAHzM-o4PsF}3G!lFdhTHJzQ!R*k3&=P_CjB1UzP!C*J2LzyUn_J6gX
za0RfV%>hSH1A3gLlD+vJ0b(P*kyJH4t_ioWBK@Df){!?TW7QdV)(mDN@&SJx(B(wh
zm+i8hPlN3m_?UBEHrg81sV<qO9bTEO4%q5xXb^^Wn=b<C%iS727-WCT^JZDh*#G1Q
zswdarE%2Nr%}mIaVEcEclpnt*eu`IM{9gGAzojgw+d(*A%)dELI}#8Ed}t1R$oVP$
zbF<mb4)P<;*lby_ZMukz7at+?S9@lUA6~Bb6YJ>cz%vz^cdLk%2K8th92`OiUnf5M
zlS72Np%|y^B}fGIAa^p{IK=lkPB5FU3kora-^d?(-rCwq0C3%0m$|^MRsbvK3jUm(
zrE0Hgh@r%QSMJ|p@HJAM+=&B})%cUid4yYqDPk=_z5@pap{Uph2!OQ{db~S>;~{_q
zV`#)PE>pxJtlkZ3@vBbca3Z~uUrtR)$!l+4OkHUMRL3KIb!jPgmMhbvE!8C}DiY$o
z^5Q`xJ|GoJHrFR~>ov2C|KXLs=5fK9T!aL?xrU=lX(J>7t7}uhmj(U(+jk+u)uNQP
z7-8A?QXV73dnJJM<&p<6Gc-G9xx>%u`gBvj&5MIR5d)mDk27XL)V4VRjc>VvpMN|(
zGRl08xSp27?mg6?{+xvGM{oMMx}W%lw9JDF$4DY&>%zi<Rd%jlt*IEXDJzsOb$6hX
zmOMuP9Iuc#;R`#nr`lnK0>8yzRc7lZ+PmPRDr=a)U&oC1=+FEK;>9n>vIB4IKL<Tr
z-;9cih<tIFVuf25OOAh?7FWpgy>vu6XkrobIQU!zwhZpl5{Sv3e{>D-D@^gP|Eb?=
z^2Yh4V1H?E45Bxw3&)V^W;g$87?;j2B(U)C>lD{<cl+<@>1n#W<x6gx;SFTmZx6sf
ze-(BKhN|3ue(CpyTEqiG&&2bI#O*#``^CH?4sWUXaU!aZAtxU#DAKI#X#jhCnZY*}
z&gPRYxOe`A;?kipD4zg{=)dq=CHaYyjg2iG8JfQtRfmoi6~J}nFX7?sUGnw5<9AsQ
zDRKyc@A>ZL&}FKQm0w+A&D?{rI$y#fLXL4A``F)+%(36UgKMq&Ex`xao~lR#@n!=9
zxSY4v!o&#R03k~0O3qGIzH)QpVP$3gC#x)C5U;X#5{p~u@V=Lg*ptyfrnl#-aK$B(
zYoT{8xRYwCP2uxK15-@&GWMLElOFy=1T5rdGS3%*y~}qXHFsO`oUyprpKCtpqckDE
zBH9F(Pn}W;AlSco6IpD(z0~Hd`>+6z0U6G&0lq;OR~WwV0i(t=gOR|519F8P647C~
zLVWG#(WA?)4T9Y#5w7XKbiulM_<v{)KN~?>0Ko<x)#hw-e&js-V*xDI3deTWN9|@n
z$7L6!n!K(+VKbO!pos0`kG2urON^R&z$8x?=#2KqRG9SjSDRM+DMa{F|F4Buc*X(!
zi^AZ-gd0*`d#9IiA*r`9?J%Yx0Z;#WgRXvHqOs(BHgLW`zr$A`PV(VhWOed@52)Zx
z)DZ^jpZli27LT#_64<ktDNNT!CPLte#4Ewi>+7NO)?rK?nk?4F%X5&*M`JP-wJXO*
z!m16DK48$e9M=b2(CrKB>w*P}vn6KWjO*(s+(&;!J`0R`{~ql#=v;&Eqcb@!bi;<d
zn|>4%B$nan)zzHfqZgQU1=>G)|G#h<_)R(aTz$A|qib(^30s%RDEKX6^v%W>0{qAw
z7!yQXoT?l#+8@g3$C@3#vBsTq**7cQNh*`4{&P#F(;Av(ZrlshH|xyf9UzHnzW$pD
zM_7O69cwl_7?AjAH{3p1OIiPw$uj+O0Wub%*9!S#rB&oa=fjzltx3aPn29&q#t)ZM
zC4bgjh5nW`Q(iy|duElb#i%sdb4E*o@RF00bH9GoeukQTM7i)*nnQ!n&y@%(!u9r`
zX&s=GZ$LXws;#F-x)DInyI!3{sF7Nxx?1LcTwwsX1z0+A_CyCw3(Q3dRwWxr#b8v*
z#IsYk9$3VCI0f7=5B4^{Kt-Qv{zl~uDJc2<R5-R4mvdR1&Bi_r#H9NpUxxfG%`7l_
zn6nKB<VY2IHGR3#VN|V}^xl@3PkrJHbxE)yelIV_bq3rr`&&AuoF_;sv+S)k5aIi5
zs}DC8s!)8Iy4}3I)i5<ZEW1r{b|`3}JWQ*7^R1a`Y^Q_2@Ae5$F6k)KlVH8}=!y}S
z>qF-^eX+);d6Z_DflbxIQlVS*eX>mJ1W5DWyxY!;oRyTd1~N0BIEr9En}Cy-2MS}U
zX5~B9;xSO4sgm*Q<n!(5Y!pV*g^mQr-TOKRe(wx7jp8J$TrK;bC2wBe2DaE+`Eiyq
z1=w6FV}zh$KV8-9x_LqaCOe5E%Ro**;MyuhRA40P2I%#7oA*(3Ngw>j%~DNAL6UJr
zGU>1qyl8iLb}iQCN0+&w<xkt?rNd|xl~$2W@9EE1BFJP5GX5zT)Cx0af!C%?Q`u=m
z=?3KTZs2p}7Zst##>RqR)*hy(yb&1TV@VIy>LQt>AxOD99VwhCAgHH_){nhR?v83D
z8{azHLC?~@b;nZe<g!Q(>eCHCc7SRKdL@#PZn`B%yAR6u*(9wT#bctyB030BO$Ys*
z?(_0=E1EM;9F=^ef`f{#!22)0=MGc>|HjmenKeu;(89%DE!|YuQ{d46yfvf*>}y*q
zd_zGEL{J=w^z&_DpWleVU=~vIYV~gly`&f6eLLDpSDFEmUkbUz2bH2ew-!4&IgL&=
zFwG-zg-a<bV{>-hxdR}Ee7X1*zl%6Z6DVQCUfmNB%dCD2LB7$dfU~+^zMIA`YB&zZ
zLQjYfvR!Y!t&vz+1=0zqu{c=krZG@Zeu95I?fk#V)=UtWwbVYZbH1PZd}b5;VxUU(
zpRN(`mMhKM`5%paB%>GqtJD_T`x=^|(Sj=&Fu?I%<4iKef2l${gx<&GGw3j!gvw~T
z{0sN^MQ7sNdzW(r3JjEgyxC6sGA~dZBsDr;yJXoprG+V0-b4e^#|_XOR`)mO@6M<6
zcx-;A(>!KFBAA0DBBFR{yGp;WaVse+TRlyzondI3JdKm04F=&VIr~pZJ<g~jV)Tjm
zK*`lfo5F@H5_wFiIP$4@iQ0Vr8}cTPEh9<7gn1W=>+At}kkuIQ?Zsa&tUSDbp;c@?
z`^7@P!Qsixzn!EAJw@$zAcl^P@#BJcc65#q3@A2jU0o#eX$DgFxweG;=ywu_6a$y?
zZ&^ya4Gj#)gWn8u#3rSj`b;`DZYU4x2r_22rz%;s-L$;9@CGm$++g+CGo@!13uik(
zEV&a^o<xQ#CjffDqY?@RF#w>ya2Y!B+iuhNea)ON^HXi+{`A>z5Ebe8Gi<e3yEAGb
zua@Ls&;rwXxquHjZ<zGv5?rj>2#~U}vTnTP?d|OsFTW$>iv|0f{au3x6a2vAu+VqU
z2sbq?0$4S=31o7CHsCoY=X6d&`J*!6VksLGYA4M{_LW8TtbglkY7*(eEZFDVwWKWl
z<K2v%+{@4$!O)bzE=ikzH9FC8z6hQZr87?I&oD^bgJLo;k^fLG$O6|Zi~UKOZ|@S$
zl2v!!IP5RB5VPr4egg&2g8&`~%?}*lA?}D$(LVzA3^6!+#2ln*r+JVOA-Px1h5{0V
zO?$r9!Z<5dWtgD0wYxi{WTP5DpKTu<w@}uDv+qxv5N|C7ob12IC*y$YeNQdwmD}1%
z(NFD7<<fp4uku-l@N8>zX82un>DHpB1MLtJx%q%5b<Suc(AR=2?!GKtr6aP(!bW^N
zfeYQf<?biDEPYP!kXhPZFyv4s(a-jZMp=bAD5J{ci4C-OiirA?076|o)tE~kj17-Z
zLCr?#<A)FbP8Y#90@PZiN!!z)n~5CSa4Eny6GlVeH=q?%SL2mr`#MP5q(s%l0^+-{
zClbG(zl6fER|k~nUcY|r^LExXltAVU(#5sEPnUPMkX*l0r7aX0RyL#e_5ZX0##DtJ
zY8)-RY>|U0+TMxH&ZmZTay}%`5t2_S=W_Ifce{FB-znd;hg4LzVxi+veS<M6!?P3y
z5ORU*oo#kwK*OimU#*}W3i2m;m3`KY{>*&sdqgRy7VGGwm^yOZN>fl$x-I)`?Rc&O
z3KP_wQ5U{s#N0$2q=bF_nV_I+CC-Pd0E`lV^?zX=btcave(gMv0ckkgzFa6?UOXv_
zopAXj-x_^aiEf?|yV&&RSur#x?$<9h{+BPU1`_CQMnyzLYXKw$g#e;3G~BxL*68*q
zuqU74tWgjk%72fKhjJJ-nlCoHl~z>rL=iF^U#-Vm9n3de4VyiL4koee6&czvPzjIo
zdHu8Rv=tw>pt^Mw{5!=?%Oi+b(u!S_u-QZ3r!6nNFsD~&^FiM|g<?+2I^~5SfaAp~
zpF-b`ij&5J#y-FaQ@3j7e}}77Kzs?fy|5U|mwUI7`09)QwL1=#5W#wXe1BO<_6N{!
z4Gf8=+8_y0(dnAtQoHbJ0!T8%_qv%^Crdl9ub-|K8}`{!`Hm3-$~H%&Z;m~TE&SWr
zsARlBz0DdR`~={(|IooyhX~N&SFpTQpq9tftu}@eoPR)1!V5SHQ(Vma_3Kyj{y1tO
zFjP>(%q$-M<aBPjYG(Z<VT(qqZ;nRkn|bOO-KTWex`>zv0ao=hVUyX%f0}2lr5&Kk
zqEnGbE+Zr@u_F9E`}O<akfjSaK3xvHb#A{q9WE^|cRzNB#sb$woU+rQ^u<y}0+c0o
zKB-k&T^)COO={*J0$u@psR9uEKtSy?u+549uP`|^WpQ_P;uc%d$Mmi&b*nk+rtmjg
z<ac8&{D`{9Z>$QhrQE}PpZ{|&em)k4ZDpoWH*=)@Hc&#jF~;rkaH)B}ZWvB4q?!6+
z0^Z5v4{tQqQ}{k~fo=gZv3IUY=j}2W`yjVgRaFHEos6R+8}Nc*<>jp6vb<jVv&~$*
z6)ZV(pG%V+9VYmj-j?*GztF^w7_7625ZRVGc)|J=cS*M>3&Ao}d!4`2=}`6s%HYss
zI03xgQ>>c#UYWB6rR#d+&1w0+_yIO>4VM9d41AY)5C;LxuVG>m+uABJIWx2Vbwcf<
zC;=mu!e{-^v?xQY9%MBK-WSi5ru8<?gCQf;wI$ma>zptJ#{5T>&Sx8=?Yse;Gq990
zP!E-SIEhc^v&EwUmpJhA^B-pu34=K%8Xv+5u7C@uwHu@Fu&w|!L_<RZ+^jbSlX!gp
z*-ZrZmK#hBx@GU8z2i#zZmrdLXpkSwxO7U@q%|CwZc@*TOeB1vH<D1AIlt8Sv6llp
z%+u3zZA<b|IF-vR4o*{9^~ZT`D@9{~k&DSZsuN#4_jJK0+#`np50{Za0oG-8x-m#5
z>~3vrzUCoix!9s^YHF%p#QjN~l;<zKR-pv$rKZ&P)4YEXJr*8bUf+S30aOQ>^#uS8
zv%yK!lyLcTH=xm~9}ld{8~_FVW>>m6k=^GYJ&0li_k|hc+fxX>S+0E9(A3xn|KE*`
zx<_=vJB3N?$V5H^+}H|rz{Z>r{xW7w3e}d-uOJh-gg|p@RXE6S`f$4x8b22PMj~pU
z8sE+4F#L+aXY=!0&|QDfYIWP8C4tHY5xCLI9Wuf1H2Z~)nf3(-M=Yj-8mLEwBSufr
zF)Rf`0z!|@uCb~Sg|u7@dhMPhAVgE9XMKxS-O+l6Q)pPs93DmKl@O22&0FLW`jRbj
zk{|c&)~_Fc8$a1dg8uoA&E;`-ghm(uMjt1t0_QeeV^)LC2r5w0uU{7}|27535D?B1
zhsxuuv>L<QM5v6r#QqYoYOxEM#l!;KqE8aYly8BT$#7Cbt7kp>H!rkpow>R&Kks_b
zFk=ab@OR*EZuboZ_<mIaBGFv|j3I6H-1`@7y)#}Uv|CzWnK2(Pdf1#i%>9!>XXzA0
zGH1hT!lMo>?ajuY1|f~&)T5Y9*%l1<w;fbtn1eST8eNt?faEELLI4Y1-2_FjgWoAc
z8uVbRjHOZ)stoGoSH?R57lK2^EoJsWAFKgH$7z0A4gBoWny;f<VFaPj0<DcAmaAO|
zlh0Lh9d$w}_6pfXF>?~E!b#!vO@@)ub<AXCp#ifVMo81LFZPV_TU?jUYnz`L0>e56
z%Jd%qNjVU)Vet--jlS=#L6HW2lvSj7&A$aVbFIN$`0re4w@1|m<Kb9z92TE@%s^n;
z6PI$#x9ZQY-@Ax1?Qslw{$<x0L~69-|8RTQh7T1xL4j8rKl>fFWZ&E>Z^(fKod%2T
zSLxTv;4A3Xn8p&T2_b?L06#<Fn8e(my74nd@)<ztLPzaqfcO7X>_o2Ixaj)LRg*5Y
ztn5?03KG$|)BrVa*<fTaN6KcDSCA<q_ADG$u~j!C`CaTwTUjvxdx98Me~j`4cYJEl
zz^v5Lz~KVE@NV?+4qwoDJ{-WvXOCC?S7|}lq`)1*%R_bSb3@Zf$Nrla@^3}SS4g9g
zrT~H3P`o3PXES^H!d!VXAhBXIYIj4TBG|d-J#O>Aab>u#0SGJ)K!u7KDM$;rNbGq$
zIChKRM%HGaa(2DhDOm>@>5oR|1$ZEAaGbH6EY<MU@^g0PfTK@$sXxU;F;2t%;F(M!
z?~s0{Z9E|*J?mszE28-X5x1RNl*>O%<Y9{@ky%aAGwlv~*b1`z^;x+2;bsr}VhvHF
z6X@8I`#IZ>W!mpA7OpnAf}$5(+EgVMG3emg;g64`d#VR}Gd2A$N<Qab#!pU8<^pBY
z;YV?ev%1VoZLPH;@Z`Ztz-_=)1$bDrm(cHLDL1BX9S3!wMRimfSzlAI8rnyl;#Gi^
zgX?9n;ZlI2G*cN96Y7<>5FgN2szKxE>@3Gk)(gHM(SrybN+(Fh(?L0^FIMRO-^BTk
zsHv%`;>_@*gUo?6p6q?2ftb132?+GN))I96^)5_|3r|3y*8vnfg!4zYKDMzqG?s!L
zJzV>iE$+u1aN_|IPI$l*0^@V{eKa7+pnD_=D}ey%_VVbJqhke4f<GA8E3_WN_m7&<
zHp?^Ben+-0cF_e#VS_HiXjQ8uwVTn%d@y*!XfSo-Nrtvp;=pn~L{wur=rKU$6EEsk
zPLKpV4NXo>6*+N_daevY3Sl_GEI>PkhKAr+FEl!*b6bR()(N=%>20nI(UHnOeBuiF
zQMU+GS2Xd*>fX?0-KHwp$IT!SS^1dB$wh1_Ro6)$(kJGS?#ki-u>lYa?k}~G_S5*W
zfB5i$HXQW0N_<>+(?yLhaYb9*{0D2a-8Kq~99(q%fg$`G2m=v+|JvLK<pezlGpl`$
zf(LEj_u)Z%WWfbX!CE+yfWAQ1gKkJCUZosS5t^Bn1*^2r_$gPQ0R(MQ7<uTF0P&rl
zpU;Ad1o#CEfRNAkX4V^KOl3Vi8z$0JAO@g8HVQ^0v-**%^S5^b9I1fdF0@x!|1WqW
zhVNIDUw+nefWb;__P(X!-#%R&`1LxbF4`okGqnaWn#ai6DTDd#nkyN|NZdlI36QF)
zs^FCU@5#w<AiKc;Wx#H80`TE)wg3#kRelkD3c51sl`;vCM7<9oT3W;z85v*Ro&2{i
zfM&iaz5p~|A2q>7o|dx{dMB>Jq^vD-q?urbZ69udWvo@u)Zf=lXpA<vH_7QY!MUZS
z10Kge049bTh)R~o6TsNM1in6gdRq5uW25ECpC1?zGVnYlOgvu<l=h;0RWL#h7Mk`D
z6#*atCQ1p&8*n(_P<|0{q@2bWS3#3}t&Eeys=#M=@iT|g(}Y6NWr56DqhxY{VVlOF
zS0E*ZHyDD0&HcZWyW6DzAg*sO4#)(Z(wEYp5Cf2cN7boegg7rW!ttk_y}dL@)B#Ec
zQlEk}b--kcSV&n4^7OOQu32#NU8iy=@h?yP3N9?cy@^n|@I&AZkor9~7Ls$nZ>wK#
zgA0f7^Y$geZl565#2`8+zs)GN*qg56A&29ZBigc>nto83UB`>pWf0p`Lj=YK(h<<$
zel$T_FBg1m!TqeT^$mR{VzkZ*5)-rXRArr0PT<W83JMBt>j45WKpkw))a1zODM7qI
zg$2?$!H=6M=scb|Z)6150jL`pRm9(m{}=rM1fiN}-v3Xmzi(%tCIQmSwl>N<^#S)4
znNFprl^2Y(Jj*YwR**hv9c3DLfm|5QqwAUaFu63Zq1~OPby=&w(xwIHD{&pj58&gv
zh!fxyLu-Zic+^<n-A6KCuJ=*C=>@~PMDKHiO;VHa|9>?>b?>Nyc@C)^LfI1gNNrHA
zi*d<e1Y3HkI&lmV`=@ZRSGMo&3uD)A)b<((D!^-mEoy*I0@P&b5v@U;H6|+?8;{d0
z_kYO>I)2~x2yD0AvXZ_#zrh=nMIpsOyxp(GfL>^4Z`m!=Mj~2A5#&lP7naT{xA9uR
zqZSFTsOhADZ=R{&fxu_~(?LM&5DEhiNZ7J`xYXu0$dqOd><*6hYXfe*w6(Q!tE<&=
z>IB{$uQs!kmc@9Tj|;jm2yPoG3m!{7EILngFRGYW+k37wI?aW`*;B81d;Krwh;T;y
z_Y_fL5`3pBI>7B=&A|VG{PT->Hxhh(fOP@J5o=o40CF4uEzy-|Oa(MNs&sDJTEoDt
zl;H4*jBy6!H~<-4g(@G<j5$aC!K@cnBMAQExB=s`-nih>qxy2GH#e>{7|0dq_1ZxE
zsP#H96ay}sOvu%IadcB_21qA%5C?}+*xVM~`gd>-K&U3;HkTSgRHP+{vN!%!f?RS#
zk^vrOWeje}{WKRxfyTHp&W{Un&D4cehD40R&RZAwNOd^K1YN^NPwvDSSUrIM9+(2y
z3=p3~Q02{ytQ0oAUBtn1oeFM`ZFO<lvZ|^dW3Qf3UxrW5vgD#M?&eKYMY7PLfMWzU
zK_y2fsYUxK?+{M3e^4o6wA+v_-O`zLL3G!{jU${61KF|^90driq1bXJM;U_wHH1DH
z?6?MvPM>HJ;Ds!B^@Jt?9u47)G&&>ibGM!lEr||Nng40;O1z<L`|#Kb5i=rEni*@W
zL`9OpLzME;MhRIel^B&uG!z-7OxjcuQ?HaFL@A-g+UlicN-2~=CR(OKzU!Xe?>pc3
zFMQ{94(1%kGxOZfZ@GTg@A}<0Q(HxdHo=E^b;CL$x;JS}Rc86W(v8;=WZRv1yrxCX
zwqIRURn@YutFWi9c~I}B-=oYNL-rNoItq@8v0&ULuDCXWR6fDMsK}#Y+)q1w2;i+-
zx32#Bme2-ZK;`b;F$RavN&+8vhCe=-aZw=4nAWgL>#`E-Q_Eh{mUkyYt8ggruFl-+
zlf}N~^dKCX#2&xKUJ;z@N?6y=Bi|6o18|>!5s0j?=}ExluP<9OYrlMm&YTytY~l1B
zYK1ZKoO>0!qJ}v)z8SLWSn(ewaiYT2tq*y>%zYw17~<C39(U@qw*$bkJrr4CpI+P9
zFj4p4M1@6QOQ+#-L}jj_aaD`~Bk{fMT}RHXF`CD?=iQ7B$56v<!#ZPEyfhAYBJM%A
z-XB<o7DEa`Z~);b2u;+;b?0&vV+Z;qr2M>p|2|m|Xp6X2I;-Ey)S?y*ni<gu46ggU
z!6vMZpvX96eME!8w(*X&>91b6s5-X|&f9=Hkfe!b>#*cWQ61cn^=6LkPr6>@`^>8@
zo)m>MnkXt56(DIJZrJd|+<<uKgIS+lZJ1WDy!o9iZ!MP=V#(#sqJz`bMOln;1dBI!
z?%dI~Y47A(*R<@}qx@7Mq*gY?Rqjdp+;8dNPiwI?*qAL%u_YVZ$<0N%G2A_1wfsoj
zQeN55xoS?wPn^i#*uEfb^AoSgg?)fNEm|^Jz4G24>4LYWvri=(E)8}rW<I=_qHrNX
z-eB#)GST*21WVFX#J~2`O~~`R6{loXGMwG3M{dE(Ag%&Iq`tB7%#0!)L?Wpc!Hka1
zPQn{-TdOKRYPdA&0}(Uh{Aw2;PkQnCbIC*#cTUP)J$;=|%cgx0gNj?2<!T4s8QLJz
z*U)No_lRAncpTF}cOS0Ax`Z+bs0_;l=HM^6lm<6jo9+(<>ayD<-DE);lY>q7Z@wQN
zsM_#j@=^!2l%S+I^)(?Om3YAWPxZGMUu{p{p1Mz8c$vhDf&=QZ&{$Oq9ru38`i7w_
zWK+V79W-{Ug7S~lVubtL+FjBfhVNTq$;f-@iRpi+JzAyO;B$Ci{-d?!30_ss$f@6+
z1(00@h~CY6Q@x`CH#E@%eQb5Q75*9pTa4uV1ylig4-;3d3AX(#Qn)1P&X}AlZ&cr6
zT{I^HJ`_-UIA5G4IKvRs3`%}PnGNVc!p_w>4UhEZ4SfCn>9$%=M^2tU9H+S4h~0~W
zBPnWfzJH~%cX+MgvvQto)`F9JZfutbi;6UlXCc@DK53iz#|3kQ(=Y217yGi7!@xX}
zNc2IJl9TQs$M~Ulz<02qx~Ah5(mle0swH%lJ@?EXW`s7EEhv>|_o_?Zg}{aWFK2Bl
zVJ6at6Q{9&7G~g35Zqv7F>P9ANr?{NlcKyVC=9SXQh&A>=G@jGjvDTqwB0>v7`ry@
zUIwDCUBl(nQ!BWYn%c9YEZ44EM<1ThQ+R4X#eDUHJ6)neD_ShSfE)#VVAw>|-asmo
z-`zWQFw8DJE7-Gq#|=HfM4wV^G}lpZHAn<1V`f^c<S-JcFAitTn6+o+j3mOq2Uic}
zXJ;#73I1F%@95-8piu>yB>({O&dAqk8z0!Acrj;EjRC|VdOiQnu)k|;(4=2u0trDe
zHPE@~=B#Ep!drlsqP?>oU$d+rya=M0mh(Bex6tMJv}<fX4Pe6_2r1&@a}4w(C}|k^
zi3);B+o?m?e3hc~u;3&sMVqxi9$60+H-PYI=#Pk<p|{ye`pl@ep&*j<wPvlia#4&0
z{%OT3HtkI@YQ14TVd;y@!5UCEdZcZ8cP!PGHyaW2D^CIh?yzyAd2wj-esp9ldl>v0
zePVRm4ynrm){4sq%hFk1pq!9b2>`zlxDI4Vu!2F=xaAS<^K#Laa9_5pE+uCK=px)8
z#iloGH;%(mnf<Gq!fvZxZQXDE>GP?<wAeB9w<?I?LRnTE$7IAUWYhYlrX&D-f<p$e
zj|Ck5DLA=`7?~`Y#F=)r+Jfwb#o=~iH14M4e1Vn-ke(JZJKa+=kqT=+xy2tW@76c*
zBf}n8wNJUEFC?R5)S`?w`E!~Rmb}+Vd7QJf^^S&HW1&9SsrSw0<v|?wt`u(tgSDbR
znChG<+HD7aSY+$InJ|8QOZQ7O5W{}FJxVk+d8NI|yH#{r2=~Au8I7Wqy=0ed`ru*X
zP_XL$C<1vRIUdlg?-=-LsN&o>Mq@39UHvO-yq8GbYhks<@{qVeY=SW*sBT2mU|(~!
ze)6oHy}jWHxkvcQ*vpp9o@yA!nKNhTk%?Zy8AW+e2A_7?w@1mDsD8sa3`Kc7$&0ee
zMQwE#*7GkI>bgygs~{i^3cZ$p)2>4zsK5!tG!8LX@&IK(kWo0oRNH)&t1HI?VuUTP
zr|lmoQg)8f4{~_<@+EqLoC2R?$S#Duy?Eu-t-Ge>NgyMQCm-Z1k0zh0j;}5|R85bP
z){^UAmFj?ldn^$wTD-V&-*!#=JM%?TcW83&Y=~@&iKSv+5MmW+C=*jMre)vR`)$5?
zO79R`)ZtpnEctKZ=Hw$sv}~B7YnvxjC6CQBbxgD5vg_*VT4Ym5e-Q1iy~C)+E@!5?
z-#LMj=(@)61a}V)|481_3!dcKZ-4}<uyagJ%5--42<df|aNX{gdto-XZdtUok!xy3
za{tLGD`l};tXEUSixBDVz2d=NV$2tC7Zx|BO!^NPdt{*sr#BflF0t>2hTE*z>)qFe
z5xfin-Q$xaiJyi^AU=ORD~>ds&$&HYc6kjQ9vKfU7J-{ryRUT;l{6*|P8?Zp#!S>i
z7?e)nB|wzVVXA-r3=g=MW~Bjs%mL5~WOO^Mh+audlPZ>a;I;Scu4nmj-~)*XXqoyX
znU?o6Aj3wjQN2%h<cvvowFU5^fz=&HBJu$GHXITPg;(`-X^;vJB~-tR+z`^?ASUgw
zbRIP#8G(*2>OeD|1N|r}h&@$RPD!J!7!_A;F#n5u#^b>ACV``@Oz$itXib4<17#9b
zoJ!RN;7LTH;1XS6HxNVwJrcDGf~Lk%9Tq!$6j}R%$GcT^Kh@|b`*k;!WM^kz>fday
zNG3Wwu}q){2(|M&1EQV`s7gw-_8Z26{3#%7y~~f$0^-3D5%vuhM9}UOj&cRLsC7$u
z-x#flo-w?Er@Z?)Pv`tj)bzIs;Cw?~Oef*{Z_xDSb^@6bK@4(DpQS@QM%Q+N97AO5
zuO+~o!%eZ|`nL)_X<||;Zw`Gfu7C1m_qx0rT9_rM!tdRHpjMmyq^MewxP_%>f=8j=
z5W0!hhJ?gt1*2i~cwVtnlJ<cq%ma2xwC_=6saX*-<PGc<A=%I$p-c=N*28`J1*j+^
zGgAhN?FyiWeGo=jZz`gwbJAFmo&nx?Llgh+o&&Nd+QUd_n7-*XVxqkYrhxBs8l!Ly
zE-o(K<<;On_N~h;8A2e@s{6+)EB8@{Mv@SGz!*Hn1>bU?v3=LU5MbGWTi!p)<k~q3
zW2ln_hwl<hno(G&K~skUWR&=Hy$RVlyT2!ycL)8T1QbjUd)BMc-OrQDx!dt6w&#@%
zClSxDZ5jtSw+jt|=q|B<<=`sEDB|7G7Sv=IESGV?NCE^M80v$-&U;BKn&?-1>n0F{
zE6mTQgY@eph$;{p{7}>)0eaI<a#!RGH6TVm_-o+nq>m&kL8mbllJf49O|pnkL<lE!
z?wth%1y^}ITx89fHJrtZH6d;~8n#E9et?{N#1`%}j^lu0>v)`l1Pc1q60nG1zvs`F
zk;mSQld!<QLpG*hkoeQQaYr%&D<lMW2K2y|rdn+@f{#uniKk1#6qbU`<@0Me8+3h(
z)*eR9gG}SFTc!SP5frdt_psBP3TgaV%$M^UBQWTcpq`ZGtKPe&t(7Bvspx2%J{!JC
z9{FVKx%T^SuPBQu(G!$HZAh4oR}LP_m+$>G&{wT6^mtaKMNy7SVzYRT#i)4vbaZs|
zUJ~#y-6nN*^i2Cqb)Ri-HAio0;uhan#g<N=A#@u_YdNCsvtfFnkq7I`T<wfmj7L{O
zeH&uPK&{zLet-LxSKExwpFbZ4UuzZf2CYdHARK7bZ8h5@5SiC~jp*(RNUg!{V6o&O
zC&&-1Wkio_0UtjF4Vmb!P=FPb)A6Jq33z;S(v_y2g!%>8f>nI)|NM|!?>O?!LATj5
zG>mGk{OgPhiFVR=gdUXpeKCU`G?YXE;Y8F0k1`!ROShVVyCP1nsyqEPx|!9$?Ph)z
zG*UVf8q$BNL4w-YD&pDX!UU3E;<h^3Fd@iNDeHSF_D3A)l3@u;2`|Mrsm7y#1H>m_
zzsS>blARnK6$q0BKCDBZhY!Bhd)wQ?7~a_m?e!GA53M41#q}>idwh8V(|mXXn-hFt
z>JO!>P?w2<E1>SsxUB#R(#_qUEDw!;*nXSn!H|<~Bm_2O%}+H87^27P=OUh4M~f)R
z3dA)&M*_;T1hhHFx?$tU2R0J#24Ap|h<#)!A)`H!C*ark5Y)#HC*4R0kvIRG%IAAy
ze<M2i5>=h*KdLT$`FqXS2%gf!|IZ6p(ElvLjJ(6%G@4AC3uU+HpC$S~wv+pVL*2~Y
UhFJDaYoy`J#c=_ZYwx%7UuCfH$p8QV

literal 0
HcmV?d00001

diff --git a/lib/mpl_toolkits/tests/test_mplot3d.py b/lib/mpl_toolkits/tests/test_mplot3d.py
index 42718166d32d..5b8485f802b8 100644
--- a/lib/mpl_toolkits/tests/test_mplot3d.py
+++ b/lib/mpl_toolkits/tests/test_mplot3d.py
@@ -110,6 +110,17 @@ def test_scatter3d():
                c='b', marker='^')
 
 
+@image_comparison(baseline_images=['scatter3d_color'], remove_text=True,
+                  extensions=['png'])
+def test_scatter3d_color():
+    fig = plt.figure()
+    ax = fig.add_subplot(111, projection='3d')
+    ax.scatter(np.arange(10), np.arange(10), np.arange(10),
+               color='r', marker='o')
+    ax.scatter(np.arange(10, 20), np.arange(10, 20), np.arange(10, 20),
+               color='b', marker='s')
+
+
 @image_comparison(baseline_images=['surface3d'], remove_text=True)
 def test_surface3d():
     fig = plt.figure()

From 7d1f71791fbee330d377d8ef42b6cb8d79bd1f3d Mon Sep 17 00:00:00 2001
From: Ryan May <rmay31@gmail.com>
Date: Tue, 17 May 2016 11:05:30 -0600
Subject: [PATCH 2/9] MNT: Raise an error for invalid args in animation.save.

When passed arguments that are only used when creating a new MovieWriter
instance, raise an error if writer is an existing instance rather than
silently ignoring the arguments. Also document the fact that these
arguments are invalid in that case.
---
 lib/matplotlib/animation.py | 36 ++++++++++++++++++++++++++----------
 1 file changed, 26 insertions(+), 10 deletions(-)

diff --git a/lib/matplotlib/animation.py b/lib/matplotlib/animation.py
index b45d3ca3561a..256d31bb808b 100644
--- a/lib/matplotlib/animation.py
+++ b/lib/matplotlib/animation.py
@@ -686,28 +686,35 @@ def save(self, filename, writer=None, fps=None, dpi=None, codec=None,
 
         *fps* is the frames per second in the movie. Defaults to None,
         which will use the animation's specified interval to set the frames
-        per second.
+        per second. This argument is only valid if *writer* is not a
+        :class:`MovieWriter` instance.
 
         *dpi* controls the dots per inch for the movie frames. This combined
         with the figure's size in inches controls the size of the movie.
 
         *codec* is the video codec to be used. Not all codecs are supported
         by a given :class:`MovieWriter`. If none is given, this defaults to the
-        value specified by the rcparam `animation.codec`.
+        value specified by the rcparam `animation.codec`. This argument is only
+        valid if *writer* is not a :class:`MovieWriter` instance.
 
         *bitrate* specifies the amount of bits used per second in the
         compressed movie, in kilobits per second. A higher number means a
         higher quality movie, but at the cost of increased file size. If no
         value is given, this defaults to the value given by the rcparam
-        `animation.bitrate`.
+        `animation.bitrate`. This argument is only valid if *writer* is not a
+        :class:`MovieWriter` instance.
 
         *extra_args* is a list of extra string arguments to be passed to the
         underlying movie utility. The default is None, which passes the
-        additional arguments in the 'animation.extra_args' rcParam.
+        additional arguments in the 'animation.extra_args' rcParam. This
+        argument is only valid if *writer* is not a :class:`MovieWriter`
+        instance.
 
         *metadata* is a dictionary of keys and values for metadata to include
         in the output file. Some keys that may be of use include:
-        title, artist, genre, subject, copyright, srcform, comment.
+        title, artist, genre, subject, copyright, srcform, comment. This
+        argument is only valid if *writer* is not a :class:`MovieWriter`
+        instance.
 
         *extra_anim* is a list of additional `Animation` objects that should
         be included in the saved movie file. These need to be from the same
@@ -720,6 +727,20 @@ def save(self, filename, writer=None, fps=None, dpi=None, codec=None,
         the individual frames. This can be used to set tight bounding boxes,
         for example.
         '''
+        # If the writer is None, use the rc param to find the name of the one
+        # to use
+        if writer is None:
+            writer = rcParams['animation.writer']
+        elif (not is_string_like(writer) and
+                any(arg is not None
+                    for arg in (fps, codec, bitrate, extra_args, metadata))):
+            raise RuntimeError('Passing in values for arguments for arguments '
+                               'fps, codec, bitrate, extra_args, or metadata '
+                               'is not supported when writer is an existing '
+                               'MovieWriter instance. These should instead be '
+                               'passed as arguments when creating the '
+                               'MovieWriter instance.')
+
         if savefig_kwargs is None:
             savefig_kwargs = {}
 
@@ -735,11 +756,6 @@ def save(self, filename, writer=None, fps=None, dpi=None, codec=None,
             # Convert interval in ms to frames per second
             fps = 1000. / self._interval
 
-        # If the writer is None, use the rc param to find the name of the one
-        # to use
-        if writer is None:
-            writer = rcParams['animation.writer']
-
         # Re-use the savefig DPI for ours if none is given
         if dpi is None:
             dpi = rcParams['savefig.dpi']

From ecfc1e724173e2ac7bb8e3d1fe9a368e976b3df5 Mon Sep 17 00:00:00 2001
From: Ryan May <rmay31@gmail.com>
Date: Tue, 17 May 2016 11:07:15 -0600
Subject: [PATCH 3/9] MNT: Improve error message when saving animation fails.

---
 lib/matplotlib/animation.py | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/lib/matplotlib/animation.py b/lib/matplotlib/animation.py
index 256d31bb808b..6a9232f1d834 100644
--- a/lib/matplotlib/animation.py
+++ b/lib/matplotlib/animation.py
@@ -234,12 +234,14 @@ def grab_frame(self, **savefig_kwargs):
             # frame format and dpi.
             self.fig.savefig(self._frame_sink(), format=self.frame_format,
                              dpi=self.dpi, **savefig_kwargs)
-        except RuntimeError:
+        except (RuntimeError, IOError):
             out, err = self._proc.communicate()
             verbose.report('MovieWriter -- Error '
                            'running proc:\n%s\n%s' % (out,
                                                       err), level='helpful')
-            raise
+            raise IOError('Error saving animation to file. '
+                          'Stdout: {0} StdError: {1}. It may help to re-run '
+                          'with --verbose-debug.'.format(out, err))
 
     def _frame_sink(self):
         'Returns the place to which frames should be written.'

From 5e5999f8a3b26a8d34d7cb69d5ab5b697af3c7a1 Mon Sep 17 00:00:00 2001
From: Ryan May <rmay31@gmail.com>
Date: Tue, 17 May 2016 13:10:11 -0600
Subject: [PATCH 4/9] MNT: Add original exception to raised IOError.

---
 lib/matplotlib/animation.py | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/lib/matplotlib/animation.py b/lib/matplotlib/animation.py
index 6a9232f1d834..52cafea7fc0f 100644
--- a/lib/matplotlib/animation.py
+++ b/lib/matplotlib/animation.py
@@ -234,14 +234,14 @@ def grab_frame(self, **savefig_kwargs):
             # frame format and dpi.
             self.fig.savefig(self._frame_sink(), format=self.frame_format,
                              dpi=self.dpi, **savefig_kwargs)
-        except (RuntimeError, IOError):
+        except (RuntimeError, IOError) as e:
             out, err = self._proc.communicate()
             verbose.report('MovieWriter -- Error '
                            'running proc:\n%s\n%s' % (out,
                                                       err), level='helpful')
-            raise IOError('Error saving animation to file. '
-                          'Stdout: {0} StdError: {1}. It may help to re-run '
-                          'with --verbose-debug.'.format(out, err))
+            raise IOError('Error saving animation to file (cause: {0}) '
+                          'Stdout: {1} StdError: {2}. It may help to re-run '
+                          'with --verbose-debug.'.format(e, out, err))
 
     def _frame_sink(self):
         'Returns the place to which frames should be written.'

From 7c98a9886e5be6186edcfb701e761776a6f7d40f Mon Sep 17 00:00:00 2001
From: Ryan May <rmay31@gmail.com>
Date: Tue, 17 May 2016 13:10:41 -0600
Subject: [PATCH 5/9] MNT: Re-arrange docs to avoid repeating save verbiage.

---
 lib/matplotlib/animation.py | 49 +++++++++++++++++--------------------
 1 file changed, 23 insertions(+), 26 deletions(-)

diff --git a/lib/matplotlib/animation.py b/lib/matplotlib/animation.py
index 52cafea7fc0f..a86871fe255d 100644
--- a/lib/matplotlib/animation.py
+++ b/lib/matplotlib/animation.py
@@ -686,48 +686,45 @@ def save(self, filename, writer=None, fps=None, dpi=None, codec=None,
         If nothing is passed, the value of the rcparam `animation.writer` is
         used.
 
-        *fps* is the frames per second in the movie. Defaults to None,
-        which will use the animation's specified interval to set the frames
-        per second. This argument is only valid if *writer* is not a
-        :class:`MovieWriter` instance.
-
         *dpi* controls the dots per inch for the movie frames. This combined
         with the figure's size in inches controls the size of the movie.
 
+        *savefig_kwargs* is a dictionary containing keyword arguments to be
+        passed on to the 'savefig' command which is called repeatedly to save
+        the individual frames. This can be used to set tight bounding boxes,
+        for example.
+
+        *extra_anim* is a list of additional `Animation` objects that should
+        be included in the saved movie file. These need to be from the same
+        `matplotlib.Figure` instance. Also, animation frames will just be
+        simply combined, so there should be a 1:1 correspondence between
+        the frames from the different animations.
+
+        These remaining arguments are used to construct a :class:`MovieWriter`
+        instance when necessary and are only considered valid if *writer* is
+        not a :class:`MovieWriter` instance.
+
+        *fps* is the frames per second in the movie. Defaults to None,
+        which will use the animation's specified interval to set the frames
+        per second.
+
         *codec* is the video codec to be used. Not all codecs are supported
         by a given :class:`MovieWriter`. If none is given, this defaults to the
-        value specified by the rcparam `animation.codec`. This argument is only
-        valid if *writer* is not a :class:`MovieWriter` instance.
+        value specified by the rcparam `animation.codec`.
 
         *bitrate* specifies the amount of bits used per second in the
         compressed movie, in kilobits per second. A higher number means a
         higher quality movie, but at the cost of increased file size. If no
         value is given, this defaults to the value given by the rcparam
-        `animation.bitrate`. This argument is only valid if *writer* is not a
-        :class:`MovieWriter` instance.
+        `animation.bitrate`.
 
         *extra_args* is a list of extra string arguments to be passed to the
         underlying movie utility. The default is None, which passes the
-        additional arguments in the 'animation.extra_args' rcParam. This
-        argument is only valid if *writer* is not a :class:`MovieWriter`
-        instance.
+        additional arguments in the 'animation.extra_args' rcParam.
 
         *metadata* is a dictionary of keys and values for metadata to include
         in the output file. Some keys that may be of use include:
-        title, artist, genre, subject, copyright, srcform, comment. This
-        argument is only valid if *writer* is not a :class:`MovieWriter`
-        instance.
-
-        *extra_anim* is a list of additional `Animation` objects that should
-        be included in the saved movie file. These need to be from the same
-        `matplotlib.Figure` instance. Also, animation frames will just be
-        simply combined, so there should be a 1:1 correspondence between
-        the frames from the different animations.
-
-        *savefig_kwargs* is a dictionary containing keyword arguments to be
-        passed on to the 'savefig' command which is called repeatedly to save
-        the individual frames. This can be used to set tight bounding boxes,
-        for example.
+        title, artist, genre, subject, copyright, srcform, comment.
         '''
         # If the writer is None, use the rc param to find the name of the one
         # to use

From 27d119900dd30846548f21d141f2af3740dbeb47 Mon Sep 17 00:00:00 2001
From: Michael Droettboom <mdboom@gmail.com>
Date: Thu, 19 May 2016 08:15:33 -0400
Subject: [PATCH 6/9] Revert 2c1bc6b for performance

---
 lib/matplotlib/tests/test_path.py |  1 -
 src/_path.h                       | 16 ++++++++--------
 src/_path_wrapper.cpp             |  4 ++--
 3 files changed, 10 insertions(+), 11 deletions(-)

diff --git a/lib/matplotlib/tests/test_path.py b/lib/matplotlib/tests/test_path.py
index 22035b3dc501..6543137c0e1f 100644
--- a/lib/matplotlib/tests/test_path.py
+++ b/lib/matplotlib/tests/test_path.py
@@ -39,7 +39,6 @@ def test_contains_points_negative_radius():
     expected = [True, False, False]
     result = path.contains_points(points, radius=-0.5)
 
-    assert result.dtype == np.bool
     assert np.all(result == expected)
 
 
diff --git a/src/_path.h b/src/_path.h
index 6b1b3d1f7053..8ba732509366 100644
--- a/src/_path.h
+++ b/src/_path.h
@@ -79,7 +79,7 @@ struct XY
 template <class PathIterator, class PointArray, class ResultArray>
 void point_in_path_impl(PointArray &points, PathIterator &path, ResultArray &inside_flag)
 {
-    bool yflag1;
+    uint8_t yflag1;
     double vtx0, vty0, vtx1, vty1;
     double tx, ty;
     double sx, sy;
@@ -89,13 +89,13 @@ void point_in_path_impl(PointArray &points, PathIterator &path, ResultArray &ins
 
     size_t n = points.size();
 
-    std::vector<bool> yflag0(n);
-    std::vector<bool> subpath_flag(n);
+    std::vector<uint8_t> yflag0(n);
+    std::vector<uint8_t> subpath_flag(n);
 
     path.rewind(0);
 
     for (i = 0; i < n; ++i) {
-        inside_flag[i] = false;
+        inside_flag[i] = 0;
     }
 
     unsigned code = 0;
@@ -118,7 +118,7 @@ void point_in_path_impl(PointArray &points, PathIterator &path, ResultArray &ins
                 // get test bit for above/below X axis
                 yflag0[i] = (vty0 >= ty);
 
-                subpath_flag[i] = false;
+                subpath_flag[i] = 0;
             }
         }
 
@@ -164,7 +164,7 @@ void point_in_path_impl(PointArray &points, PathIterator &path, ResultArray &ins
                     // Haigh-Hutchinson's different polygon inclusion
                     // tests.
                     if (((vty1 - ty) * (vtx0 - vtx1) >= (vtx1 - tx) * (vty0 - vty1)) == yflag1) {
-                        subpath_flag[i] = subpath_flag[i] ^ true;
+                        subpath_flag[i] ^= 1;
                     }
                 }
 
@@ -196,8 +196,8 @@ void point_in_path_impl(PointArray &points, PathIterator &path, ResultArray &ins
                     subpath_flag[i] = subpath_flag[i] ^ true;
                 }
             }
-            inside_flag[i] = inside_flag[i] || subpath_flag[i];
-            if (inside_flag[i] == false) {
+            inside_flag[i] |= subpath_flag[i];
+            if (inside_flag[i] == 0) {
                 all_done = false;
             }
         }
diff --git a/src/_path_wrapper.cpp b/src/_path_wrapper.cpp
index 4a4fc2db1eb3..45c1e288e86c 100644
--- a/src/_path_wrapper.cpp
+++ b/src/_path_wrapper.cpp
@@ -93,7 +93,7 @@ static PyObject *Py_points_in_path(PyObject *self, PyObject *args, PyObject *kwd
     }
 
     npy_intp dims[] = { (npy_intp)points.size() };
-    numpy::array_view<bool, 1> results(dims);
+    numpy::array_view<uint8_t, 1> results(dims);
 
     CALL_CPP("points_in_path", (points_in_path(points, r, path, trans, results)));
 
@@ -152,7 +152,7 @@ static PyObject *Py_points_on_path(PyObject *self, PyObject *args, PyObject *kwd
     }
 
     npy_intp dims[] = { (npy_intp)points.size() };
-    numpy::array_view<bool, 1> results(dims);
+    numpy::array_view<uint8_t, 1> results(dims);
 
     CALL_CPP("points_on_path", (points_on_path(points, r, path, trans, results)));
 

From 2c5d2927688870f1f3cdfdcb06f85dd3e24a82ba Mon Sep 17 00:00:00 2001
From: Michael Droettboom <mdboom@gmail.com>
Date: Thu, 19 May 2016 08:16:44 -0400
Subject: [PATCH 7/9] Use ( , ) instead of [] for indexing

Fix #6444
---
 src/_path.h | 29 ++++++++++++++---------------
 1 file changed, 14 insertions(+), 15 deletions(-)

diff --git a/src/_path.h b/src/_path.h
index 8ba732509366..b8cee3629410 100644
--- a/src/_path.h
+++ b/src/_path.h
@@ -18,6 +18,7 @@
 
 #include "path_converters.h"
 #include "_backend_agg_basic_types.h"
+#include "numpy_cpp.h"
 
 struct XY
 {
@@ -112,7 +113,7 @@ void point_in_path_impl(PointArray &points, PathIterator &path, ResultArray &ins
         sy = vty0 = vty1 = y;
 
         for (i = 0; i < n; ++i) {
-            ty = points[i][1];
+            ty = points(i, 1);
 
             if (std::isfinite(ty)) {
                 // get test bit for above/below X axis
@@ -135,8 +136,8 @@ void point_in_path_impl(PointArray &points, PathIterator &path, ResultArray &ins
             }
 
             for (i = 0; i < n; ++i) {
-                tx = points[i][0];
-                ty = points[i][1];
+                tx = points(i, 0);
+                ty = points(i, 1);
 
                 if (!(std::isfinite(tx) && std::isfinite(ty))) {
                     continue;
@@ -183,8 +184,8 @@ void point_in_path_impl(PointArray &points, PathIterator &path, ResultArray &ins
 
         all_done = true;
         for (i = 0; i < n; ++i) {
-            tx = points[i][0];
-            ty = points[i][1];
+            tx = points(i, 0);
+            ty = points(i, 1);
 
             if (!(std::isfinite(tx) && std::isfinite(ty))) {
                 continue;
@@ -242,11 +243,10 @@ template <class PathIterator>
 inline bool point_in_path(
     double x, double y, const double r, PathIterator &path, agg::trans_affine &trans)
 {
-    std::vector<double> point;
-    std::vector<std::vector<double> > points;
-    point.push_back(x);
-    point.push_back(y);
-    points.push_back(point);
+    npy_intp shape[] = {1, 2};
+    numpy::array_view<double, 2> points(shape);
+    points(0, 0) = x;
+    points(0, 1) = y;
 
     int result[1];
     result[0] = 0;
@@ -285,11 +285,10 @@ template <class PathIterator>
 inline bool point_on_path(
     double x, double y, const double r, PathIterator &path, agg::trans_affine &trans)
 {
-    std::vector<double> point;
-    std::vector<std::vector<double> > points;
-    point.push_back(x);
-    point.push_back(y);
-    points.push_back(point);
+    npy_intp shape[] = {1, 2};
+    numpy::array_view<double, 2> points(shape);
+    points(0, 0) = x;
+    points(0, 1) = y;
 
     int result[1];
     result[0] = 0;

From 7976970830a10e5572e114fb1da2d4458ebd01c9 Mon Sep 17 00:00:00 2001
From: Michael Droettboom <mdboom@gmail.com>
Date: Thu, 19 May 2016 08:18:40 -0400
Subject: [PATCH 8/9] Rename operator[] to subarray...

Since we now know the performance is so bad, this should discourage its
use and flag it as something special.

Also, use () in more places
---
 src/_backend_agg.h | 26 +++++++++++++-------------
 src/_image.h       | 28 ++++++++++++----------------
 src/_path.h        | 31 +++++++++++++++----------------
 src/_png.cpp       |  2 +-
 src/numpy_cpp.h    |  4 ++--
 5 files changed, 43 insertions(+), 48 deletions(-)

diff --git a/src/_backend_agg.h b/src/_backend_agg.h
index 4816cd826f60..2467084eb5cc 100644
--- a/src/_backend_agg.h
+++ b/src/_backend_agg.h
@@ -961,13 +961,13 @@ inline void RendererAgg::_draw_path_collection_generic(GCAgg &gc,
         typename PathGenerator::path_iterator path = path_generator(i);
 
         if (Ntransforms) {
-            typename TransformArray::sub_t subtrans = transforms[i % Ntransforms];
-            trans = agg::trans_affine(subtrans(0, 0),
-                                      subtrans(1, 0),
-                                      subtrans(0, 1),
-                                      subtrans(1, 1),
-                                      subtrans(0, 2),
-                                      subtrans(1, 2));
+            int it = i % Ntransforms;
+            trans = agg::trans_affine(transforms(it, 0, 0),
+                                      transforms(it, 1, 0),
+                                      transforms(it, 0, 1),
+                                      transforms(it, 1, 1),
+                                      transforms(it, 0, 2),
+                                      transforms(it, 1, 2));
             trans *= master_transform;
         } else {
             trans = master_transform;
@@ -989,13 +989,13 @@ inline void RendererAgg::_draw_path_collection_generic(GCAgg &gc,
         trans *= agg::trans_affine_translation(0.0, (double)height);
 
         if (Nfacecolors) {
-            typename ColorArray::sub_t facecolor = facecolors[i % Nfacecolors];
-            face.second = agg::rgba(facecolor(0), facecolor(1), facecolor(2), facecolor(3));
+            int ic = i % Nfacecolors;
+            face.second = agg::rgba(facecolors(ic, 0), facecolors(ic, 1), facecolors(ic, 2), facecolors(ic, 3));
         }
 
         if (Nedgecolors) {
-            typename ColorArray::sub_t edgecolor = edgecolors[i % Nedgecolors];
-            gc.color = agg::rgba(edgecolor(0), edgecolor(1), edgecolor(2), edgecolor(3));
+            int ic = i % Nedgecolors;
+            gc.color = agg::rgba(edgecolors(ic, 0), edgecolors(ic, 1), edgecolors(ic, 2), edgecolors(ic, 3));
 
             if (Nlinewidths) {
                 gc.linewidth = linewidths(i % Nlinewidths);
@@ -1274,8 +1274,8 @@ inline void RendererAgg::draw_gouraud_triangles(GCAgg &gc,
     bool has_clippath = render_clippath(gc.clippath.path, gc.clippath.trans);
 
     for (int i = 0; i < points.dim(0); ++i) {
-        typename PointArray::sub_t point = points[i];
-        typename ColorArray::sub_t color = colors[i];
+        typename PointArray::sub_t point = points.subarray(i);
+        typename ColorArray::sub_t color = colors.subarray(i);
 
         _draw_gouraud_triangle(point, color, trans, has_clippath);
     }
diff --git a/src/_image.h b/src/_image.h
index 0385f45fd9aa..8a9bb8626ad4 100644
--- a/src/_image.h
+++ b/src/_image.h
@@ -126,14 +126,12 @@ Image *from_color_array(ArrayType &array, bool isoutput)
 
     for (size_t rownum = 0; rownum < (size_t)array.dim(0); rownum++) {
         for (size_t colnum = 0; colnum < (size_t)array.dim(1); colnum++) {
-            typename ArrayType::sub_t::sub_t color = array[rownum][colnum];
-
-            r = color(0);
-            g = color(1);
-            b = color(2);
+            r = array(rownum, colnum, 0);
+            g = array(rownum, colnum, 1);
+            b = array(rownum, colnum, 2);
 
             if (rgba) {
-                alpha = color(3);
+                alpha = array(rownum, colnum, 3);
             }
 
             *buffer++ = int(255 * r);     // red
@@ -164,13 +162,12 @@ Image *frombyte(ArrayType &array, bool isoutput)
 
     for (size_t rownum = 0; rownum < (size_t)array.dim(0); rownum++) {
         for (size_t colnum = 0; colnum < (size_t)array.dim(1); colnum++) {
-            typename ArrayType::sub_t::sub_t color = array[rownum][colnum];
-            r = color(0);
-            g = color(1);
-            b = color(2);
+            r = array(rownum, colnum, 0);
+            g = array(rownum, colnum, 1);
+            b = array(rownum, colnum, 2);
 
             if (rgba) {
-                alpha = color(3);
+                alpha = array(rownum, colnum, 3);
             }
 
             *buffer++ = r;     // red
@@ -295,13 +292,12 @@ Image *pcolor(CoordinateArray &x,
                 a10 = (1.0 - alpha) * beta;
                 a11 = 1.0 - a00 - a01 - a10;
 
-                typename ColorArray::sub_t::sub_t start00 = d[rowstart[i]][colstart[j]];
-                typename ColorArray::sub_t::sub_t start01 = d[rowstart[i]][colstart[j] + 1];
-                typename ColorArray::sub_t::sub_t start10 = d[rowstart[i] + 1][colstart[j]];
-                typename ColorArray::sub_t::sub_t start11 = d[rowstart[i] + 1][colstart[j] + 1];
                 for (size_t k = 0; k < 4; ++k) {
                     position[k] =
-                        start00(k) * a00 + start01(k) * a01 + start10(k) * a10 + start11(k) * a11;
+                        d(rowstart[i], colstart[j], k) * a00 +
+                        d(rowstart[i], colstart[j] + 1, k) * a01 +
+                        d(rowstart[i] + 1, colstart[j], k) * a10 +
+                        d(rowstart[i] + 1, colstart[j] + 1, k) * a11;
                 }
                 position += 4;
             }
diff --git a/src/_path.h b/src/_path.h
index b8cee3629410..71da86b0f42c 100644
--- a/src/_path.h
+++ b/src/_path.h
@@ -384,13 +384,13 @@ void get_path_collection_extents(agg::trans_affine &master_transform,
     for (i = 0; i < N; ++i) {
         typename PathGenerator::path_iterator path(paths(i % Npaths));
         if (Ntransforms) {
-            typename TransformArray::sub_t subtrans = transforms[i % Ntransforms];
-            trans = agg::trans_affine(subtrans(0, 0),
-                                      subtrans(1, 0),
-                                      subtrans(0, 1),
-                                      subtrans(1, 1),
-                                      subtrans(0, 2),
-                                      subtrans(1, 2));
+            size_t ti = i % Ntransforms;
+            trans = agg::trans_affine(transforms(ti, 0, 0),
+                                      transforms(ti, 1, 0),
+                                      transforms(ti, 0, 1),
+                                      transforms(ti, 1, 1),
+                                      transforms(ti, 0, 2),
+                                      transforms(ti, 1, 2));
         } else {
             trans = master_transform;
         }
@@ -436,13 +436,13 @@ void point_in_path_collection(double x,
         typename PathGenerator::path_iterator path = paths(i % Npaths);
 
         if (Ntransforms) {
-            typename TransformArray::sub_t subtrans = transforms[i % Ntransforms];
-            trans = agg::trans_affine(subtrans(0, 0),
-                                      subtrans(1, 0),
-                                      subtrans(0, 1),
-                                      subtrans(1, 1),
-                                      subtrans(0, 2),
-                                      subtrans(1, 2));
+            size_t ti = i % Ntransforms;
+            trans = agg::trans_affine(transforms(ti, 0, 0),
+                                      transforms(ti, 1, 0),
+                                      transforms(ti, 0, 1),
+                                      transforms(ti, 1, 1),
+                                      transforms(ti, 0, 2),
+                                      transforms(ti, 1, 2));
             trans *= master_transform;
         } else {
             trans = master_transform;
@@ -770,8 +770,7 @@ int count_bboxes_overlapping_bbox(agg::rect_d &a, BBoxArray &bboxes)
 
     size_t num_bboxes = bboxes.size();
     for (size_t i = 0; i < num_bboxes; ++i) {
-        typename BBoxArray::sub_t bbox_b = bboxes[i];
-        b = agg::rect_d(bbox_b(0, 0), bbox_b(0, 1), bbox_b(1, 0), bbox_b(1, 1));
+        b = agg::rect_d(bboxes(i, 0, 0), bboxes(i, 0, 1), bboxes(i, 1, 0), bboxes(i, 1, 1));
 
         if (b.x2 < b.x1) {
             std::swap(b.x1, b.x2);
diff --git a/src/_png.cpp b/src/_png.cpp
index 05ac9b28960c..974b47d5e56a 100644
--- a/src/_png.cpp
+++ b/src/_png.cpp
@@ -91,7 +91,7 @@ static PyObject *Py_write_png(PyObject *self, PyObject *args, PyObject *kwds)
     int channels = buffer.dim(2);
     std::vector<png_bytep> row_pointers(height);
     for (png_uint_32 row = 0; row < (png_uint_32)height; ++row) {
-        row_pointers[row] = (png_bytep)buffer[row].data();
+        row_pointers[row] = (png_bytep)&buffer(row, 0, 0);
     }
 
     FILE *fp = NULL;
diff --git a/src/numpy_cpp.h b/src/numpy_cpp.h
index cbf8907f7cf0..887ffe48aee7 100644
--- a/src/numpy_cpp.h
+++ b/src/numpy_cpp.h
@@ -284,7 +284,7 @@ class array_view_accessors<AV, T, 2>
                                             self->m_strides[1] * j);
     }
 
-    sub_t operator[](npy_intp i) const
+    sub_t subarray(npy_intp i) const
     {
         const AVC *self = static_cast<const AVC *>(this);
 
@@ -318,7 +318,7 @@ class array_view_accessors<AV, T, 3>
                                             self->m_strides[1] * j + self->m_strides[2] * k);
     }
 
-    sub_t operator[](npy_intp i) const
+    sub_t subarray(npy_intp i) const
     {
         const AVC *self = static_cast<const AVC *>(this);
 

From aef86ada619fe81e3b2d665ef3c286b509a9a5be Mon Sep 17 00:00:00 2001
From: Michael Droettboom <mdboom@gmail.com>
Date: Thu, 19 May 2016 08:19:19 -0400
Subject: [PATCH 9/9] Don't use contouring on the path unless necessary

---
 src/_path.h | 11 +++++++----
 1 file changed, 7 insertions(+), 4 deletions(-)

diff --git a/src/_path.h b/src/_path.h
index 71da86b0f42c..bb730be71ff2 100644
--- a/src/_path.h
+++ b/src/_path.h
@@ -233,10 +233,13 @@ inline void points_in_path(PointArray &points,
     transformed_path_t trans_path(path, trans);
     no_nans_t no_nans_path(trans_path, true, path.has_curves());
     curve_t curved_path(no_nans_path);
-    contour_t contoured_path(curved_path);
-    contoured_path.width(r);
-
-    point_in_path_impl(points, contoured_path, result);
+    if (r != 0.0) {
+        contour_t contoured_path(curved_path);
+        contoured_path.width(r);
+        point_in_path_impl(points, contoured_path, result);
+    } else {
+        point_in_path_impl(points, curved_path, result);
+    }
 }
 
 template <class PathIterator>