[Insight-developers] correction

Samson Timoner samson at bwh.harvard.edu
Wed, 14 Jan 2004 09:55:11 -0500


Content-Type: multipart/mixed;
	boundary="++----------20040113203926-074758809----------++"

This is a multi-part message in MIME format.
--++----------20040113203926-074758809----------++
Content-Type: text/plain; charset=us-ascii
Content-Transfer-Encoding: 7bit


When I submitted itkKLCompareHistogramImageToImageMetric, its name was
changed to itkKullbackLieblerCompareHistogramImageToImageMetric. Alas,
"Liebler" is a mis-spelling of "Leibler". I noticed this last week,
and we should probably fix it before the next release. All the changes
are in files that are attached. I submitted an experimental build
from horus0.ai.mit.edu with the new files.

The following files should be removed from CVS:
Code/Algorithms/itkKullbackLieblerCompareHistogramImageToImageMetric.h
Code/Algorithms/itkKullbackLieblerCompareHistogramImageToImageMetric.txx
Testing/Code/Algorithms/itkKullbackLieblerCompareHistogramImageToImageMetricTes
t.cxx

The following files should be added to CVS:
Code/Algorithms/itkKullbackLeiblerCompareHistogramImageToImageMetric.h
Code/Algorithms/itkKullbackLeiblerCompareHistogramImageToImageMetric.txx
Testing/Code/Algorithms/itkKullbackLeiblerCompareHistogramImageToImageMetricTes
t.cxx

The following files should be changed, and are attached:
Testing/Code/Algorithms/CMakeLists.txt
Testing/Code/Algorithms/itkAlgorithmsHeaderTest.cxx
Testing/Code/Algorithms/itkAlgorithmsTests2.cxx
Testing/Code/Algorithms/itkCompareHistogramImageToImageMetricTest.cxx

Note that I also made the following changes:
 1) I added some comments to itkKullbackLeiblerCompareHistogramImageToImageMetr
ic.
 2) While my implementation of the algorithm was OK, it in fact did not
    correspond exactly to the paper I referenced. I therefore changed
    the code to make it do so. This changed 2 lines of code in  
    itkKullbackLeiblerCompareHistogramImageToImageMetric.txx
 3) I changed itkKullbackLeiblerCompareHistogramImageToImageMetricTest.cxx
    so that it would have a minimum at the appropriate location using
    the changed algorithm.

 I'm sorry to have caused problems,

-- Samson

P.S. If  you do not like attachments, I put all the files in:
    http://www.ai.mit.edu/people/samson/upload/itk/

--++----------20040113203926-074758809----------++
Content-Type: application/x-gunzip; name="submit.tar.gz"
Content-Transfer-Encoding: base64

H4sICGiVBEAAA3N1Ym1pdC50YXIA7F19c9s20s/f+hS4tHMjp471YjtOnbQzsl5sTSTbI8lxO70b
DUxBFi4UqSMpvzSX734L8J0EQJB2097zRDOxZXJ/u4vFYrEAlkx3jD+REXU9d8978F78IZ9mq9l8
8+bgRbPZOjg6OoLf8Gnt+7+bB/sHzfaLFrvUarUPWi2gb785OHyBmn+MOunP1vWwg9ALF69d21LR
Eedr6POVP8NBfTj7MD8Znncmv857w8lODSG4OOuO5vBvesb+RqjT681n/emsfulQy5sSc/m6Y97a
DvVWaxd9/zki/8LJo8/3nxn36cXVpNtn3L80rjxqUo8St9HDjteI+HVXxPi05xkmQkUsuvaCNGLx
TMP+eS+ttH8h27TatD+rd0anF5Ph7Gw85U2agoTuL7/w7/PLzuzsS4N6n2L2MwKDY0eIbOtA2zns
Wb/T60+0hJ8RvCAO4wPKR52QIuEWZBTAKqfgFySm3Sli1i7DrV3Ibr8Mu/0MO4u6tufYG2oM7K3j
rS4csMmI3BFzSrzhGt+SATU930wCOb6YMkyy8j3aMSl2T6iFnUdtgSpURsLWs9fYo8bM3timffs4
Ju5qCmoaRClCDUvL8NUYU2uMH7pb5w57W4cMTPtesz3aeKFY7H7a7+m1SoHIsrZPiVXAKyYRgs+J
d287nzR4JCnTrLrYsh6n5HZNLA96w7ZKuqYuPiPVdhxicvquTZZLalDAn8FMbt86eM2xM5v/GhPP
ASdRqVCFWVaf9QY7pLwC7UADHXhGJLghcctaWwVK8++RtW25E3ILSjncOMXsVZg09/6DRywXCD4S
0zb4lKjZhGJkTpKDDY8Npq5tWcTwyIJpCIoWiFHB0jIG2PW6K7xeEqcHTceWQTRbU4zMSxpjx1hR
67aIb5JOziWyZwmNi7AZaf1xOTeSAtJ8TwmkQi41OoZH7wj0kgfBsuSAKMNDQ/p0hTcEpnL7WfQo
5PYcGs1Fyc4TdJpn8qFTBy9YPP0IQ8h2Ssy5hcC0nCldb0wy2P7++2M0YC3iulMDm9jRlFmKSVp+
FLzH2OODQ1NkEU4ipdwUV2pm8W/abEZWhsgsnYDLZbfDHWYMnmT2XY9CxmYrzVEAE8hIhgpo0spe
SPkjVABC2vznLWULZCB9/tJhqQTp8xctSgr57+vzP6jC/0Cf/2EV/of6/N9U4f9Gn/9RFf5H+vzf
VuH/Vp//j1X4/1hifDUrDbBmCQnVhnCJMdyqNIhbJUZxq9IwbpUYx61KA7lVYiS3Kg3lVomx3Ko0
mFvZ0WzBlLyxYaWoXHZnyDI8XDtIpEouFAqBOTlP3TQqwSIt+wM219gaUYtgR2vOlwIyfNcEW65+
KiGjz3DdmuYNNj6NCL0xifPElXwlbmmFRnhjYoNi6wm7KmV4ZKQHROeE3q5ubGdl24tgLay2diFQ
LCdYxT9WllfIIC0XsmyPuOOtt8Xm0Frazprbplw6XYpJRj745PTfW3AKV1eYFCHlXHGVoM8hI7ny
vmoxMiNpMtDlnKPMchqeUOwOKDEXfAPQ0NubKMJlpGxNj06Ia5tbRnb5CBalC90m6ICz8jIuWdUT
SvJJazEhxtZx6R15UvPLcVHbvdfvEd4WCMrJ2b6M+VU8lNLLr1Jr5ZigquIli9iyCgCb6iqIU+TS
KrTzKmRcuILXyxz8HDuOfX+CrUWlXXhNuEzmbAUBemWbiyckCNW4ZTQi90oR/u0MhpnWpL+TReLU
RTMo6WA1pF3a/NDd0w6G5bhIJ+XScnWwsvY+10RQnWNaswvP3Uaexim72DS2ZlEOX4TLSIF1in1r
W9ic3mNv1e5dws/iwVAAy/coTE1049gGAHj/9OhySRwSLsoKzIrKckkrkPGF7FQki+k6yOw0zo67
TiFOtHvqiTpNl+VCLepR7kQlg5QGNHMGwDaNe8Tzk7KS0rTAAnn8MGTcuezarjfYWkZRXqEEyfiL
gnM5YUUcMpI3mBn+4uZfYA6Jv0hlamAz0rYb4txRlwRD3MSuS5e0oL9UoAx/Nqy7K0wtVj2lFxRU
mDT355iPq8/Cp/TmRnlMExFkz/fDZFYFTlOlObR7ejyydFlPW0HQ06uOEdGmuV1Z7go7G1ZAU7IH
ipE7KCXqo2UOBjPlblBMslMbDupX0/4c/r7eCWv8UqfkcKOIX4ImLPKLeYYSpt3pSCiB3SiSkKDZ
CWsQE5J83ulivPEwdVivDkq+FDkk3ZvB0W/uOLa4SwP76+EzQm3Htmzao5glN9KpL5AgIBayS45r
bfULoWlR1xguuyuiscL2BcgAsoQj/KYxE7SzKYcMy9yJlYjO+uNLv9yVXaXWbaJ29UtjRtYb28HO
Y74IKto5ko1uPhYkxEwqrzrN3dy78z6FNbOpmtb5dNKdstVpqarKPePhwQepKiMjKnVxY0imW5aY
oRdUF8YUUc1g5lKiBDC8o1u7F9Fr7YbH1BXK8iJwrjuDG6oStZCmuMQsQamqEgvJstVYouuyWqqI
VlYUFRKUKWVSYgpLfJ6GnrcjfFGFT0hYqiwnBBUV1uTo1B7lX46rXcLrlU58ooYJEpqUQHklTIpM
vCmmQTJv6RC1dYj2dYgOdIgOdYje6BAd6RC91SH6UcuYTS0qLZu3tIze0rJ6S8vsLS27t2LDpw+6
o8vlJ0XpAXBEIDnKDe+XOe+MMEUHllnCwpPGEFDqeDACSY/5BBSa4UU/M8gfl8V31AdeEZ3OKU1I
XO5sRyZCdSIjwWhEyzLnHhUw7RhTbkNVisuTBFvv0Z8aG+dKWun2s8A1C2mrbymHHIo2hSM69b5u
UnntvdgQVLiXGrt6enc0vl64uRllCTo7kyli2fZinqhoWzDOVIo39SJa1c5cRKTYXwtpyuyKhZji
fZwo/Qw3xuK1QXKbKrya3b5KzHPqiqxklq9+xiO5apEeB0bjR+/0ME9e2prC5W87Xv/mH16M1q+C
xyQV96IkRi97fnLaLd9/yg9v2b5FSKm5x5Qmr+LYoq2mzL1L7HiUsZKvdor3lEJK2eZQ5Nb+NmHg
KWxXpv9Lv3s165yM+ukHYePneJHkOmOyU5t1Jqf92Xw0PP8AP04mncmwP5WyGs4+JJ7Chr+GF5mH
pws/ADrBsID1G8eZwDJb9jxvQnZqlyn5RPOX4EllmSkY2kX5S8wAOtqD4Pxm1Bct0/miJTJKGFNs
NPZryvwJ0nPD1bFCW2AGHiRyjWyXbWVb2B49xUX79Ehf8HLp3e/7P5eae/QluLuGi3by+/HBw/69
zqwzn1xczDKPyxdHB9GZcZttnL5+bfiRNXgTQSTiSwNsSUxYsiVeRtBIShLK2dtYtym/8tkygfOL
q9nl1czf79VlpNU6gcBEO4bWZus13OXBs+qGWqU7QfzWgefvhPZzNbT9F+4ErhtqZjuh9Oz7XH1S
ftrXNURVzpXModOLH6nLcrJre42t/iOZmhQSXn+tE3KVNE6L1R9rmPSsd7jXRIfN1CV2XhW/XyUZ
csMw3J1eD0/Zq1hYhFeeVs0Ms+2rfz2cniWmy9SLZFKHYopXwQiPsZgI9mKZ5EthEgr+2W/g+XM/
imz02WSo3//UOmrCvfT7nw4P9pvf3v/0NT6NVz8916cGI/bS4WvPYxizQ8ultysPJaMO+jtKLiLR
zLbNT9RDbEyyUDG2F1uTHMO3EbZutxCagFP3hx/gQg/WYuzGR8hbAXoMwaZrbx4dLqNu7KA2eBA6
53yxiUb0hp2VI3uJxmRBDZiMah3TRJzeRRARiXNHFnuwvicEdS8uf50MT89m7CVoaGk7yIh4L4iH
qekC5bNZ6lWj9h1dAuclqLWoz8fT7vxjf7JT+27j4Ns1RvfYYRM6qqMFddmuCzpGB0dv30Do/Y5Y
C7qsAQPLMLcLgt5TG+xJ8PrnWqOBZivqoiUFxAq76IYQC90Si4C9yQLdPKKTLTUXiXEOcZHD2Dzv
2ShkSrCxQitO5nNjNgl6NCH7ZZkaALDtw8ssVlYKICKWVwQIqE+wadq2NbDhtpSluBpATgiehM1z
mHP8DRU5Yf7UIdxiKYNR2yOoSthbia8H1Qq529KSBYWGOmUOIli07ckSQhFB4Z6VEJQuM2AdJ1fB
dkgH/CbgfMUszbf13Y1tLdju9hiDnAcIXAZxXTGHckUYIh6anlCi99M7nwU8hVukAtp4m5cZdb+n
IBTXkOiRypVV150IAIP+eMBOERLHFSH7kY0XYoSW3oPBjLmnSR7A4f61vWVnvPaEYLNAocGMEc3s
LLgQdl1R3HUlecodeQlAVahTAFFTKgt3FO6iW/CjC81X7pQVrqz9ETEJCoAyp21FAUVZNySgLxeu
VHVDAvLMEZck1nOqU7wFKmyl6xhk1KKaBxltmkqlQ1DBFB+YygjFlUcy6vxht4xSwFSlb1FvBUSp
48lSxHqaq876xPQVk8Khx9JVNqRM24W8lB8/SRv0YTSOz5fVjIU1Nrn8qPSxmkCSshhHEVJ0i3hE
UFUtj4JeWdIjwIkPDkWEk87UxDdDSLI8HheENAN1q1SlNyJy3bIjEba4BuJpLIRlIRIuPlbX0lpV
UWpgWXr9xoRLpktMnXTmLc+5Y9TMoZDVeNrAUguv0ksuvSKnYqC0bksDmqsPE2L0y4p04Fqo825X
02E16jeUqLL0ujvgIibSwrFCc2iXkSnBT+xKRQGYiFpRViYgl7qxVIAcEXpLdjrWtZiiGk06KvMl
LhrsRdUwUgmT0xNeY8VXBFKrTPBjF9ZKiWpi3ynlI0i/elQI1suXCurzBAh2GF9lFctwVVax5cez
uJBQYenCykMZRlSAWIZWsiTQq1ssDSxokOTpkwRKpq4ECQOjmsjcAy8ibOaZEhGJ52wN7wy7K2nk
kZVnCrgpJxhFnzxhYrrS3JBWF4MKABrVfFJUVTuUO6gWMeCPGVcJPT6ySvDJViNqkJz6RzHCyVRd
PyIHSOoYSwMkA1gAKhy7Ud3kib21FuCnOjR8ErtT8+v/e0vvsMnWABPIqG6IqabXVBPWAfSSbnjJ
Cj/BBHcHX1hD6piLChFIS4HAbDOW7ujQOUSXTOVKWfq8jneWOV8uvfk+b2ANMhq0hkCH6oh93UUG
xIxXCO3UPtf4/8fjEFggWaj5rval9v+8YuL/1oeKC1CfVYa6/qPZbh00M/UfB4eHrW/1H1/j02gg
LyoZ8AsSXLhC+P98kygV5gUCrAiA3fPYHwsHlh8OKx6AMImwaQIEkYcNTOUuoowJ9tCjvUUrfEcQ
Rstg+keQwJlkERSC+GUIbi0dmxosPkGPwGycC8CMfgyxCm6gWq12Z9MMr7oftWIWc1btQRb1t0dH
h82j5s47uDvpnw6ns6BmXfYfaiFt0nYJ2n0ZreaLd2Xwaq/xFfNSvNFEDNB4Q4YYWKEeVcWo+E0p
ErTolSdKUtXrTJRARYWwBCh9q4lETvzOnXe1ZKFRsqo/h0q+WeddWGaURYdV+zl08q05EfpbrvCX
/2iHnSdkBWz+h9Arq/88aO5n6j/brWb74Nv8/zU+f4n6z7jykxd7T7pTlo4c6/+/cMw5d+/Q99LC
0aCM3P/ebjb3G612o/0jar05Pnh73GxxaFhaynUgd9T/o7XH7sLtdLlp1Dx76d2zhyBg0nFtiOzb
9R4SVpoixEpNh7MPEaM9Vm0KOdXK8zbHjcb9/f0eNHnPdm4bZ7PxqBETrrw1z77iYlTeIF7t6YYa
wPcFsy292bJqz+vh7OziaoY657+i685k0jmf/foO3UMyYm8hU7sjVlDpz9M9SJQogICRgy3PL5zt
T7pngOqcDEfD2a9M08Fwdt6fTtHgYoI66LIzmQ27V6POJOB0eTW5vJj29/ymMr74xobsL66mtWyP
GpBissasbaZ0fLLxFyuxTRfcBIUrfgGHv6uYzUz5vdxVv95Aurf/9PKDXIIMHej6x038K7MvX+U3
Xr3yPYbn7xDQ4Nu9jdq912HrXFRntR0L7CzAdHfUH6X+ecGU/k4a7Z29WuR4nA17BgjczV80rPlh
VbJP0Q3x7lkVMl83gDTKNAcHhg5iWw3aQ7ye3Jj47Z9+ju85j+gzVwiWIV2HgHW5EJdvHkeyGIEB
wxMazUYtOATf5uASenTt17Khn1D7HSf1HjcEHAct7C3zi0v6QMwZXBPe7dq2s6AWSJ6QDRvsQaiL
AYF+EBNsgynIxbopXmCF42N+/X0kbTet3s8QIu/CjDGvjRaHAdxZZBiE1usRw2RBJHQF5Du5QE/B
UHif0S2lbPBeJ7HKImZpNZNqy1gJpR0f86UIcfK3gxtpJiIxMY/s3RSL0IjnNptfmKsvbejuexZg
Ftyw3CdcxCxs8Eo5WAUjFmsffQtHjmnagFnGwmDU/fZP8M3PkCHtIgQ/v7wTINaJJkoh/MfStGFh
npCwwQZAAcEhe83lLmI/A4xPnmQf0YvIhSIuIPZTy0c0OaLJEUgsQU3OfxR0+Dp7G/0UPE4mAZ6T
+/qO3wK1GyQN53OOWIuBCc78R0611z/DEpf1WT3ThwFKDPBtVM9ZTQkKuq6e700lLC5GqC+x6aoV
Y+dzpN4+ZF2WbHfWbnGz084eYITkQaOzjqWChE3O+bsKJGmwWICkvQLrXG0WEP/BG1io4A+m+KGC
r+fZ0zfO1pIIiqFF2OzwkA4MGF0CHU/BylsPZvT8cJAOBBY48irnOIVBEoyGths/X/AzmPy0UDoL
el/yvRfRJ92+3apsMhavyib8wATqZz75qS66fnwcpXfsTxT9xU7y1EDutC5HxV/T5Dlu0PUwjUGP
sd1ghkz/mejjpKTI8fg1cJTkzWxsTGVn1gmFCROSssM373JsI08I2gNjlesEaSkPKj4i/Ou3JptL
OMPMjVbmhq8mH9iRCB6kQgSoG/rxpbm9DdZOLKFjSts5r04wjD0tEY12cnITnpQM1AnBwQDCkP6G
6X1+BIlWAe9VyWro+tmkMeVkic5KXY87OtIJDJshSXR3or0RUT1uDsSOgF99Rxw9IEGl0ZLKdgQG
UK+7shlr2Ho/q/8ZDRPMMw3P3orbntQImp8nTFggScvtwB/3z3W7yBQJ4yVF1FMGUZiwxxfF3FUd
vrRDsEZ30P2KGquEA6N7apqwfguXeAuFR/tLxIRf8/h/suW1z8GDelyP7CQQecsmiia+mEywiW/H
bsJlnG/XN8S5WMb0sRy2yVBPBZVP0C3Nd/DrfcbP4doPP4Tvo4mF/faJBQg2t6dVd3nUDLZkHHrH
n1dgJ1u8WJLvZTDqRJz1Ec+lfDGbVIN84X5jWgH/RD/2oiZMPbIZEevWW/mq131kouPO7XvkxjM4
jF7K91GmHnS2YGzzu5pjXIDUDJp5cm2PlHCSh2EJQHsohpYcRsWV0lkjJqknJx52kojYJl4wUO+w
uSX8GDZ2xr3cnDkm2N06/jy59r+zHDAQBUp+ZGzqsfPv5Ofd2FM4m1hc2iynSZ9KcNxNIMIGud7i
+NhgrXn/Hr0cJ1v0E3rJLobKwldOTKyFmUmwGDSWGAITQzOJBV1DU/YfiGNQl/g2re+wFrBK3lRr
/HuRnmLF+XncBlb2ECdfZsUh9jYTgz1Mh+p8cuo/GGTD4oT/BNbfEXlgY/VzgjFxHM44IgUG29uV
97c89wyEPAgoggqe/i/D2XzQGY6uJn1fsVr63vSq2+1Pp6zQ588+m/j2+eM/1Xacy8koeP/L4WE7
rv+Bv14025Dzfzv/+yqfv/r5XwXnLHUQeNBothqtfdRuHR8cHDePlAeB7W8Hgf+jB4EWWxbO51Uc
ar6qfecfI1ZnkD6g0zvTq1mQNLkbbPBaPH7cxQ7w/sGPD1BpLfyOeYXQP24cCsboJs/tPrBFr4NG
lNyY2FjVP4x2wtwyeX63CtmHvMA30ud6bCfyxl9Wsno88CNeHMjWK/4yIWLhp1ivQk7cgf2m8ZNF
8EL+yiK+MuUyWL4ZyONJPefsZ+Yhk5QOwcsJF3JBwbrWCWwweh02+TEAM4VT7XazDfdfSRcTMHde
+29aYveZ9ibeMHKW/0OkYG80ge8hn3j7aJepLTQUGnrMJjfYZfbgPRHC41OmDd4Q5zjd1I4Jq0IP
dfeme6i72lq3u+galvQU1BzvoWtimi4aDoe7QEgewJpg6XPbATmrXWbckM31Huozs4xgNJ86lE2T
u4g/d/Z6bS+w6a+k04H95jFEs9r6NURQ0DH02deB06LwTQK7EET913RF3Hz3TPQt6+3AZ53XHddl
/rXwt1ju2BPmIPQ1Gg+73c6QxfX2Lhqdd6eoffD2x92Qx2azhw7bh0B4uN/O+UW07wHGvqUsKkIv
fhjVL+esm//zH/jNOmcnGkn/be9qm9u4jfD3+xWIPniOFkWKkpWkeklHcRxHrS17RLrTad1qKPJI
XUPymLtjZMaj/vbuLl4OwOFeSCptx3PoNKYOwAJY7AG7i8VzF6y/mt9+Dnvt8OiR8WI+/dVis2jK
fONRV9QXf7YklQe8qyRq0yhFOTorXMbR3fAunHG54b2ar2HxlAcUkkomojR3LOzxf46IDImZEMA1
u0PHami90hY7PoL447J++fbV4NVNX4351TIJZ8Br4NFwPIYZAKmFHSQmqkj0Lkrvtd50+Au3jHGb
SdX7Mw75roos/i2IIxroLMDysoUEdqfZWKCXzWZrfKfR8YCviSSzIMcHS+YYUozrxHABb0hA5MfZ
a5hGKQiWKCxGr7+B2rvWj9q0JQ3R74z18R+V37Z7d6dexfkKjLtZkCS8F73g4KTDOOA0BUBDhYcF
D4oe3Ue09cGqEAhyOElqUFHKiKMhfvANXhAQi4Sf8g4ROBt5TyNmwSICi7BDsgtb1BDeScyEtntS
vEEPI58mD8Umj94i4NM2gn5OcXOxp364As0gZn3SiNkgnANLY4d49D+8f//uZmCsrAg5BroD9HC1
xN5nkyA9r+IVxwXhBlQixAfSpqM/ug9/i5b3cbAIh8S/OyDh/3h10x+02PdhNJfLhFQXRomiI4kI
6DPmf391c93WtarFXRgvOosgbbWFvxNj31eLTFgUDYkk+JIu0JBQqP7iXSWKjmAwnTg4WVoxfQEb
Q0r7K8zyT8FwBm+Ff331U6ujMwrfolkSFfQgWMURbQ7sEkivYcmTnfHfv+ixm5ve8VHv25Y9fcDZ
OFotmXVBOUZO8UJdz5M7LDvne+Eg85y1xfY40Fxg33n8GWKnvvorzvoWGsiptwQ9A37UOM4z+mP2
BDQhTudUqER9GbLE+yhOAWApgYFmhwKb97cfzCZnGoWdus3okyDUQ51mfz6MpXfwHBv8jtVJWvSJ
kxKPdqpDD62VVJETHOX32TkQJkZVhaR2xLjacL2GvFZsQiAxa8HokL6683Y4iiMfW24pcjerxUEa
zoUSp0eH+bhJETgCiD/3vSWtjB46GTnBjWePb/olBbLuYSuJcKkHMSqtkzia00BR6RJiZQgT/ov6
uTarp6d0iRRHaCeZcVZFwTzQ1VLuEK4eGXVYb5NxyI+Tki4dzNkhU35qd8w43fHMvuXOlevR/NNw
FN2Fw4XknTOjkpxEVlNRdHLEuQySn1JausfdSFpGZYcsl3uWzIxKOmagg0YnHxVYPijz4FQbVD4+
sWaHLClzZ2zSMb2eZ/Utt+CVEjUCDfQpNDIq+2YEE+ikchnVw+SyM1fof+oNymD7XEU2GOuPcfDL
KliM1qKn7oz6g+aIbrwb5qqsZdQnJyZPsdGVUUnNPpxXE+s6+a9NyVgqHRlqz8FDZKH8a5poZlIh
/rJuUakxqr0RQ0Voa5Q2SVuGJGc722u9FcPSgS1OTCU2t4JdLhS+FfUc+nOoWnuNGFyLpLrJG36W
hKSEwQUKcHYCyJB8GGdmgVqleVPGgbf7aBtvfLZEPDfU+CyPr+a3gxqn4o90+Ic9/REKJhZXpF9H
emEyjxJYsPStC7J3h7MVhZpDUeo1khPnt+ztCvqFhipaaNrAJ2EMGVSc7qzqZ7ukXD0w3zoWRK6C
YZzSXUap7dI0IHgJKGdkWYt86tmCDkrJyvw4YhRukl28hdKrhApyJy3X7VAh424YqTFvrGvxmJZf
w5ii//+9Rf3Pj5ItdNiK6iM/cBXY5s9YlKCXCMEEob/4jxABJXg4nTQj2asi7yurIIXxMB1qU4BV
ZJ6qVTobvKlXQgAKLz0o/VHzLPBmBY90zUBSE898Y2d5lpHQByxeu1ya34o3kyQH1QMCzceggRUs
RUkwW5N/ge5HEDQE+Uc3nzGcId6fZ3zykZ/Rki/nF2Z2RfOP0NdHOowHW8Bwenc85btHw/Pt5fWH
yze3V9f9weX14OpycPXu+mmwOyWmPv2LXdna0/+/Ps/bNG3LtE3aqMB/+Obrk2Pr/Pe4903z/Y//
Svryzn/TDa+C6ifAvdPet81V0OYE2FrWQaR2PQNGEk990VOV5bc7u/A6jaNT5d8fjhAkl31Igtv3
YFiQOrtAjRTeRlBThz8H9LfUH5GArDoJp6hHonA8oC4Jj8YRnlQt2F6CgZ57KGG0ky+0U2Gg8HMQ
L4KZ0jn/mD/G3tjzvPldkGIv7Onp5uRAIaTTd6Vb2UrXBT9xIT1HnQuZh1K9LhhDeM3gn0eIirEx
D0i9elJOeGi4ahZIkdJLQzfNXT0mdZvBKCt6x/F89PIKcO1kiYWmkO/K5m26ZU/NhjYBjw6BBeAh
YmBszimGHi+Kkc3JIxXIOsKLlpwkt11gKYhhQc7iGeI5r7n0P7XX/Bhb/DzQn/1CP3HrsXycFu/g
3QAzPMD4YbCz0uRcKw5c/1sQR2TPuF1FoqtjWHLpNcusqe8DvAhIlkdFVVDpzapgZYhgandNafvL
RsmtYJqIB99Vtq+oiPadVLKuPNwj/pavD/irC2MQPF5AXXZ/pc7cpwGoFuocWTqQQv6dSmEnqvD4
OPgFOqMNEYPTlVfPb+1rBqRWX8zaWNTXullS33PJw/6FQew5CJOvP+jqXZWR1vv7Wotn4pE2CP6I
QpnBfPP1GfzqwpgKzkWCVCL4B+7O2sMDdNkLiptQHgM1W+IYfyzg/FPNu4URBZ094R7oktMtizrI
/Gmo5/DvqukxFwVOV6psTZtThl4j5haU1SbhrIKsNZtZd5y0xMqCCtyIfBfZwO0RogcTuJGy2TCe
kued4OWGtjtwFK0WKX6dbjREz+NDgHKr7a96TEnm7Lgc/2uFQT8DB2vy7NqnmZ7fGg72vx/+47n9
qIePNKF3tWZxLM/FLVvLu16AuvWkm6DtNLodDZPUWD79wi5yKT/Albq4ULeQn3zShVfV6s12Oz6p
L0+sx23qKMwpMhkBrCLK0tCjBC9Q8Cf8DgafsFNxH0UJqnnNAjhT4MjSXUtbmw+/i/2/LW7gJnBg
5fhfR18fa/hfx4eH6P/BYo3/57+Qvjj/z65wYCffnr74Q+MD+iJ9QE8HB7a784bf7nY+5ZdTpTlB
zh2j0OUER2Fge+2CODYI54h6fhckL6PZjOIFCAv+fq8UlDfL1L5AjBhjngoG3QVnjKi8W5AYJ/fh
RESSnrAlxt8k2RFbBP/B+/92uyPzgoOpZemHda7bs4nHw+tjyI3DaCX6wMtS8LR/0Dts9w5bouWu
J+HMtlqx8MJwBm5GuOtet3uwQ/IyRLRJtIpNSLRdaXtuQDMDf4GenqvzfBxb+0iPRiyJF6pHoDhw
qV59x1Xv3QjZHQJCAVg97HMeai6H0GOVeKTauVIqfCjBa+MX7DNhfCHCF+Fv5amCDvtpwKMu4RfV
OGwXleYvIBXnwBBnZAfgr45Ea+Itt8wcasYXbbTcXVcgGXPBcBcXFECGKkV37N+gPYlRW7B6iA9H
RQtfolc4Klj3/EtK3qDxmqSVRaW0C5O4CBYJalAWDM4uoo+N8jYYWla+YmRZwaqBiZLWuMqmzXxf
qqbPLL3BWB0VKwbtqFE1ertK/ek1XvaKaTbKbsGC2tOer1CXAW4xIN8xQvrz+CfUEvFKzA9sKnZw
bY3UFklTc8mhQpJWeBPIT2oWBP8V07OAIYncABmabk6rYPEXREXmDn117wkmeVffPUGTJO6cO6Nw
txkFMgiW/+KYVj4v0BJL8WuxSLcgt3vUOdSK90qL90RxTwKlCi9YUrMNvQ6MVasF1pWjzCIKk8As
iMVQ7gK6CrNasjV+9BUXmjtU58BwWs3wQxEF/GHLM5nFY161vLGU6b5wzYKOjx/apRCk4lpaIXIM
an9z7p/kHhOXOR+dsgNmoK82lTbnKL2ReVFgKRUl+TFKlskmm/IGTOl2VXc0OE2NulbLuOGH+sEL
nVn4X8HDq+QypbMM/ZBiyXkE2a+lgoAicyYye/nMnsjEVXV5kAk8Ptm/yM2HJVKo3IxVA0bWmrIk
eWgU5MDXwJVQo2uB6YcC+Ry/YeKzA/jz0/NPbJ+tn69Zq+snz5MWo//JA4g41M4dUhdr0nLWpGWs
SWuyZjsepE/Eg1TnQf6LK0fHxy+Ojr/h0jN1is+0Qn6mpQI0rStB8PI/tRxNCwWpJhfFqYHBOH6z
14fFsK0WyVZL8XtqCN3UKXXTCrGblsrd9PcVvGmh5P2OTNOkdCfD13OCL9JHkHahapmblnuHI5uK
bck2JB3wjJKWjc9YDt7pVQM68msNxaiOBfVzGHAlJLS2OLwhfTw+Xfvm1pWF+VtlnsBnouZXh3Db
le5m0JS2XaVuWXDtMX9BxXsaUMpiGi6MuzJaTzcPGTLeE8/BlgjDFlywPVUWcK+3ORguV1RH/AtW
CrhBCWIx+qc50U5K6vXyChBYjZfd3Rl1IUabmSKY28wbQvJVAHKreYRUk4lQ1PNe4oQKONuT98B/
kGV9dtI55G0bVXQQ+MI6Zkey4AW+5yTi2hQvpdGWt476yyGGyfV5aR9R+AVVJ8rx1y/OvM0hjgsA
jgvgjWuCG28CcSpfqVrg0iVoqLjel+Kh1kNDhRmqA4iah0PdHAxVkw0HiG4kjx8qXhHpnMn8cQ6Q
UvZ/uqLmVpE8Dqp757BXAwfYa97TVlTJXkosL2V1W9ocWE6xgrmogLh1aDxKZEIX5iu7CyZ49smX
kwLs1yfUbIYTjv/swF7efZM1V7fce8ijDytfemRyMQK0i27urdSuXordKhSaomd/PIXi0ZTjhrlm
akYqW3c+hMn5hLRTvij55oAxhJqWJfyRx5iGp/v7TLfJ8iRG0YyTwB95EvBUkZBEDMRqGoxwQ/HM
cOLzjl0QTVk3q+2u3ztT2Y/i1/4+5cqMx8wUdTItmkxgNrq/kmetmGnloNwGu0rGafRtN0vPM0CW
7bNjfgaMZSZ8VLJHL+iLNfxc+HGn9r1K5GZbWygEZi55Wfh3Atyn/2TFi7+591b80QGNKU79vTdR
tBQxryYQMjEEePExtRj3MR3boY1jWTiPaSwlRZg+tGgAe4G7Z+KP8wv2Qv7eR5D2kyJRoYmhglw6
bMzry8VYg6vWarYlu3XgauYEgBZ9AgZ8TG3Uavkoo4FdOj/XBiydVIEEpFahDeouTeLsu95d6Zdz
w3Ab43IM59Ga5miZzbKecxMgsNdum5KnD3UeJiPoMynY2WC3p21L5DXeNcGwY4r7pBhOm01Y5N3k
JRaAPd6URXQdWgSjjm4J0KV6ATnucUejTtptEFS3wtUeMrgwB00g0Uq+kTIjqLoloSrVa6rUeLLb
8mwUda1lLkBiGfoz3RuTLhDN+3Qf8CwTqt4sjrKIlDS1zCqg0RFl1YOSnlkI7/EaIdgJynII0ith
2e/WbD78GZ9K1qVRxK9/7Tm4bc1EXUv2sHN42FOdj9diYj5bE2RqjTzPHMZfTGj98tXEQcDqP64b
+C+h2zvB7dmzII5bVn/NPr0kZHt+I0lWzqHoOyoKwHtnl+gn/B+PIlKHI2Ej7mc+gU1lYqK9xdcf
3rzZO7OJGNtegSMFazLhIl1rm5x72nec9LIpVxsFTTcrme+sk1vMdY2ZVqqn8+sFzbfHm9SkJjWp
SU1qUpOa1KQmNalJTWpSk5rUpCY1qUlNalKTmtSkJjWpSU1qUpOa1KQmNalJTWrSl5r+A965eYUA
GAEA
--++----------20040113203926-074758809----------++--