From 3210b268aa8355b858b438f14b3275eb770a197d Mon Sep 17 00:00:00 2001
From: Oli <olifer97@gmail.com>
Date: Wed, 1 Dec 2021 19:56:03 -0300
Subject: [PATCH] rm refine module

---
 fftma_module/ref/FINALrefine.c                | 200 ------------------
 fftma_module/ref/PyArray_NewFromDescr.txt     |  27 ---
 .../build/temp.linux-x86_64-2.7/FINALrefine.o | Bin 37560 -> 0 bytes
 fftma_module/ref/exmodule.c                   |  96 ---------
 fftma_module/ref/refine.c                     |  96 ---------
 fftma_module/ref/refine.so                    | Bin 31168 -> 0 bytes
 fftma_module/ref/refine2.c                    |  84 --------
 fftma_module/ref/setup.py                     |   7 -
 fftma_module/ref/test.py                      |  12 --
 fftma_module/testRes.py                       |  56 -----
 10 files changed, 578 deletions(-)
 delete mode 100644 fftma_module/ref/FINALrefine.c
 delete mode 100644 fftma_module/ref/PyArray_NewFromDescr.txt
 delete mode 100644 fftma_module/ref/build/temp.linux-x86_64-2.7/FINALrefine.o
 delete mode 100644 fftma_module/ref/exmodule.c
 delete mode 100644 fftma_module/ref/refine.c
 delete mode 100755 fftma_module/ref/refine.so
 delete mode 100644 fftma_module/ref/refine2.c
 delete mode 100644 fftma_module/ref/setup.py
 delete mode 100644 fftma_module/ref/test.py
 delete mode 100644 fftma_module/testRes.py

diff --git a/fftma_module/ref/FINALrefine.c b/fftma_module/ref/FINALrefine.c
deleted file mode 100644
index 27fbcb9..0000000
--- a/fftma_module/ref/FINALrefine.c
+++ /dev/null
@@ -1,200 +0,0 @@
-/* Wrap an array refinement function using the Numpy-C-API. */
-
-#include <Python.h>
-#include <numpy/arrayobject.h>
-#define NDIM 3
-
-
-/////
-/*
-int refine(double** in_dataptr, NpyIter* in_iter, NpyIter_IterNextFunc* in_iternext, long nz, long ny, long nx, long rz, long ry, long rx, double* r)
-{
-	long x, y, z, cur, curnextlayerx, curlayery, curnextlayery, curlayerz, curnextlayerz;
-	cur=0;
-
-	clock_t ti, tf;
-	ti=clock();
-
-	for(z=0;z<nz;z++)
-	{
-		curlayerz=cur;
-		curnextlayerz=cur+rz*ny*ry*nx*rx;
-		for(y=0;y<ny;y++)
-		{
-			curlayery=cur;
-			curnextlayery=cur+ry*nx*rx;
-			for(x=0;x<nx;x++)
-			{
-				curnextlayerx=cur+rx;
-				r[cur++]=**in_dataptr;
-				in_iternext(in_iter);
-				while(cur<curnextlayerx)
-					r[cur++]=r[cur-2];
-			}
-			while(cur<curnextlayery)
-				r[cur++]=r[curlayery++];
-		}
-		while(cur<curnextlayerz)
-			r[cur++]=r[curlayerz++];
-	}
-
-	tf=clock();
-	printf("%lf\n",((double)(tf-ti))/CLOCKS_PER_SEC);
-
-	return 0;
-}
-
-*/
-/////
-
-
-
-
-/////
-
-int refine(double** in_dataptr, NpyIter* in_iter, NpyIter_IterNextFunc* in_iternext, long n0, long n1, long n2, long r0, long r1, long r2, double* r)
-{
-	long cur, i0, i1, i2, curlayer0, curnextlayer0, curlayer1, curnextlayer1, curnextlayer2;
-	cur=0;
-
-	//clock_t ti, tf;
-	//ti=clock();
-
-	for(i0=0;i0<n0;i0++)
-	{
-		curlayer0=cur;
-		curnextlayer0=cur+r0*n1*r1*n2*r2;
-		for(i1=0;i1<n1;i1++)
-		{
-			curlayer1=cur;
-			curnextlayer1=cur+r1*n2*r2;
-			for(i2=0;i2<n2;i2++)
-			{
-				curnextlayer2=cur+r2;
-				r[cur++]=**in_dataptr;
-				in_iternext(in_iter);
-				while(cur<curnextlayer2)
-					r[cur++]=r[cur-2];
-			}
-			while(cur<curnextlayer1)
-				r[cur++]=r[curlayer1++];
-		}
-		while(cur<curnextlayer0)
-			r[cur++]=r[curlayer0++];
-	}
-
-	//tf=clock();
-	//printf("%lf\n",((double)(tf-ti))/CLOCKS_PER_SEC);
-
-	return 0;
-}
-
-
-/////
-
-
-
-
-
-
-
-
-
-
-
-/* Wrapped refinement function */
-static PyObject* refineFunc(PyObject* self, PyObject* args)
-{
-	npy_intp r0, r1, r2, out_dims[NPY_MAXDIMS];
-	PyArrayObject* in_array;
-	PyObject* out_array;
-	NpyIter* in_iter;
-	NpyIter_IterNextFunc* in_iternext;
-	double** in_dataptr;
-	double* ref_array;
-
-	/* parse single numpy array argument */
-	if(!PyArg_ParseTuple(args, "O!lll", &PyArray_Type, &in_array, &r0, &r1, &r2)) return NULL;
-
-	out_dims[0]=PyArray_DIMS(in_array)[0]*r0;
-	out_dims[1]=PyArray_DIMS(in_array)[1]*r1;
-	out_dims[2]=PyArray_DIMS(in_array)[2]*r2;
-
-	ref_array=(double*)malloc(sizeof(double)*(long)out_dims[0]*(long)out_dims[1]*(long)out_dims[2]);
-	if(ref_array==(double*)NULL) return NULL;
-
-
-
-
-
-	/* create the iterators */
-	in_iter=NpyIter_New(in_array, NPY_ITER_READONLY/*|NPY_ITER_EXTERNAL_LOOP*/, NPY_KEEPORDER, NPY_NO_CASTING, NULL);
-	if(in_iter==NULL) goto fail;
-
-
-	in_iternext=NpyIter_GetIterNext(in_iter, NULL);
-	if(in_iternext==NULL)
-	{
-		NpyIter_Deallocate(in_iter);
-		goto fail;
-	}
-	in_dataptr=(double**)NpyIter_GetDataPtrArray(in_iter);
-
-
-
-	/* refine the array */
-	refine	(
-		in_dataptr,
-		in_iter,
-		in_iternext,
-		(long)PyArray_DIMS(in_array)[0],
-		(long)PyArray_DIMS(in_array)[1],
-		(long)PyArray_DIMS(in_array)[2],
-		(long)r0,
-		(long)r1,
-		(long)r2,
-		ref_array
-		);
-
-
-	out_array = PyArray_SimpleNewFromData(NDIM,out_dims,NPY_DOUBLE,ref_array);
-
-	if (out_array == NULL)
-		return NULL;
-
-	/* clean up and return the result */
-	NpyIter_Deallocate(in_iter);
-	//Py_INCREF(out_array);
-	PyArray_ENABLEFLAGS(out_array, NPY_ARRAY_OWNDATA);
-	return out_array;
-
-	/* in case bad things happen */
-	fail:
-		Py_XDECREF(out_array);
-		return NULL;
-}
-
-
-
-
-/* define functions in module */
-static PyMethodDef RefineMethods[]=
-{
-	{"refine", refineFunc, METH_VARARGS, "Refines a 3D numpy-array."},
-	{NULL, NULL, 0, NULL}
-};
-
-static struct PyModuleDef crefineDef =
-{
-    PyModuleDef_HEAD_INIT,
-    "refine", "",
-    -1,
-    RefineMethods
-};
-
-PyMODINIT_FUNC
-PyInit_refine(void)
-{
-	import_array();
-	return PyModule_Create(&crefineDef);
-}
\ No newline at end of file
diff --git a/fftma_module/ref/PyArray_NewFromDescr.txt b/fftma_module/ref/PyArray_NewFromDescr.txt
deleted file mode 100644
index 1c6afa6..0000000
--- a/fftma_module/ref/PyArray_NewFromDescr.txt
+++ /dev/null
@@ -1,27 +0,0 @@
-PyArray_NewFromDescr(<PyTypeObject *> np.ndarray, np.dtype('<H'), 2,dims, strides,pBuf, np.NPY_C_CONTIGUOUS, None)
-
-PyArray_NewFromDescr( &PyArray_Type, PyArray_DescrFromType(m_table.datatypeid()), m_table.ndim(), shape, strides, data, NPY_ARRAY_WRITEABLE, nullptr));
-
-
-
-
-/******INFOS******/
-PyObject* PyArray_NewFromDescr(PyTypeObject* subtype, PyArray_Descr* descr, int nd, npy_intp* dims, npy_intp* strides, void* data, int flags, PyObject* obj)
-
-
-nd dimensions
-described by dims
-data-type descriptor of the new array is descr.
-
-    If subtype is of an array subclass instead of the base &PyArray_Type, then obj is the object to pass to the __array_finalize__ method of the subclass.
-
-    If data is NULL, then new memory will be allocated and flags can be non-zero to indicate a Fortran-style contiguous array. If data is not NULL, then it is assumed to point to the memory to be used for the array and the flags argument is used as the new flags for the array (except the state of NPY_OWNDATA, NPY_ARRAY_WRITEBACKIFCOPY and NPY_ARRAY_UPDATEIFCOPY flags of the new array will be reset).
-
-    In addition, if data is non-NULL, then strides can also be provided. If strides is NULL, then the array strides are computed as C-style contiguous (default) or Fortran-style contiguous (flags is nonzero for data = NULL or flags & NPY_ARRAY_F_CONTIGUOUS is nonzero non-NULL data). Any provided dims and strides are copied into newly allocated dimension and strides arrays for the new array object.
-
-
-
-
-
-
-https://searchcode.com/file/101698327/numpy/core/src/multiarray/compiled_base.c#l-1251
diff --git a/fftma_module/ref/build/temp.linux-x86_64-2.7/FINALrefine.o b/fftma_module/ref/build/temp.linux-x86_64-2.7/FINALrefine.o
deleted file mode 100644
index 6a621afd0f5c18dd8bc81df950e683f3756447e1..0000000000000000000000000000000000000000
GIT binary patch
literal 0
HcmV?d00001

literal 37560
zcmb__3w%`7wf;FXVHgM`K^_KQ42l{MLW0OsG(!@YKuAUs5L7x$k_j0}X2Q$_5~_$i
z)QBkfXl-kiwmxV}m1-YotHr7nw6~YqR<X5LYHf@4MNwO{)!c8bz4n=tO>%qx|KHs|
za?W1gUVH7e*Is+=ea@V7xURB#p2xN=oowq8D@$qAvQBKv=F2!+W}R*oTEWYMOM^>-
z^-JnPo8RdSZ9Q1M?R(=6fFb|Z*C405`#@Fqr`6rhhw>lUGHyTc>h9`s{_5_EaWfWm
zzaIk866&rWR~YKfj4KE|UokFB!kK#?Wq^bPIT1pwvQYO4F5OF|09J;!T{A93sF#S%
z*NmIN__F+4?}P%2x<3tdUo$R1f{IAcqV9i$w$%gZK14OjsD__wl;q!9jT)uzR(Bs=
z)csji_q(C)=R(^mt*Y){S9|YQ+HYTPS8YGzWooqGcz%Hu+Ey{H@Ob_(Z_lAcGyj!;
z>j^ZmaC`lqq3q%OF=ZC>Y~Jq&ixryr(~kG6yWgnV{G44qvvr)m<7?1{+S!&FR}z}}
za{jGPqk7f$Ay6?hE>N{S=ata-mqXphn~U?uY$5GQUe2MA5?FR3H1j|QYFD5(Oc4le
zMkj5ei$XJD;ry+`dU|@$7cw<8UdCcWGyk!nooo<71Aj+rpM}OBH|0L9<Srw*A(9(r
zxnYvqlE3xF&w6@7o1g70-GAkZ;8npD!OsOF&s2AR@>T#1(trluN>Yw2{nnYBBmGF&
z+utPxF=7DTYDQ817%SU9cv7hQA9bPa7g29R*S}MOA!!P2%Q=@uFVwvS9mE*z4{aO3
z4BM;WJDcC}hqi?aLfbJG;nMxp1XhD6I7&QQ_g}F*e@hzGLffta_qKE(w9N-6${EK>
zgtjg9hq|98hsb}VVjS{hzik@|Kw^|}d^+^>5ihw>35XLcG_-YpRrj;`Th3sytV?C-
z{?da_c+lcgmgRx2ulMY*2BOui`{79WTmFWU5N3kXy!7uqJ36Vf`<W1pj_S_V-{o(4
z9DF@HPEZ-lhnevn*IEz=)`z-(+mjhrsQ;?BRTWfk3!?Lblu{yP^Rodc6WUI04%;m1
z{=?%#NG&gEu&uM;YOBRYdsm~GqtF)v8DJ%!$%M5c-S7A82tb9+3(Iog%Asxdu;H6R
z7oNc*(J>)({j)SSH2%+I$4}MD(xv;~euj0+AF~myTZj`fwA<$pU3#$eCDJjt?;~iZ
zdZ9|Rne(?SA{$FB%*gKNRng`)%ZBwbSB1!GYqJ9pnpy78By>-t;<g!5>jFD7E41xQ
zYFr`PX2Sq}Y2eBgZ&69Jo^P1{aESMN!X<CtP4hXltsM00phd>pPvnoOfPUSti?rNO
zx47nIJ)x~Hg}R4l7QnjsTY4b3r(n(?3u>!@>Zt13akmdfJ~2p;{2T5jt{ppoQk|cI
z+I+(9XdxA^T=C2fcl+GT+uz{!d*J;%WxAh!!qq}(+l<ovq8NAiEk932&TE@)o6Q_=
zy-Sl?n@jb<ss7eM6e?|>7!0K_UlixtclBdvW#YH)&I*0|oBuOw8(3?{wzjrfsaSJ7
z5wj8<ZS7sfjmcE3xUHi#6OX1+(Jm_zT)a5AJQ567xw3&oG81S{b|jkIU_2eDSyEjc
zXivrynOMqN%yQF#XyC$%085|5GK>3eD$pE_x5k<RnPec|)}BmdEJHXo8ty8uZMZtt
zm<c5L1mhsJC7U{0V}ZtGTRV%0u8PJJ=}aIvw<@qUmP*Hyi9ku`xt)QAj!Yoa0!1~N
zY*v~DfDy;Pkr)2|lQ%C|AFQ4gI7!WPpd+z5kzAJu#1c*MXhOAaS^e|3#xoi8RIyLw
z{xrohum=p+6mM>hrDBPU6enV77(Ww8btE$JwwM;~I%5|9Kc(M|v8)`f;peIp=X#|u
zW|B@6o2{I~GjR0!5!`UH{e1E$Dn~%B<Mr?y$jKQ1=RpOgSj2$Manw@TEtC~$1N;ix
zLzvIMi%M*krSK=5YQB;uQOK<w9{9IeC}uqd9&Z3F`G0(wvI;*z`2v;ynDT`x{|V&-
zD*sQ)7peTeC|{!TpHhB?%6~@rvOzeYahr&cB5WqYitsS8Oc6PD3hJ#?gwO5-(V_^y
zP1?3AVh9r%MdaD-DCtzhDfT=N8x=9!t_QJ65hLskkhWP7h4!r=Zc@Zp`}1Jip@{MJ
zb`W<fVgeJp6mbC)yA?5!iM@&_wkM&^K1EDs-S#V@)V>q#JZn)Dvj%>F$DnZFJ=95M
zqpwm{^7m4AWemRF0Js3-8U*I^(*Wk~!hzHtasx7i(oiybzea}x{T7I!yTFjwfWuRu
z@~KP8g?!DoA?x!~k$9Q-lh6N0uoiHoPpQ)HodZ`NMl$?)_Cq+K#IgqMA~u!%AvNR8
z;mP6k^NIWFA3S|AmRj!X;F%3nP331%J_m|SLh&iY!~6xgS8vk%xvNJYtN3$QkAZ9Z
zt>Gl$COfynCYo)WTWK?&CYVtZR8pXbc{W+8P!aRFPQXnRDPo}=8VYY#M713QF+&lH
z>@R>QQ$&q@07OU;wRR7Pup+{CZk}Z=Q$)Re48%%BEV2I%qD2u)?SFx2SHu<esUSKP
zai#rv^lhUeu3~97DPo2FW3;nF5udX^0CBe>qP8dBvUVw=$&SNHyA{#Q-QKH+tL+z&
z*{6s$`vJ7GUlD8U8PMQ>A~I~9LyG9M8&Pst5nXl?#1TbYV^^Zlql&oJ4ud$Rhz<6A
zAWn>+(aOD!$I7x5wb7;^=@;sH`}2@jAk-%NMo@)9ZB}C>)E4^=lobiJ)qVj~iBLD#
z_ko%r)Q$E9Fie?HH*qgRLUnWNVWDnjYMD^mm|7{+EljltwVkPUp>Aa#=@e?Gy$Idh
zXv3Mwcj08{bn6NyyLC1?jERh69`&96@_D83`OigORk$IFtk1s+N9H0QU@pa~m~Om*
znQqZcG=NzqH|I@znezw4*#L!})11;=_aQESnaXLB{K+W4`52e~oyw{HU`)?^=Kw{i
z14^mwr=!g0AB!Vvei5-!aX!yW#dK>KGtJXXkzP!0&Yf2?XFG9H%iX;>Nw740`<=}B
zWzD&-Hz&33Hpv~&F#AEx{)WTO%Ccaee>=qc{Ipmp#lQM7*Xp5KB$3v_Z13_(h;ZIZ
zoT}1&#o+S!r{c(!LQW;hp@bwU)?Le)HKJKN9ah%=D9NWkYh~{JBQd*1Hgq&wxet(<
z)<|put=tDmbN|TJSi;JEh$^Y2J9oEM7w-5ivSfhTwsOBsdEADe{2uD4TE=s~L%OMD
zJojPBt7Sa*5z3b-{=K9E(Wv_<>9EYorDo0r%*&E3-Fl4rp?diJ*(ikn<H!Sg1v1Kl
zoFi>eOcE6T4^}h(8XQMm--}lmUb%%CZpV=YKIkyHoa69!ndyg`i2|Q$l})Fbdi7;y
zdQCHZ<S?;J<~((Vo%_2#=jw&%&86gk$^w6)ykF)2OnI^VUnozG0xKP*yjcD%$`>j8
zuQb8M+HaR2FV=pC@?!0GDKFN3kLG$<)qS7vWh#G+WUf^CzfrzL<&RUoUFAQZJPio+
z`H=D(RsJKI^*5>f-$~{UmH!9j?^gK}l;5TDA5(s}%6~%ny(<4NlDSXiKc)PBmH&+L
z2UNa?^0K5@HW9L<*u31ylH%bKSyH@A$dZ!7ge)lon9xg#k14&R3}#9%DY;DPC1nUx
z0SVkgnbJ#2flWQaT?XFh<PNjR84AO>!|e^&=OYi}jj$<>3Neza`&G#(y8yGaKoO@h
zQK*PQn*w4$5u<Gid_{^l&8FB^qKMON^7R>tIFm(}DdH@a781S2*vH|Eg@w6e?H@UM
zo$c1^9JgNMY>J8nQu$o_X4EMx$vxj5GgSAZ$bJtj(y8%wA;wEOHNidygmmfxy9R`G
zYNAbQOQ$B;D?t=V8^!kXXhRvX)SixJeEvB&j-q{)+KH;dDRw;<UO|Pl&*JGRH7~T+
zaN#-}d6MtY4EdassPQNiSJI~1U*`(DsX`c|L|Y{_;VH%m*F`ohv#Q24``29Kw^W1V
z(S9mh!>!G9o7N6hZH7(ZN9AYQw6?4KEc+a6LVW)5IOfG6jy6<oA+zlYW>`QBq#^CA
zRC`3_bC;B?i(PW(xbm0Sv;j~&m%4N+bLE5fx4D_`>Sk!Clux_jY+E8?uKhBzzNT4e
zv#MArsG~C!_k4bu06ssto6k>vL+)Y)hMY^=Kc3w~&ZE7x<j>znu)>R|Qw1tNf%o4K
z;Yays65;bNM3$u$4b7!3BXbTN!Xq|vRV-r-9r_D`N2V#C_W|++Bb$<q*3f+7Az#JC
zPQi8@0}!_q<;pr=v1u>nE#ONL(#;xj3Jtx=P9zmLJJkOwvOfRoIP#!S6UtB|RG0>R
z#!NW{;wfh$leMODYM=iBFrAjO2(|L|%mv5kBtRW)&U=)&ydMmLZmefsseLJURQ___
zx1T|kSc+%*xhObU3L{=V=J0U<AANir$VWdP2k~(*A9MLQgpWh{n8(L_J{ItC7#~mJ
z<8VHX;NwU>j^g8~d@SVSXg;3C_xCK-#gSqG&+Qd`>W*S^3D0A36rV+1y^e+SHj}cV
z=O=geDT6CW|8;Xef2>bdWpK|S5`=$fc%1zLGPscgkMT^fdq7}1qu>I2`YAXoVxm0<
zgbIV6V!IrK#8A)qb}pz<qj1Ql1Dxk6M|zQ61SpSilBOA`6`%hMAi3mKb;!tW2?c07
zYDWOE|D~xKF}Xser@_9C%E@wfQ#r>Alc~wRA7!Kd!x2C}s6=s&4<y1Zd6oStRr<iG
z#8OSJ7P|)wd5a)$I4;3Y(y!fq)1wW%#=Zb7tZS9S%WXR(Y(aqPd)mGdY@@dH=4Caw
zR_~ltT5Ip1N{={|dULsZc&+^-aIrqG5<4#+%uW#`pD>#LII&Z#^i1&_z@@74g9|+z
z2sM$tx8RxOp&+2vM$g5bU67*IM$aWoj8=8ZnK(@ml}wzji1{9}@)?S#Vj`f3g-nc7
z#ATkTg?h)}x!gk|uXYTctC*M(XomlIW`U_p)o5fqq=>kOd@-zuHqUbCyi9C(o&9U{
zlMQnm84jpwWa*+q8Ac7!Ci^qZdFn_M^)|MM+rmt$!d81S_*nEJrveYZGo(^$yS-d<
zu5&mYCighrW<TQi*iQR4Fpqk;pBnCJci2x5>#GiHZ{_o;Ez#vJ`%mCvtxp}LyUFZX
z)mgRdfufEzw$t|9XD>TdJI#Zf7Y};aPM~178h}Ua42Zpo_`ba!#6Cs*z}^UAzaoBY
z-vHu(B7SP$0^*P&4%oMYIIM{OwC@ISL=nHS?*(yG5wF?zgE*#$H|%eLP&;SOQF|{4
z**SaOvwr|ecFvyn?Vp2^-LdB%e9J>_s(3!L-#}TJ+BtiA?4zJULRp@J7#X>#;_-Oi
zLz&!E@pwHSf?6q@Ii62IwFouXLp^J^t=twu@y9`N1g0ZJ1D~Jbq|Z+g(dVbQ=krr=
z^Z6+TVJE`JT0T-B@cC)p`ur5je139ipP#%h_h*nye=OYtl$<#CH6ZlopyvSf`N<1?
ze(I3VznSyoWj;SSkaF%J9`X(j76qPmo}S9kYdqDMg);6bo`rtY90Jhi$Kwl9KjEQZ
z<MWpRVJTtHM$aEWC@I&wL(cUcSIF7yp@>xwz|T&etsEmH;@s#V2MwtDw=hwpi2K|D
z<7=$n3{~<C7Eq>$2UN%q1t0V*c9eaTHS_u3b&%)0x!An%{!Q8AfcFegJZ&yOCZ7)I
zUJq*I(o6}0W+RUQS^(l~lv|!NeBL$?=ioOnB0dKDI4NkL^SN}er*odNNaYw6O7=b9
zJI%7r<Lq?GQWBC|&yv(QrSlUd5EfSisM2246c<&L(wv^1R#8ebdUpDYQku`Rv!E!2
zI!~Z&&(6Z}jUAw7YU$(Skl>;q!bL%bLroxk{%^<&q-$2r&Z8TA6YV8AJCAI@V&U_X
zl0N@=*m6$d=B}nSM5kBfH!%C2!`OH#8#)XISxlB`r6Kk1+<R>?nMSC6pwKU&_Fp@h
z8)_|u@x40_3E4(N?A>|bTGi`>Rzs0VJ#dO_*8``>HF}g4nW8>XDxrhDhmMp==gV*t
z!CQZv?cM*09pjXGDVnlU^IS+=l-e;32c!G2n)d^{exhAZ8+4T>+KWl&iS{z;^hA5P
zcjvxqt*N^B^Qd}}CW1<aCk;vU-m*F4Cfb+TQy{@EL4CU<%Xx)xt|d-)LBc7p9QA0L
z+`@Rhhbml#2C|}lpi1luQ9#C~<}T*upseQ_7v~FE*^d4<QUsbt9q{fRKy%F=itYUL
zG=78p7Sl7h?Arrqv{Zf%=dGDkVABs%oC40=`~chfS>#RbWIw$f@N9xEeoM+O{!nKu
zbqSH}@(!rA5g<k+I~3XCub}t=Rm>J~jaQcGAuv75OzvW5uns%*v$lLkFxc0z{-Y3+
z1KpvD$8Gn*Eaq7AN$@%qsr?FOa<|_<Z~IBSbF>9XR4=bYfz!?b6uzL8$=a$HH@P7^
zP`ROZFdlw!wD=V&U-CaRb19X1GakCnot2*Z3@&v?qbkRB+fwuW`-6aFnaZfDZ#Vs@
zsj}a6viTR#wKAgx8_Ba&_J(YBJ7t3qz-id&b9MHVm+W$#cglwsIOVyo*bZD!<7D&a
zI@y9}o$Rm<Cp%)EWe>B@%^RNU%{x7Bn0?B?0&EUEqbzUUtlZhTw7`1DS-H8Fpd1b%
zS=@gZkPj3cr4w-t9BzB_a&xJY4bnf<kuU>ngM=JVHJr%7T0YiF+pZc;HHM6}s)kQS
z;n4Gt!OA+)K7X_-${)Z;0R`q^QWcW$3+i+(7|v}EK-&XI(g;^$c%Wb;NWba~wLOaD
z65FZ6?OdwE`WAv9<ypY!;ZXi!grC!LfL0CHOs5wCKLaw($l=arOU~pHG<FtehGL83
z*n>p|L<up#SZSNf&Mt;DYWw`r0aoUm0ZL=;<2VppN*zA8*uHeM8X7VVjaOOW7w2++
z&J#u`?7<kGFG7SYa*<e&#!q#OLd+9D%77DjLH8ybZgIz4a8A5=6Mo$&+)4RKqG@in
z@h_;xi(T@-IN2o+WXWI=Kg`Bzb*5d)nlH;;&fT3dz$#^xFI3HzvOi1(VGRQ-Lb<|e
z(57W=RW)2$Yq~Ve)|w#=GcnU`EoCcFcNXce*j{p~)&X9y1RcM4s;DwYtAaS8vUklL
zO}bo~ZNP(`m1R46>1cA9V6jD+M4qyzWzFpcz(vZlfS%Df7fIf>mgwv^9B75ke$Rm-
zI{SMETA{N?94M@_l`uRd+ghr#|8$@gIy=OpxN@uy>?QYUYJ{MutFz9Q)LNlxoP69Q
zkZslLtZv3@VUwl_-3ql?@3g5KVAhj`D@4n+t;=*)OGUk-eKv7fUrij{@b8;g+(#4p
zuG8(gd$G8$CiY$T@0(cGM-w;o*~GHGnz-qI(ZupTnmF92p_cd6#NmHyD34WgpKe3D
zC<f#1>ojph>;UlC22Fe*AVyMW>8W2zwiTHn*=7g2>~6^(kTXj=CWIaoZdZ?tc3sHR
z*>gO)iOY1hNFbJ{iKz}>z0PX?#B^0MPCjZYxGnN_R^46Knd;hBSSzc=LycW0jp~=B
zEp1V1K<DqsHgIBtHqti*bTy#6h6YY-Fgry8&HnscfFp7y9-Y?_#kN-H4rmhtI!mml
zo9s)vwv}r*@4imUnj!*tI89yq7y#2g=(aYG%QJR6^0i3F-+i4XW=f-|th25`*?rQW
z=oYk(LfEklx(B1PJvg>O+e5n?d|qd@iMX;D+BxYFctsFsrn7p0QE#KuE!_@-n$Rtv
zvLjeB@5?zmnpbK4v!i*HSt~4ETk6<riO!CfOHs9&^jGn7PQ`ki{kj7!(OI{w`jedV
z7joVDik!(blKqPV)tMTMm!g>tOIT;i9B7HoW(6(vNykRY89u19qa0|d*<6LltZ`(A
zbyn+dTXkA?R?ta$(=%3-ZpjK#K=$~VSTHxsqO;k4Ty=6=Ck^ir5s<BQk~*5@O=0y|
zQ%@S8E3xmPl$0|D=#XUpD+`^AiyunGa_+-PBVyk!$IfA$ebRxJz9QL#BYUZC`TJQ&
zXJ2q2M=81tql7Jw%4CmPgH6%~g)^Ps;;0bT*{oJ|TB|p+Tu#<eq0TGWwk)J|-<yTB
z?%Il^h0X^Y<?FOX7iOVsPgHIqPRd*OS7oPIo4^4iQzzQ%jvW!YMU-oobAR<MxQWO5
zEtVwP{%77$9QoBaq9=ExQ<LL~#@3Fe*yQ%EOiMB`rFi;eeu=>1!q#}h<W-H0lRIZj
zi%gq3sWqPH=$y1F(J}cVT~G08At&)Q#53t$EWL<qqrKsi#H72@cwWD%CRnXrqbP19
z>m-`^O@~ZZdn{et0_U+Jk?czqk!X8dq3KLhQ>+<z*j+QjDuCljpfD)<sk0r4I0_Wp
z-bHUnAP*bU!~Mu%VAVSj3T=zQYE5pWDVA<bxq(D`R|M}&v?UX$0v#LOh3#FjR4SPw
z_0a}iT}XGdXWN7#&<tfzQN699_|#b8jRQsEdSBuX?jB6HfA8(OxZ4vSc`&xze+L~d
zt@_U9%NJh0aI$~l<CZ7KbFKe0|5(rY!zk~~yy?yN@A;Ged8XRl^q%ft=1uue+3M}Q
z*PDB<cg+*tBL5!m;EZ=|CgUCKFZNtC!oR#~&MUh%Ua@oboHqX}v**P9i|)D0o7?PN
zzpT}}et|d7f7~;DMCy_$kbh$Eaa_G`>sT<?d!=W*?HzsGJ7j@3|B1=|F$;6N-ggG=
zdGb#G*vsck_72@~(7znKpLHDAT-7lNg~tcI>@7I%U$<q}+8=#k-I}R)tU)8Iyt$je
z^>W8DaNIC>6LiW<e9gPu|AhAv|GnPT{wJEek>lP;p0jOl%MQ=k@c2Y*U2{jG(aJ<u
zg}Z|E#!W51*Ai)t$6A}xRzwYnm02HYN@w^rj!3!%Lr+gZ*0gt3;YFJW9cp5onRz5B
z5^0Y%U1&AL6VX%`%Wh1@QjK^8#)?$cMp~1Nt0Nhz8tmw-h-RXZ^44S`##PgqR2)uX
zg|!m%W0_zkld1!rSY_3Omq)5<>aFHvsx6waaL+%1Mq8tBloFN=og#QM1@%#}E|#f?
zHCuRYB?=Kq)<n7(sYoha(~Z&AXlivFbx8Q)U{zfuYP2?6&8^W@X)DtnX^KT#p|i!#
zYb6>Y@kD!Tv@sUxOs1?@YdXF@h8LKkP02)S7YkWfSsAWfTv550YDd$NwrG31*1j%O
zySSdY>X(Nr>!2u`JQ7T;s)J3H{a~bZsdy&FT}e{|e850NQ>)rj(52T%Q8S4sUa!Fr
zwIK>+V(E1;+&-ob&8@Aiwi=ovsaR_)nvTIfctOXCsLm0?lIq%;`Shn{M(P@)3D!g2
z$_jU($C>s_3WGMUs=AVz(>AplY3QV%BAlCDUSDa&6Y-45Koz1WJm4R9A(CuuiX@wx
z>9HbKq^7F8wt}T+Qc-$4D(3bPxTl5ZYk3fP48W1G*bxKtzD_(5;ehAv21a%vMwUC%
z98G81J2Gw2^lI%9N>)?6jr3m^i>_`(bh9FQNGd2;S?M;6ekR(`8nf^)6gd)CQ*V|n
z;s}?vczDRCVbIH#RFyBRqfwq;vm{VnS{k^hcxrJ;U`k2pw32BhGXs;_G98K7oK>+z
zEER7A8BH~|%t2V6G;L~N(klK}J0&owd0i^nz82?1auUsp#>}K>Ydo5UD*;SrqK&I3
z!IJdqRWd~!NiGOimG@hWsp$+Q+&Po#pt{*vOfStPqOG;$$u677C^gcmVZjCAWLA1j
z1Rlz!6T4WE_9Rc3#%N<p%vh0J(<a)~M75*obR^vxZ;Yvl1vg;B!uc98RGB0j^YX>P
zurs2VBRrI3<zUTnHN_euSZ1P42qT!WNdyH~OYO;ZaCUe<G~@{-9!f)sNkmRbZwWSF
zVo)_$0^YnTz=w`RI=(7_ciPCgXkmokDnFe|z8&tW>}-rsXo};FzE~wKR01wmFS{wI
zv8`RXS%X@A+m&U-WJ#VQJRO?jcn`8Anr>lsQH`P-so9R`m$95ERl~6(f_Doe<+Y2#
z!TPGXYFes(k*jB-t7BHOBO{Z7;%);RrlTDv^p*v3OWND#rIKxkvTR^IEnyQ3X*?5O
z3;R%j=@ms-Z5oTZ;JnI6d9bd&3Q+{h2Bv9JUg>0U6daZ1xeEG2{Zz}uMb7BXgVtd+
zLAe%%Jxr-|OSC;^A%><_BP5V9FcVtnJ<Dh+=2TwFlVAxgI9%D~1n`PvD$*2Ri-m>P
z2vSE`fx`~1C)pSa?{20w4C{i5HDtLo!V$+}s!&s#?JJKD#*)0KITb^|0z>WG1(oIX
z<P51)qz=oFTD!EN%aiS0b?c(-tY<@4CYDS!p*lL2?r3Oi#kAn{uO2Q}t6nfAq2=Kv
zk-1g##hKHwRgHM}RRx@iV0~qM)gp?w5t;{x)Rr=LJkgva!O)PRY*qs{9xK2m;PBtr
z5+x5-%MsMZlHSHF5TMN+!6qz5lUGe#mzBl7y8e7ueYh@16*f+Jbyau~g?X~88dNJy
zZ`!gVHN9MkHizUX@B*5hjfl;t#KT8@)NRef%(*<CX~DXMfX^yZaFRv3F^PDRu{c0r
zF_QTJd5PGnD0{Xx?~<B@HMN)5L@H}4s)9Atq>3NT3%`+IdTBD=6h_-rr>=f6IlNY*
zvL=XGIIlW5Kf9{eBopgnsU+*<GFB?x5^v5}H09SpakU4rToH=rNkpO!T-&t}TJRFJ
z(8=D*!U7gUuu3M<tX3r!prwwB!V+Qu6$WXPB?^);)i#4e0(B`Ate)qJffRI9psV87
zp>>oirw(OkKL-xAF~C+#wqZI~%A>_ohvl!0Ba*mc#-hr|!m64IM^jo&crV4CkliP-
zbHWGPVr>l+0pf`@9r08bjZZ4phPSL~_UQehGa3|4>JbCA>BZ#Xu0?7yrdMuAf$PGp
zEm?TfU{I*V5~NBgaPa02^Qf7E1I0oPMcgqCx(d(k4Y3+DVxp%p(Gyl%djv6YRi=d|
zBQ}J%h(uQ`@^3PriiEHi?VG0%4?SMZ)KsF>byV7FN=ow$ktp`7E{9?g(i)P<R;tm<
zmd*xGZNBst5)pyIR%9)2Zgj|Z?aJpNpt1S%?58afrV#Iv(0EI%(`svF?Wj%IBHqXY
zh!)s7@`ihawr*;p)@xCPh+eT-IcRpVP>$G??4Z3Dw53j2S^w2*36BhgINF(TD1o^#
zNw9f~w_{&O!3WDI4-j3wtfvht+NM>W^j2G8dh;K<(+GC}lSamiV?#yl5?Wtf_J(-u
zQQZNFE;^iz2o!0Hc3Ko`aP>gDNQ7psh?P)5f#mUZ5pSOnyISFbblDBDyaS=x2m7j>
zJ@lmYg%?h!iD78u$*ICMe;XmfTf#_V%jyXChUPc#*Qpa~%V5#gps_XHE(4~#n43gP
z<aHFyRz<>5?7r(e+F^KX@Y5YwX-UYU%Ej|5d25S##9LpwzG54Nf{XRsQ>&J~eAM9)
zijp&_7&Q;u<Hf-mwS!Y3Dns8(Ahp_3NxePf5vP9h;6);On~f_HiD9W&yV{S#nA)xw
ze%@O;F3`*C%quC{BO!vOX*iq>BsRt|y=CBG$I?TWtroclg=IBDy##}gE(qD!wX{LA
zB57#JcTy6tO<inFM=a48(~*Z5=T+AR>$Tk?&IO>_I)r)Cr<Zysa?%Q_K3~E9g4&BX
zGSX16YT>RLjxV|)_f>LfuWRNhmeT`X{kf~^+O}0KQ4~wkkyQwJ6xy0F<K?aZ3J}QB
z(kDw_v8jd2Sy*V91}8a~Ttmo2U^6YLtHf3EJkEzIgJJP3wbn$`X%T-{gIOq!$yDSG
zh1^@9K2t=chHyngQ$|f`%q0x2HMu3(7ULUs)ya-{Ivh<$@iB_%<nZ$PP;Jem@=2IH
zlQExYVO2_|^_8BwWdOO2xOh%kyz!L%e0Eb%Pwt?-yn<(YgrX^SE%4JST=l9xOvOz?
zLbMZ;U5J~mrDa<uPK^@Wl!34D^#_j(&oOKPtE=kkRVb1HQm!&j-4DSaMH|w|R$P!c
z8qw~Eqj)nm_+%fgs@{jzaH(s@ggF(%jw#$Vmv0t?alau=p^=>naaL)`BTWtBYL2_o
ziY<RSmRcL@wKAPFlazb$ok+SC;*O~}Adx;2=e>q()NsYSDuMk?uj{((x=+=Y)>c(m
z2zv-1>;vj5oyL&k2eF9Tpp()$<swtat({n|4MSUpa>PisO9ZnAVY;=W4VSz$`w`nY
z8g)1~x3~)~=Omi;;}w<F^+7Ark?2Uru&h%!hUwImhu#~I!W6?qTNTS#pb8D_Nz<mn
zm#q2a<+B1sOB%54?@+fCCQa*Ls&qr?wBnMf#ib?Vje8Av{z&eW+v^uwc4vWo`jCPC
z9k|h24O|}mQ_u9y%`3RgGdHjBR`1-r!1kQEc}2Gj2<DaC;0xx>*gUWzZ?EnBn>{#h
zMtNRIFs}$a;6<IedHx_CLb6)$*D7uC1Fa~##S_dcyun+M7ucLrlXu9|n0MH-An%B$
z9_Juv&-2L42T+lB6o1D&R|2dCRtKyS*j!)>fOQSd3qVwmn_HMW1MM7^wgPl>7X>#^
zTbsR0^8%jh2j>-r@(O}^{+e9REg(<~VlJuKZ-Pm1Bl&MVeKC0H&U7EiOFtj<k-T(E
zx{u_g+t__1FWopd5}wV&9eModBsZ`1L&EgXt?iQU{+-0D-s9^f={)Aoa+SGCz3N>*
z7oxIW^{tdX@F6E*1b8hxW7bFVT6j#akL1ORTzw_4diedMBCq=9NFO*}2<RhuExa6H
zq^V4wY-fIi7WOww)B~BmlX_Mza?Ua2hnODdJUXR|7YVLueO@fM`OEh*!2^cNeNOOQ
ze9j)4IP0^Dvk>RhgYW~bV+NiSe$|iOg<dCGWbq#&dy32d^>YUZNPUjt&HdnC=?DLM
zKlo$7sofHpN<7bjJ=G8Yp?>f`_k(}h4?YBT@5?*K^n*|C2M_jx)3@#Vk{|5{PxONu
zw^y-$!3WfM#nQ_dT@x&3SZ@T2tIJC*94@58R60zf!wfoHM2FH6K5}^}vrcDV3LR%!
z5&TXtH&`91oj0$pvOZEz_d6oIpXhbPAa_@>U8h^r+3OM4&jEdHjQZO4>C5`a{9Zo;
zowWI)d$GN4yPkaWjaXC@`Ys0Do>P15ti0?Uo?aU(5GUX5g$G&I6%wE6MA7J9_-Y=d
z&%1DPVU6D>IAP;l_=hh1To*o2mKZ8_%ekl@d_g~Wx8Nk>Y?u5U{qX;)AN*VrfJ7%+
zTXerJ5}fp;g;wJiyKsui8n1Q9DRSXUUHqrI@GD*X<6Zb_7k{A(U+dzh7glt;I|L`}
z0vG;O7fw$hX#V?M^6^BigYR+4r+0cZ|Klz>6q7anl#73|3;(%GKCSJV|I>ct43RAa
z6_>d9M+t8F$;^K6h5g{y_JePC;q>yi*7GZZlRe#bdq8l~v(&}E$Hniq&%S=}p9oI<
za@*&5mwdPVf92w*{jJvXw=Vf^J>PW6aqIc1OU{`tIX>AaP=BYm@S%c}o@cr6k%E&Q
z+<KnjlH;~RkxR}Pmz+`;zgzzqF8+Xvzf5qm-T5v#?sk7gj7TQGLTS&5_bcR&FukKB
zocgTysWb;EX@0eDQFRF)jK86jku*fkgyW@A>k5fG{9YqyjV}}Z^o~#CM`SUjxk*Xu
zqbR&vAQ|G<G6|x-Q_^ziv#gZ(^(dS*K8-fwNORT6!9Nwe%D|hYLtz7dBToUB8~B$*
zDQzd+?vrFxq*iru+PzbZx5mKhM9wt^{?*|M=r-_4qW=y9|EB2q4Fgvn6Tw^52L7&$
z7wy+5Y5gaO{<L4Cr14ip{}&YO;2+Y!A-$|l4!%bm?F|E`&)HCVPn{h8cV!$uG4PKE
zDZnS~Xnhuo|D0mrw+nrifwzb~Cm8raX>Xc=UqKBal^gh<B))_Vd<<UlqO`)mACP`s
zZQwtnNrSY`z-xtni-CV6{5uVt-V392uYsQ>{r$Fq_sG@sJ_El{iLstD@bksrUNZ3E
zqUY-dzFFkIW8gOn{;`2CmUey84?Vtj3;%Ef?~-=)b%o}iDE_AXSK|?}!)!xNo3uOE
zz}JbK8Uvp#<E!_LTK*McpJs!9fyBL^7`T3b_ay@#MgoyuHSnS0Z~9?e-R>OWA1w1&
z<71`%NCR(^ICHjvzaV<*=a{veZt=HbgZ~|gld}x`MH$z5241APWqr=T1A<>|;Ezgt
zSZm<ti9VYR{7Vuab{e>rbFYDGIo~$$TxoZofj=Ytdd|Qr#SSkS_%Fq8UpH_akKQrx
zugUm+Y~Yv5xah~L^|;Vi94HN!a*dA_zYQ4p2J!z%2Hqxm&NA?65?ALN_*G(udIMi0
z^D%1Ro5h|913xJ8KX2eEvBOOUuHW0Z%fKtd&%bWqpUSvAV&FTazdttc9U|wTfsd8>
z^0I+HEdBkXffovW%)mbpT3;9F@w!X;J6NAFrZlwoiJqq#_ybfC32$<8(sGVSTq-s2
znIh*>1E;U;Qd(r-Pm102L-tyJr`REG@UN6{TxZ}hvFBC;KPdDa27bH9*YiNP`+(^A
zh`~=^4x;p=f&X3jUo`OTvM+ko!1Z<PUkv<4vFAqyo+I<qE90;Aq|aMX8fM^^i#^XX
z@aL!^lD;m`a;D2Xz1-kmuGp<>4E%dyhX)OOugL$Ifqz@Z<<AD*B=UzzzqLLSM4xE}
zuD-l(St||veX;Y+2EI$;#^VP5SMksH4ZK108!hon>(ec9UVVHuyDz0L-%+YH_#YPg
zcNzF^q`hw%_)TKB?;E&&k?Cmz?-6_chk<`jsbjrr;P*+t{$$|Wh5v5`K1k-prv^So
z{BV%iL)#}ReCHYXV=^wK23{<F8#VBWqNn=ud)Cfv(%**+{uZ&@&kej%{IkcvFP3(5
zrJmORQ^8L)@Mf8Z=NLG>K~HJ2f$tM~wt;^p^n3$Xp8>*m<qZ5Iv1i1<PZ#^FHt-i^
zyw(}`&&2+l4g4{&=WPc5J(2Sj1OKDY4;lEQVu!~JJS6ktCk8%H?D?XBcS-zw#lQn1
zN4HPKG*lnS{CeNur>{d&`lo?EAod(6`f54TM9y#nZ<cj)jDe3AJufiutHpn&8~AK#
zx7@&Y3xC+a4@tW#3_Mr-`Dz1i5&vFi;ExNv#lY!z2THda_(Gv~8F-i2=V1dMCU*X@
zf&WJQ^92L1ll~qy@bjc!e=+c1iu`{W_;9iRAkj;Y<9~@AMj80q((XA1K3(iwV&GF`
zyyh7AoidLW7<h-s(d(_&=X&99GWdTj^4A!+PxQRbz?;PmcN%!3*yU>m{*uso{n7gT
zQuO(;!T&BXApO$7yM_KQ1OJ2c_iY1zSk|=@2Clwphfl4@xM_X1N}M^x!2c-ziEoW?
z()@*DpNR&(OZru1;3s6AYBq4aulk~aPZ0hm4Sa;?dC<T=koftbf&W3NjZeSy*8g)7
z_l6kwlR_67_$|Ue)4(T)K1~LGgY1{u4g7T(k82G4Jn^5K3|zkjs{LKte}m|w?X7X$
z-4_ix{BQ>f2Z{e`{)>cfq=A1`>^RoI^*q$~9krZ85;tZV`~yXPwSn&!`f>yRjdY;V
zz|WUBc87tF5%~`q_+8@v-!t$pi~S!r@KWjs(z6Dx?em`o9u)Z>8u$qr#}cuJwvR3G
zP+w1L`~mTwiw*v-i~Mo}UoY}&4E%uLR~Y!0CB8Kqct+;+S_7wdXenJ|;Frn#+GOAp
z#Sd>XaD9LCRs-KDcD~cV_o!yAFB|wjW!!fe_&p-$K?DCx?6$|ie=qiV)WG$}8}=D^
zr;O`U27W~1-m?b&fY2`(`0Wx0^!(7{ex}T?mks`)@V{o@i=<y~82ClvS4Rzei;UNM
z2L7bz|DlPC{7(%0GSNr>=%V%ai~W7#zZ!o-`ZdJBYsH?Y7<f?TZ=r!}zZ_%W`ucy0
zfmh3VF~h*Mzg}YC=S#bn8Tg}m|6t%7#BM7LT*voD1FzBZ%fKg!ov$%)eZ92Fz%LR1
zzsbN;V*gtW{6?AYUo-GMV*iH>{7cf`?-=+X(O+MG=<&T;<p0Fre?a_DU-xMKd&Cb9
z8~h=$!=DWN65&5?;0wh5e6fL8q2=dEJU>&)HQp)lq~5?6iT)V_A1e5l4E&hr@e2cQ
z5&b_j@OrVs2$?rppT7!!z`*tW%CLdIC;D$S@F&F|_8Pc;IsS-&7s$ML*TD7t@J|d}
z-wz)s51eZK^?k+>2L3IH|6>jO=i+b027ayBVYY!6iob;ne80q_B?f+p@HZIvy)qt2
z16Lo4#=2(UlSO{FfnO-=#n%jcfyB3m4g4oEkA7s}-x2<U2Hqof`@q0o5q~R`AIP=+
zH;FyR8TbaVLydtyF8<tU;D=;>Z87jc;s=iy_(t)!-x~NYNHL_d<q=J-=Z{3s$p(H{
z#^quIe^BOcm4Po8`!6-{4$-sGz|R&r?FL>a_P^G^za{c-Ht<p6pI<caqcSe{8#w*f
zDwMu!;5W+nK4swRB|iMzz&FWw{no%=kbWIC@NLrW-wpgTrH++D4_zVYaq&xkPciV@
z#BT!z{;tR`Ht@fQ|6FX~FG+tF7`T4!)0PKHv_278e}@`4{SP*jPBU=ze|q8zlm@Qv
zqXrGUQpT&!z{SkgN(0yNu+6}AJX~+!e-=S^82HmNzPk*3yXuzp9Rt_>dcwfR%e;Qx
zz;&GbwSi9%KY7Ez>3_(f^qzr#U*xN=k7xblanXN>Jcy&m%PW2|%)sX;@m9#d)z{E1
zD{A0c|Ac{`CVtf^IH@NPw{V-m|ANfZFS&5(D)Zm#!YTA5B+l=4;l%$-=C-WwxNzd9
z|3rz>4_!F%PZhg8XW(BH{8tA4iQs=Q@D<XzbLe3>B(gKf-!FKnf$QIuFE#KBWn2~+
z_-?@$yKw5)JZcE(Di==e`egn_T{!Xoh?tR9yKv&K6aID=PW)k!zutut{{i8@&V>{I
zA7s7gcHzW7LDrpHT{!W-FZ_4AaN@sN_`fQ+_Mcp_|I;phlB0iDd(q(6&%eE4@DCR`
z@3?SkcfYjzHy2Lw7t4J7rwb?kkA;5%J$#9z=T{FEBF!}L=fxjF2L3y-=VBL5A}VEG
zT<*e2AMO897f$@27yf1!PW;;c+g&*E>wRp83nzX(ALYLwa{T#1iIZCl{<j6c)rFJ%
zr8=&<aFXwn`S>Nl_4w*f+C6CCM`Yjov`Y@jzfa`9VBkL${E#8%8)E-=T(}G&q~RH9
zdHG7~b3*tF4g4I@XRd(<1;5h3_4oLa2EJYRw;TAgg8$IKwH;nGaQ*GWR|VI8xPy$1
z^pSyoU*?@n578oN{s|Ii@(f%*7d1j~Z6CcKm~7zs3vkn2awKGPdvp827rJm7U%f8Z
z8~9{d*H#!f{nwV18U)vR>M!!GbKx}Z;)d4sE}VvZz1ZO%7f$@=;E&RSE}Z!H3jdE>
zIH^)Z%t+6<aN_@~@V{!{`gxLf4ZJ|&!+^6`obx>Gdcj8<_%{VV&%pKb6SG~on2ZHi
zx^U9-sE$hpK3XQtDi=;ru~N@k<HAW!h4623;RKZn|E(^Z__hCh!-W&Sh_Uv%aM=p8
zq!(Q{$@!_+`LKcO?<~IM!rgZGP;fmj=)X*)WYfdnbjF{EX8JV2HUBi3zvsF5-Fi+o
z_*cq&oMqsz%A~0@@M*Fi4x6})*Odl7h7?0;HgNsE!5Rb4qsfQV(GPx&fqx*z*w_z#
zlY!HJ5l-ore(;?Jew*0gu72=)4E!!7-rCg<zT3c;$-23xAN(-`A1luC!+!9m44nQS
zJW4<92R~@wR|@}s^n?G}!1eDahx@_xydgJ`tqkj<W3rBu4C3$?0+38Rz^K~G)BMt@
zT{48xpWb)9X5f0gu!_@NZTP|o&Y6@tx9A+cW)&^Y&?mNv`SU=<4e7MSKygQ+HNHBA
z_tT40Nh+h-rNww9&Q0RQITpbT$Q9!g5z$r+;#+&G@XZ8#fvLEO^U;QeRBWxL=mTU@
zU6rPyiB<R>hcIhOHNYJBz-UYt(~JDY^x;f;Grl;1kH8h<YjgMr9NxGSv9T6>GK4<b
z^?#BAkV0dm$5#t2lZ^R=RXyICkG4THb;~p$=kOAhIU@6#E_bLM*;dhi*~RHR90`|j
zr$8Kcm<mY=@w(GQ5G1p|`no+`$D~V9&{zNex9yin{VBQuIqUkJ2Cntj=h}Yi`z4;N
z)r=xX>#yNPsc-gwqcE09J$(x3jGRA*b6?|Ea)n|D2_KCgiPY`u@uU7}U4)fS`-Sm!
zi6^?u<V~ZXul75o{iBmrO=_QFkJ<hfV5Ur}abFP~Eo*~e7?byQ@R%6YJChs>>3UiB
z6m>H1*HO?{{d=dW8YQ}-oK4=Z0`IH;D{t;ypZZTrs@ebV0yAY&&0k3SWkk4=$@?tI
zOpMs`rmIT&xmHD;%=;4*^rim}(O>%s=}$|NssABh)V^t_Q)U4~dPwf?Dk|H)*7<iR
z=&SzYi&PCUgRbuMRo8!0>N^aQd0*<AxUQ%1kC7pa)+~2ACiM%93cCNA_d{g*vfpm8
zUq~3qespVB*VlF>n;I_8HI7MruhG6i_NhAypumK5fRsa`+ZJcz&upKzF@3dvpbiB{
Kn~eqy()w>PlrI?o

diff --git a/fftma_module/ref/exmodule.c b/fftma_module/ref/exmodule.c
deleted file mode 100644
index 94587bc..0000000
--- a/fftma_module/ref/exmodule.c
+++ /dev/null
@@ -1,96 +0,0 @@
-#include <Python.h>
-#include <numpy/arrayobject.h>
-
-#define NBDIM 3
-
-/*  wrapped refine function */
-static PyObject* refine(PyObject* self, PyObject* args)
-{
-	PyArrayObject* in_array;
-	int rx, ry, rz;
-	PyObject* out_array;
-	NpyIter* in_iter;
-	NpyIter* out_iter;
-	NpyIter_IterNextFunc* in_iternext;
-	NpyIter_IterNextFunc* out_iternext;
-
-	/*  parse single numpy array argument */
-	if (!PyArg_ParseTuple(args, "O!iii", &PyArray_Type, &in_array, &rz, &ry, &rx))
-		return NULL;
-	if (PyArray_NDIM(in_array) != NBDIM)
-		return NULL;
-
-	PyArray_Descr out_descr;
-	out_descr=*PyArray_DESCR(in_array);
-
-
-
-	/*  construct the output array, like the input array */
-	out_array = PyArray_NewFromDescr(&PyArray_Type, out_descr, NBDIM, dims, NULL, data, NPY_ARRAY_OWNDATA, NULL);
-	if (out_array == NULL)
-		return NULL;
-
-	/*  create the iterators */
-	in_iter = NpyIter_New(in_array, NPY_ITER_READONLY, NPY_KEEPORDER, NPY_NO_CASTING, NULL);
-	if (in_iter == NULL)
-		goto fail;
-
-	out_iter = NpyIter_New((PyArrayObject*)out_array, NPY_ITER_READWRITE, NPY_KEEPORDER, NPY_NO_CASTING, NULL);
-	if (out_iter == NULL) {
-		NpyIter_Deallocate(in_iter);
-		goto fail;
-	}
-
-	in_iternext = NpyIter_GetIterNext(in_iter, NULL);
-	out_iternext = NpyIter_GetIterNext(out_iter, NULL);
-	if (in_iternext == NULL || out_iternext == NULL) {
-		NpyIter_Deallocate(in_iter);
-		NpyIter_Deallocate(out_iter);
-		goto fail;
-	}
-
-
-
-    double** in_dataptr = (double**) NpyIter_GetDataPtrArray(in_iter);
-    double** out_dataptr = (double**) NpyIter_GetDataPtrArray(out_iter);
-
-    /*  iterate over the arrays */
-    do {
-        **out_dataptr = cos(**in_dataptr);
-    } while(in_iternext(in_iter) && out_iternext(out_iter));
-
-    /*  clean up and return the result */
-    NpyIter_Deallocate(in_iter);
-    NpyIter_Deallocate(out_iter);
-    Py_INCREF(out_array);
-    return out_array;
-
-    /*  in case bad things happen */
-    fail:
-        Py_XDECREF(out_array);
-        return NULL;
-}
-
-/*  define functions in module */
-static PyMethodDef RefineMethods[] =
-{
-     {"refine", refine, METH_VARARGS, "Refines a numpy array"},
-     {NULL, NULL, 0, NULL}
-};
-
-static struct PyModuleDef crefineDef =
-{
-    PyModuleDef_HEAD_INIT,
-    "refine", "",
-    -1,
-    RefineMethods
-};
-
-
-import_array();
-PyMODINIT_FUNC
-PyInit_refine(void)
-{
-	return PyModule_Create(&crefineDef);
-}
-
diff --git a/fftma_module/ref/refine.c b/fftma_module/ref/refine.c
deleted file mode 100644
index 123ee5f..0000000
--- a/fftma_module/ref/refine.c
+++ /dev/null
@@ -1,96 +0,0 @@
-#include <stdio.h>
-#include <stdlib.h>
-#include <stdint.h>
-#include <time.h>
-
-
-int refine(double* array, long nx, long ny, long nz, int rx, int ry, int rz, double* r)//////, FILE* pf)
-{
-	long x, y, z, curr=0, cura=0, curraux=0, currnext=0, curraux2=0, currnext2=0;
-	int i, j, k;
-	///////double* r=(double*)malloc(sizeof(double)*nz*nx*ny*rz*rx*ry);
-	///////if(r==(double*)NULL) puts("not enough memory output");
-	/////////////////npywheader(rx*nx, ry*ny, rz*nz, pf);
-
-
-
-	clock_t t1, t2;
-	t1=clock();
-
-
-
-	for(z=0;z<nz;z++)
-	{
-		curraux2=curr;
-		currnext2=curr+rz*ny*ry*nx*rx;
-		////////////////curr=0;
-		for(y=0;y<ny;y++)
-		{
-			curraux=curr;
-			currnext=curr+ry*nx*rx;
-			for(x=0;x<nx;x++)
-			{
-				for(i=0;i<rx;i++)
-					r[curr++]=array[cura];
-				cura++;
-			}
-			for(curr;curr<currnext;curr++)
-				r[curr]=r[curraux++];
-		}
-		for(curr;curr<currnext2;curr++)
-			r[curr]=r[curraux2++];
-			//////////////////fwrite(r, sizeof(double), nx*ny*rx*ry, pf);
-
-	}
-	t2=clock();
-	printf("%lf seconds.\n", ((double)(t2-t1))/CLOCKS_PER_SEC);
-
-
-
-
-	free(r);
-	return 0;
-}
-
-
-
-
-
-
-
-int main(int argc, char** argv)
-{
-	long i, j, k, nx, ny, nz, rx, ry, rz;
-	double* a;
-	double* r;
-	if(argc!=7) return -1;
-	nx=atol(argv[1]);
-	ny=atol(argv[2]);
-	nz=atol(argv[3]);
-	rx=atol(argv[4]);
-	ry=atol(argv[5]);
-	rz=atol(argv[6]);
-	a=(double*)malloc(sizeof(double)*nx*ny*nz);
-	if(a==(double*)NULL) return -1;
-	for(i=0;i<nx*ny*nz;i++)
-		a[i]=(double)i;
-	r=(double*)malloc(sizeof(double)*nx*rx*ny*ry*nz*rz);
-	if(r==(double*)NULL) return -1;
-	refine(a, nx, ny, nz, rx, ry, rz, r);
-/*
-	for(i=0;i<nx*rx*ny*ry*nz*rz;i++)
-		printf("%lf  ", r[i]);
-	puts("\n");
-	for(k=0;k<rz*nz;k++){
-		for(j=0;j<ry*ny;j++){
-			for(i=0;i<rx*nx;i++)
-				printf("%lf  ", r[k*ny*ry*nx*rx+j*nx*rx+i]);
-			puts("\n");}
-		puts("\n\n\n");}
-*/
-	return 0;
-}
-
-
-
-
diff --git a/fftma_module/ref/refine.so b/fftma_module/ref/refine.so
deleted file mode 100755
index 7e257c8f7953347c2bf30e0986c8f4eef5778081..0000000000000000000000000000000000000000
GIT binary patch
literal 0
HcmV?d00001

literal 31168
zcmeHwdwf*Yz3*Oo^4Ob!JS0FcT1E{T5JCb$0*YqxfC+>ULqJgJI7ueSNHP;<COoVn
z0-`C#g0_~nRynm=Z0(W0ETu23wxaFnp*^kQsXe8w+ECky*eY!;x!>>FYt2k1$!+ib
z>)y}D4cUABe(U#Izx7+cwf1vvuWMLh3SsEYWn66#sv75Ds#L|%URA)M%BVE5u`f0*
zlC+u0su8943LKO^8iq$As3Qx1$1jv1jviUi$2^yV_o$=hGl3@GpyC@;eEN7?HANpS
zM~snyBma)6`t)P5qSHsxa%*}m-p=ygi4q)foaLYmebju{f{%9ot6QETY*zZl@=F;V
zfK4A&b7b4-YPxG_<7(s3=}-Kkd%^SXy$}k&P`cke_(I!d`Q*#S-z5BL|KEMwBaE!t
zde(7c-DLbt!QWK;72)rV|9Uegwrs}ttLuLLhnF|5z2NGdB~!mK@K(+@e|z7ps2+23
z8Md+39$tL`K+kZLgdWFU`yjRYCu86TfIlylP?g^ugMJ(uJevGjW8mo@kCK^s_$GkS
z<P461-wpg|cJ_@yUpa<70nmH!m%g1p2K~A*@C{?&D<RmEE*x+({w^K`&low+m|*Oc
zQI#ZkNaDvEy}1s~qvG{$sdtKT(&I@Iqx2DpFC2mY4C>A^#v4^`hrC9~q3_QdC9XcD
z<eW~EZ}|P4J&~|K77s+@e!t-lg+p<}-+_bSuWwo5Zx2R;ouODf7;RZm(;W#1TLP`!
zL8or?8h_iC00jiPL$`o%*(QIpl2Fqfh{b|2qiL%*+Uai!L}S60zFtV{33PWy+F%RJ
zXxrd#>)PP&2!y%~8WRmdRuDY?dO#~8?S0+B3jY-qz|{BjMqssbaB6^E*?M!ZE$&|$
zjC<qp=&E=$6z&ANM&Vnw_EK8imNtKLUpOA>3D!lUktk8osDD*3p5P#fBRiF%3G(78
z^1QJ%kx+Y6TuWFIiS`6YKz3RXYVC_F9o?bU-gs9u7-%nzMM}#Iz-=sDVfcgXfp`EO
zXpP03&HyUhZY*u6udeaWEuC9BKT)o5x?Jv6w<}CyGDyI}Nd3pWW*XaMEI1^6=ZBDi
zT*d)~D+q-q7@&CqFgHEl=y@Vf#oGxruW0;nH6LzQ@ul$^K9GigOpS|!Y52X+&h2m-
z{+yyeJk2BPId`PC{b&k4bzV7^f=^Ag$5ZgCNmA~K6ucUu5`QuUpPJtWQ}C%n{d5Xm
zO^s5{*%Z8*`X#<lwU_fxYA<?H@OthfLrDrgb$wBmf~PFKRi@x|->|MK1+ThLviMT)
z<FybRn^N$qt7Prk6nst!{Y@$O)I9Io`3K9lCo?$NWBBeo8qXYj228%4&ssdDY(w*R
zfMjU;F6?usRiQ+<ixj6$4q=;q3t>U{xP+sG85*aKO892Nbm^%>5^g0-SDrc`;p++0
zg{SsPxS23R=+t%zFC$D>p6ZowHDS8&)J+mzM3|v<s!784P4XBWdxxeUkw@Ra;!l9U
zH;~E0YTt~rzW&#IJJ059nuO-B0S%jPuOB+ukvpx{FxEfAR(6!;PTK{MaLK^pb&`KR
z?E?Qf@;`mf<?BD~dwTF1pZJ0AC+FgY5aw0F#yi3s{&acU7u%3y^j+rLxp)c~AhQ3p
z_;}x*#k6~{4ow*JVYmJVnJ0k|Z;!O+`SW1(&TEW5_o)wp(y<%*A1`nkbBZE|j;}ve
zg<9`*-ZkFU-j>y?(BoTtcRt^+=c~oX(6HQlUqPoe^dGD5|3^drbH3ascNHH6zM;RN
z*lOsnEv{VA|E3Q>m#@F2xX{-hFV6ElS6kc!8beYaW&s$V!oah&#Z|ukv$FOOYk|1Q
zx98SkAF(b1J8vznl=!u|_kI)xtmyxTum9F!4+WJ_(2D+deS2Df=s!-5Dsotoqb&Dc
zt`5ro*w8<?qW@fd{~vw*-}mjQGwS<)*5LZAWBtG_qJH3l7pXMwOm3dx+f!RycqTXB
zHFSK%!uN9TJqs1f2U=c-u_xHml4j>o3nYeb;rII9Z0P?@{m$=;hK1e5R^PwD8tT~-
zFD~;fd@=Xlr@>x7Fdinxi#_!N887)}yy)vclZ+S0%Q)_H1Xi8(Ej-o-?pkoe6&~MC
zG*S>Q@hyahbMMSY7ojb5H7s5$#rhV$yRDZ#@Ik>l(00x@<4n5TryaR#Dc482O;T<X
z<#y%XdDpq2A&llN<wtK=@4eBx-upqX|CxsVzn$_xAv@G_D#8*yed<D4qJ0$h#+PUy
zdJKqD9jMBkmNFg?NxuGfSNR6!n;65sJN{6HE@}7e$(YXG^Y!mS1JOrE5sxz_dN6z%
z5IZv0H0AjQ&=*bRM;nN20FXC$;d?_vcOJcdUGA<J*nE3#1ofVn$G0a7L6lJ}jqvSR
zWBK}@W<caVSzC-U{kLaZ9zgUm_D}nsK8XQ&E8Lv><k`W)e@DyiJX+uXZ0@cLq}a?)
zW6ac*A1#0W2s6`nICuH_e>Jq1=NlkJAm#3QOE%y+SRMg&`8z{<x3ISV86SHm^Lh@u
zJAaY8>l>gO+IyCDa34m-pL`fUdzN{;Ex!K$9*P$i>c56P^?40@yhyb6U<@sVlma<;
zpY^~P-vGlLzFE=#>u)f}e}`Gpw`WV<%}A)oh&MycAnZj)dZdxh#G4@EC;2d>|IMMj
z9+<Fmc~u5N*|+CG>G*cvybEMc^v(3$@htm>{r^ldJ{_(kTz>S8X9k!4jlH}b+70B1
zf%%tF!@%N8p}YL~@}JU1@8Q3|y1{?^6|ERpJP|d6UU0#z++8c^<AM1FG!-MV|GB}}
zNw%}c5b$37M(MRpPDi|#3(*T}taunL`4xeI`D@9zr*iPP<hoqM7y0&F_*>Ks{+O6O
zS8`qgodvoZ)}MMDNF+U#pZC-uAl@8mDtlu;$Fpxw4d7n^hKz?H<BcOSuLC!DR7!LE
z`ZMp%qvC5}{!e}V6XVNdGltOop}fVIN-&f%H-MzCerWH0Ohr4-+Q6&>CU^IKQtjOf
zESvs4Ks(QhzAl=0!}@3Trs~hU3><GfBJ~d;`saTfKInh?h|=)u92l{uvi#`a0%_QV
zr-3+S$-#KxtE?QDKS_G+%E7Ol9U40IM-J+N`4>rdWqu~XewLNbqtdbV2-gu9F9Y*8
zgZqt-!Ssh5X}6(*j_w2l+_Rnf_B%o!B<bTIUG9)}D$=VQ($ysGkx2#X2DDW#6INc-
z-Q8`dbA)hTPw&>!wn#Ks+SAt^4+Wypz*fWWZEp6i^Lv}>Q;ME&B<|^m^kIpetoMXs
zp2pP;4W8Zz*7(7w(JaqBVxE9!Uadz;pDksU8l$T6@T4W!?ukb{A$e+IqzR|e#F>aE
zB2QG%P;5Q&yvEZO>FJds0-b?SI2QMKtLr_Rg3(wg684mBnZCt?lL=3}3x;ZJq{Fd{
z#*^Zo$ZPt)lefg%;%!*u8DVD3(-+<lj%*HlaJCW(gdOdMF*bd7C?4+)D)L#fJ?+6b
z&e-6%_E1L$&NstxRUHn-;QY8J>YS)SxF~X6$!keS{>JIQAdGN4^zP75FOCKh!ciPf
zBL6b=a4T>QhPb8K5L@y@(fDj@FY1bc=e3ptkin&bk(0M1r*K)Wd$YCOxaRyTXU;3W
znCJ%l)#9%bI^;Qabxz(#%<7!NdtKEzo`H<&oRWJoy*XvOv%EQ#JF{zZ4hh#=Vq8vT
zO-`9Nrvx;h1y6O3<>h%sH~v0Kc^Gmp(?-<YV|sH6ce`qHJUcTQbB>#BIVa3zIVa5)
z9KC?e=TKM*L~YI>{!W`W0NDWCD&XpXs|Icva9hXacp%Cv)h(B&H5E{i0X+xd8xIDq
zfy&*~x6`#I$79|xE~n6!ljqH`8r|kS08kB}I`xzYOz&^+d*FQ!yzhbcJ@CE<-uJ-!
z9{B&)1N!>4zFw_cqLzzBSQRwPYvtU=JBP&mI(_X^-{0ZAKyLbaIqwT{vz&u*zB=mb
z=3IMn)AfHjHxwa!TIti*bWd`ji%nnCuHqUFo4!{2l%m%QwnC>mc@5d)kr>5hoL2RE
zU5rzBxlK_=d1Vu2jn@Tjr@YdNdVMWj>%*y;wEws=6eoicR^}O}75spz$MRR!>vk#i
z|IZ!tee-m~9#9T=OdX$6#~-NUi|Y8gI=-WhS*qhEtK+4*BhwWqc)RMqF{%xHA4T6c
z(f-!=OUB~U_@y;9i##Q(TXBNb=b2wxQCc?piarUIZ!5o|w5+1EyljS3;h9@jenr_8
zWeerad?%osu+TRlo5twPUrcjO*bwew9xmgIFGp0rk&C+zDr{33%`+HV87C@Xq07o<
z{JVse3u2d<{Tk?98QHH%WLCaJW<CT{GBPrGSI%W+J_$#gO?ifqxibV2<}JAEkhzOR
zNL!Nwz%t*2hRi33%d!rzMr0L2miYtNmU%?7Re~g|?Axf!u0Zib_FIdf#`-15vd^>v
zZCK0D`0RH-$)fdj<cI9DEazFDgcq{^%5tG~pN;b0SoT=!$@_PfORU}Kr|kDwF0<-j
zVfH^*uC$)VKKmTYRrc@6YX|~9`!OVhY$1Ur`$yDZN?@(sOdvx<Vct#lp?m;YVheyS
z+r=77(Ar-6SF9N?fw;X6_1QUsC(T=Ii~7$KO8{)Qcf-Tk6GaPv9rirzvkSyFQ17&x
z;JfTXaW8<o?e!FJk@ye@_S)mfGeZmj*l*_(m??n+b~6QBE`fvgG}g?Lz#)4!Z7UVC
zL2%gqPh^|yInuVH_Ej{eTznK7o;AkxP{IS*=fl^9*$=Xj)*x9w&7xs#q)`vCxt4Va
zq9Xe<tw86YU%kL(S?oaf0qkh@_}wTF%dV8z&uZ*MU_TFF!T}KEv|{%iK)Gy^VZJp1
zOdTIVG0S2VqU`hN8^qHupfbz)ElBeu(?7`c%fwEy=O05<Svlf!6qSuS`vA$zn-Qb-
z*Qv&pkx_wmxU7s1FpM@}XN_CXOt$YLPR1>41*>6|U^i|N%NbBQ8`bBLCffr{p8Mt<
zq@Ms&-5Ux}G^{D;arcI4@RnsvoC|~R7VcU>*p_5<Lekliu1;21_MK4fc53qMuS2_g
ziJ+GX?Pt+_?xm8)W3NpDC3YdnmJ1(rm)V^NRCj|20;seL;UD)3@c@7-`|Gr)Q5*x{
zv!6wmyH|=K08RD{=yG?HaKi~}?LQ+Ux?99)05{oQj2ZW8@eY74yNP<%i1z^W+SAbv
z_w`~5fGzePX!Q-^L%?mfe@31grL?>4N}97?d<%N^+PrM+{-Agpz<zr=fq*bUeZanz
z0@_6gemiIvvT+@<*@x_nB)eJs5dJ%C&qANMd&HyAbJV_p7Hkxiu;7^e7F_O*OYa=F
z|AWomBHB=M!oGnPY!wjzC+$lq`c_c~rGxe?w)8g91mLuNHyyG~JPhD$!BN)SF1y4K
z1=qt~_jbV~X^DcXXzCr}Ly(sz3i@c#4sjQtLQ${*-0q!DABln`G<}!2A9W?7;KQ`z
zPVoYuGEwjm^6eH61F94SnbdHXxEzNnQ80&6?v}0eiGpic*Dtj<iGu$k^kE6D6$K@1
z=pG5(1U@$59tm}cg4d{hKtjEupn{!$uZ)o`qTofkVV_t5cW=kzRTg8{Y+$3i*Rxo_
zfVlHYG#CRJ`y95Nf4Op|$g<e`<~p#ucLz|+vO2MsREz;h<&c`+Wn=EzEs5^YMC`yM
zk)m8c{dXUdluwbe74pBYDGl=;HmLu&tbdX996qll>)UD0hfmA;cUVvUaqvlQVh2On
z0F&a~14XFIvM$12T0Wnotj?9=()<#M?yZtUOEi&x7?Gm<1P8{xn<Zs0DXF=CI3)!e
zCU1GV@B65v{G_HlJe-o+l}k3VpYDG~l0UD>f0H1W#!A6i)&Rt3S)8oQ$u#`af0bNA
z<f6nvWbkC`aw)_xFUF{HKSWCNDMX_C;Zjg#SrynzCSQVyCDc%oIh!Nkz&c6l*Q8q#
zq|*LD%Fn(H`<vXKodh+LTKhT-_oFn^n1n0uhWjyEZcXYAh7I@U$dqT@!1f+=%w0%V
zeV(55SSQ)_Utqb|`YQJBC)m&utAQ>3B5f;^4sicF%av9EN7s`qS6N(&yARO>pOpcO
z`^&Upt>LDc=|IbA3D4#Hbpk5C!giR`I9k4TB`OjB*vk&R9tFpPvP5f8O$p{nj_s#5
zNcxS~PrhRqtzkB@0WaMn3AiSa0=af|1UkWS7O?xomn6~GH4zh^LnK{iGF!P{=JGbn
zdPNicB|#)*O3Epf!u^Za;ZSQfTl4A~2B7sjmS1DpvS!ih*ICZ99^v5o9m|E*H7pOZ
z?6K}-`4r0~)|cqL-*bSMS@*L1Mj6VL*3T&C4=h(%-=gRL$g<BG&ye~P$9j{+<%Ii9
z;@4V_lIApJ-elcO`ERk@Wi4dypJBPz@{{*%mbX|xXO{aj%iFErF!uh!QGd5pPnvfq
zbFURA{#};$Tjj)`W%+<rLp%S<@<Hosl>axD4_Vhx=6jTR*m{J#|6uv3bp!o#j^$(4
z{~&tZLo6S+?<3w21WwqGkXgv-?xcMoJH(VVgSJ3NxLp!BZF6<)&XB;_f=d|AnG(R1
z@>4cDOG21Z&XRDPgfOKnqWNwKVM^HnjqdRh@`!?3nlV8_m{KaKI!~}gxZjfhTQcSg
z218+!d!pEecG&}LT!CO7HSEnaZ<1uU?4MKCWRZtf=Gl2{(i904+COK*3k4IR$Np~)
zi>ZQ%uf%R7+4+LGt;~LoHARB4UTH_^=nJLjD*Ijb@CT$cpOJqX?V2Xepp}J%?u*1X
zY1i-Rfs2)0*V4#KlwDj>xQhj|VxEy7C*5@MVNev7xi1sbpwzZVUn2el62tx~ZJHsN
z;4FI$Q0|%H5&(HNEpuNk8UYmAJ807^Q3jmHeu=#?TigVo#K^bliBj<#q~^cQMwN>N
zIAmFiv7gL!m9rAf{|R%+T+t#cxh$B>b(S2SM*dSA2J^&5S-BZ|Immfp=!{=kQd#56
zFg)K)X%*r#lHnj3n$SyJDme_yaYp_$+C5)zn$7<iN6{7HXOe>*BUNxcmE=&?{2ILO
z1%h)&{z8hV6ih$)E_!#N;M|_yOp6zZOMuF<W?-Ka0v8ulO2~XNUMXrN!7>uiLawVE
zeS~sllk(fRWVp(a`$2ZkVyFB?_TSZl3xNC^*i+XyHq}yhl~ZOr-Cpqpspd;s4OdFJ
zTrMYd8TrR3u3Ee(Nng>VT&y~zOzLP1^In$40gz=e+_EgjO*Xs0Gk!Xkf7U;M9)Bs<
z)`s;hkdMFYFwxd8upeK-M&((ppdUX|uD>C|Le12dq~_%)X8#;CB@^6SGFo`l*qAV0
z_SmG(VBDB6;m1Tzim{yYHp+RE+9Peogj~`TTDOosZvbTvB(Fhz{$hsirv%qxuDpz?
zFwcAeJvRP4cD+f{$Iqe(=AQshuwF(n%lZ}evQwzSaTJ9b=Iv<G1pAyM%0M^DNVyP&
zq&J<Q&axf_(fJuGz?JhvH7JTci2-92J%FMpS_Z33!xRvg^OOzS@>jB0F0KJhekCY#
zu9NHb3wUWtN--Bq$ABvOF<_%;E*vb1{)V16i#`a9tEdoNm{IfuDl>}?qA6KLmjIJp
z^cXZ*MeUGo7Y%}PT+v71Yj@EL5H`N(5R9Bq6hV~b6nz}N$u0T-%6Ub9frR{`AH$CG
zih|HJvFLN4FDME^+N7dV&`d7+H@I<15eIW&Q3qsBEqV~;^Dn#uGGSaVLUEGvfO#KN
zY|(UvQPBYEMbQq>ny6+-m|qE@`U88=oZE{cC>ljip;{DG!(bAvAklX}N*OQ1GNUL6
z-J<9WRG393+UY8q4qG!c>5slb(q~|lQ51z#QDjR=<^)ig$0;cPYKSn41^qDXHZ;*=
z&J;rcY>Sp%E*7AXhCK~6<}9%ofMs8aYO_?-&_zc6$KX@*GU28(CS#Y&9fIdOl+Fhh
z%o0%oR1Wc!rU|Ge%env<_W;ONp`cbv`KK8d0pSNq#$SwVBvPsJ5bImT?X0Kg_Oo8*
zixi??BgD*h@mbVO{*QzJ#-JlgmgJ5ilw?j=W~X?WOm8Qcq|_9uKhl3)Vh98|D<E*<
zL@kt*M*e-2(kott*YZEa=-enS2Z^+`K0zyWCqVdui2Of6fEg1vfNb)v;k42W$u+!W
z<ll^NFgFRls5AM=1k-RT)xvie0=J28f=b%+GRfuiA<3CRa*5;k$6?xAXv9tCT=N*N
zD$8pN<~&m-LV0ZgfV>D{pJ6i>nM?wCHXlkbuQCrnN}+u|!{urTOtq)eGc^)8-~JyA
z0kcj5MfMK~EH&xn3+%ZRP%i<GeK!Rxmq4*y&v?DotiZMD5}Oz1&Ff6|e3{Khlgt|>
zQ0eJF{FsYCRAqmXmbFQ|&)!cUWHJ_;>=Jsd$6N=S*BbdpIgD->KSMjEW6q!;6RHy~
zJn3*8lYa~h<_>XAQ%=E)Ln$#uKN|U;z_2rSNuvBHn{%g_13D>sMS?+gf1*o`{OcH5
z17e+~+?=3H5UGCr68gpbh<K7Q_BK7VPkaQ#lmC4T4#if%Hs3Ffko4sQX<DREmyy4k
z4*i&T4OG(VDR}S2X{IDknl4RC?wECqi?#^!VX+ps9_*je6OYNVVSkUN{<{dH;-Ed9
z^U#wb4&ab|lIp%DZUJ!EZbk>0Pl@dSj@r*LXuc(O12|?sOP=qEdjK4_%b9$RiTePY
zu%BVm{!{D+aMFI6ee`qj5P(5DmlnJtJ`3Qqy@LJx8}WGnXA34WoCd`q09ZLUQ1+k1
zQ-H8?t|#=S_z@tij{i*PU3tj^7gf%I-aIFMgSskQaG~i#Vi1r|6z~Cb%xcf0XK+!4
zD?HOQ|AabRR5`|8bD4hzbdxB!pQdD(e*@Gd3T*O?Gug6UVYs`9#h*-y1sIOZ23Z#K
zWR}H@m}N2VWm!yaSr&5;Rw8o0Qtp`uvMi3<EQ`4;%VJ1pS&TmS4<MO;Qo0E&gV_BF
zF#JnkF99~oViaasY*3c9Q<fQJSr!8*|4}w@yvgX8_Cq)+&+L`M(|&~gvC(Y6C^T%I
zO`B0U3N3Js2Qtetu+Kk3mxN6wn=Gpem}wXTMw7YS{51e4<=l~)a_&$m=OOxQr^$?!
z=fSgn=AAM}6xuw6Ht#YSpdOpCipx6$N*t3ORte)cd-UI=eU(lwdPEASvL9zyJ?f+k
zETt*oF|(P*KF*Q+WocQK^+#fwsqFhMH*(R`SGa^(O!8x{3xJTr=5iErxkK}Yz~Sad
z@dB1ko(60lfQw(|$iE=V)dS!XJcvl(uOKg$S@MFcRtPYrUjt`}1uW~^EBfUy@!(g_
zH;hZ87%O4{i`*{V%2mb!YX1qr7tv&~fXXY-kjuhlfI_I?!=sml%N6uD3C&f|dZ;Xk
zYD#k-08MD_vr3|x&fJ$*5(Uo@=r;Ei&S>icw2;z|Lc)v?BuMC|aMFwrL`Y~JI&nq_
zG9+{m&`e%`aKJ2#kC}Yy&;eHwJejCHM(}M6ZKrl7B5bB$d!?uk1H61W+RH6t-{7{a
zSz>j@zLVR4&$4J_mUSsine7Hb+{`(IvmiJ5>>iHd!4tTZd+9{D!yz~dqDy~*cyPeK
z0xYFxy4j_!eTQx<rE3&)7!cZ}prg0VDWR#|<#4`h-|^e#9AsDau!~*$j@>r%L9|AM
zX+|BoaHj0I2s`dW1ZBtdvg3YFQ1<dh_VPT0&CFfci74B&2Boq&Xoq-$>wa$KOTpz)
z0Q@UAaqvHJ7kh3CsvZW*+?!DN9ty_1-zq2c;=bNEG#_4by)9a1i5A;+fSDcJEYVDx
zXNk3J^enN?weRq4Mg>pqxP1t$^LZl2&Fjc8W8i9WRV^-_C9V~7AwiUZUz8d1?@*N2
zE6PozWW}ddML1Dm%=-_8<1o2L;#*8ITnh!p{0!C1r<@uw4;6Ga)m<gk!C3QFP-=&~
zkd$2!QrvgJe&dQ;&{)^OOpY}%;X<rkW4H(Tc?nnk2v)ehkjdVvbit>tCuG@JcqKda
z+r3JV2M{;$;M*TCj75y5`&E@hKRp2Sl`LL$3U#le1Frfr3Ud(SR~gXB>vi*Z4#;gU
z%-Rkldc<~pQFVne?})1Y6sjL}s-;KJU(Y0nPAj6%f#_LDgz73x3UZr2lM_C+6A6A}
z{vw6@gTleT^LX!-+m|6G1Gd8y>GvZFUzo&IoKm=NgLa-<Swj7_k_h?>wf>3{RV``P
zrO@_j+8^ttEZ&t;Eu-?cx(a%<dyhfS3yv|iOW#L#V%X!6O<))z0-Z2)jl&xV6{fQI
zzrcL;6a~;?$(XR#xMm8*h%0XLI=Af_R07@4z{gWTE%UlC=1f6$5nr=`!N^fWRhtnx
zSqFJjv+7mQXFV>9UMS1Goaf3_8OQ^9+ZkC^htd0ayJc}Iu>&l6`3@)l@{6i}fY9UL
zc|R33mjJ~-?<XOJe-rbnKZcg_&;63{;&wAmgdwHTdjxUDzuao?lTggRyl1^@ao}Hm
zpZD`Hl79tD43RIU=S+0Fa*A^D#d+CzSR9y>4Oh-0_mytWsjgyNPrMrS2nbb#<jh=P
z(K;STg}ctK9JiZH0wB&UwS-EL*$T|8pGa_=mJh>)sGmrV@fR8O6X&3E!lfu+W}PH1
zo9a~MW=c%nrD$KiVud7lWr3`iD0OE-cP1qjq$q<^Cjqpa#!&ZU$|c#9iNZ}DX>TC_
z8ZQM*oe1NvLi#yB1K9eBny9D*_zNK8f(+T%WX*-L2FgAl3+SY234chD9%Y1bz(q>8
zth=}r(y05gsUB&}C7F)JvW>+6WGx#!y;NK?)p0DH!~Uu&{CKr&&!vhG2Ak-^%ajlW
zmZX4~kY-dQoI>K609A(*yu5#hKv>*g4a!*`-hsOth5J~Ztt@pXl|SGpFHMmL!Z|7O
z0G5qY;`1@lmR~5!rR8hg>tu81W*X(v<av(Tav2X50E~Q)B9%LMcG?w5U)4`^ytP27
zmfosViY2fx)mx6Ql(~y&L$g>tMcaTVSdE5XRiR8-tW80laEc$Ro=RJ;Nh&bWvZ`c5
zubIk_@s=7agmWB!;thde2!u#Y5`lQ;{Uy`*Z}>%szvoQAUZ89cp3@hL&IyIvy8GIL
zb9%SNyCUJar3>cBcUKHqi5Jl4bhfq4*;080Ua6nm9SZktncW%gn=@bYIdodc2)fo#
zJT{DE7?4yt96drzY-<efjMO)J8=O~pO5517aJzilCXP4uW2IeelHZ?vvBw|i4LR6Y
zyuCfxfihC2CWcpl4n=@P_Cq^$u`i6T8MHg--mUzA0Lt(xzXX61(x>xYj)Uz9!fWlR
znD$_-Et-l5_ipv$y&Zf-0W7ewEw!?DYmi?Tp!v{&*JomVy-A%g1eT!=COU5nIdo26
z;jI@3MEW!JTHdul{L;1cs(v#x>G|L~>wfO8ssG~6>y}@)e2%sJ8$4OO%{t$@$ixRt
zP<F*%b>&)5yk<Qop`KS=Mb=tZ)H?4@*OrG|?uT3(kGM*#CtTy=u1)c{Yn)YT&M&am
z)h~YO!1n9+UAef&dg;o=A#25hA9J}oT(_+4cHOefm1CVT7ZgOVo(uVB$DKhk+tatK
z+I53DL%61%agAT*$~`j2nzlT{<@$sD#J4|cU3A^zIj#wNpSRY5_FUhwo%MaQQF+FG
z(Uo_`+PrJgrf)v5d1J->8=<7r<=z3R7yH(NV)wWmuqi(471uiJi0f+WA=d`$NW05_
z#x>i#Sh%|OninJD!@<oRec?7E9_Zw^CGghQO6MaG9id=%d(7}VU1G#<@wdm~>}kUv
z>q6J_rf_5L)_T0`<>#(3xFx=XlKlSOK>IwS6`!$)Zk4jzBEe`IUc<ttFIM{TwTcb?
zI9a`YTWay;2EY6^hGdQ5lMe_Id;vlkG5kvujZN$P^^Gk?hx45XoC1ZRv^#)6DU-5c
zll=5U(^fFzdla?74g;6*1Nf9gL|USn=&y}9w#M25-GS(a5YBHYyxCj7stz389Y(wt
z-|uj~tr2d;r#*VR18qV77JTg@*d627Iru4xNVt2eWL#cX*R-;^wyv3cftbH1(A&#z
zXHe)W-^%6|S!!9=RJRI^lWz9oLmsQ(N5^My%;sn)9+c+Cs6g)E{6Ms`H;UU2!+M8m
z!U6OJx}^u1CmxJ#4q^pOebp;hHW;lP{%Ei}7>EVo6TBEEYgRX`Y+TAeEzG~FEfAIl
z;2R-$3a~dGMb|*KcAC+~@0O@WrAVq?*HV}Ii4m~ygCiXO_%H4CN4nenk&X`DMxj}&
z8|!OU*49y0JR0CDc)`>lK>Qhao<R1JbT}yDhKw{1pHB>h{W7DanvSkpj;@mp=m^B(
zy?ybXKx~8Vaz|EssE0Z?2Ll_rk-iMS_J1w2krCquRO0f(EV%i@NRe#L>x(O73X4hm
zuU%bVvwRi15+6bG)RdQ}J=*D+-4pK%2N!n+!@+2%4PYSJ*0mV<diE6+p4px9|H`?Z
z*&Umsf!<9xh9k2X@NMzg_})t(hIj!|e#T`sJjnO@B2iLAWJObb&6w2~bBW%=_=!jR
z$eSFVd?OPdm08JfP4R{5oi<Jf2I35nWW+Z55r)!b${&WmHxeUTTcE8gm^Rg<jEz8h
zJGt?l8h@-i)D|>kx<~g*ha%A0&_(fxbmn!<-ljxPVj##aWR!Rt*E!>;&5!vbAiwOA
zqN6vm8L^DuhHY}ZsF+k=!GV!5K5Bpi%}3~VdgQJz91C@ZgY7bsIL#qQouHI~Jn9E$
zM)}SS#L79e*KworPecv`8S@<>yiwZ~h;_*(D{Gje=)PY0^&2*_Q6>UEUQqSdtX$FL
zZK<zrkORZXHw@)?U_;P|^u<-sGjX=!`#ycWIAGvNSyG9P<|WZc&#HJ-IudrG%R=$c
zCZ?&BEFg`f1*^PE>ijj{RW0>M2KeF-MqWg{3M?fsxjC+tEU=!fa#H*JM8_?Gp-s-1
zlU%W`06q|e<QLt5;mtN9y1V#dZXn7J2RV6fjSTG7oE@~aI=yyuM<P*wduS780h!Qg
zhtp!2LYVxM$q^CT!3hUR+{qjCQVgj@CQ<3=mB~iRZa|kYK3KQ1dRbjf3#>~2c9COc
zO{910s?C93X=3ZvcrX&h_ng?KSYK;f;%iL&B+^RfWn>DiX<F^Cu3xGmG#2b^!+X+B
z{;T!2)V0*FVA}LEtdVOC$Hk#=M}&f59FwK?Q;p*{X@G3~@XQ^IW7h4FwcyoG_G6tC
z<S=xGV9L-@Ue(f<rtl+o$owGvQqxf1w1SD6-gLUs87hs#f`&^s#suaIj<z-=R?y4-
zV#~CeB^VLch2mY9XppU?wahSTI&F&}Da09gD$XG)9PA9p$kqN_-MGAQ<#mnzy2jdi
zZzHuiSs?K?xIg(FDe|mpX=YSw3+fuZ7+p&myi1c)Z(}5UOE4Od=A^hN8tV#m#0`$(
z&9K>792hBiMh*tQN+1a1lnKm0SU3}uS`-<Wc!J0~k+3s=)L{ZylL$vRK{eh<W^G-8
zsLEyZqq7!}L!C|Xc^j6bq%+QkPI?;tjjb9sD9(izmPTq>q)X+pf8^+pVmF15AX4(d
ziaP)D`o`LXfimKfOBWe?h_jwxPb<fED7>*R6y3_+hz5J`!aWCvUbZEC&jiteOr{;J
z9Ne^Zg|j}=^DCqvJQ2EvnkF3Ab6tiB!>MITkTVK~Ne72L*FQ4J$o9z8>EM&gHB16+
z7{@V;<FL`w>qnaFjCXOYVU31!IW)y!Aan9`k{jH_<?CG1@SzV0RS<GcXJSTsM5%A}
z2e52R2_FeUT5BZIO%BeZa;ZpPCYD~#5=pP>FyNVx{F~$wB(YxVO)QR=M1zj|xPZYx
zk!u(z?Fw!&db*`m)CC`e+N3f_mR?aGIq-8~<t(Fy&4LgyEZaI^n9Tc5xV1<6_{ln0
zDHjFFDAm(~>=9-zt{-Ap&eF{o3Rqf&da-I_V!`|;JBB9@w{STL-JF(bcSm;s=0;lm
z(sqmr)kSn!P3_9noETGF49QqLrY0MH=?$uT0$U8`5G>lcszZ*|s>v=-K;$wOX{Z}P
z#`9Z9fdidYGz@{CX$M@~CWn_Z4aiSlq<>qFBXmhUETfgqWNvWMZtD*9s?KmCQ!8O4
zF&_>4avzq~u|CX;6tbeOd1;+oTw-j<MJG>Pq!WDJ=EQh{O&m&2LRo^fXtTG`SuQ#G
zB+f5i(1M<5M6bnU$FX&?5AdBpdKy_9D}I%5qQPEg(F;$7oy^Eylnc2;PjW4Sj25Fi
z5je=&C}5!()YCOU<|TQ4XHW;zB#C&=BptPq%O0FPLan^O5QYt_f*bpS;kKa81SDM2
zu+rO-XmjEm$(gC0ZwTs{HIcwLYdK#Qk-ZD<6*5!NjZIq{<xtl;QqLR9nJN=!k#mj2
z7-6W>dOcNjJ97hvb<E$1?8B;d3_@U#Ok%1XhCa7(riw%==bQ*i283F5s}Pm0S-q+b
zCyz^H*;nUnQjy?f_u(IDq=T4C)Yy;{w>r(H!<p5n2<HHdZfAHQ+R=T+oUTYuP+q%g
zi1dYGO@UYdA8`!KX<FCfTiH0fW;TL-4&t9vnPX(kz{pOXJwt9AR+CXfF6(-FWke;X
zVXlRcFv04C*NBf!CSs4C;5g-ZrdlhhSVRr3NH~#?wNyiWON)~Ml>eO|P#=m!9|T%s
zk#3y+Cc=da4vrC{0}En$L9YY!g04~4rYu4b{z3WG$!d9~+=P>DnUQ4tqWz9#vOlR1
zXF1aRZY*13!RV&ousLAFh;XhOBu*RTkamWw-??n1Rz5fd><nYUAjg#sy5!VKzBMcB
zYYlwi6N6lKvC}o|B00uMqM`z=Q>}6$U@pTokvp|VI6K!MWl4wlF%FO?yZd^=7zFyt
zjWjnU+z4|KTU*!A;x+tz;l5Z9lO+=(+~lm(^^8j!nOc;APKuENodGT*<%$Z#irz|8
zK^Uub@sR{-y)J$=f%5C(FB7Qsx_B~yYSP6zUfjhdj5WIW_XKLaE{<1s6*G*t#p;JO
zRG{xQ8f$bh@o0dt-jO?UzeAyf(V~l5jmyBzat&xL)Y+2IsTD}lk-{B}mMe^Fbx}(N
z-{7b!t{qLq;I@CRqInb*hi}*VQd`kHnu^1>|8o^<M^SP2s4CWursD2@LB+aJRGb*q
zQR_xiapIpk$~4Lz);e^EqBHj2u7Q)v4?v#Yrh&H=5<RJl`~)61Vfgv=W^Nq`)V2Fn
z@t8VFY5FaA-6>M9?in4rkf)26n2ua%)x{EplJYcAk)UhQMIE0Qu8xe6d+mi(kGzpI
z_rUGx<_e=p8>_{`i~}Qj^^;1M_9zwT^8HB#XSZo5@zEk~DGGGcpy2Gb^hQyj!;#$o
zNEJ`2BWZM5XB1(q*A37v2I>RKds^ja(zUNrx^w^SS{9$h<|ex{-PTb$Al*M`wsucS
zWE@P$*CHYR;O!b%sFZ?P7gIXr;D}Dq8gz_8*y(M$1yhqPIK55#Lx<d0(dRTexz|1L
zl0v~UUDO>6rP~wD(t0FP1G*EyoRsX8-c*QWZ{DcwPxj`G>0C{!X=@VxTCIySaC?fI
zGn+W&<X$s=Fu~ZOi=RoLR_kJ_td<eM`D5I+<Mxs|(lx61y98>LLzCRkP*n>PBu%<l
zl|Zf5#iXD$S*l@^)DbbLi<1+mHAzKIsaDBsOvrB1MQy(@R%zKuK_l$uyRO{sNeXfh
z$^ALyLB6ZQ4Gz)8WIJvgscS^{o>T&oTqDfUB(Ew`)2JBHL42E;TSOgU&v8}!pCoEJ
z?xJ!lRmV|+)UW)zHR0zbUHo<ewdN&N3@2o-(VD-OMCsxS2~@%;J`TeT%O$6f+-nPV
zC>=~^y1Xl4LX$2gty-n6dNoOvC?-r;<x<6-Bud+UD2dXxYcJ9YUG^l5U!^^|Jc&xS
z#3`*()@VadDp-Ep8YI0+sdh;^|38IY@S@J^sl+bach^(kaJ<W?Qt>++pJ8a@hT}7h
z2UPop<Fkz8qvG+b5+7Mg+gx}ZstZpTrEe}{<nL{e#)YrxY{zp{>6;60T;jo#5u1y@
zyD%z#<nKL@*=6Jy`mxAVVJ;&#@p~!5@p;D4iK!T;>@xBb-xo@)b>SBu^mm1Z<0l#?
zFZ86MWZ7jD;P+uh!%s2}Pa6%7=e9kg;inj*eV@jKXSPp{N{`1|M%i46uOkhwbje>=
z8IB|5l3(r|4UfkVN7-Bkp7I-Ia~ab{!BW0w1fC184CAe#^a?&aER2zUW=qXc;sj6B
zO8oY;_$gKR^zTyDO8H5D(hhD+-eQz~zXf>e)!(@yEw|nH!_!bB;^$G|M{CztfTx`F
z=bgVJ<rEtFiAFulUXb`H#sSsuRe-sT<+l^EE)TC)NV3H4b%`J8PZEq~=NU;qG7i~M
zVw5<yqytBeyCUEPe$^*EZZDJgk@3G6cn@@>Z|i`^1Ck@Q4P)T(h#Kfe>g@-9wBL{6
z@4@ItVpElTVGQ}-83WIkBS*9IJml-qepiOimw8eZs`9BZ<P430&r$NzAAjb;dbD=U
z20jn@B>i`5cq%cP{ao6UKD{4TjzQlB{AlgkJO+LL82E?Az<*&3`QHcr;!*thYv4!w
z{T`z&8jr<sZLgHKbC%UK<ErAS7W`ljZds{c*YU@DaAgy>Ch)uyZu|H<yCbcEZa;25
zMPvRz-xl0O>cNlg;5*Bu3+6AZ7@bJ|@{e;h4>zFWxH;Al#Vsm-dtXlvKBbL23<(IA
zjzj9!z9Jp|03R;G;)i?Wbvb>pjyEEjC0bv`_WPGKdvV(i_eJ^ZL$YPM+#ZYI?rgXn
z_ip^P>l(c)>T5uO2Uq-cKE>v%Z3YUzOr!<i$)IX)gMa0cC9CRM{4Km|3bFs2A2Pxn
zgb|l`oqOSFzaXT4p$Cn^Pa~z>O{uEVKh#r!YqkH%Zwuirk8>9w?FX2~`c<Q0KVL-8
zJMJ}{hene6`9neS(SZ@SYEpk&D0#g*dE0sPUuqiVN03JTuF<dy+S1!PD2M-O)Ucn1
zO0$`5!WlzsYY*-s;uw!Q$1YvM$J2wQcs8}PwJ+4&K0DNINJu`@QQE#W3{FEHalKsD
zsGEnw06*&RG)aJjs?^&ZH%g`UQWQ%&BLL%k+@(}@b!jvrAK)nscB$y^YKIUFIjrjX
zx8%^}o7+GJHG!T`8*GTgDIOU9&QmFV=E=a9>U-qH^3rxZ;MNJgK)4g5SAp`03SHIO
z8VzpJHF!)nsLQ-NposVrTB@idjgnyg4^dpNaEXnyBsabO*T-Uo*9UdZp+SS|ns@N6
zN6oL#6ZBDEyili{Sh)}1Dgoko`Go!|B}g9+D42R#)ckm|T5eU?Q?}Nx&pY%{e?L<a
zV~1&L*yaPoZ&_-7eSV^kw2zy%U(@L$FDmlem>SmSG5T00xw*D+bZh$++Gs~Xw@u>o
z`HnhvDTB3s&9B=ZR{YgUpgu3s$I}W=eynmv_-_kvRH0A5^zRkwBi<()w$awqEe0T!
zU!P~`<Dpbyl=Bp<kGxEs%CFDI^ik7kF}nR4zAKGipSS5_llG@#NS8N|#;?Eks*mZX
zjp^;bU-8dX^7VP3KI(CxsuTWK-D-SF@$2|W|DDZQno%9o^*^4*ug@#>@hYuI@u&Nb
z=K*+cClTNJ`?mV`H`D3U;qQVMFB&EI^?54aBIA}$pAP>JG&E7iw?3ae`7iPx2QOn+
z^Xv0o{rjlLRr|Hyw0<22JJR&)^W($n)X|6f)J^m2<IAW_wcqe~9CCH&a-hUs^XcPn
zP?O5PUGZ;M{5K^yoq|Fqu!AU3mTtdZ7oSvwx}BOq>(%fdQK0@ZjZjC;KR*pmT5kII
zPD)u;sNJNCoI<)>&8KNIfk?-(u1@j)O7(-LPp9SkLg_ft7JD6zZQ3zPKsxQTG59}I
z>o8>e6aLE;e|kF~Tk7!dN;~P%`t|XOG4#*8-r>I^?K*<yQ_V5FWAJ~W-Jz(`gzA{i
z&-W(NnkXqF5r<)A8bcZ=eXnzuQXvngmHz!I|6T&@0{rRt)AJYe>uCD#?v`SV-=zsp
H*Z;o&n+5;d

diff --git a/fftma_module/ref/refine2.c b/fftma_module/ref/refine2.c
deleted file mode 100644
index 6367273..0000000
--- a/fftma_module/ref/refine2.c
+++ /dev/null
@@ -1,84 +0,0 @@
-#include <stdio.h>
-#include <stdlib.h>
-#include <stdint.h>
-#include <time.h>
-
-
-int refine(double* array, long nx, long ny, long nz, long rx, long ry, long rz, double* r)
-{
-	long x, y, z, cur, cura, curnextlayerx, curlayery, curnextlayery, curlayerz, curnextlayerz;
-//	double* r=(double*)malloc(sizeof(double)*nx*rx*ny*ry*nz*rz);
-//	if(r==(double*)NULL) return NULL;
-	cur=0;
-	cura=0;
-	int i, rxint=(int)rx;
-
-	clock_t t1, t2;
-	t1=clock();
-
-
-	for(z=0;z<nz;z++)
-	{
-		curlayerz=cur;
-		curnextlayerz=cur+rz*ny*ry*nx*rx;
-		for(y=0;y<ny;y++)
-		{
-			curlayery=cur;
-			curnextlayery=cur+ry*nx*rx;
-			for(x=0;x<nx;x++)
-			{
-				/*for(i=0;i<rxint;i++)
-					r[cur++]=array[cura];
-				cura++;*/
-				curnextlayerx=cur+rx;
-				r[cur++]=array[cura++];
-				while(cur<curnextlayerx)
-					r[cur++]=r[cur-2];
-			}
-			while(cur<curnextlayery)
-				r[cur++]=r[curlayery++];
-		}
-		while(cur<curnextlayerz)
-			r[cur++]=r[curlayerz++];
-	}
-
-
-	t2=clock();
-	printf("%lf seconds.\n", ((double)(t2-t1))/CLOCKS_PER_SEC);
-
-
-	return 0;
-}
-
-int main(int argc, char** argv)
-{
-	long i, j, k, nx, ny, nz, rx, ry, rz;
-	double* a;
-	double* r;
-	if(argc!=7) return -1;
-	nx=atol(argv[1]);
-	ny=atol(argv[2]);
-	nz=atol(argv[3]);
-	rx=atol(argv[4]);
-	ry=atol(argv[5]);
-	rz=atol(argv[6]);
-	a=(double*)malloc(sizeof(double)*nx*ny*nz);
-	if(a==(double*)NULL) return -1;
-	for(i=0;i<nx*ny*nz;i++)
-		a[i]=(double)i;
-	r=(double*)malloc(sizeof(double)*nx*rx*ny*ry*nz*rz);
-	if(r==(double*)NULL) return -1;
-	refine(a, nx, ny, nz, rx, ry, rz, r);
-/*
-	for(i=0;i<nx*rx*ny*ry*nz*rz;i++)
-		printf("%lf  ", r[i]);
-	puts("\n");
-	for(k=0;k<rz*nz;k++){
-		for(j=0;j<ry*ny;j++){
-			for(i=0;i<rx*nx;i++)
-				printf("%lf  ", r[k*ny*ry*nx*rx+j*nx*rx+i]);
-			puts("\n");}
-		puts("\n\n\n");}
-*/
-	return 0;
-}
diff --git a/fftma_module/ref/setup.py b/fftma_module/ref/setup.py
deleted file mode 100644
index 0393449..0000000
--- a/fftma_module/ref/setup.py
+++ /dev/null
@@ -1,7 +0,0 @@
-from distutils.core import setup, Extension
-
-
-module = Extension("refine", sources=["FINALrefine.c"])
-
-
-setup(ext_modules=[module])
diff --git a/fftma_module/ref/test.py b/fftma_module/ref/test.py
deleted file mode 100644
index 3827d09..0000000
--- a/fftma_module/ref/test.py
+++ /dev/null
@@ -1,12 +0,0 @@
-from time import time
-import numpy as np
-import refine
-
-size = 420
-a = np.arange(size ** 3).astype("f8").reshape((size, size, size))
-ti = time()
-b = refine.refine(a, 2, 2, 2)
-tf = time()
-dt = tf - ti
-
-raw_input("")
diff --git a/fftma_module/testRes.py b/fftma_module/testRes.py
deleted file mode 100644
index 226da4e..0000000
--- a/fftma_module/testRes.py
+++ /dev/null
@@ -1,56 +0,0 @@
-import numpy as np
-import sys
-from refine import refine as ref
-
-
-def get_p(pn, pdir, pprefix):
-
-    p = np.load(pdir + pprefix + "0" + ".npy")
-    for i in range(1, pn):
-        p = np.concatenate((p, np.load(pdir + pprefix + str(i) + ".npy")), axis=0)
-
-    return p
-
-
-def get_k(pn, kdir, kprefix):
-
-    k = (np.load(kdir + kprefix + "0" + ".npy"))[1:-1, :, :]
-    for i in range(1, pn):
-        k = np.concatenate(
-            (k, (np.load(kdir + kprefix + str(i) + ".npy"))[1:-1, :, :]), axis=0
-        )
-    return ref(k, 2, 2, 2)
-
-
-def kef(P, K, i, j, k, pbc):
-    # tx=2*K[:,:,i]*K[:,:,i+1]/(K[:,:,i]+K[:,:,i+1])
-    # ty=2*K[:,j,:]*K[:,j+1,:]/(K[:,j,:]+K[:,j+1,:])
-    tz = 2 * K[k, :, :] * K[k + 1, :, :] / (K[k, :, :] + K[k + 1, :, :])
-
-    # qx=tx*(P[:,:,i+1]-P[:,:,i])
-    # qy=ty*(P[:,j+1,:]-P[:,j,:])
-    qz = -tz * (P[k + 1, :, :] - P[k, :, :])
-
-    kz = qz.sum() * (K.shape[0] + 1) / (pbc * K.shape[1] * K.shape[2])
-    return kz
-
-
-def test(pn, kdir, pdir, kprefix, pprefix):
-
-    K = get_k(pn, kdir, kprefix)
-    print(K.shape)
-    P = get_p(pn, pdir, pprefix)
-    print(P)
-    print(P.shape)
-    print(kef(P, K, 1, 1, 1, 1000))
-    return
-
-
-pn = 4
-kdir = "./test/"
-pdir = "./test/"
-kprefix = "k"
-pprefix = "P"
-
-
-test(pn, kdir, pdir, kprefix, pprefix)