From 82f888f89cc33a8446321ee23d5af0a11b0e0400 Mon Sep 17 00:00:00 2001 From: cdguarnizo Date: Tue, 5 Jan 2016 15:15:09 -0500 Subject: [PATCH 1/8] Changed add_paraters to link_paramters --- GPy/kern/src/ODE_UYC.py | 2 +- GPy/kern/src/ODE_st.py | 2 +- GPy/kern/src/ODE_t.py | 2 +- 3 files changed, 3 insertions(+), 3 deletions(-) diff --git a/GPy/kern/src/ODE_UYC.py b/GPy/kern/src/ODE_UYC.py index d02eb1d9..57c41767 100644 --- a/GPy/kern/src/ODE_UYC.py +++ b/GPy/kern/src/ODE_UYC.py @@ -18,7 +18,7 @@ class ODE_UYC(Kern): self.lengthscale_U = Param('lengthscale_U', lengthscale_U, Logexp()) self.ubias = Param('ubias', ubias, Logexp()) - self.add_parameters(self.variance_Y, self.variance_U, self.lengthscale_Y, self.lengthscale_U, self.ubias) + self.link_parameters(self.variance_Y, self.variance_U, self.lengthscale_Y, self.lengthscale_U, self.ubias) def K(self, X, X2=None): # model : a * dy/dt + b * y = U diff --git a/GPy/kern/src/ODE_st.py b/GPy/kern/src/ODE_st.py index f9d4e684..0b4fecae 100644 --- a/GPy/kern/src/ODE_st.py +++ b/GPy/kern/src/ODE_st.py @@ -38,7 +38,7 @@ class ODE_st(Kern): self.b = Param('b', b, Logexp()) self.c = Param('c', c, Logexp()) - self.add_parameters(self.a, self.b, self.c, self.variance_Yt, self.variance_Yx, self.lengthscale_Yt,self.lengthscale_Yx) + self.link_parameters(self.a, self.b, self.c, self.variance_Yt, self.variance_Yx, self.lengthscale_Yt,self.lengthscale_Yx) def K(self, X, X2=None): diff --git a/GPy/kern/src/ODE_t.py b/GPy/kern/src/ODE_t.py index ffd349ec..d5dae665 100644 --- a/GPy/kern/src/ODE_t.py +++ b/GPy/kern/src/ODE_t.py @@ -17,7 +17,7 @@ class ODE_t(Kern): self.a= Param('a', a, Logexp()) self.c = Param('c', c, Logexp()) self.ubias = Param('ubias', ubias, Logexp()) - self.add_parameters(self.a, self.c, self.variance_Yt, self.lengthscale_Yt,self.ubias) + self.link_parameters(self.a, self.c, self.variance_Yt, self.lengthscale_Yt,self.ubias) def K(self, X, X2=None): """Compute the covariance matrix between X and X2.""" From d8b5a72ea824f4061344eaba3bff4cafe8e1f3d1 Mon Sep 17 00:00:00 2001 From: Max Zwiessele Date: Mon, 11 Jan 2016 16:13:13 +0000 Subject: [PATCH 2/8] [plotting] was failing on some 3 dimensional plots (latent) --- GPy/kern/__init__.py | 2 +- GPy/plotting/gpy_plot/latent_plots.py | 32 +++++---- GPy/plotting/matplot_dep/plot_definitions.py | 70 ++++++++++---------- 3 files changed, 51 insertions(+), 53 deletions(-) diff --git a/GPy/kern/__init__.py b/GPy/kern/__init__.py index f8f7d016..e2990f99 100644 --- a/GPy/kern/__init__.py +++ b/GPy/kern/__init__.py @@ -28,4 +28,4 @@ from .src.trunclinear import TruncLinear,TruncLinear_inf from .src.splitKern import SplitKern,DEtime from .src.splitKern import DEtime as DiffGenomeKern from .src.spline import Spline -from .src.basis_funcs import LinearSlopeBasisFuncKernel, BasisFuncKernel, ChangePointBasisFuncKernel, DomainKernel \ No newline at end of file +from .src.basis_funcs import LogisticBasisFuncKernel, LinearSlopeBasisFuncKernel, BasisFuncKernel, ChangePointBasisFuncKernel, DomainKernel \ No newline at end of file diff --git a/GPy/plotting/gpy_plot/latent_plots.py b/GPy/plotting/gpy_plot/latent_plots.py index 240f35ae..976641b2 100644 --- a/GPy/plotting/gpy_plot/latent_plots.py +++ b/GPy/plotting/gpy_plot/latent_plots.py @@ -50,6 +50,19 @@ def _wait_for_updates(view, updates): # No updateable view: pass +def _new_canvas(self, projection, kwargs, which_indices): + input_1, input_2, input_3 = sig_dims = self.get_most_significant_input_dimensions(which_indices) + + if input_3 is None: + zlabel = None + else: + zlabel = 'latent dimension %i' % input_3 + if 'color' not in kwargs: + kwargs['color'] = 'white' + canvas, kwargs = pl().new_canvas(projection=projection, xlabel='latent dimension %i' % input_1, + ylabel='latent dimension %i' % input_2, + zlabel=zlabel, **kwargs) + return canvas, projection, kwargs, sig_dims def _plot_latent_scatter(canvas, X, visible_dims, labels, marker, num_samples, projection='2d', **kwargs): from .. import Tango @@ -85,12 +98,8 @@ def plot_latent_scatter(self, labels=None, :param str marker: markers to use - cycle if more labels then markers are given :param kwargs: the kwargs for the scatter plots """ - input_1, input_2, input_3 = sig_dims = self.get_most_significant_input_dimensions(which_indices) + canvas, projection, kwargs, sig_dims = _new_canvas(self, projection, kwargs, which_indices) - canvas, kwargs = pl().new_canvas(projection=projection, - xlabel='latent dimension %i' % input_1, - ylabel='latent dimension %i' % input_2, - zlabel='latent dimension %i' % input_3, **kwargs) X, _, _ = get_x_y_var(self) if labels is None: labels = np.ones(self.num_data) @@ -101,8 +110,6 @@ def plot_latent_scatter(self, labels=None, return pl().add_to_canvas(canvas, dict(scatter=scatters), legend=legend) - - def plot_latent_inducing(self, which_indices=None, legend=False, @@ -122,17 +129,8 @@ def plot_latent_inducing(self, :param str marker: markers to use - cycle if more labels then markers are given :param kwargs: the kwargs for the scatter plots """ - input_1, input_2, input_3 = sig_dims = self.get_most_significant_input_dimensions(which_indices) - if input_3 is None: zlabel=None - else: zlabel = 'latent dimension %i' % input_3 - + canvas, projection, kwargs, sig_dims = _new_canvas(self, projection, kwargs, which_indices) - if 'color' not in kwargs: - kwargs['color'] = 'white' - canvas, kwargs = pl().new_canvas(projection=projection, - xlabel='latent dimension %i' % input_1, - ylabel='latent dimension %i' % input_2, - zlabel=zlabel, **kwargs) Z = self.Z.values labels = np.array(['inducing'] * Z.shape[0]) scatters = _plot_latent_scatter(canvas, Z, sig_dims, labels, marker, num_samples, projection=projection, **kwargs) diff --git a/GPy/plotting/matplot_dep/plot_definitions.py b/GPy/plotting/matplot_dep/plot_definitions.py index a33e6bbe..9eb9efb0 100644 --- a/GPy/plotting/matplot_dep/plot_definitions.py +++ b/GPy/plotting/matplot_dep/plot_definitions.py @@ -1,21 +1,21 @@ #=============================================================================== # Copyright (c) 2015, Max Zwiessele # All rights reserved. -# +# # Redistribution and use in source and binary forms, with or without # modification, are permitted provided that the following conditions are met: -# +# # * Redistributions of source code must retain the above copyright notice, this # list of conditions and the following disclaimer. -# +# # * Redistributions in binary form must reproduce the above copyright notice, # this list of conditions and the following disclaimer in the documentation # and/or other materials provided with the distribution. -# +# # * Neither the name of GPy.plotting.matplot_dep.plot_definitions nor the names of its # contributors may be used to endorse or promote products derived from # this software without specific prior written permission. -# +# # THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" # AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE # IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE @@ -41,14 +41,14 @@ class MatplotlibPlots(AbstractPlottingLibrary): def __init__(self): super(MatplotlibPlots, self).__init__() self._defaults = defaults.__dict__ - + def figure(self, rows=1, cols=1, **kwargs): fig = plt.figure(**kwargs) fig.rows = rows fig.cols = cols return fig - - def new_canvas(self, figure=None, col=1, row=1, projection='2d', xlabel=None, ylabel=None, zlabel=None, title=None, xlim=None, ylim=None, zlim=None, **kwargs): + + def new_canvas(self, figure=None, row=1, col=1, projection='2d', xlabel=None, ylabel=None, zlabel=None, title=None, xlim=None, ylim=None, zlim=None, **kwargs): if projection == '3d': from mpl_toolkits.mplot3d import Axes3D elif projection == '2d': @@ -64,10 +64,10 @@ class MatplotlibPlots(AbstractPlottingLibrary): fig = self.figure(figsize=kwargs.pop('figsize')) else: fig = self.figure() - + #if hasattr(fig, 'rows') and hasattr(fig, 'cols'): ax = fig.add_subplot(fig.rows, fig.cols, (col,row), projection=projection) - + if xlim is not None: ax.set_xlim(xlim) if ylim is not None: ax.set_ylim(ylim) if xlabel is not None: ax.set_xlabel(xlabel) @@ -77,7 +77,7 @@ class MatplotlibPlots(AbstractPlottingLibrary): if zlim is not None: ax.set_zlim(zlim) if zlabel is not None: ax.set_zlabel(zlabel) return ax, kwargs - + def add_to_canvas(self, ax, plots, legend=False, title=None, **kwargs): ax.autoscale_view() fontdict=dict(family='sans-serif', weight='light', size=9) @@ -88,18 +88,18 @@ class MatplotlibPlots(AbstractPlottingLibrary): legend_ontop(ax, ncol=legend, fontdict=fontdict) if title is not None: ax.figure.suptitle(title) return ax - + def show_canvas(self, ax, tight_layout=False, **kwargs): if tight_layout: ax.figure.tight_layout() ax.figure.canvas.draw() return ax.figure - + def scatter(self, ax, X, Y, Z=None, color=Tango.colorsHex['mediumBlue'], label=None, marker='o', **kwargs): if Z is not None: return ax.scatter(X, Y, c=color, zs=Z, label=label, marker=marker, **kwargs) return ax.scatter(X, Y, c=color, label=label, marker=marker, **kwargs) - + def plot(self, ax, X, Y, Z=None, color=None, label=None, **kwargs): if Z is not None: return ax.plot(X, Y, color=color, zs=Z, label=label, **kwargs) @@ -122,23 +122,23 @@ class MatplotlibPlots(AbstractPlottingLibrary): if 'align' not in kwargs: kwargs['align'] = 'center' return ax.bar(left=x, height=height, width=width, - bottom=bottom, label=label, color=color, + bottom=bottom, label=label, color=color, **kwargs) - + def xerrorbar(self, ax, X, Y, error, color=Tango.colorsHex['mediumBlue'], label=None, **kwargs): if not('linestyle' in kwargs or 'ls' in kwargs): kwargs['ls'] = 'none' #if Z is not None: # return ax.errorbar(X, Y, Z, xerr=error, ecolor=color, label=label, **kwargs) return ax.errorbar(X, Y, xerr=error, ecolor=color, label=label, **kwargs) - + def yerrorbar(self, ax, X, Y, error, color=Tango.colorsHex['mediumBlue'], label=None, **kwargs): if not('linestyle' in kwargs or 'ls' in kwargs): kwargs['ls'] = 'none' #if Z is not None: # return ax.errorbar(X, Y, Z, yerr=error, ecolor=color, label=label, **kwargs) return ax.errorbar(X, Y, yerr=error, ecolor=color, label=label, **kwargs) - + def imshow(self, ax, X, extent=None, label=None, vmin=None, vmax=None, **imshow_kwargs): if 'origin' not in imshow_kwargs: imshow_kwargs['origin'] = 'lower' @@ -178,7 +178,7 @@ class MatplotlibPlots(AbstractPlottingLibrary): if 'origin' not in imshow_kwargs: imshow_kwargs['origin'] = 'lower' return ImAnnotateController(ax, plot_function, extent, resolution=resolution, imshow_kwargs=imshow_kwargs or {}, **annotation_kwargs) - + def contour(self, ax, X, Y, C, levels=20, label=None, **kwargs): return ax.contour(X, Y, C, levels=np.linspace(C.min(), C.max(), levels), label=label, **kwargs) @@ -191,13 +191,13 @@ class MatplotlibPlots(AbstractPlottingLibrary): def fill_gradient(self, canvas, X, percentiles, color=Tango.colorsHex['mediumBlue'], label=None, **kwargs): ax = canvas plots = [] - + if 'edgecolors' not in kwargs: kwargs['edgecolors'] = 'none' - + if 'facecolors' in kwargs: color = kwargs.pop('facecolors') - + if 'array' in kwargs: array = kwargs.pop('array') else: @@ -231,8 +231,8 @@ class MatplotlibPlots(AbstractPlottingLibrary): # pass a, b = tee(iterable) next(b, None) - return zip(a, b) - + return zip(a, b) + polycol = [] for y1, y2 in pairwise(percentiles): import matplotlib.mlab as mlab @@ -244,51 +244,51 @@ class MatplotlibPlots(AbstractPlottingLibrary): x = ma.masked_invalid(ax.convert_xunits(X)) y1 = ma.masked_invalid(ax.convert_yunits(y1)) y2 = ma.masked_invalid(ax.convert_yunits(y2)) - + if y1.ndim == 0: y1 = np.ones_like(x) * y1 if y2.ndim == 0: y2 = np.ones_like(x) * y2 - + if where is None: where = np.ones(len(x), np.bool) else: where = np.asarray(where, np.bool) - + if not (x.shape == y1.shape == y2.shape == where.shape): raise ValueError("Argument dimensions are incompatible") - + from functools import reduce mask = reduce(ma.mask_or, [ma.getmask(a) for a in (x, y1, y2)]) if mask is not ma.nomask: where &= ~mask - + polys = [] for ind0, ind1 in mlab.contiguous_regions(where): xslice = x[ind0:ind1] y1slice = y1[ind0:ind1] y2slice = y2[ind0:ind1] - + if not len(xslice): continue - + N = len(xslice) p = np.zeros((2 * N + 2, 2), np.float) - + # the purpose of the next two lines is for when y2 is a # scalar like 0 and we want the fill to go all the way # down to 0 even if none of the y1 sample points do start = xslice[0], y2slice[0] end = xslice[-1], y2slice[-1] - + p[0] = start p[N + 1] = end - + p[1:N + 1, 0] = xslice p[1:N + 1, 1] = y1slice p[N + 2:, 0] = xslice[::-1] p[N + 2:, 1] = y2slice[::-1] - + polys.append(p) polycol.extend(polys) from matplotlib.collections import PolyCollection From b2ec4ae765a35cd498b7c96523a61bc3ab4173ea Mon Sep 17 00:00:00 2001 From: Max Zwiessele Date: Mon, 11 Jan 2016 16:44:13 +0000 Subject: [PATCH 3/8] [plotting] baseline updates for 3d plotting --- .../baseline/bayesian_gplvm_latent_3d.png | Bin 8832 -> 9214 bytes GPy/testing/baseline/gplvm_latent_3d.png | Bin 8832 -> 9214 bytes 2 files changed, 0 insertions(+), 0 deletions(-) diff --git a/GPy/testing/baseline/bayesian_gplvm_latent_3d.png b/GPy/testing/baseline/bayesian_gplvm_latent_3d.png index 795e89f94e15b3f436a01cdf2ab06d52bb4ea053..ea0009f26d9de58c086270fb6eaf3a5d7f324437 100644 GIT binary patch literal 9214 zcmd^lWmJ^y*YBVxC`#xfNH>boAYI~+f~cfYA|O5F(B*>&hyv0fpi)Zr5E25?At5<{ z0z(ZU-Ej8&&pPkdv(`E1%e&63#Ukgvuj{(^{?#V*;REH%7wIk{5QxjFD#*tO1W_b> z@Q|K|pR8bTjqr)sML|`Y6ka~0<{|L?1xJ-9E(ipLDdB_Yi(Ixf0>Oe%Mc&o+_^>kJ zsjt0#(y-o3ghW~lIGZ$`d#EMi%o}}&=n?JrpPvgpy|6jxT&gnaHrVCPO{P`jvVM^{ zQJqkg6_7Rg7F(+sXT=vg7gL#>P|~%3{Jx3NX{1gdY|l@o*OSR|t}AiPBS36Q!t%Pt zFq#PQF~^WW27zeE4-G&dVkjvQ0awnPK|EBDMq3Sy>PWPGVxj$5#RW zn?FK@Gs$cG6zSGNN=kV#AfBBiK^xUiFWlAFx_yStzr`#forDM>`cn(JBk4ssG#D2U z64KV6d7qPorJ=705n$8V-5uMPt}Nws9haj+bb*vK^7rp&$`of2C^7~l$^TkOv#NTM zNPFKsdCsD=56p*cem1`ZBa<@l)pirrI_EE4XDKKs2xkh{=t{bCzre7HjVoI9X}JyO zOlw48-?nWdh1?Mf6BAQ^HQF{=E>2j|_n?f&s9KwxL7XK;$~`VS&9A!gQ2e=E?9PPy zlbT+x-RxBJ~uN{L~`7@qf@s(N543bt@$p` zgtq?lq-GX}Yd3PA(3z7|N`JIkIYZvw)fE#E5HQZiak#q{-`tFZ3&r3%F>Rb{Ye?J;lI^I3lT{AI@z~}uN7|6-VnK`?$*(=9i>%AK%YCooL zG{tn-W{j#KEOxZszh1BvQl52D_X>EgHdOlapH<$2h!eo=n^s_;@!EiJYBJ^Dh`2u`f+r}T`>%n8@}2-yRbsOV_Un>W8S7F#8`HHj6ZrKK5u^C)f^ zD$EOKlv3-AbS$zPD^rC-)?0}bAMne$Ww3Ii6xR9%-nmfC zH{X+TRZgy+9wnBj9+j5HZ)|MbpRLK<-rim;Nj6$)8G&`_l__ywn@sjY=)D$ROP{E5 zpAul_?jAa9OSSQp&kJf=ZhOz}yfStt(=q0x<*SOHU)|o`p7!23ezxhkIW7qai4yN! zl#`Rw`|0WD``g1N_v;o4jcaew@EdcYjSQpVR~yLEu;@w9Fd{Mb_bp*e4Tepunictw|D&J%;msTu6O=JyJ>OR_ndlUWW>R58)qGsh zUuQRpv3^CapvY)?{{ntRUdmoq)n}8xDoh;^HT%ytG za7j;3Z+&Yk25L8V*ZcQN67!X@a?bd8UIVW6^>qu`2Rg%gI@;PC1{HSiyOYK4uT9o+ z3J4Uo^H0|K)c}=8$rKGhSemr+Eyv zDNROI<+Nk^A~i@mu5=|GjUQ@Xov6WO#vA7h?yYRsnp2eTaPqvDGdX1tcl`LzwQKDa z_Tw2vMGpzI1%)k^Z%f`@ZaZu~USZFZIdiVmJ!x4sw>q%J07k3(M3<7VWQb zUyB-t{oh($UFEEQq+8~Pj>SX}%7Q>n)Vxn)Jm(@f)BA2t`5op1uarH;$Iq!BE&@~v zxvv^lIsW+L)L&>KM^!xKS0cS*%FV%n6cQHByc>AlB302mHRvB2nocN{ z481~QG3+Fi%P7XdBsM%qyi>NLvorduR(8iiUpnpwjR{W`pw7+)a7I|jaee%}&4KL5;7kC4_J*oGrZEkK3Ngh9y^j`SMy1(Aad|#o& zqBBki^Kyd`SO$y5&VmHECAL(@Py_3fvjppZ9N z2OKX*_YRjD;=I3M;pE~H4RR&7mH$DCBxn86Mp!@Xlk(Uf#h+5E3pVGmjQLw{Grpa5 zf})81bT4?))fccbLo53Mr*PD_Z_)S`1pdt2{Kc6PR5ozH&OIDWX0|Hgw41mTpq8`L0( z%B2nAB|s^i7|L$AwfWr$NslJ5d#$&75FSW4EWoh!(yV9$1z)Qfv$&kb07^WNto&ZS z3>g{XnZbXlzQG(b7F?|=DFVfbY3wb1il&ZBZv`#)yUfz`S6qO{8x>3Ch_eX0-yaL` zv8Py4shn#d9Y>f+`iTCIl;nTi`2Sjy{sWHXW3`yF@|*t@0W$B)Za;}!We>8>eS=mqMy);JYK{`ubj&P={4{3>s$t`#pUdWvum5#rq z6*j(e{ZefW8b8O)7L}ZgFDf2OvQUAt{vZGnOPS)r&G;;Qb$V;SxHhb|^6%fj?=zd6 zpvf`p%He|!6V>J3^70K&#Kgp+_bW8zjzz=-1qBCL&z^TYOa8gZln>PN9-viys;qCd z3UU^F$rz}xg?xOqgR5(qvb;RPBrV|T3@49xSB_3@-7<*IcKGb7{Z-Y|MtHJXc4g195C+XlPhJly9)PeNmhWcEg$^=CHc^orKx$A_awm^TurZ^g8$- zpaqRMoaCoZpQudy_ASB-@bQJpxOKUcz4X!2URkN1cbh7zYF4XrSBlMA!|~Oq91r`3 zVp39u=|sy^&&Bd#V{hyaT(91vf&TszVW;Wl_tz%)Sy)&E#Kd+3x=9vT_`m&almhU# zPDxCx5@Bbb(NZ{L*&;|eVzfEeeGKyZdG|Uc6|Jb9e*MVENXVO2gX9@jwrM=xK;T7J z0@LE+Vkr&p(=skr)<3Ko1ADSTTGM1#-NbALzi`8%#eAlXp-v6i+1XF0zX(dugZU8* zo`v%>reS21`E-kwaOJY^&jhTJS`zJwXwRIfDG@ys$X(`{b%UGSe)Q%lMYMN!zw)boL?;;=ml=&C?Mdt|x%zIx|XKtlmAq^H-x z)}li^+8cB^qsD+Vc;e^1P7(xZx&WN72(SG=fBxt-H8wV4?8k>(THex{_?+y|J5XH_ ze_Md(;la$p>c<7(1}%LvU)piqVpHnB&=X(q)$7Znr8fx6&oZg^y9xrqwZo^@lberE z6m;8m4X{T&(?9NOXV@mTozDVJ?SL2+C(QJr8o29(1YM%`!I-W!QtmqoD?6F;~=Ai^stDJjIuyTb43 z>8Z#TWnR}kJUnbWT4Hez^yau7h;WL~wvr?-b)#X%o&KK5TJH$l?5w&`osZ{hPllmz zf3Umzz5XXIZqRAAU^DlJ07*1sF-Fnt9Lb{g9i1H=PiviV9fLv``|&!cF< zQ_mC^#TX-zVd3H7YWMFSV*0G4KjEtdWn_Fqagsp9DUYp9{P@w5mzQ_?AX#*U z)8GG8X>L4QGwoy@Tje764({lS6&pDft3_)TQr&)r}}KPa}U861*iUW z{$NqAzIvsIcG^Se>vxa^qc~Z#veZ4lOCJjG>!bw5kOzB8J)n!^f|l0cWef-XgMos) zYPBK@GWyIB3I#5eV-fW_u@dw;g zo;h3Fp`oF^;0yGQX}-V%<7M%b$b-OgI>KlqP%wPSKgK= zkPC;jv`<-<^=W^nI6KrEMd0_M>33q23x9Apo!gUV`ydWE01M1a$H>Tt;^ol!^;wZ@ z7*5>}p6139wrcZzAPEIiZu|OYH&zzMi`4~qwvxvKWjdpyqhha+!)Sm5YeohJpLxit zykk>+6Ldt#`q0NYS=RSzYinD{=|pat?*0990?jSzqr0PZm6aeHx!KvxdRyMn%N;>< zqH8#HLq3s|^uwF-$mW(7*`J_6B>^8w0B*kJ!|B|_pMYqOzdhpjG9<$~X=gP-(98y{ z0?<(Exjy~lYN6b`5Ec0sHWo)UYSz)ehimm~_EUAfonWf30iMsezoi#Lx$EiaEg$T}1bYlE8}(=Fr2E4KMp8b9E12M0#Rd1KLXEGb z`Je3ffB`AzVI)_MAj43+U~zmp)ABA^+`KKaBD1h?iHCz5UIGG+UqY_mdb+C{8XVI0 z@+?Ui5PM#3*W0s2=VTvVe84Ixd4Sd!P@9bT1PjpDA8~YdchAKWsu6;+zWUEAHc4|= zsHsH{-YYX*fO@=;`<+z^OTz5)_or~LbaZ3h+3?0o#N=X*zZ+Iv=DI>1`ab_qrBb~{ z1u)t>5lq26c|4cwuzC~kS65fZ*xbRBVSqv11;ewl%jVxB`x^eBdwkpg)O%}4n-aau z<(k7k!OY>}UY9DnB@Xic^6(Va?u^w`7tqV?&qfy_Gf?OXPBtcC%O30r_}wi^s`wva zG!8!tPpiV=88JtMhm*T(ZCaRqp7od41h{+?-vn_!C*SU+#f>CBq_P>n8CU4&Y>K3*A74w#@_ZE>OaVt2pOBDn*QnanvvGX<@HP51IpbJi z1tS>|Lr_42>3MR75g9J77m+yhyMY(kk9|!_x~w!Hf`OlEIuLpXQb=-}ks{N;t2Z9# z8;RY#X%&-YP3CG;uN%7JWo~XB^UbokyF0(FyL;u|khb6711G>~iJx}z*JCX`({7}C zOux3*f1M7hV(1jhA3#q@InvD@A-$pR?CjhMf|pa~^k_3$SV+jY@HbzqmtG|-ysV+2 z!R!6|_g9e|KgL#e$b+-2w*h25c0n_`*d6N!_2I^T&Q`1CeXhkW#PiH_x(`dLbp%8@mp%+k+jF|t)(WG0d)1J)n)<7C$1M{6A z(X%ZrEPTycThH74L-f)}dZUd!;@x#YHhE+QxNt8I509uUYv#pjW5Wr1Nci^l-wdg0 z)UG>r=$ZI#8_CPd_W)5RiveOAJ@FU%K@HI?9%uy#NJjnv%9&c>UxRCwrr+m8G+gLK ztbVF!_MWCVS_oXP=Hld(S^grpSnYQQX!ilw=!Fc`=##|v@7K-l-mOZj(Vl^eO#l$n z3Bt<#yCDcrg=fcHLyq{u^@62H0?OuzMz4Z4TtWvqXO*JMQ%U+_o5G7+F z{~4_VeJ-O2%54$OJ*FILw6sfs=gDhA8bKudqN~z61yElTG+*&B*n=R8&;St(BFRqoyJ)Ba@tPYR}664-$1*=4`{xo{GW8H!pfVNNhs#!El4^ zv7sRkxP7@QrEs^gET>YFnZu6{g3ZU4|+dhJ8?bpQgHyFkKc=Xz2Ob%1diA;3bR8-5a5bjb(MJr|)1 zn*7=J6Er;=YzP@=YBT= zlTG)^o%$TYl!v+y>{1|VUf4iHBUcU<(SBR6fWZyd5o^n~Q6a<3j3Fo!H8nN56&&u@ zW2ABMwi!@%JJhohq_RqNTmARxLBmMtadLptMfdmjSNzIyL@Fy|Tp=p*KRujMUY0nI z9U~>RD*)P;fq_8{W{W~aW55ss-VZP-na{R{X{w%zii(0Zq_1i!Am9JO zFXC;r#1Yqga37UG0F{Am@;rnv2=&kYCT9+zRm-8)_$c4w;_%RS?|{bF<^Vx7W%p+} zs62o5nCV;$g^TDgyMQ|j*3HvjuM2E^q|g23ShYEVT3gc z&i-tL2Sea@fjr9<1(H3hShI${;()?G!X~#nz9Tf1fns7^yoz`h_an!#c@jXaP zObooEVNef>{L$p(>9hK0s5 z-j|8WFA?2+2l^n9;mVb1PPQwwwEGFc!NG;!pH&(lGhk%(6cCB9{m_sw-10pnkV1Ds zDW&ZnXt0WWW?ec_yDxGF3+YDFGc%SF0vtgKDi6<+lG?oT+%c__wi_w(-rU$2^YZqN z&ofI60<`sb+9TVwEJ)_0Ew1J=6bUoKLtn~^C{JSIJB`jeV36nMhd;vP;GG87fG177 z+*IM%mq+ZpA|fKAO-C-Ciuvn2 z_h684`Um{)j*kbnk(4wie`^^Mi9U3#ahj|z76vAZlZ2SJ_YJ|pPl0|gfffn+B{yZP zt(~E$q?8rXW^hNQ2Q-QtFWdB`KSfZC91yb0o`5*n?F9Fu5l+{-9qF_5rL{FbIN;-$ zTCXhN>rz!s&2ZIbkCX+J__hX@)=N)nA)md!w!mj3a7Mu~Xs*W!3jd}kt?Gjr-T-|2 z(ALvqoER4;syDoSZ)h zuCA^$XZ|%f$8SpH1{1Q?bhcA6WiPJ-cCZLC^mlBUTW|f2L&_thSQ)<9laf9oNy@4M znmWGC%bV{H3olO^nGB&~zcMi~v8G?`x`-Xx9r8jT&c+b_7a-|6`$L_h^TZ=B3-DRM zEu6i8{oRn*pXx&jx1N#*lXgx>N{c$L_o`GASl3k9grQN=p=;+fO8sLn0OcI=m-Es!(U`34)l|F`=IoLU4MNr?;b{ zVG8O?#A<*$c{Y|*dp9>>-#QP0Nw;C zcNFs3GnPLA-4JN1E>s&oTq;pT%Jzw@OhKorkzT%>@2Vo|FnNUaUaPgdbouhTu(xmR z=9ZRhN}XnQ6y*^&vadwp)AIJB0HT(_G^Mg}6rDSB2FsQgGz(rFm)fYKtFkQoe9-NP zm*+P%fhy|j{d2*2o*IrG!qn3vf>wBr8)kQ=K=;=(@CTj-;Z)`HL`0~2Nm+OPd4GH7 yh}39M_CL9@`{y(N5B_c9zfDB_|8bNbbxIOM!WRE6*ZDRAUaE=@kOd0Q{{0_JR&-nd literal 8832 zcmdsdXH-*7yLLbxK~OT?ATP-Wy)DOBB8elU9YJ17>PXhJ#C0~%$3jl!m*`JH1QWfq70Pq38 zTIyyYgFUN^Rs6A4!%fn2{uSEW6EzSXw@ZJc(l}+W142X z@$EAbGmyqK_c{~z2}b**ySeo7{%;t?;$tZ3(k%Yxhs^9;i_6C@52a%;2lu5927EgA z*!J^>PgWnKE-?V2;`|>&XaG^ttTZ%$C;WT>fH)l;;Aw0W08q=q0*JavOAGj*p#}h0 zFfahVGX9?iDP`EJ-DP)L6beiwC;;AFdbtfW^n#x@Kk*)B7jV*l;ytia1ur%U3TnR2 z${L6-&vhql&+*Tx+_*LTwrp{6J>;TOYuadcER(F`cZTwc3NZcyz&FC1VVs8AyDccEC!Y9f!Wwy z=xOq9pv1Dtgpqwa9)?54V+3vgY`Eu&+s;x zhCT(JJqV5uvTzYa<~~`9;6jaV9^AsKMp94^++l`*P5{=cMN~vYM^sdF1A{0D`aQK* zx7F-DS3O#$Q)5u@$bGalRveVcYT<cSPwz9Gs*-T1GYN}s8H_JN8P&p0(r%*ye=Npb+j=4^5CR8jD4oc2# z!;c9W*eweemx}t1L+g(7Gm31zG8anJb|^oJF3#y_`-XXMOZE*Vu5C!nQatHIUD_U3 zt*kRgwC~9hHEj$%aQTCFJTvcP$r(YU7m2*skwX0jpXm($%*e!alb094>GNf8oR}$M zG-($+S1zEFx?Y3mL^omV#leunY1gKpEo-ltMsWO9ZnQY!ODPpRk#%d8qh(g&@9Bea z4E3I|wsst~BQxF8fmbqDZ!T!ZVxHA#{Py)0DO}7y-xJ0((Ebx*TU@I44Hf!s0iaSm zyR4ZhSJ}qA?E#y!?I+9oLzEXIo99wHenJ7*n(vo&d=M0m5wHG%0pda@8v&X4i&B~P zT*Z)TH7031%I;rO{y>mFv8ncI+N9MWG91rcH)UJfPrfXwvN9ehj(J(prYHB@#;<-G0t7Grb@uwB}O*5BX1}axH%?P?d`p{pj93~TX=;<}rP8`4Q`0Kv$yWhV( zon2k^pkM3}+ihz+SGa(?<3xLeCFdj0k7Z0#9YHM+#o7{*Vi*)t>Wus>8Cc z_+Tsgb$BO!PdD5Qo_6-<(z(5|(w`@1q{t>Y4j0hAdtOms#P`hLDuE(+Zbjm)Y171%u;_WMDUE{r4ejlU5d*2wx^V$$QJg}oA;L}n!)8v6r z!4cjn1g}rN3)DLBJEj^|ZMj~qJ8@$|)_ZImC0dxjMYmW2^bT7$m%tVI-{GO5I$K5= zDv%ITt}n4tsp98}_>l$Axl}G38ymZS^s>LdzxoB1>K<<0F_+wNogLIilJT5+ke!`f zbo^J?O132d4=<^Hu{N|Z_qoOX<=$g3S{7qlCTd!-zFOd8u{l{6I2VX(0O`UAv@-|DJw8Y>BVH@6YoH z<=(8WUL>emlnbE(=F%ePDt*aROQ~sWG{d&|AYNcU?{wM-2?<$uK$@=ds?<^=#(*D} ze5wF=W1=S49)cIz!gV?PSyfm+?;RDpzOA6m?aw=t`LZ635vLHU4Vp76?W&Dn*XKJZ zZY1d8@9DC+&%rhC-th@y<>Cg|owsMoOS8qzZ!Gg}#f_VVt2}#kv4fa{RNYUQ-S$yU z&NH?M5)(XsvDH*vSy_1yRQBIib52-Tlb4Sr566<>%q?pKSJfjpiz>CZMf~XWVhPh5 zlcvK&wS6igX2!b1f#1a!07G8(c zf8JJ5K0i4)UN1`v_&p^x7fzbvPo+ni^bd9)DdLie6~l zUx?2ht+%VA%;Z5!HJ6b14SJq*lqspKX*4LA3uO!j_f@H)%gf8fL`AFoP7OknY2nOW zUqxZZgatNWES({5#3NHsBRAZt!3CeHxIHr4_QFG%w7^y$skboVdEpM2y5CGxS`#za zITk`(uh2PL5LI_P=^Oo)w5iLw zt&hCtkaih#8}=xX^~f1gCfN#!h-@BWeK5H5v*14?hKrT+A3S@7Mx(QGbK&FTW^<9} zr|Ou8R1#oNqGGp%OfiK(+W%zgx6-b~?GdYkZ{c&9)gBZi<7E`g(tb33KvIz?jg4_S zCGTyJRtwrUg1byDyuCXx*9=h4_ZE8-)6y`#Y!0>5z<*a#((uxyOYnMTj8wHsIQa;k zt{H#zuOgL`6H^co-!tsO_@-hhV5fsbwJ%Ul+!a$36DYMA76^Br`iK@LM`kHz-gbiM zS-1JG^jqWfKbDqiNIO6@eZpI<929T2`-LP%JUE2D!{+4VY~a`I$5Tv^B|$rZD;q`x zl3dREm>pESvuTEMxKZVR!e>QIWkP;W`E&q&|Ad!VAV$kP zQjj>vq1kXpSd)~{EA#;;4GXOs!uUSjw&?~zTh{_5qhb4MZP)|DX4%gDx* z9i_Udb0Zl&>RO`X&-3xLq8{-P3gKU%lj=UFB zc8N{vhU?mTyDXn2WfuB2vUNPa=_PlbeaY?Kckgr+Km_A2VNXSn8r(Pj7Kwv zd3fg&!2tlZF@C=PCZ6S1MHj5PBEn(0006(A@QJdA#c+4xzj&_JR>OyUJ+=SEd(9&9 ztsAJUnePPn=IZ|oiT=N0+W(=0wD^rYO&~A~V7!+V{}+I)jh0_qfj91idKA=`?F(H= zAwZ&LSXlTN$vvYHyID+&HqA!B51~hKAjG9gCK{SMl@Q>a4xV@X3=AwZF&BQ)h>S)* zkVby*t!VB*JdTgfWnnJtV_+zSN(l=~wBWdrd5s{1&^Q-LzH{dw*79cO`Jil=N|9n$ z<8ucq7=~`@7ymu2n@;-r125v|(-luv#h1{Zh7=oD%5sIcvtY2~p*g3;v3i9MZupLC z_AVf<`NO?1=)9gA#_754=WSK}pWD*ox#KgO_&^%sV>OnWv^^O)*d3)@RLxR+kgmeD z)T|(7Z!gu5{YvO-MJ9nS!jdpRN3)ieZ8Dpd2cl8{c2NC>NS=n{Hn zTTU`K`@H=^K~5b`=KQ1aHLC|vW4_n9ubEm|e0g8m7-IPtY~1?>kye1{dU?M$BQs~F z)6&z{+9pCZ`YuxQNkZOh8w<%Fzi(&$P5vel2_CxFsv-d&PH%0jU2@hry#kTKbM2!AY+MmLet|lB7wKLc@l4PktdV+;a`j>}`{R9CF!s6#{pJK{kbXLU`|eK#YwBQIr@ucY|W}(eo6jZn3=!JH{!>t+nRFsA&YPc zh09;=NqioD>f@*v;Shbd8yuje^=HMTbMerKAvGV!XmNjkaD-jJX~pC61KJkw(C?sU zqeFpu)K{JVujwqlT`|Aud>flP%)#^CS+4D#lYbD%h8tB?9p=IJsZ}RFVg1i$T$vNV zxiDBji=u0Ka5Zd<*8d-rE4=g2NW{P+4mF2kNrFu|qpg8GzD zb5P3X7=L=XT{%AA?h8ShSR7-dEiV{Mu?P7T+qEyN(U{ zT{HgyUnrl?)s^q&3t*mHwHBU%g$-C|RO1=J0b21)3qwdk5Vz6+64 zjDN0=ZSNE@f0(NmzAYs9{7#{k&J*rAo3@=(owU=B6`M)S7QIv*&y5njy4GH5m9b2Q zCF9=yWu-1+OimE@SuA@tf~wTxI0%$sd%P`+&|#CFoOo{CJXB7+3F3B7s))nkasQpkd-hc%7sn_khmYiTK`eo0Ehd`^P}Kk@kbGv9l;#Rc~} z?pAeVY$H00GR95ajYaKD?q-9IhuaUgK*t(0Eq2x&=ErvlHK^&BC_eYwPWrl(Vzwq; znNau8z6?e%>j|b;_%jrB^62Q(V6oI~o7`MG2{)HYZ)qs=)=Pdbd+|AW2gqnw>QlG` z7g)UfRIePv2JXTw0>j*~2m!pN+s;G$L@YY_dSA=1-=SRZa8v_5jr#bd%R27fxafJZ z8L5-VWEht5?9BON@bc9|7xrA^L%@!p(_k%9pC6ewWl^LvJ6nB_z$BF7Mqqx1zk+&e zj3uUAK{aQNwtma&I-g}Hnc(m5$Pay)Iiwl?^AvPUIv=(iGwu4EDnuMyY(=);26Jp} z8n+I5Gw#I;_$+(n!dD4+-;mAx%9HP!I_SiUaZaJ>zT|Q}=Oq}E4(N4+_uzF6T24B; zYnC`)av&!sWoedw{_;7dQ_%7?Jbh!qx6O%z;v`cJT`P~Qg+{QlrzRyIol4@3Jh<0v z5JL|4f7;wpwp8YtXv&zi1Z!TpG^76x71iB>X^z0{}}EQpdfhY>8~)BREOyJmE-kL>Ce38qRHC-}t7T@&>5brj!X zD(_6&H2cBa5`9aW!O&m{GzZh=Lrf2ONw4TiBFzUHMz0M}{RN?7njj_oh#JznW{4#t z16wH`w-C`lBVC|TyXOZPg(o16>XhY95XTn>-61KLrCX?Ju6|-^C_i5UqezLL43=Gb z?MJGS5!0I`pzd=kw2JGymY)mU*=v?1DkmHH8dPa4g|uBgM4msW^1_IZ`5y-bgqB9q zR;B$Nvz%@w>tNkxdyJ*UJA4^o38YRcrmhLaI%)aa^pdaJ&74}{=e>AuDqD?_EOYSB zxheNE3G_~bC5c0v^R=1}C-Y(vx-CT~Z#}etPm`kf#Quu;#~ZR%Yd%kg)n?rbb|(4x z)#1*}M-wIVba6h*XRqS~WopuD5wN9!^2cE?5sxJL1vRH@ELu(5mfznfZTH$H76e1C|tGRHRa>J{XG5!bHDT;z32 zs&r0%b2=~3$E1*w`t#D4likDqPc4z=GF6>afRPg$wxl*nXDzvu2dhoUUca*UL5epu z{S(Lat9)5`YAj9d=X{6FS9`+cC+_>K-#~{Z`)CX(4Ld{YD?PB%LOjq$tA|*a90UEG z*x?K6WrJT_d#z5HF+JGD7%Fs4Mzy9YexW$)SrYZ@g?31%IZ)l6i`G=vyr6pX{|PuE>Tv zcz%+%RIf;%EK-^~47b{Q8{l{hB+)tqJCy3lzy5;JgTUWJB!9qe067IZE*8!)Fx-*! zsSLI*7?SdvJ{3niKq}#ZhWH`bl>wm%_OyI|GfJj)v`JAmBOd>*JM~HM zDUtr6_)|5UjL=g+ZBgYUfPsdZ_xm|FiY7n37O{uU+L!%4n-`}Yj+N2+6`J|w*?dNw zLc7nN5pTqqs34sl{Sbue!K}b=vaAau{?Dx3+GKwr(Yhk|<0M<7;SKDASipNtH5RGp zvG|g$8`pv1z*GU0v29H%&(7J;n>dAG%MwvR+0Ud{YDz6cr^^CdISoq%YG)1;+MZ%& z5@W5_P&aISeb7IiV(Ikxn<5k(-ckLnDoe-2+$8R^Ke2Z46ssaenUrpN2^MxxFV4u$ zpQCgy55C%Lpf|eR7&(TvBA$JVYG9#JpJI1XD%&DT(Jh|dr9xmfQT&ROKDu<%r}1>04ulFZPwYt$A?6{!2*!2;SBqN<)#u%y-j{>ITYH`* zgVzO!PZofB`Nz#s71TU>3&}p*+hxeX%jz01zdrS z8en+mY*rp2zbmGi{_eUbT6S<$H66a)&!gn5mM4;WD-#xZ@;&b_Js3PkehB?w{t=4_ z9o@50V%SGaAO|=0G1PG}B~y7s@f-U=;S%s0Z;yMFk%N4MQ4yy(C84V9a~E0WQPUQ* zEKIK1AVse)u&sw=^uFkZv|0PEB6_!ky&)m}mws=-Tv=$uBk64HP!uzAL`#1Ef?LjBQv>88YLfX0-P?! z1=6SqNJYQv`->{Z-xq6m%BPy%!AGV@88vPwE#@o_o`{=`P_XUw_{E()Pd z2jq2LNdsuqGZ~dg#&dly_%Wbkk0~|_gAeJz?3HKXk1X1k4RZIbb;AXzfk>J5mGDi^ z!P|BPL%T+o$ zgX-$;(L~pzbK5cE3X0pv-^MZyw)M=`qR8;<8uRf?#y!z3ZG1%pdg>B#&MMpdzJe32>flabu%rI_VT_R&&HgfQ^lTh>%9Vd&Hj=A9ausHTFnOAz~ zzeNuW-DX~+UtL0;=zNj_%X;Aq6*^7Dc7)4=(ua8O5Sf?!qtoy3sWH&SEeE;zMN6^$ zg_Q)0@d|MdWH=)A1a_cjhtf{PyzYY^cI(-nD-A^h_&*s{KjD+G5k6gEMTOR=;4=Ll z3iiO>t4gnSIHOB&FA94NIxP5ZFfa^pAla#!m_M@V=~Li*ohz_`Ec)>YSLcBktz0o>tXD5jm0mj3-uCSW=`qJ|oaaS@6v ze01H;3K}ss-ZG(%Hb*mG)McCV9GoS=bi6eD)e=M1AzkvZ0_AqY_m=!Ozpk!u_Ox9r zoMvFiw}oKxRKhTpwbf_zX!fPf^>46CA4ThLA@gphE^#rTAnoA3IhZ&-m<=1V?2K=r zb~bgKKEF)=C#%CA?mO*YpK0BHBzZSnM`4d;VuV-ZBv{vmlnJ$GoNYc}$9*>}3!hW| zsI81bor2IlJ#7})bQxhDaoGJ^y2Dy`xjHB9?vpPqs8#OD*vFX7KfrBKgSi~1Q?g)D z@JAnwH1VZ&7v_D#LkHoaJe5B=B!t!0zb)kxJp^^z>H3|ItKKI2L=laBje04&f{h^U z+Xp_`=!)R-(1lbKINN-soL+1ol{c`ubg{_EH%@?$;cfCs-WMuEN*d-~*~C3-4W{za zY~`P}CbrcBg@qdSIm_*T(4UX-+3uo-=(N}+TTnL48C}1E zE|f0`i~!JeVlj%9M@5Pl?)pxQG9~b${#pU&rhf#c(2i)Ywg|nxDIvCjhnpE(@WBSGLnO- z5x)bdR~jZFa65H_RugTR4;;xuvv>ks<*yj6jPm;0D9+QSP&85lU~jtR!*^}yg{$3>ggMoPmz9m=ntTA710IqIoK z?2+3oBmI{tvMQ{2_<~h1%v|?aZKhd2|I`R?jPgs}!*IG@r+uYer4kN;F8wgKdcOoM zv*1=9DveKz$Qws?^hSpFEr>ao#EI*nn|%>-94uF5*j~Q7!66hD)e@##5$5ZK<-$5N zd5@K%8_u|0BA*U+F+=MwUyc`~>J)2loL_xj>*itP(Kd;yOcM##(Z_szYIy;&6+IBQ zV4|;}>#C9Az!rwGK5wA&sFUccLi6qItN0I553JUJ1u&oP4CQ||n%vsLDIt!TA?o*s zKy-8rtliP8aZ!M0Pm}%1T}Z_f5hKg|QhQYu+dGwEk#3=M$V!_9@u<$KTM5yLh&f#m3v|2Gm@54HWuiPLR>iKhbZACgA_R;u4f}R(rI2d*4 zE#@_AL_;kJIabGiZ=HpAFq>=C6dp7jr<^0Gy6Q%Z7?;}JYb<=-v>N+QetD=RlGy6^ z%l}bYI@{1ZCm8Q!?zqXXsZK|XG_fJeov&Y&xlK!}b%|;O)!}wRcwS^)B+vBTbavm{ zxeab;Y0>8P@o~~)#Z2h~l(6f~f7pBQS7Uf+Xi5tWKw|xg%zBWKrIeD<>VsWHhD-Uh zG&Dmj-!x3DAhnasQPrY?3szG3it|iNK+W4uww=p_AffnRwn@V+}! z8vvla_2<6;he9XVp_}(p0=U+mn$}boAYI~+f~cfYA|O5F(B*>&hyv0fpi)Zr5E25?At5<{ z0z(ZU-Ej8&&pPkdv(`E1%e&63#Ukgvuj{(^{?#V*;REH%7wIk{5QxjFD#*tO1W_b> z@Q|K|pR8bTjqr)sML|`Y6ka~0<{|L?1xJ-9E(ipLDdB_Yi(Ixf0>Oe%Mc&o+_^>kJ zsjt0#(y-o3ghW~lIGZ$`d#EMi%o}}&=n?JrpPvgpy|6jxT&gnaHrVCPO{P`jvVM^{ zQJqkg6_7Rg7F(+sXT=vg7gL#>P|~%3{Jx3NX{1gdY|l@o*OSR|t}AiPBS36Q!t%Pt zFq#PQF~^WW27zeE4-G&dVkjvQ0awnPK|EBDMq3Sy>PWPGVxj$5#RW zn?FK@Gs$cG6zSGNN=kV#AfBBiK^xUiFWlAFx_yStzr`#forDM>`cn(JBk4ssG#D2U z64KV6d7qPorJ=705n$8V-5uMPt}Nws9haj+bb*vK^7rp&$`of2C^7~l$^TkOv#NTM zNPFKsdCsD=56p*cem1`ZBa<@l)pirrI_EE4XDKKs2xkh{=t{bCzre7HjVoI9X}JyO zOlw48-?nWdh1?Mf6BAQ^HQF{=E>2j|_n?f&s9KwxL7XK;$~`VS&9A!gQ2e=E?9PPy zlbT+x-RxBJ~uN{L~`7@qf@s(N543bt@$p` zgtq?lq-GX}Yd3PA(3z7|N`JIkIYZvw)fE#E5HQZiak#q{-`tFZ3&r3%F>Rb{Ye?J;lI^I3lT{AI@z~}uN7|6-VnK`?$*(=9i>%AK%YCooL zG{tn-W{j#KEOxZszh1BvQl52D_X>EgHdOlapH<$2h!eo=n^s_;@!EiJYBJ^Dh`2u`f+r}T`>%n8@}2-yRbsOV_Un>W8S7F#8`HHj6ZrKK5u^C)f^ zD$EOKlv3-AbS$zPD^rC-)?0}bAMne$Ww3Ii6xR9%-nmfC zH{X+TRZgy+9wnBj9+j5HZ)|MbpRLK<-rim;Nj6$)8G&`_l__ywn@sjY=)D$ROP{E5 zpAul_?jAa9OSSQp&kJf=ZhOz}yfStt(=q0x<*SOHU)|o`p7!23ezxhkIW7qai4yN! zl#`Rw`|0WD``g1N_v;o4jcaew@EdcYjSQpVR~yLEu;@w9Fd{Mb_bp*e4Tepunictw|D&J%;msTu6O=JyJ>OR_ndlUWW>R58)qGsh zUuQRpv3^CapvY)?{{ntRUdmoq)n}8xDoh;^HT%ytG za7j;3Z+&Yk25L8V*ZcQN67!X@a?bd8UIVW6^>qu`2Rg%gI@;PC1{HSiyOYK4uT9o+ z3J4Uo^H0|K)c}=8$rKGhSemr+Eyv zDNROI<+Nk^A~i@mu5=|GjUQ@Xov6WO#vA7h?yYRsnp2eTaPqvDGdX1tcl`LzwQKDa z_Tw2vMGpzI1%)k^Z%f`@ZaZu~USZFZIdiVmJ!x4sw>q%J07k3(M3<7VWQb zUyB-t{oh($UFEEQq+8~Pj>SX}%7Q>n)Vxn)Jm(@f)BA2t`5op1uarH;$Iq!BE&@~v zxvv^lIsW+L)L&>KM^!xKS0cS*%FV%n6cQHByc>AlB302mHRvB2nocN{ z481~QG3+Fi%P7XdBsM%qyi>NLvorduR(8iiUpnpwjR{W`pw7+)a7I|jaee%}&4KL5;7kC4_J*oGrZEkK3Ngh9y^j`SMy1(Aad|#o& zqBBki^Kyd`SO$y5&VmHECAL(@Py_3fvjppZ9N z2OKX*_YRjD;=I3M;pE~H4RR&7mH$DCBxn86Mp!@Xlk(Uf#h+5E3pVGmjQLw{Grpa5 zf})81bT4?))fccbLo53Mr*PD_Z_)S`1pdt2{Kc6PR5ozH&OIDWX0|Hgw41mTpq8`L0( z%B2nAB|s^i7|L$AwfWr$NslJ5d#$&75FSW4EWoh!(yV9$1z)Qfv$&kb07^WNto&ZS z3>g{XnZbXlzQG(b7F?|=DFVfbY3wb1il&ZBZv`#)yUfz`S6qO{8x>3Ch_eX0-yaL` zv8Py4shn#d9Y>f+`iTCIl;nTi`2Sjy{sWHXW3`yF@|*t@0W$B)Za;}!We>8>eS=mqMy);JYK{`ubj&P={4{3>s$t`#pUdWvum5#rq z6*j(e{ZefW8b8O)7L}ZgFDf2OvQUAt{vZGnOPS)r&G;;Qb$V;SxHhb|^6%fj?=zd6 zpvf`p%He|!6V>J3^70K&#Kgp+_bW8zjzz=-1qBCL&z^TYOa8gZln>PN9-viys;qCd z3UU^F$rz}xg?xOqgR5(qvb;RPBrV|T3@49xSB_3@-7<*IcKGb7{Z-Y|MtHJXc4g195C+XlPhJly9)PeNmhWcEg$^=CHc^orKx$A_awm^TurZ^g8$- zpaqRMoaCoZpQudy_ASB-@bQJpxOKUcz4X!2URkN1cbh7zYF4XrSBlMA!|~Oq91r`3 zVp39u=|sy^&&Bd#V{hyaT(91vf&TszVW;Wl_tz%)Sy)&E#Kd+3x=9vT_`m&almhU# zPDxCx5@Bbb(NZ{L*&;|eVzfEeeGKyZdG|Uc6|Jb9e*MVENXVO2gX9@jwrM=xK;T7J z0@LE+Vkr&p(=skr)<3Ko1ADSTTGM1#-NbALzi`8%#eAlXp-v6i+1XF0zX(dugZU8* zo`v%>reS21`E-kwaOJY^&jhTJS`zJwXwRIfDG@ys$X(`{b%UGSe)Q%lMYMN!zw)boL?;;=ml=&C?Mdt|x%zIx|XKtlmAq^H-x z)}li^+8cB^qsD+Vc;e^1P7(xZx&WN72(SG=fBxt-H8wV4?8k>(THex{_?+y|J5XH_ ze_Md(;la$p>c<7(1}%LvU)piqVpHnB&=X(q)$7Znr8fx6&oZg^y9xrqwZo^@lberE z6m;8m4X{T&(?9NOXV@mTozDVJ?SL2+C(QJr8o29(1YM%`!I-W!QtmqoD?6F;~=Ai^stDJjIuyTb43 z>8Z#TWnR}kJUnbWT4Hez^yau7h;WL~wvr?-b)#X%o&KK5TJH$l?5w&`osZ{hPllmz zf3Umzz5XXIZqRAAU^DlJ07*1sF-Fnt9Lb{g9i1H=PiviV9fLv``|&!cF< zQ_mC^#TX-zVd3H7YWMFSV*0G4KjEtdWn_Fqagsp9DUYp9{P@w5mzQ_?AX#*U z)8GG8X>L4QGwoy@Tje764({lS6&pDft3_)TQr&)r}}KPa}U861*iUW z{$NqAzIvsIcG^Se>vxa^qc~Z#veZ4lOCJjG>!bw5kOzB8J)n!^f|l0cWef-XgMos) zYPBK@GWyIB3I#5eV-fW_u@dw;g zo;h3Fp`oF^;0yGQX}-V%<7M%b$b-OgI>KlqP%wPSKgK= zkPC;jv`<-<^=W^nI6KrEMd0_M>33q23x9Apo!gUV`ydWE01M1a$H>Tt;^ol!^;wZ@ z7*5>}p6139wrcZzAPEIiZu|OYH&zzMi`4~qwvxvKWjdpyqhha+!)Sm5YeohJpLxit zykk>+6Ldt#`q0NYS=RSzYinD{=|pat?*0990?jSzqr0PZm6aeHx!KvxdRyMn%N;>< zqH8#HLq3s|^uwF-$mW(7*`J_6B>^8w0B*kJ!|B|_pMYqOzdhpjG9<$~X=gP-(98y{ z0?<(Exjy~lYN6b`5Ec0sHWo)UYSz)ehimm~_EUAfonWf30iMsezoi#Lx$EiaEg$T}1bYlE8}(=Fr2E4KMp8b9E12M0#Rd1KLXEGb z`Je3ffB`AzVI)_MAj43+U~zmp)ABA^+`KKaBD1h?iHCz5UIGG+UqY_mdb+C{8XVI0 z@+?Ui5PM#3*W0s2=VTvVe84Ixd4Sd!P@9bT1PjpDA8~YdchAKWsu6;+zWUEAHc4|= zsHsH{-YYX*fO@=;`<+z^OTz5)_or~LbaZ3h+3?0o#N=X*zZ+Iv=DI>1`ab_qrBb~{ z1u)t>5lq26c|4cwuzC~kS65fZ*xbRBVSqv11;ewl%jVxB`x^eBdwkpg)O%}4n-aau z<(k7k!OY>}UY9DnB@Xic^6(Va?u^w`7tqV?&qfy_Gf?OXPBtcC%O30r_}wi^s`wva zG!8!tPpiV=88JtMhm*T(ZCaRqp7od41h{+?-vn_!C*SU+#f>CBq_P>n8CU4&Y>K3*A74w#@_ZE>OaVt2pOBDn*QnanvvGX<@HP51IpbJi z1tS>|Lr_42>3MR75g9J77m+yhyMY(kk9|!_x~w!Hf`OlEIuLpXQb=-}ks{N;t2Z9# z8;RY#X%&-YP3CG;uN%7JWo~XB^UbokyF0(FyL;u|khb6711G>~iJx}z*JCX`({7}C zOux3*f1M7hV(1jhA3#q@InvD@A-$pR?CjhMf|pa~^k_3$SV+jY@HbzqmtG|-ysV+2 z!R!6|_g9e|KgL#e$b+-2w*h25c0n_`*d6N!_2I^T&Q`1CeXhkW#PiH_x(`dLbp%8@mp%+k+jF|t)(WG0d)1J)n)<7C$1M{6A z(X%ZrEPTycThH74L-f)}dZUd!;@x#YHhE+QxNt8I509uUYv#pjW5Wr1Nci^l-wdg0 z)UG>r=$ZI#8_CPd_W)5RiveOAJ@FU%K@HI?9%uy#NJjnv%9&c>UxRCwrr+m8G+gLK ztbVF!_MWCVS_oXP=Hld(S^grpSnYQQX!ilw=!Fc`=##|v@7K-l-mOZj(Vl^eO#l$n z3Bt<#yCDcrg=fcHLyq{u^@62H0?OuzMz4Z4TtWvqXO*JMQ%U+_o5G7+F z{~4_VeJ-O2%54$OJ*FILw6sfs=gDhA8bKudqN~z61yElTG+*&B*n=R8&;St(BFRqoyJ)Ba@tPYR}664-$1*=4`{xo{GW8H!pfVNNhs#!El4^ zv7sRkxP7@QrEs^gET>YFnZu6{g3ZU4|+dhJ8?bpQgHyFkKc=Xz2Ob%1diA;3bR8-5a5bjb(MJr|)1 zn*7=J6Er;=YzP@=YBT= zlTG)^o%$TYl!v+y>{1|VUf4iHBUcU<(SBR6fWZyd5o^n~Q6a<3j3Fo!H8nN56&&u@ zW2ABMwi!@%JJhohq_RqNTmARxLBmMtadLptMfdmjSNzIyL@Fy|Tp=p*KRujMUY0nI z9U~>RD*)P;fq_8{W{W~aW55ss-VZP-na{R{X{w%zii(0Zq_1i!Am9JO zFXC;r#1Yqga37UG0F{Am@;rnv2=&kYCT9+zRm-8)_$c4w;_%RS?|{bF<^Vx7W%p+} zs62o5nCV;$g^TDgyMQ|j*3HvjuM2E^q|g23ShYEVT3gc z&i-tL2Sea@fjr9<1(H3hShI${;()?G!X~#nz9Tf1fns7^yoz`h_an!#c@jXaP zObooEVNef>{L$p(>9hK0s5 z-j|8WFA?2+2l^n9;mVb1PPQwwwEGFc!NG;!pH&(lGhk%(6cCB9{m_sw-10pnkV1Ds zDW&ZnXt0WWW?ec_yDxGF3+YDFGc%SF0vtgKDi6<+lG?oT+%c__wi_w(-rU$2^YZqN z&ofI60<`sb+9TVwEJ)_0Ew1J=6bUoKLtn~^C{JSIJB`jeV36nMhd;vP;GG87fG177 z+*IM%mq+ZpA|fKAO-C-Ciuvn2 z_h684`Um{)j*kbnk(4wie`^^Mi9U3#ahj|z76vAZlZ2SJ_YJ|pPl0|gfffn+B{yZP zt(~E$q?8rXW^hNQ2Q-QtFWdB`KSfZC91yb0o`5*n?F9Fu5l+{-9qF_5rL{FbIN;-$ zTCXhN>rz!s&2ZIbkCX+J__hX@)=N)nA)md!w!mj3a7Mu~Xs*W!3jd}kt?Gjr-T-|2 z(ALvqoER4;syDoSZ)h zuCA^$XZ|%f$8SpH1{1Q?bhcA6WiPJ-cCZLC^mlBUTW|f2L&_thSQ)<9laf9oNy@4M znmWGC%bV{H3olO^nGB&~zcMi~v8G?`x`-Xx9r8jT&c+b_7a-|6`$L_h^TZ=B3-DRM zEu6i8{oRn*pXx&jx1N#*lXgx>N{c$L_o`GASl3k9grQN=p=;+fO8sLn0OcI=m-Es!(U`34)l|F`=IoLU4MNr?;b{ zVG8O?#A<*$c{Y|*dp9>>-#QP0Nw;C zcNFs3GnPLA-4JN1E>s&oTq;pT%Jzw@OhKorkzT%>@2Vo|FnNUaUaPgdbouhTu(xmR z=9ZRhN}XnQ6y*^&vadwp)AIJB0HT(_G^Mg}6rDSB2FsQgGz(rFm)fYKtFkQoe9-NP zm*+P%fhy|j{d2*2o*IrG!qn3vf>wBr8)kQ=K=;=(@CTj-;Z)`HL`0~2Nm+OPd4GH7 yh}39M_CL9@`{y(N5B_c9zfDB_|8bNbbxIOM!WRE6*ZDRAUaE=@kOd0Q{{0_JR&-nd literal 8832 zcmdsdXH-*7yLLbxK~OT?ATP-Wy)DOBB8elU9YJ17>PXhJ#C0~%$3jl!m*`JH1QWfq70Pq38 zTIyyYgFUN^Rs6A4!%fn2{uSEW6EzSXw@ZJc(l}+W142X z@$EAbGmyqK_c{~z2}b**ySeo7{%;t?;$tZ3(k%Yxhs^9;i_6C@52a%;2lu5927EgA z*!J^>PgWnKE-?V2;`|>&XaG^ttTZ%$C;WT>fH)l;;Aw0W08q=q0*JavOAGj*p#}h0 zFfahVGX9?iDP`EJ-DP)L6beiwC;;AFdbtfW^n#x@Kk*)B7jV*l;ytia1ur%U3TnR2 z${L6-&vhql&+*Tx+_*LTwrp{6J>;TOYuadcER(F`cZTwc3NZcyz&FC1VVs8AyDccEC!Y9f!Wwy z=xOq9pv1Dtgpqwa9)?54V+3vgY`Eu&+s;x zhCT(JJqV5uvTzYa<~~`9;6jaV9^AsKMp94^++l`*P5{=cMN~vYM^sdF1A{0D`aQK* zx7F-DS3O#$Q)5u@$bGalRveVcYT<cSPwz9Gs*-T1GYN}s8H_JN8P&p0(r%*ye=Npb+j=4^5CR8jD4oc2# z!;c9W*eweemx}t1L+g(7Gm31zG8anJb|^oJF3#y_`-XXMOZE*Vu5C!nQatHIUD_U3 zt*kRgwC~9hHEj$%aQTCFJTvcP$r(YU7m2*skwX0jpXm($%*e!alb094>GNf8oR}$M zG-($+S1zEFx?Y3mL^omV#leunY1gKpEo-ltMsWO9ZnQY!ODPpRk#%d8qh(g&@9Bea z4E3I|wsst~BQxF8fmbqDZ!T!ZVxHA#{Py)0DO}7y-xJ0((Ebx*TU@I44Hf!s0iaSm zyR4ZhSJ}qA?E#y!?I+9oLzEXIo99wHenJ7*n(vo&d=M0m5wHG%0pda@8v&X4i&B~P zT*Z)TH7031%I;rO{y>mFv8ncI+N9MWG91rcH)UJfPrfXwvN9ehj(J(prYHB@#;<-G0t7Grb@uwB}O*5BX1}axH%?P?d`p{pj93~TX=;<}rP8`4Q`0Kv$yWhV( zon2k^pkM3}+ihz+SGa(?<3xLeCFdj0k7Z0#9YHM+#o7{*Vi*)t>Wus>8Cc z_+Tsgb$BO!PdD5Qo_6-<(z(5|(w`@1q{t>Y4j0hAdtOms#P`hLDuE(+Zbjm)Y171%u;_WMDUE{r4ejlU5d*2wx^V$$QJg}oA;L}n!)8v6r z!4cjn1g}rN3)DLBJEj^|ZMj~qJ8@$|)_ZImC0dxjMYmW2^bT7$m%tVI-{GO5I$K5= zDv%ITt}n4tsp98}_>l$Axl}G38ymZS^s>LdzxoB1>K<<0F_+wNogLIilJT5+ke!`f zbo^J?O132d4=<^Hu{N|Z_qoOX<=$g3S{7qlCTd!-zFOd8u{l{6I2VX(0O`UAv@-|DJw8Y>BVH@6YoH z<=(8WUL>emlnbE(=F%ePDt*aROQ~sWG{d&|AYNcU?{wM-2?<$uK$@=ds?<^=#(*D} ze5wF=W1=S49)cIz!gV?PSyfm+?;RDpzOA6m?aw=t`LZ635vLHU4Vp76?W&Dn*XKJZ zZY1d8@9DC+&%rhC-th@y<>Cg|owsMoOS8qzZ!Gg}#f_VVt2}#kv4fa{RNYUQ-S$yU z&NH?M5)(XsvDH*vSy_1yRQBIib52-Tlb4Sr566<>%q?pKSJfjpiz>CZMf~XWVhPh5 zlcvK&wS6igX2!b1f#1a!07G8(c zf8JJ5K0i4)UN1`v_&p^x7fzbvPo+ni^bd9)DdLie6~l zUx?2ht+%VA%;Z5!HJ6b14SJq*lqspKX*4LA3uO!j_f@H)%gf8fL`AFoP7OknY2nOW zUqxZZgatNWES({5#3NHsBRAZt!3CeHxIHr4_QFG%w7^y$skboVdEpM2y5CGxS`#za zITk`(uh2PL5LI_P=^Oo)w5iLw zt&hCtkaih#8}=xX^~f1gCfN#!h-@BWeK5H5v*14?hKrT+A3S@7Mx(QGbK&FTW^<9} zr|Ou8R1#oNqGGp%OfiK(+W%zgx6-b~?GdYkZ{c&9)gBZi<7E`g(tb33KvIz?jg4_S zCGTyJRtwrUg1byDyuCXx*9=h4_ZE8-)6y`#Y!0>5z<*a#((uxyOYnMTj8wHsIQa;k zt{H#zuOgL`6H^co-!tsO_@-hhV5fsbwJ%Ul+!a$36DYMA76^Br`iK@LM`kHz-gbiM zS-1JG^jqWfKbDqiNIO6@eZpI<929T2`-LP%JUE2D!{+4VY~a`I$5Tv^B|$rZD;q`x zl3dREm>pESvuTEMxKZVR!e>QIWkP;W`E&q&|Ad!VAV$kP zQjj>vq1kXpSd)~{EA#;;4GXOs!uUSjw&?~zTh{_5qhb4MZP)|DX4%gDx* z9i_Udb0Zl&>RO`X&-3xLq8{-P3gKU%lj=UFB zc8N{vhU?mTyDXn2WfuB2vUNPa=_PlbeaY?Kckgr+Km_A2VNXSn8r(Pj7Kwv zd3fg&!2tlZF@C=PCZ6S1MHj5PBEn(0006(A@QJdA#c+4xzj&_JR>OyUJ+=SEd(9&9 ztsAJUnePPn=IZ|oiT=N0+W(=0wD^rYO&~A~V7!+V{}+I)jh0_qfj91idKA=`?F(H= zAwZ&LSXlTN$vvYHyID+&HqA!B51~hKAjG9gCK{SMl@Q>a4xV@X3=AwZF&BQ)h>S)* zkVby*t!VB*JdTgfWnnJtV_+zSN(l=~wBWdrd5s{1&^Q-LzH{dw*79cO`Jil=N|9n$ z<8ucq7=~`@7ymu2n@;-r125v|(-luv#h1{Zh7=oD%5sIcvtY2~p*g3;v3i9MZupLC z_AVf<`NO?1=)9gA#_754=WSK}pWD*ox#KgO_&^%sV>OnWv^^O)*d3)@RLxR+kgmeD z)T|(7Z!gu5{YvO-MJ9nS!jdpRN3)ieZ8Dpd2cl8{c2NC>NS=n{Hn zTTU`K`@H=^K~5b`=KQ1aHLC|vW4_n9ubEm|e0g8m7-IPtY~1?>kye1{dU?M$BQs~F z)6&z{+9pCZ`YuxQNkZOh8w<%Fzi(&$P5vel2_CxFsv-d&PH%0jU2@hry#kTKbM2!AY+MmLet|lB7wKLc@l4PktdV+;a`j>}`{R9CF!s6#{pJK{kbXLU`|eK#YwBQIr@ucY|W}(eo6jZn3=!JH{!>t+nRFsA&YPc zh09;=NqioD>f@*v;Shbd8yuje^=HMTbMerKAvGV!XmNjkaD-jJX~pC61KJkw(C?sU zqeFpu)K{JVujwqlT`|Aud>flP%)#^CS+4D#lYbD%h8tB?9p=IJsZ}RFVg1i$T$vNV zxiDBji=u0Ka5Zd<*8d-rE4=g2NW{P+4mF2kNrFu|qpg8GzD zb5P3X7=L=XT{%AA?h8ShSR7-dEiV{Mu?P7T+qEyN(U{ zT{HgyUnrl?)s^q&3t*mHwHBU%g$-C|RO1=J0b21)3qwdk5Vz6+64 zjDN0=ZSNE@f0(NmzAYs9{7#{k&J*rAo3@=(owU=B6`M)S7QIv*&y5njy4GH5m9b2Q zCF9=yWu-1+OimE@SuA@tf~wTxI0%$sd%P`+&|#CFoOo{CJXB7+3F3B7s))nkasQpkd-hc%7sn_khmYiTK`eo0Ehd`^P}Kk@kbGv9l;#Rc~} z?pAeVY$H00GR95ajYaKD?q-9IhuaUgK*t(0Eq2x&=ErvlHK^&BC_eYwPWrl(Vzwq; znNau8z6?e%>j|b;_%jrB^62Q(V6oI~o7`MG2{)HYZ)qs=)=Pdbd+|AW2gqnw>QlG` z7g)UfRIePv2JXTw0>j*~2m!pN+s;G$L@YY_dSA=1-=SRZa8v_5jr#bd%R27fxafJZ z8L5-VWEht5?9BON@bc9|7xrA^L%@!p(_k%9pC6ewWl^LvJ6nB_z$BF7Mqqx1zk+&e zj3uUAK{aQNwtma&I-g}Hnc(m5$Pay)Iiwl?^AvPUIv=(iGwu4EDnuMyY(=);26Jp} z8n+I5Gw#I;_$+(n!dD4+-;mAx%9HP!I_SiUaZaJ>zT|Q}=Oq}E4(N4+_uzF6T24B; zYnC`)av&!sWoedw{_;7dQ_%7?Jbh!qx6O%z;v`cJT`P~Qg+{QlrzRyIol4@3Jh<0v z5JL|4f7;wpwp8YtXv&zi1Z!TpG^76x71iB>X^z0{}}EQpdfhY>8~)BREOyJmE-kL>Ce38qRHC-}t7T@&>5brj!X zD(_6&H2cBa5`9aW!O&m{GzZh=Lrf2ONw4TiBFzUHMz0M}{RN?7njj_oh#JznW{4#t z16wH`w-C`lBVC|TyXOZPg(o16>XhY95XTn>-61KLrCX?Ju6|-^C_i5UqezLL43=Gb z?MJGS5!0I`pzd=kw2JGymY)mU*=v?1DkmHH8dPa4g|uBgM4msW^1_IZ`5y-bgqB9q zR;B$Nvz%@w>tNkxdyJ*UJA4^o38YRcrmhLaI%)aa^pdaJ&74}{=e>AuDqD?_EOYSB zxheNE3G_~bC5c0v^R=1}C-Y(vx-CT~Z#}etPm`kf#Quu;#~ZR%Yd%kg)n?rbb|(4x z)#1*}M-wIVba6h*XRqS~WopuD5wN9!^2cE?5sxJL1vRH@ELu(5mfznfZTH$H76e1C|tGRHRa>J{XG5!bHDT;z32 zs&r0%b2=~3$E1*w`t#D4likDqPc4z=GF6>afRPg$wxl*nXDzvu2dhoUUca*UL5epu z{S(Lat9)5`YAj9d=X{6FS9`+cC+_>K-#~{Z`)CX(4Ld{YD?PB%LOjq$tA|*a90UEG z*x?K6WrJT_d#z5HF+JGD7%Fs4Mzy9YexW$)SrYZ@g?31%IZ)l6i`G=vyr6pX{|PuE>Tv zcz%+%RIf;%EK-^~47b{Q8{l{hB+)tqJCy3lzy5;JgTUWJB!9qe067IZE*8!)Fx-*! zsSLI*7?SdvJ{3niKq}#ZhWH`bl>wm%_OyI|GfJj)v`JAmBOd>*JM~HM zDUtr6_)|5UjL=g+ZBgYUfPsdZ_xm|FiY7n37O{uU+L!%4n-`}Yj+N2+6`J|w*?dNw zLc7nN5pTqqs34sl{Sbue!K}b=vaAau{?Dx3+GKwr(Yhk|<0M<7;SKDASipNtH5RGp zvG|g$8`pv1z*GU0v29H%&(7J;n>dAG%MwvR+0Ud{YDz6cr^^CdISoq%YG)1;+MZ%& z5@W5_P&aISeb7IiV(Ikxn<5k(-ckLnDoe-2+$8R^Ke2Z46ssaenUrpN2^MxxFV4u$ zpQCgy55C%Lpf|eR7&(TvBA$JVYG9#JpJI1XD%&DT(Jh|dr9xmfQT&ROKDu<%r}1>04ulFZPwYt$A?6{!2*!2;SBqN<)#u%y-j{>ITYH`* zgVzO!PZofB`Nz#s71TU>3&}p*+hxeX%jz01zdrS z8en+mY*rp2zbmGi{_eUbT6S<$H66a)&!gn5mM4;WD-#xZ@;&b_Js3PkehB?w{t=4_ z9o@50V%SGaAO|=0G1PG}B~y7s@f-U=;S%s0Z;yMFk%N4MQ4yy(C84V9a~E0WQPUQ* zEKIK1AVse)u&sw=^uFkZv|0PEB6_!ky&)m}mws=-Tv=$uBk64HP!uzAL`#1Ef?LjBQv>88YLfX0-P?! z1=6SqNJYQv`->{Z-xq6m%BPy%!AGV@88vPwE#@o_o`{=`P_XUw_{E()Pd z2jq2LNdsuqGZ~dg#&dly_%Wbkk0~|_gAeJz?3HKXk1X1k4RZIbb;AXzfk>J5mGDi^ z!P|BPL%T+o$ zgX-$;(L~pzbK5cE3X0pv-^MZyw)M=`qR8;<8uRf?#y!z3ZG1%pdg>B#&MMpdzJe32>flabu%rI_VT_R&&HgfQ^lTh>%9Vd&Hj=A9ausHTFnOAz~ zzeNuW-DX~+UtL0;=zNj_%X;Aq6*^7Dc7)4=(ua8O5Sf?!qtoy3sWH&SEeE;zMN6^$ zg_Q)0@d|MdWH=)A1a_cjhtf{PyzYY^cI(-nD-A^h_&*s{KjD+G5k6gEMTOR=;4=Ll z3iiO>t4gnSIHOB&FA94NIxP5ZFfa^pAla#!m_M@V=~Li*ohz_`Ec)>YSLcBktz0o>tXD5jm0mj3-uCSW=`qJ|oaaS@6v ze01H;3K}ss-ZG(%Hb*mG)McCV9GoS=bi6eD)e=M1AzkvZ0_AqY_m=!Ozpk!u_Ox9r zoMvFiw}oKxRKhTpwbf_zX!fPf^>46CA4ThLA@gphE^#rTAnoA3IhZ&-m<=1V?2K=r zb~bgKKEF)=C#%CA?mO*YpK0BHBzZSnM`4d;VuV-ZBv{vmlnJ$GoNYc}$9*>}3!hW| zsI81bor2IlJ#7})bQxhDaoGJ^y2Dy`xjHB9?vpPqs8#OD*vFX7KfrBKgSi~1Q?g)D z@JAnwH1VZ&7v_D#LkHoaJe5B=B!t!0zb)kxJp^^z>H3|ItKKI2L=laBje04&f{h^U z+Xp_`=!)R-(1lbKINN-soL+1ol{c`ubg{_EH%@?$;cfCs-WMuEN*d-~*~C3-4W{za zY~`P}CbrcBg@qdSIm_*T(4UX-+3uo-=(N}+TTnL48C}1E zE|f0`i~!JeVlj%9M@5Pl?)pxQG9~b${#pU&rhf#c(2i)Ywg|nxDIvCjhnpE(@WBSGLnO- z5x)bdR~jZFa65H_RugTR4;;xuvv>ks<*yj6jPm;0D9+QSP&85lU~jtR!*^}yg{$3>ggMoPmz9m=ntTA710IqIoK z?2+3oBmI{tvMQ{2_<~h1%v|?aZKhd2|I`R?jPgs}!*IG@r+uYer4kN;F8wgKdcOoM zv*1=9DveKz$Qws?^hSpFEr>ao#EI*nn|%>-94uF5*j~Q7!66hD)e@##5$5ZK<-$5N zd5@K%8_u|0BA*U+F+=MwUyc`~>J)2loL_xj>*itP(Kd;yOcM##(Z_szYIy;&6+IBQ zV4|;}>#C9Az!rwGK5wA&sFUccLi6qItN0I553JUJ1u&oP4CQ||n%vsLDIt!TA?o*s zKy-8rtliP8aZ!M0Pm}%1T}Z_f5hKg|QhQYu+dGwEk#3=M$V!_9@u<$KTM5yLh&f#m3v|2Gm@54HWuiPLR>iKhbZACgA_R;u4f}R(rI2d*4 zE#@_AL_;kJIabGiZ=HpAFq>=C6dp7jr<^0Gy6Q%Z7?;}JYb<=-v>N+QetD=RlGy6^ z%l}bYI@{1ZCm8Q!?zqXXsZK|XG_fJeov&Y&xlK!}b%|;O)!}wRcwS^)B+vBTbavm{ zxeab;Y0>8P@o~~)#Z2h~l(6f~f7pBQS7Uf+Xi5tWKw|xg%zBWKrIeD<>VsWHhD-Uh zG&Dmj-!x3DAhnasQPrY?3szG3it|iNK+W4uww=p_AffnRwn@V+}! z8vvla_2<6;he9XVp_}(p0=U+mn$} Date: Mon, 11 Jan 2016 16:44:56 +0000 Subject: [PATCH 4/8] [plotting] skipping plotting tests, as they are inconsistent across platforms --- GPy/testing/plotting_tests.py | 54 ++++++++++++++++++++--------------- 1 file changed, 31 insertions(+), 23 deletions(-) diff --git a/GPy/testing/plotting_tests.py b/GPy/testing/plotting_tests.py index f833faf0..441854d4 100644 --- a/GPy/testing/plotting_tests.py +++ b/GPy/testing/plotting_tests.py @@ -27,13 +27,21 @@ # OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE # OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. #=============================================================================== + + +#=============================================================================== +# SKIPPING PLOTTING BECAUSE IT BEHAVES DIFFERENTLY ON DIFFERENT +# SYSTEMS, AND WILL MISBEHAVE +from nose import SkipTest +raise SkipTest("Skipping Matplotlib testing") +#=============================================================================== + import matplotlib from unittest.case import TestCase matplotlib.use('agg') import numpy as np import GPy, os -from nose import SkipTest from GPy.util.config import config from GPy.plotting import change_plotting_library, plotting_library @@ -41,7 +49,7 @@ from GPy.plotting import change_plotting_library, plotting_library class ConfigTest(TestCase): def tearDown(self): change_plotting_library('matplotlib') - + def test_change_plotting(self): self.assertRaises(ValueError, change_plotting_library, 'not+in9names') change_plotting_library('none') @@ -115,12 +123,12 @@ def test_figure(): import warnings with warnings.catch_warnings(): warnings.simplefilter("ignore") - + ax, _ = pl().new_canvas(num=1) def test_func(x): return x[:, 0].reshape(3,3) pl().imshow_interact(ax, test_func, extent=(-1,1,-1,1), resolution=3) - + ax, _ = pl().new_canvas() def test_func_2(x): y = x[:, 0].reshape(3,3) @@ -129,21 +137,21 @@ def test_figure(): pl().annotation_heatmap_interact(ax, test_func_2, extent=(-1,1,-1,1), resolution=3) pl().annotation_heatmap_interact(ax, test_func_2, extent=(-1,1,-1,1), resolution=3, imshow_kwargs=dict(interpolation='nearest')) - + ax, _ = pl().new_canvas(figsize=(4,3)) x = np.linspace(0,1,100) y = [0,1,2] array = np.array([.4,.5]) cmap = matplotlib.colors.LinearSegmentedColormap.from_list('WhToColor', ('r', 'b'), N=array.size) - pl().fill_gradient(ax, x, y, facecolors=['r', 'g'], array=array, cmap=cmap) - + pl().fill_gradient(ax, x, y, facecolors=['r', 'g'], array=array, cmap=cmap) + ax, _ = pl().new_canvas(num=4, figsize=(4,3), projection='3d', xlabel='x', ylabel='y', zlabel='z', title='awsome title', xlim=(-1,1), ylim=(-1,1), zlim=(-3,3)) z = 2-np.abs(np.linspace(-2,2,(100)))+1 x, y = z*np.sin(np.linspace(-2*np.pi,2*np.pi,(100))), z*np.cos(np.linspace(-np.pi,np.pi,(100))) - + pl().plot(ax, x, y, z, linewidth=2) - + for do_test in _image_comparison( baseline_images=['coverage_{}'.format(sub) for sub in ["imshow_interact",'annotation_interact','gradient','3d_plot',]], extensions=extensions): @@ -194,9 +202,9 @@ def test_plot(): m.plot_errorbars_trainset() m.plot_samples() m.plot_data_error() - for do_test in _image_comparison(baseline_images=['gp_{}'.format(sub) for sub in ["data", "mean", 'conf', - 'density', - 'out_error', + for do_test in _image_comparison(baseline_images=['gp_{}'.format(sub) for sub in ["data", "mean", 'conf', + 'density', + 'out_error', 'samples', 'in_error']], extensions=extensions): yield (do_test, ) @@ -216,9 +224,9 @@ def test_twod(): m.plot_inducing() #m.plot_errorbars_trainset() m.plot_data_error() - for do_test in _image_comparison(baseline_images=['gp_2d_{}'.format(sub) for sub in ["data", "mean", - 'inducing', - #'out_error', + for do_test in _image_comparison(baseline_images=['gp_2d_{}'.format(sub) for sub in ["data", "mean", + 'inducing', + #'out_error', 'in_error', ]], extensions=extensions): yield (do_test, ) @@ -242,7 +250,7 @@ def test_threed(): m.plot_mean(projection='3d') m.plot_inducing(projection='3d') #m.plot_errorbars_trainset(projection='3d') - for do_test in _image_comparison(baseline_images=['gp_3d_{}'.format(sub) for sub in ["data", "mean", 'inducing', + for do_test in _image_comparison(baseline_images=['gp_3d_{}'.format(sub) for sub in ["data", "mean", 'inducing', #'error', #"samples", "samples_lik" ]], extensions=extensions): @@ -316,7 +324,7 @@ def test_gplvm(): matplotlib.rcParams[u'figure.figsize'] = (4,3) matplotlib.rcParams[u'text.usetex'] = False Q = 3 - # Define dataset + # Define dataset N = 10 k1 = GPy.kern.RBF(5, variance=1, lengthscale=1./np.random.dirichlet(np.r_[10,10,10,0.1,0.1]), ARD=True) k2 = GPy.kern.RBF(5, variance=1, lengthscale=1./np.random.dirichlet(np.r_[10,0.1,10,0.1,10]), ARD=True) @@ -325,10 +333,10 @@ def test_gplvm(): A = np.random.multivariate_normal(np.zeros(N), k1.K(X), Q).T B = np.random.multivariate_normal(np.zeros(N), k2.K(X), Q).T C = np.random.multivariate_normal(np.zeros(N), k3.K(X), Q).T - + Y = np.vstack((A,B,C)) labels = np.hstack((np.zeros(A.shape[0]), np.ones(B.shape[0]), np.ones(C.shape[0])*2)) - + k = RBF(Q, ARD=True, lengthscale=2) # + kern.white(Q, _np.exp(-2)) # + kern.bias(Q) m = GPLVM(Y, Q, init="PCA", kernel=k) m.kern.lengthscale[:] = [1./.3, 1./.1, 1./.7] @@ -341,7 +349,7 @@ def test_gplvm(): np.random.seed(111) m.plot_magnification(labels=labels) m.plot_steepest_gradient_map(resolution=10, data_labels=labels) - for do_test in _image_comparison(baseline_images=['gplvm_{}'.format(sub) for sub in ["latent", "latent_3d", "magnification", 'gradient']], + for do_test in _image_comparison(baseline_images=['gplvm_{}'.format(sub) for sub in ["latent", "latent_3d", "magnification", 'gradient']], extensions=extensions, tol=12): yield (do_test, ) @@ -355,7 +363,7 @@ def test_bayesian_gplvm(): matplotlib.rcParams[u'figure.figsize'] = (4,3) matplotlib.rcParams[u'text.usetex'] = False Q = 3 - # Define dataset + # Define dataset N = 10 k1 = GPy.kern.RBF(5, variance=1, lengthscale=1./np.random.dirichlet(np.r_[10,10,10,0.1,0.1]), ARD=True) k2 = GPy.kern.RBF(5, variance=1, lengthscale=1./np.random.dirichlet(np.r_[10,0.1,10,0.1,10]), ARD=True) @@ -364,10 +372,10 @@ def test_bayesian_gplvm(): A = np.random.multivariate_normal(np.zeros(N), k1.K(X), Q).T B = np.random.multivariate_normal(np.zeros(N), k2.K(X), Q).T C = np.random.multivariate_normal(np.zeros(N), k3.K(X), Q).T - + Y = np.vstack((A,B,C)) labels = np.hstack((np.zeros(A.shape[0]), np.ones(B.shape[0]), np.ones(C.shape[0])*2)) - + k = RBF(Q, ARD=True, lengthscale=2) # + kern.white(Q, _np.exp(-2)) # + kern.bias(Q) m = BayesianGPLVM(Y, Q, init="PCA", kernel=k) m.kern.lengthscale[:] = [1./.3, 1./.1, 1./.7] From e75f78d7023cbc5ce772cfb7c29bd23e6fd4e4ad Mon Sep 17 00:00:00 2001 From: Max Zwiessele Date: Mon, 11 Jan 2016 16:52:40 +0000 Subject: [PATCH 5/8] =?UTF-8?q?Bump=20version:=200.9.4=20=E2=86=92=200.9.5?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- GPy/__version__.py | 2 +- setup.cfg | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/GPy/__version__.py b/GPy/__version__.py index e94731c0..f8c6ac7f 100644 --- a/GPy/__version__.py +++ b/GPy/__version__.py @@ -1 +1 @@ -__version__ = "0.9.4" +__version__ = "0.9.5" diff --git a/setup.cfg b/setup.cfg index 0b515a11..c8c61b7b 100644 --- a/setup.cfg +++ b/setup.cfg @@ -1,5 +1,5 @@ [bumpversion] -current_version = 0.9.4 +current_version = 0.9.5 tag = True commit = True From e917da8cceb871653c8fded1ecdcac963a74d79c Mon Sep 17 00:00:00 2001 From: Max Zwiessele Date: Tue, 12 Jan 2016 12:33:13 +0000 Subject: [PATCH 6/8] [plotting] catching 3d error for plotting latent space in other then 2 dimensions --- GPy/plotting/gpy_plot/latent_plots.py | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/GPy/plotting/gpy_plot/latent_plots.py b/GPy/plotting/gpy_plot/latent_plots.py index 976641b2..85a98f49 100644 --- a/GPy/plotting/gpy_plot/latent_plots.py +++ b/GPy/plotting/gpy_plot/latent_plots.py @@ -229,7 +229,7 @@ def plot_latent(self, labels=None, which_indices=None, plot_limits=None, updates=False, kern=None, marker='<>^vsd', - num_samples=1000, + num_samples=1000, projection='2d', scatter_kwargs=None, **imshow_kwargs): """ Plot the latent space of the GP on the inputs. This is the @@ -249,6 +249,8 @@ def plot_latent(self, labels=None, which_indices=None, :param imshow_kwargs: the kwargs for the imshow (magnification factor) :param scatter_kwargs: the kwargs for the scatter plots """ + if projection != '2d': + raise ValueError('Cannot plot latent in other then 2 dimensions, consider plot_scatter') input_1, input_2 = which_indices = self.get_most_significant_input_dimensions(which_indices)[:2] X = get_x_y_var(self)[0] _, _, Xgrid, _, _, xmin, xmax, resolution = helper_for_plot_data(self, X, plot_limits, which_indices, None, resolution) From f68001a84e281d437b5dc05550528ae8e4056eb0 Mon Sep 17 00:00:00 2001 From: Max Zwiessele Date: Tue, 12 Jan 2016 12:33:42 +0000 Subject: [PATCH 7/8] [vardtc] these two lines are overridden by the next two lines... --- GPy/inference/latent_function_inference/var_dtc.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/GPy/inference/latent_function_inference/var_dtc.py b/GPy/inference/latent_function_inference/var_dtc.py index e05dbaf9..ec055120 100644 --- a/GPy/inference/latent_function_inference/var_dtc.py +++ b/GPy/inference/latent_function_inference/var_dtc.py @@ -190,8 +190,8 @@ class VarDTC(LatentFunctionInference): tmp, _ = dtrtrs(Lm, psi1V, lower=1, trans=0) tmp, _ = dpotrs(LB, tmp, lower=1) woodbury_vector, _ = dtrtrs(Lm, tmp, lower=1, trans=1) - Bi, _ = dpotri(LB, lower=1) - symmetrify(Bi) + #Bi, _ = dpotri(LB, lower=1) + #symmetrify(Bi) Bi = -dpotri(LB, lower=1)[0] diag.add(Bi, 1) From 4461e35c9e3bf8cae609d4d96ba9552d6b9dd904 Mon Sep 17 00:00:00 2001 From: Max Zwiessele Date: Tue, 12 Jan 2016 13:42:41 +0000 Subject: [PATCH 8/8] [plotly] scatter plotting was defaulting to color='white' --- GPy/plotting/gpy_plot/latent_plots.py | 2 -- GPy/plotting/plotly_dep/defaults.py | 2 +- GPy/plotting/plotly_dep/plot_definitions.py | 29 +++++++++++---------- 3 files changed, 16 insertions(+), 17 deletions(-) diff --git a/GPy/plotting/gpy_plot/latent_plots.py b/GPy/plotting/gpy_plot/latent_plots.py index 85a98f49..ed12ad9a 100644 --- a/GPy/plotting/gpy_plot/latent_plots.py +++ b/GPy/plotting/gpy_plot/latent_plots.py @@ -57,8 +57,6 @@ def _new_canvas(self, projection, kwargs, which_indices): zlabel = None else: zlabel = 'latent dimension %i' % input_3 - if 'color' not in kwargs: - kwargs['color'] = 'white' canvas, kwargs = pl().new_canvas(projection=projection, xlabel='latent dimension %i' % input_1, ylabel='latent dimension %i' % input_2, zlabel=zlabel, **kwargs) diff --git a/GPy/plotting/plotly_dep/defaults.py b/GPy/plotting/plotly_dep/defaults.py index faf343b0..24170b95 100644 --- a/GPy/plotting/plotly_dep/defaults.py +++ b/GPy/plotting/plotly_dep/defaults.py @@ -72,5 +72,5 @@ ard = dict(linewidth=1.2, barmode='stack') latent = dict(colorscale='Greys', reversescale=True, zsmooth='best') gradient = dict(colorscale='RdBu', opacity=.7) magnification = dict(colorscale='Greys', zsmooth='best', reversescale=True) -latent_scatter = dict(marker_kwargs=dict(size='15', opacity=.7)) +latent_scatter = dict(marker_kwargs=dict(size='5', opacity=.7)) # annotation = dict(fontdict=dict(family='sans-serif', weight='light', fontsize=9), zorder=.3, alpha=.7) \ No newline at end of file diff --git a/GPy/plotting/plotly_dep/plot_definitions.py b/GPy/plotting/plotly_dep/plot_definitions.py index 613cdf4c..54f04a75 100644 --- a/GPy/plotting/plotly_dep/plot_definitions.py +++ b/GPy/plotting/plotly_dep/plot_definitions.py @@ -130,14 +130,15 @@ class PlotlyPlots(AbstractPlottingLibrary): except: #not matplotlib marker pass + marker_kwargs = marker_kwargs or {} marker_kwargs.setdefault('symbol', marker) if Z is not None: - return Scatter3d(x=X, y=Y, z=Z, mode='markers', - showlegend=label is not None, - marker=Marker(color=color, colorscale=cmap, **marker_kwargs or {}), + return Scatter3d(x=X, y=Y, z=Z, mode='markers', + showlegend=label is not None, + marker=Marker(color=color, colorscale=cmap, **marker_kwargs), name=label, **kwargs) - return Scatter(x=X, y=Y, mode='markers', showlegend=label is not None, - marker=Marker(color=color, colorscale=cmap, **marker_kwargs or {}), + return Scatter(x=X, y=Y, mode='markers', showlegend=label is not None, + marker=Marker(color=color, colorscale=cmap, **marker_kwargs or {}), name=label, **kwargs) def plot(self, ax, X, Y, Z=None, color=None, label=None, line_kwargs=None, **kwargs): @@ -169,10 +170,10 @@ class PlotlyPlots(AbstractPlottingLibrary): elif X.shape[1] == 2: marker_kwargs.setdefault('symbol', 'diamond') opacity = kwargs.pop('opacity', .8) - return Scatter3d(x=X[:, 0], y=X[:, 1], z=np.zeros(X.shape[0]), + return Scatter3d(x=X[:, 0], y=X[:, 1], z=np.zeros(X.shape[0]), mode='markers', - projection=dict(z=dict(show=True, opacity=opacity)), - marker=Marker(color=color, **marker_kwargs or {}), + projection=dict(z=dict(show=True, opacity=opacity)), + marker=Marker(color=color, **marker_kwargs or {}), opacity=0, name=label, showlegend=label is not None, **kwargs) @@ -284,11 +285,11 @@ class PlotlyPlots(AbstractPlottingLibrary): if color.startswith('#'): colarray = Tango.hex2rgb(color) opacity = .9 - else: + else: colarray = map(float(color.strip(')').split('(')[1])) if len(colarray) == 4: colarray, opacity = colarray[:3] ,colarray[3] - + alpha = opacity*(1.-np.abs(np.linspace(-1,1,len(percentiles)-1))) def pairwise(iterable): @@ -302,11 +303,11 @@ class PlotlyPlots(AbstractPlottingLibrary): for i, y1, a in zip(range(len(percentiles)), percentiles, alpha): fcolor = 'rgba({}, {}, {}, {alpha})'.format(*colarray, alpha=a) if i == len(percentiles)/2: - polycol.append(Scatter(x=X, y=y1, fillcolor=fcolor, showlegend=True, - name=label, line=Line(width=0, smoothing=0), mode='none', fill='tonextx', + polycol.append(Scatter(x=X, y=y1, fillcolor=fcolor, showlegend=True, + name=label, line=Line(width=0, smoothing=0), mode='none', fill='tonextx', legendgroup='density', hoverinfo='none', **kwargs)) else: - polycol.append(Scatter(x=X, y=y1, fillcolor=fcolor, showlegend=False, - name=None, line=Line(width=1, smoothing=0, color=fcolor), mode='none', fill='tonextx', + polycol.append(Scatter(x=X, y=y1, fillcolor=fcolor, showlegend=False, + name=None, line=Line(width=1, smoothing=0, color=fcolor), mode='none', fill='tonextx', legendgroup='density', hoverinfo='none', **kwargs)) return polycol