From 83f86b8dcc92c549f61984f0c03d2ca44198e218 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Jos=C3=A9=20Carlos=20Cuevas?= <reset.reboot@gmail.com>
Date: Sun, 1 Sep 2013 02:25:41 +0200
Subject: [PATCH] Made the julia mode a command line option

---
 main.c               | 119 ++++++++++++++++++++++++-------------------
 mandel_classic.c     |  50 ++++++++++++------
 mandelbrot_kernel.cl |   2 +-
 mandelclassic        | Bin 12111 -> 12502 bytes
 4 files changed, 103 insertions(+), 68 deletions(-)

diff --git a/main.c b/main.c
index 96f079f..4f39e02 100644
--- a/main.c
+++ b/main.c
@@ -25,14 +25,23 @@ int main(int argn, char **argv) {
     int res_x = 800;
     int res_y = 600;
     int current_line = 0;
-    int total_res = res_x * res_y;
+    int julia_mode = 0;
+
+    if (argn == 1)
+    {
+        julia_mode = 0;
+    }
+    else if ((argn == 2) && (strcmp(argv[1], "-julia") == 0))
+    {
+        julia_mode = 1;
+        printf("Julia mode activated...\n");
+    }
 
     screen = SDL_SetVideoMode(res_x, res_y, 0, SDL_DOUBLEBUF);
     if(!screen)
 	    fprintf(stderr,"Could not set video mode: %s\n",SDL_GetError());
 
     // Prepare the resolution and sizes and colors...
-    int i;
     const int ITERATIONS = 256;
 
     // Load the kernel source code into the array source_str
@@ -40,7 +49,11 @@ int main(int argn, char **argv) {
     char *source_str;
     size_t source_size;
 
-    fp = fopen("mandelbrot_kernel.cl", "r");
+    if (julia_mode == 0)
+        fp = fopen("mandelbrot_kernel.cl", "r");
+    else
+        fp = fopen("julia_kernel.cl", "r");
+
     if (!fp) {
         fprintf(stderr, "Failed to load kernel.\n");
         exit(1);
@@ -101,19 +114,23 @@ int main(int argn, char **argv) {
     printf("program built\n");
 
     // Create the OpenCL kernel
-    cl_kernel kernel = clCreateKernel(program, "mandelbrot_point", &ret);
+    cl_kernel kernel = clCreateKernel(program, "fractal_point", &ret);
 
     // Common kernel params
     ret = clSetKernelArg(kernel, 0, sizeof(cl_mem), (void *) &kernel_res_x);
     ret = clSetKernelArg(kernel, 1, sizeof(cl_mem), (void *) &kernel_res_y);
     ret = clSetKernelArg(kernel, 4, sizeof(cl_mem), (void *) &graph_mem_obj);
 
-    // Our screen in a linear array
-    int graph_dots[total_res]; // = (int*)malloc(total_res * sizeof(int)); 
     int graph_line[res_x];  // (int*)malloc(res_x * sizeof(int));
-    float zoom;             // Our current zoom level
+    float zoom = 1.0;             // Our current zoom level
+    float stop_point;
 
-    for(zoom = 1.0; zoom > 0.00001; zoom = zoom * 0.98)
+    if (julia_mode == 0)
+        stop_point = 0.00001;
+    else
+        stop_point = -2.5;
+
+    while(zoom > stop_point) 
     {
         for (current_line = 0; current_line < res_y; current_line++)
         {
@@ -137,8 +154,8 @@ int main(int argn, char **argv) {
 
             if (ret != CL_SUCCESS)
             {
-                printf("Error while executing kernel\n");
-                printf("Error code %d\n", ret);
+                // printf("Error while executing kernel\n");
+                // printf("Error code %d\n", ret);
             }
 
             // Wait for the computation to finish
@@ -154,53 +171,51 @@ int main(int argn, char **argv) {
             clFinish(command_queue);
 
             int line_count;
+            Uint32 *pixel;
+            // Lock surface
+            // SDL_LockSurface(screen);
+            // rank = screen->pitch/sizeof(Uint32);
+            pixel = (Uint32*)screen->pixels;
+            int iteration;
+
             for (line_count = 0; line_count < res_x; line_count++)
             {
-                int temp_val = graph_line[line_count];
-                graph_dots[(current_line * res_x) + line_count] = temp_val;
+                // int temp_val = graph_line[line_count];
+                // graph_dots[(current_line * res_x) + line_count] = temp_val;
+                // Get the iterations for the point
+                // printf("Point %d\n", i);
+                iteration = graph_line[line_count];
+                if ((iteration < 128) && (iteration > 0)) {
+                    pixel[(current_line * res_x) + line_count] = SDL_MapRGBA(screen->format,
+                                           0,
+                                           20 + iteration,
+                                           0,
+                                           255);
+                }
+                else if ((iteration >= 128) && (iteration < ITERATIONS))
+                {
+                    pixel[(current_line * res_x) + line_count] = SDL_MapRGBA(screen->format,
+                                           iteration,
+                                           148,
+                                           iteration,
+                                           255);
+                }
+                else
+                {
+                    pixel[(current_line * res_x) + line_count] = SDL_MapRGBA(screen->format,
+                                                       0,
+                                                       0,
+                                                       0,
+                                                       255);
+                } 
             }
         }
+        if (julia_mode == 0)
+            zoom = zoom * 0.98;
+        else
+            zoom -= 0.01;
 
-        int iteration;
-        Uint32 *pixel;
-        // Lock surface
-        // SDL_LockSurface(screen);
-        // rank = screen->pitch/sizeof(Uint32);
-        pixel = (Uint32*)screen->pixels;
-        /* Draw all dots */
-        for(i = 0;i < total_res;i++)
-        {
-            // Get the iterations for the point
-            // printf("Point %d\n", i);
-            iteration = graph_dots[i];
-            if ((iteration < 128) && (iteration > 0)) {
-                pixel[i] = SDL_MapRGBA(screen->format,
-                                       0,
-                                       20 + iteration,
-                                       0,
-                                       255);
-            }
-            else if ((iteration >= 128) && (iteration < ITERATIONS))
-            {
-                pixel[i] = SDL_MapRGBA(screen->format,
-                                       iteration,
-                                       148,
-                                       iteration,
-                                       255);
-            }
-            else
-            {
-                pixel[i] = SDL_MapRGBA(screen->format,
-                                                   0,
-                                                   0,
-                                                   0,
-                                                   255);
-            } 
-        }
-        // Unlock surface
-        // SDL_UnlockSurface(screen);
-
-        // Draw to the scree
+        // Draw to the screen
         SDL_Flip(screen);
     }
 
diff --git a/mandel_classic.c b/mandel_classic.c
index 29c6e54..a92c5d8 100644
--- a/mandel_classic.c
+++ b/mandel_classic.c
@@ -38,6 +38,7 @@ struct piece_args
     int max_iteration;
     int total_threads;
     int thread_number;
+    int julia_mode;
 };
 
 
@@ -51,13 +52,15 @@ int get_y (int linear_point, int height)
     return floor(linear_point / height);
 }
 
-float map_x(int x, int width, float zoom)
+float map_x_mandelbrot(int x, int width, float zoom)
 {
-#ifndef JULIA
     return (((float)x / (float)width) * (3.5 * zoom)) - (2.5 - (1.0 - zoom));
-#else
     return (((float)x / (float)width) * (3.5 * zoom)) - (1.75 - (1.0 - zoom));
-#endif
+}
+
+float map_x_julia(int x, int width, float zoom)
+{
+    return (((float)x / (float)width) * (3.5 * zoom)) - (1.75 - (1.0 - zoom));
 }
 
 float map_y(int y, int height, float zoom)
@@ -124,7 +127,7 @@ void store_iteration(float pos_x, float pos_y, int iteration, float x, float y)
 int mandelbrot_point(int res_x, int res_y, int image_x, int image_y, float zoom, int max_iteration)
 {
     // Get the index of the current element
-    float pos_x = map_x(image_x, res_x, zoom);
+    float pos_x = map_x_mandelbrot(image_x, res_x, zoom);
     float pos_y = map_y(image_y, res_y, zoom);
     float x = 0.0;
     float y = 0.0;
@@ -191,7 +194,7 @@ int mandelbrot_point(int res_x, int res_y, int image_x, int image_y, float zoom,
 int julia_point(int res_x, int res_y, int image_x, int image_y, float zoom, int max_iteration)
 {
     // Get the index of the current element
-    float pos_x = map_x(image_x, res_x, 1.0);
+    float pos_x = map_x_julia(image_x, res_x, 1.0);
     float pos_y = map_y(image_y, res_y, 1.0);
     float x = pos_x;
     float y = pos_y;
@@ -288,11 +291,10 @@ void *thread_launcher(void *arguments)
     {
         for (x = init_x; x < limit_x; x++)
         {
-#ifndef JULIA
+            if(args->julia_mode == 0)
                 iteration_pixels[x + (y * args->res_x)] = mandelbrot_point(args->res_x, args->res_y, x, y, args->zoom, args->max_iteration);
-#else
+            else
                 iteration_pixels[x + (y * args->res_x)] = julia_point(args->res_x, args->res_y, x, y, args->zoom, args->max_iteration);
-#endif
         }
     }
 }
@@ -316,6 +318,16 @@ int main(int argn, char **argv)
     SDL_Surface *screen;
     int res_x = 800;
     int res_y = 600;
+    int julia_mode = 0;
+
+    if(argn == 1)
+    {
+        julia_mode = 0;
+    } 
+    else if ((argn == 2) && (strcmp(argv[1], "-julia") == 0))
+    {
+        julia_mode = 1;
+    }
 
     int number_cores = get_cpus();
     int number_threads = number_cores * number_cores;
@@ -365,13 +377,15 @@ int main(int argn, char **argv)
 
     printf("Rendering...\n");
 
-    float zoom;
+    float zoom = 1.0;
+    float stop_point;
 
-#ifndef JULIA
-    for (zoom = 1.0; zoom > 0.0001 ; zoom = zoom * 0.98)
-#else
-    for (zoom = 1.0; zoom > -2.5 ; zoom -= 0.01)
-#endif
+    if (julia_mode == 0)
+        stop_point = 0.00001;
+    else
+        stop_point = -2.5;
+
+    while(zoom > stop_point)
     {
         int iteration, max_iteration, x, y, res;
         if((zoom < -0.02) && (zoom > -1.0))
@@ -393,6 +407,7 @@ int main(int argn, char **argv)
             arguments[thread_count].max_iteration = max_iteration;
             arguments[thread_count].total_threads = number_threads;
             arguments[thread_count].thread_number = thread_count;
+            arguments[thread_count].julia_mode = julia_mode;
             pthread_create( &threads[thread_count], NULL, thread_launcher, (void*) &arguments[thread_count]);
         }
 
@@ -441,6 +456,11 @@ int main(int argn, char **argv)
             }
         }
 
+        if(julia_mode == 0)
+            zoom = zoom * 0.99;
+        else
+            zoom -= 0.01; 
+
         SDL_Flip(screen);
     }
 
diff --git a/mandelbrot_kernel.cl b/mandelbrot_kernel.cl
index d42f85d..9a18e26 100644
--- a/mandelbrot_kernel.cl
+++ b/mandelbrot_kernel.cl
@@ -8,7 +8,7 @@ float map_y(int y, int height, float zoom)
     return (((float)y / (float)height) * (2.0 * zoom)) - (1.00001 - (1.0 - zoom));
 }
 
-__kernel void mandelbrot_point(__global const int *res_x, 
+__kernel void fractal_point(__global const int *res_x, 
                                __global const int *res_y, 
                                __global const int *line, 
                                __global const float *zoom, 
diff --git a/mandelclassic b/mandelclassic
index 4c61ca940241435362b07d52600dfa4257e3ab69..ba3552fc7ea36e2a45d9dbf029c32c4922bd0bbb 100755
GIT binary patch
delta 4692
zcmai2eQ;FO6~Avc3rXDMz1{3)HydD=n2_+bva}!{Vs>FyUNVUjfoK{NNYEfKfRt{4
zQDmdbDz9s3=%r(7r!_DwjU6j#aTKj%mOxC}+G()XMs%DQweuFF#w2`5(ARVBzNb&}
zhj-?^bAP{c?z!ild*9uA_n+(EZESN|pUCzTEU7mkgVPSrIlk~MCnJd_3Gpz3gZC78
z22}lkvw`sgPHNz0S!I=$Ra$n^W!(CiHMz`))a3eEMU(qdG}%XA<*2Rp*x%<qzFYn{
z@#vDfK6RBpvhC%b%##g<ord!M{07xBC8{6{fjCLmP5ga~&zr;-18*V7Q<uqW2-;7U
zn`RgyvMc2l2V~TBpt&^x=hE@%YIjqC!$fFB(D4S11nJcAdjFmjA0g-nJJCvst^o$6
zl?6H;+h}Epj!)M}kbWJXq2omzpGhL$YN~XKtZ!*jqvNx6e1neXb^K->FX;H4I^Lr3
zO;FjZQ&=?x$QB)+qvH?ic$<!I)A4p4-=0s)Gp0u#vv`OYON34Ga!`uVv*~A2M#`dH
zpi+v$@aElqs3G2r%H)s&rD!W+obBWQ!y6IfEGH8T*CNK*O?EN725}nV4u;DCM>gXu
zC)*gY5)E*6lPwG{LyWVU+{y4l#5kMD28QP###u~OF<gK+6S1FR7vftGFJRb;*o@f2
zaHho*A)9ellL8|QXn?bsBn)2{0OKqsN4^H4XcV!4c!=Q(h;i1E0}Kx%#>q-109T9B
zr-Azd_Xes1>&m3@wbDel2Xtcev%N9~3|Y^Q$j?AcjIT&QQ&H-9W1Zbkd<pp|u-&_$
zF24evq7z(S=m@kWVRv-G(CjLUhYQoYp8)-<n4huZ{E-CzazY9nfwt9<5~)Xev%4Sk
z{D}cyLaOeROTmxt?N`Ct|3(eWQnv@DMvgoWG|YdQ)Z3gdn^UnYQe%-)vE*Vbq+&Xe
z??7Y8RD}8mLY>{e2LrLy){C2%7#$uFXTJk528f8KhfJc>)FRFv5TgmJ7##xbSyB4k
zUhpP={FS1_n~#j~5y6BrHpW{*ZLurOS!2AAwKhdmwgE7$1i;%8fj9mDPWjsYQQlH{
za1^XEK(UGURvz3f@Lp3fpi>WmNswPX60htG`92Zjzs9~|f8`Io2#R>9GhPinT2`e8
zr3Y+%MjXGd@k-ySG)#xIN+^5athClCwwhA|ka~jBPt4-!VG}fp;ff1H&vN4Izlrfy
zbSy^COcS|tV!XOVY`vrBwGJ!H);1VeaCuWFgwC2YOEh41+I9gA(=AIM`V!qQLDuAs
z8+5x_%-H1R-3K5#LSHe@C@#GLok&y5qJ3!<ODDp9F>b2J_K*s0Xl49|OIq(zqSyBG
zPd-m?nDZjagx<C}AsXvDDvV9-(d{pA%`q5i>}@{wG<b$&nKr^#Gs^!cbZ`^oauZ*h
zsE|I_MkO~xuUfpkH`K;Il>ze%8N8UeiFXGHiM=0g_zQmMV(VjF^1IM+EMW^@XnCy5
z9N>c=TuiECAhazjH*#+(u%s1o1bv5Kw(3mp!7|nQPCg788RHjQn+AP}`vY}>`oQ|Y
zhWazUOJ&lz^`bN^Mn4&;T36PW@Z%EeJKbeN(U-`C)GdOR<G3ZAdIk2jYH@xS2r+ss
z8@{(ahNAj2r!b+osg2inqfQy(-7UzSLc|)-+U(gOrW;muEh`S2J%lt{@UxZoLx}Zn
zn)U<5jh(^H#OR9l3@`}if8*X#-96-_{n^V5*RRtwzMyb1Y+jJ&bC6`nG1Ph>0rnEc
z)}6yuTR4T5^CiVd{}*<?$;N|>O&Nbc%lOi)r*%v_{RQv4?QUS0|GNG$e*Fb>_JfqO
zfl<CD!P_qJwU@v()66!}%Y2D}H&M5+%wROp`-Kp#vd+y7z*Y?A@5A)iD_Om6dSx0t
zZ1qh8#6R^*NEi&n2jdV{2+TlOB&5*~t+#S%)R41-b_jWk@W6o|Z4|yHwnq2^&jA&R
z@3DnqSNKC~Q6J#vya~(ATnOQwmaDGP{|WirRhns;8@bBX-nv3zjGvc)=K_C1=H#I9
z!NQWIvJ@<w#~(TdfmFoj*Qa~nhbRc2_;OJ4hxP*vo`u!v@w&nia33rzyy)ZIu;ujz
z;(O9U@zrLC3|Q&Bl;j4ZXE?FFvNaqd@-VFU;L?u+(v<S>@<4pGGY~J0C8TTN-tTYw
z41VFPZF6t5ms`FKLX2&N<w9S#cpZO)F9(M^6-Ua3?7JjST&A4Wv2ysbqVxp|OCWAY
z*}O#ZJQ_pp{sB2UX-KXCD9^v5D4Y0_{)<0X*CqU1q<61_Ci2G1bYITQ3`ZL52(3lw
zY4mu`e8ZM;`tO|9Df72BZ`;yH8*Og$nx-8OHtpK7;~}rtOCvTr?X>yf^T%?dl*%?w
z<i*t=Jz{(D^U_CcSUXS(hYCMv;b|#-p||*UIMODd-p%D3m*YVce-TT@6r~^J=5a+K
zuwTx(rYP-zonI=-2<kzm!#G5c9+2%Imw=Q(;u=|UT~Y9`P%@zuC@oMCVGi0+fujal
z19Aa4mO<9g!}emsPp{I~>|R67D7|2x#d)YD*OBLhZr#uw5wnFg+35sNI$?M`G?44y
zYT%u7&u;;@fn6C$t)SiO@bm*etdSs*e+Qlp2-68sXr%%wpTW}!VJB&8?kw&e`ckgP
zbO%OGzsSv5v<_+wql(glC!9dG@T6f?))et)#zSeXrl+R_qf?E4GMcBR5?=@=2f*G#
zmpX1Wy#Td?w9(<P_<28QfLA4<haC>u;SjAs*!HmGu{Yi-=yk^o!=RJ8ocj$UF8ZwV
zYa`qx=sC`QvMX<a$?!xb9iHLM)EWlu)a71msZ#w^sYkzb`jOjBTik+0yJ=`QkbU&1
z`!>s9uAlWdsNPG)!Tz>P-Cf4zS)4hj-eDk`T7W|;guU=tfVvybd4%d!LUvKt%*B}_
zDuf}_0_=bd8B0e#eHI|@$XMVya(?EW2E3IN+-Z0#pVk)KnNfpjgwjEe7Tlh3ZC(m}
zc=DryJscOLM+@yno>ZZX2;W?&Qbwt_5w1Gg0KdyyA;ArLiq(xIhiGHfD~+&hv_XHD
zw-YvcTvhsFB)LSUjYW2&gN&e@=7bS9=uKGNn7WB=EL!io+5a3cFpw<LmF@pdGx&HJ
zKSaMMa%}u&A*nW-L2}^ql-e1rJXZEFe)8yk13dQP*A`^ToZ#V$BbX(+b(YuF@ARvR
zl=u%*V^Zj!XL%#q9BVdTvWj(5E;Cs;lXE{Ec#f#K)SQ5v1-yq$PMr^U-4%&$D7Gxb
za|0IGQCO*FHd#cw)T=b9Q~4GPJyyKirF}k{`h=?9$<u`;!JI5IqGqO2wHT*9qa)7%
z-B=-SZE#G%qFu~xp}yG;*t+oTqvuK-d6i&kA+S!>_|lj?A$3f<8?<<Rzj<*)(L>&I
zYA<-w#FB*Iu$*B2PcfUfRGStxPPSQs&D*Gre5QH0X%t@pZ!t{%+aLorH{?5STCBOi
z<9P8{hDBm;m?4#K6a0)BR{7h37f9a9N&i2rQXBZy2HKuomTExyZ|eAN`uH3NcOSm<
z8o(y?-$euKM*Wc^VP^9<vpK^2Odi;;n6^vxkWhQkJe+3Qe}NVUiHA-sY^Mh~c=E?{
v74DYpjXO8&-H`g-eR%$tLJDthgkRrn8+Sb%rVn}F%h<AmY}m46OPKs0%1yLi

delta 4430
zcmZWse^6Z26~4O*$+E(}-DP1}NMs`mAxH?D8XNMX+1157QsQ6?HMNkS21;6u6hZhy
z*=<)e>jquDnJLw@W?HL^?M!HxI+!}L`AKx9Epb|#s#9rdXWm9*C>VsO*?#Bld+{Z?
zGw+@EedpYB&pqe9`|f-1zQ1md7dWPJ{ES-q&RDqUxr$d-cDV$}G|8Ap5C+^A;ObTN
zy{;C)_qw>zQ09oMbX=w5F22UFeTgSe7%@+tUp#2CD^0U5{)T}&pZameOXt7)#fI1R
z*FLJu|E_(aV5@7S#n@ro)KS=?TBbz<!mzOdw(T1J9>Ke><ClWBGvaw+W|h(F9NCn4
zqcIeBr`^(oOu7Lyx0&Dyb^NeuzYrWHLTjTs-l&m~2_3KZe=5y~7&)Rrv^JsZKtXBk
zjE<*1TAS4IW{reYbbOYMpVILb7Wzgrr&Hv7Lu2|a0@0vW9dFk0l8(3Oc$<zd(D8PS
zZ^eU0r*LQr$Z{Q@tK+M6yi>=o*753M6Iy>E-;}jDRBH1uxj)gBsjLMpM~<71r8}oq
z6+)FNU4*+-?8g(~c2XuMQ)s0Q!c^_#guoWUROMtsU?X9wZgN=QOA>Gf;R6Cs1BcqF
z%E`DO&XECCHyIZAEMZD;vP0na2~#zbEdswwn5vj;5cnO!7Q%jkj}pFtaJ9gP31<`b
z2z*HLgjhROHE9#XYh++0%mn^5VX9(s>N1GZ7YI}Jl9K{IPuNa)Lg1$fQ?Zf>;Knfz
z9_8qN1NQ{J7ibJ@y<3ia+N(6e6B1se?8ie$j;((Pz2xZ8w{CIznJ=Nd0(xu<o|QM@
z1Vk&Q96cI6IW`1RIym7=L>q^che`SF8yLa=);45)%!5=b{{)M2)<sWr7b)xD3I20b
z|M#c;-$DLa{bCk)>or>GN*ncwW`1Es7!TAAjqQV-Jm5T$Ptgx#_;H<}Ny(A3Q}VLE
z$&tw^IocYQmmQa*zwL!L#R(77k}b!JVvN`BJ%vRvBiS+=0jDMV=6I}jsMhzf9Q!3%
z_a~%hUWB_?_fXvfq2n`>{dO~k^9Q33ILA#ii}BktXhwm^aYOVGn;iQibqhux$%bt|
zm0pe<E071W-;-mFVc)Dgu<q#L12!b={MWEs+d727g<bc4-v{-v*6=4;=wYa79_lNr
zi>>?6H>=FR5=Bt5qyO?H#=@AfvgB+2S+;DHius0|8*dx`DwPUE7uH3_GwN%;?DEU8
z%*S#(tlltrTkMwA19zk%bA8g2&tV?Hn%TP8<L1C<3gEHBzPEkvVt}W<;xFYaYOar7
zP)actvXglNOF8n;Uc<EX{n~h&<V>_l_a-2udryfhb=^fd{5h<sGK{rJRh#3|lXoJy
zQdQ=Q;751Tc+$qHWDUqTueRogX1^4HCM%E^B>7MDLdpY=4=c+tgZ@N!asLP1=cM8w
z3jHD&d7$--FY$vwQ=mDpJ+PzsSorbb?0^*f@RN6Z=j)=Ux69G9a^&NwhOKqu2@kD=
z@gwnE65|OAvNVq2K5rBLMt5LqS=A^nr>}w>nX_W4_ZXXxjUocEr+nT#Myk}LR2(L5
zl#sCJ=d8<Wh3>4X3$K6s1`lK1_C|P7`Y@U1lbNrBGRz#K?MIHRpD@FqtLTb*d*j$1
z7vCnWF+Tkz-zVKv@(T76B=@w)X#!72k<h$B!8V<y<nH6kY#VawbxZ^Aw)s|A!Jrt^
zQq!p!Y5O_qHiMlCE+(d>Z79-t>E81|LX*tkf47wz$1`}!w#GP;!PnYr`2mM}l{LNY
zbM8i0Q3vUXW>B@U9%pU;tn`eX^Z}{jB-U}Z7Xuwp8s_<Djv~W6x8&X&nop&$a<_VI
zGg4&&*E#8+Vo-vno!;{L*m8T8ht(+2ptrKGIUZaKvacTAY#Yt7CT}@h2EASs<_XF)
z3&b`AVm%qr_iJMt@!Hs+izcIU(!mpkGk6u`!Oa6*{Y=TkG78rGClFoG{pX1JgkkN2
z-pcNm`Q_ZI?H8b;U>8!HB4FQnW&NB85y1zg%1Kgh%n<YNLNOlnmebf2eR2^hh;eg1
z!izj&Ozwj><zX1LN#&<LdA^PJ+P{^57j#JP2S|Hfqj+vw+-T%S?7l@n(o?^HGt|<D
zy4cni2m{-y0wcFNN(@H6(y==m8>)G1E6kLi&GP-uC0QSt@Qx0YPMG*f=W=88BERSy
z;L*Grm)3T6KHTZqwcAs=)6>=7*?Qkj3)|AVduMCsuH8SZs;Z*E(fD{?fZOstp{8e#
zUMz3AumQ(DKX_bWdLA<rnby+bjJ_HzpQlm@$O9KssY%p!0`rRlCs0TQasrZ$LsO6i
zkp4@lR5fHAQicrUU21`>#s)DA*?=U`-q~?EWl7mEU-}BI=BM+Q8duNpOZipC!5LoR
zDl<&*yId~!DGXN(BPLsITdZbAUwsN~g8$UzT4lzdb);1rX<x?W2Y*&0A?tCyiK_rX
zy~@9acZquoJgwOXq!$5F;94UeUf?h`5^?qjT5We?<Vx`Cz>m&gg7n1?SZz-l<(vhZ
zpEW&`F_5`;L1dvxGFh@0ngjGb_TqXF_S^V#1tnSg@RWc@BR_4(wG3PR&>~4z{(eE*
zN-wq`YHo3-<8_O8DpIjom4&JdyLhL&&p7Dj)9%YAoVMXKZ^T?wooT#3i?8uiSu~Gc
zC-3y!>Zn!+sHVxGO)-AhQ)Sl<;@Sbdi_d#j*ay{P`k;D}`jFQwSz*y)!co6-$*q>K
z8X`hD1h_|q*sG=ibw2<rZvs@$mvJ}m1RPW$;y(l^W5c8-7xHO9+A*o&w7#<BJ6YS*
zNiX83%X0bN60fnVkRO7Al?XGohtHRMJ8RF)X*51k<Lxo{wE3sW^xecHrX}KP-o`xq
z52bk~I~yb!a$Kpv)pof|oA%YblWFt6n$KfeBCh81nU;~Oc^4C@QHyj%XTraZ%zsjr
zXTm=WN&c6zP3~)ltbn}{udIjrmbyX>>I&1GY*DrGYxs8X1}wZ5HS1WWqV)*=`U(FE
zJdL8Q;IuVk^es2AD~aP5%9h;pw&3in!{yh^g#3r#DQVnNUUk2=95kDG@!+^F*A@}m
zl_4~9wivvDU7zbR@E%HCA%CH~Y9)Q=snM20tr2PmW}^K8r&-3DRKCH%FO@%9c<o9$
zsH!Ke{OM&ueyhpho>i$193iM_e=I5PGr0^-Zm)0|N_n8d<qlwQyKWic8!B3Z(EbQo
z`mfbq2vx~?$b*|hOYzbnklL#lmD+EF|KAInglaRTR$1(vusKQE&?Jp+y#BxWQ<ZN}
z*MM+~+q~+Hjz-j?;?7tWc0Y<_Sn><ztjZUHw_z=5Z?H<Uvb22#Yf(FBTXv<|;o2X3
zAOGLYF2fqS+i!<W`sckx_%m}uWwBu=>D)b{|Me5<6I!nJY*(kH?V>M}*7rKJGzps;
w(`Gf=H&vbvg4D!T4|!zd-<2r?|Fmk3UtDoU>e$iC9=dPWZnk6B?p<B%|A-1jMF0Q*