From 703f63e55105981f23a492aa0ba1da9f6d97b63f Mon Sep 17 00:00:00 2001 From: "Himanen, Lauri (himanel1)" <lauri.himanen@aalto.fi> Date: Mon, 8 Feb 2016 22:27:04 +0200 Subject: [PATCH] Added unit tests, updated readmes, updated gitignore. --- .gitignore | 8 +- test/examples/Si_bulk8-RESTART.wfn | Bin 14132 -> 0 bytes test/examples/Si_bulk8-RESTART.wfn.bak-1 | Bin 14132 -> 0 bytes test/examples/Si_bulk_energy_force.inp | 63 -- test/examples/Si_bulk_energy_force.out | 760 ----------------------- 5 files changed, 5 insertions(+), 826 deletions(-) delete mode 100644 test/examples/Si_bulk8-RESTART.wfn delete mode 100644 test/examples/Si_bulk8-RESTART.wfn.bak-1 delete mode 100644 test/examples/Si_bulk_energy_force.inp delete mode 100644 test/examples/Si_bulk_energy_force.out diff --git a/.gitignore b/.gitignore index 7d4e30c..e79b289 100644 --- a/.gitignore +++ b/.gitignore @@ -1,4 +1,6 @@ -cp2kparser.egg-info *.pyc -nomadtoolkit/nomadtoolkit/metaviewer/frontend/static/frontend/node_modules -nomadtoolkit/nomadtoolkit/metaviewer/frontend/static/frontend/bower_components +*.inp +*.wfn +test/examples/BASIS_SET +test/examples/GTH_POTENTIALS +parser/parser-cp2k/cp2kparser.egg-info diff --git a/test/examples/Si_bulk8-RESTART.wfn b/test/examples/Si_bulk8-RESTART.wfn deleted file mode 100644 index 51352fc1f4e1c3dd50ec6357065b7910a8139073..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 14132 zcmWe(U|`^2U|?WmU|`5#U|?WkU|?WjU|;}=D?s@mF<cm=#sR7aWClntIW))&kRw2P zL716=fdPbB8BhU84unB+AUxdUL17@kz`y{*f(#4{pfCpUL2@Tr&TjD+Ilo_~{?A26 zzwP^NguWjK(UuA?!Ssgh@ryzH)lVY9^o8(aGeP`?XY==1IBmC&cB}%^^4C3%g7~Ek z_Lo7lSFF?>5S_YC1x%lqSGMZ<qLcPtch`gInp5&8L40UzfTN3nfr0)sBLjm2GXn$v zB<XgZcNgrVvJC=+RTtV<HZ`d8WlP%``CW~vnI~n}qc^3aB1qaUYwZH-|Nn*dnKIu$ zIosK6U*PR0=fcy3_xV>lG)#Z_Y44K19MTQ7pZ0b<V$3-?t8DMWiT+o;OFr#Adf-Oh zmy~APhW6<jLH4vCb=xSSy3qc^L{X+l9+Q2ieSWreIGF5<P7Jl-m5{dE!)bl<k%_Gx z_s?TeApz=ke{Q~D-sz;cFX{Z%3-0SB?6_xbd+&Af(_W{g_D}WCUbcO~aYmu^jD+37 zE7#@csXg3l@^TT_9`1@$8$k9n?+E3OK54RV(!x_&o%f{e+<!DT=4P4fV?Mg3INrx) zpGhU_(U*r5?UXDT#HWRc?UM{!YkEKC)7}logMM%Aldzj}UC=(^*W<myLT4*<6C~^+ z)Uudm4>#?7{upA<pL-|4_MD8|Ss`yAZFlY4md`!&O!l3ucjJ!8HrW@El>a-`M{(a$ zy*<7c*v<DzOKf{49xuMHPXA`+lu8M^m#1!DIngI!_v+nMnb#FpY_A39PV71GX>Uj2 zffoUq4{d$;{O5qeZRYO^Mv#9lYkqNkx75V$YSM&_MOD)KB%Xa!=9(_G&x|#gRsX!8 zUH1R)vY(@r_XVW>&e*4{U>7JN(lguk<6efpy0i3ue%#BqW#4{l&Wn4uxwTxllO$nx zKPo><Y2Mwvm&3V_gY4m3@;pXVb)h{Iqxa{a78ASB8*>({-)LetuOyk-OVY%S!AdE2 zpSgmaocN=gK2er-O$Q?`Z(U-wFM65Y<2iB?c6x0Sb~-dl*sV}^*(<a`&Gxy<?&I6c zCF~l#J-%k&Q?Whw#~bXQiPrUC|A-t~-Tbv#YTt@Bud?SCB=<f4<Y0Q9$J8#tfkSZW z6?wbQ^Dn+%U1DjsVaoRx`?fpVIj&>pU#;<JuWQ#b*1hHucF&I~yuT!Naj#kK(Ti;R zKkgOMljgrZ>#psr9!R*|T~Q6TXQt0|jqXONeI<LcDh$(2?1H7zAD);lwXZ;6#w^uL z1-syt4QnJKwf3cK^Am6tRIuCo*1>B^lZ2hmQ=uM~pC9)=kWxJ{N93WcTE<ZZ+k+qX zw)Ag!5WmBBFDMUzuzQ~P<#X4M+rLf<=2t4)WWVvmr;h3u`}Ti0_;t^WSNrzQd3keR z)2n^^r?3Za>{>8+zeq(;^}~ge_vdFz-@LhS^8QWU{wsv;cG#CD<-BLO+hN~Uas0|8 z`-^r8^$Z#Dw>#_|lRrq^5Lj)uE_xTs&8zG8pKOjzj$1Zkzu*3nqt!3=*@p^;{F?D< zpZ$+zywc6D_U-?$A?APkg30!|Og(H57EZQ5aco2O@r9H3A34_cOz3XM{=FX!oEYzR z*eAB?cb-xY-7k6}Y<c|cj{Qp|wx>^ZOtM$F-*oNp^7Z?DUZlNQ-!x<Yrl;{oK=#Zn zlvSGfYTy2vCj#G_UhT6l&-hu?wqUaT*$or--d;F)zn^W-iwg@U+qcY%{Q&aM)+`+p z#=9N+H?J|?VkVnvzh~aE74f$__RAPET)ix?df%)hWAk|r*YAJP5?ouoe#ZV+Z|}YV zg<I^oqd!3Yc`q#iwnv=lw0`q~$@^QmgO1!?IN3g~a`!usJ)9=BFF^ib=dl9YbI`-4 z=#g^hei^N)8$scg-g*_}pZ(Lr#k=Fp&fD)OUg7#!Ylr<8e;<i6#q;;8N2&CbPMB|h zbH&HJg_ZN|_goI$`RqsW{<DnoQ&qT2?C0?WKiTxB*nVfIbXHkXrhT!)`DCepOnY7P zZR#A`O6+9|Opj%HXYOZD|L|?e`$l_jCWXQcY0vk+l;srfk-M~i+Y-jdAbT!dPTNyD zVgCLt5&tGEten4p<E)xsiJ!&$4?R4#!<V~cfBw^(iO2sG@3*?rTUefyx!+B?xK=tK zbANNFQ~mU*%lGl7)v;uIXYM~6dGIUul}-D&IQy#Hjn3PzUzK-wyZR1$SD8o0&lJzM zFLil87vvw7&D+l|sGPt5%Jk#^czzb!FFq?(W6ND)A9c&6WB#Aw{jwqkp&<V}cR6V) z6_B|<=gTjql<g(<r<(GgWO`@X?{)kH^3VPY6B>A=jn3OQ@s-6)Ro`Ks6g+zq$eyf< z&u>d6%->)7qbU~@Zo2kG-0y!B+Z+4{a1!S(*`LKze`3R*V*Bq~&zO}ZW$ynrb3-a9 z+{|7x?+@SBvVU38xxXNLSnuo4<hruS4wQ#L*fE)N6Thg$ewUj)iv3HO_V3(NRPxzz zy?s&P;pB`R>-Qh9yn4&_dw~7Z=w*KbZ#LUkWiH&%rP8ur$tG|8)NWOKexslA)0-;o zzj`zlUiXyVuW{rj{~hTj`(*;NErOLR_up9)!28s`#QsK<v3yu?(|+j>pJU!9<m|sS zeQnlf$MyS{XKjzl*s<Qe&{W>v?t6g!7qzY1{BJhz7iIB!+@aE9?^!QvT+yRyU)NUV zIlZZJKcm<0eb+pt_n$I;vU&5X#{F;8l;VSxEA9DDtmf8#RkHs^iqZ>VgGcu5F&1^9 zT4(HE6<qHG*%Q1XWgW;rlXPSZoxTU`fBjYaas18Z{nx8^ub!gPvOl{s&xL)W>i!z_ zMzDW!#jjVr^pv*$yO^(}S}DU`vdDHV$etZcr8-aTOZJ}$*LcqQ?BV`nm8a8}Z$G_% z`7Ne>pl}P6D60VZ$IWiVP0Q~A`!D~0`62pdvpw@dN$nXbE%tZ6K8)dLSKTi%>DOeC ze;i{%!jF1N+b`lOPLzF~vA^)QB{<xAR^4LPe^p}d`1VlL@#gdPO&S^Dm(q6Fdvxd5 zRqZ)w|EuNJxoDO{_NMb*wj2ph+Q0tm&fnc_)9l$d_0GQhaGE`L;M{myFZ2Be-4<4D z`qg89{P&XYcN=u~%Q$)WJbM*te{^rOyJ&ll{q`3|J9iW&?cZOa`&vKG&|dO#_)p<~ zQuZ36H<yC!>2TWoI-2Fs{x$WJ*$#y#+24(BnBCAe&E98&QN`(p)Alp-|K;84ZEkNC z9dvNhub%xABh_TDH|W|6Ex)h6@mc<Uuey+YkUjj%BBz^QN!%Z~;q{&i90m51UiC|? zxDsJ+*Jn4mYR|#_|NfTW0QpDHQS1J}@Fe>sN=CVjZPWIz)cjC)=HWDZ(J5l0hrG@0 zk7v722l?mI{~I5#HR#$mI=k*Z(4B7|r*qOxw7q9Ps|&O5j>07Sq|?DQNo*JP@0@tH z=E;Q}`|a-R?*!Q+yproM$Ul=6SfciaC+$z}i)!g=o3{Tw&sDFh52x+-Osy2R^fli< zQ{7<}$Um)nHBO&z(B0qpzE+CUG<5&_1(E6?dnP~q#BY8j(H@kCK-lV=Z}6>dv;9)@ z8MjN7D(>GkDW&+x;zj#YTV!v~JF;+pmDOVY8skm&xd(YK_^c`4@AW`_=c5T_`{nDq z8W+p9+h+(cZqS{QW550pyJ5Xs?fyL5El2b1O6+y^Yi^PFook=9wki3hRMY-zn-yB~ zFF)Du;!;>Du5e+0)U3nZM;0%#ufE{-W!{m6_S1K6umahWHCOp+{F-um-{@;*-zSvW zpZcr*biHi5eJNAw4818i`}bT}s}HiL)S}K}=Gs{M$Wxo=Oa9K?pHM2a^K3@?{;KH? ze2c@b*t@iy{QJ6Yzda|zv~Z9;lP}-+1M-i0=2xK_<4yY~?+n_exu)EHt?uOy*Cv$h z-{*ZldzEba{-fPb4}tu{`QTD#y<4sQ=3j?DeSH~a&$j*cEQ#N_`~PoZtG_ALWY2wR z#;1FsSM5K3@ih`_Jz#(1=C&!IaJ#o>)&Y=zO1#sALH4}Zxn+;ln)3b2Z7ZUWOenLT zezoQm$eu|$v%7Vt<m|sZ>2^2Bo^ZJtH#LnK_kS?Bat0J`iuR(%&t|0Cmw8BCZ#aKw z|6Jh-E6o+x?!VWxLH>d4*8OH;f<^xow%FS*o{-17^vM2ci*f>IU+mgn-F-|hb3&(m zB7689X^UC=AIY)ls~9xcD@Dht2hN(XKiEu6;(mCZ{RiFR!%e>%_TM^Ywfowfko_JP z;`pjxeA~C<MtSMMJrDQ&?`jDG*%P>`vHHKlmi<fbtm|Z5dc;0Q`Fib)i(UJ}cD`7e zKcRE~MaN&yL@Z|6=N?EiRWWGTUr;~IBXHIP`}Ds~%dXr?+`ou@|22?3b+@i7+<l*K zzx%)PO#Apt_8D$bGZKsU*ncostN1{6tNjJeHgAxB(xvvtfb7}(??Bgri(U4`TeyFf zP3W}mFf=%=Vlm5pG4o0WkbmNicp8B03GX>xW3248KU=7GLDTOB`|6Knb~omP*uPCk z6nB4jaQ`=_@SXWbR_||}mhTR-=cmv*FOYw#_eSJ{!cFD*zSKz<yX;?Hy5v$gp>zMa zyR!vlEoSY%sb8uL@()*>4mjNM-#bYCy5PS5XhhaIkUid^+In~1=kEvQArM~s(fNP2 z*f+cX^?kO6QV;F!MZHpresRz~^{RqG;LU^kZxy<2I4x<vpV#JnuG7A0_SWa3)@!Vu zwx4D53F}h_h4=sH-o3bdYmfbj$9zt#%$)Y!Wv!C2GxGOuu&hyCakgha^UE5Q>9I-s zFSxC1;<@X+UxPu!M{u&~{%NZ(J4C-YxIf2B&o1!hL3`7QN7+tF+VAgMcy6=fzG?fv z6a`vote$3{@chQw$9skMCv0y#UB0zvzm$tm4+k@+eQW968>y@E?YBjE`>r_KV_*4J zus6ak$-b3Y;#%(R^ZO^>`(pEQ;*R~k?CUi__PB@c&;<FXbBVR%4@rCbt2>ex*zKFP z|3_R}qRQ%N`&GLbx9>kFY~SKyv=-!_lKb5AO_@3O|DRYV&Bh<Pf5oLU@*sOSj+f1v z7Mo;mU}mE~>)LtyNxL`Id|tN0es58U4k+BNOKwsF`N!H}<B~a&_V$KTtnS$En`Xbf zH2uEH>S^|uG-GVu?G?7argVHR$Uhd|Va_+0IQM7wY@PbyRfzpv4K008xV<aeq8DkG zwEt(M%B?`*NA{B@<)oeYd)j_pyoW!3&-(q>dEdQ1lD~exRm?Q2C$F{kdn8WXpTg0y zzo0+1spMPpeqWW0ba#2i{lCJZgLOPA_j~PDd$4BN*L}OLT@Ji>Fx~zevunhgxXS&% z{F{AVZY|k=)8tQ4U^ByhhPuilF_T{IGjiK|3uKRBF>C+P{Pp%?r{8C)zR|LmtZrVO z!O^mRa{7MZyl>6++dKFDogu+!KkdcENje^t_VQ7!+Gm!1wR2n+bne~z#{GI?`#*u~ ziKvO*&v2_`zuVk*cW>Q~+rPHqQS#I0?)!_RjH~&3*4s-mJ=Z^)zka{zZ=-iAZ?x?H z#4b%q=V-D2vGu5X&bQ|M(tlNzW=SyG&;J))0`gBX8|$Le%f8yZzr9&7>Q|#ZtM7!% zZ{jNLpRHKf{&H)Hy^wSM{rSv~_S^iNY+m>D^nPadccvhFj+cjVg8Va2>(9p2H(LA4 z{L4K+_Wb=esj%={v;Eb}R~NhyWZbX#YR4ySk4pQMpW4-5FZ;T0s%o`<?d$aYvwL3d z1lc3?u;&iLtrB>BB%*V0nZ9~}eMgzL(gtQrd#9xh_x#(p?9Xd?6m7L<%l<oCcmyUb zIb!co_dYL8sB8a~$vboAg?8<4GdNS+Wj@P3r26aUUuzodSzUXRzRaAkzeB;~6lY4( z{sPu5b8?y*?5}^iuw%i(eEZYeS3TE`x@6C8e(6?A_8$A&OI%L*w{Ni*wR`SrwP%a{ z+sAKZK=$-Gnia<gb=iw$HPrNlcG-V)F<uC=$M)0So4?jH?9ZHiDHUW-8(aP|achtL z_dVy;=QK6!-@DCR=6P<Y{mLHye9`7Z`=cE{F+NzeW`Dp}L+OC_E%qME?<Rr#6I!aE zHfhO`{WX?PwbO;V?1Pk=UiXG}?LS|q-q~e7Yk!r4-#?Il1g0uvf0;SKe#tYj4a~wG z`y0M;p9k6V+kEe-`3v*+FI2GDKOykY{-RCE!iI5c_CJ2{{5dGx*jG8u0Qsla%N1-- zGjoedtWekf3;!3qoEzF@-?4r5WRN|6eUH_CuW8s{S8fZor{CyaaesN9{b#}Vi$URb zVvE*?C%K{f%TH+9OYFX~-{`@!LOb3A`_q@Uq`I*!vhP@$>=$OfX#cG<{O%X*HrYSQ zGui3;q|82M=WVWbrSko<$K%_1b=vKVJ)eed^vT`-Pfm$(M?|gt2@9pkr+pgjGtJG! zkN?fFXZEx0trpAJ-<7cXyZBr+`_^FZgGEh(_D?pfSqZWya7JlNnE4|6TmP*lU9j7< zzxJ{JIfEx<_Ir;mIo7IFZg1E%m78CueLvID^tg>ax%M9?tTNsaQM<qX=FBORZ<Op` zKj+>PkUgC5k9og%*Jyt}EcHmzl_&O7E=bM0tZ>2p$Mrm2H?~Fl&Ao2k4mDr2U%oxi z31p9!@GlX+CuRG2`%YR;S1R9M+JE*rw@&;1V@B%lHu&W3uPVM(zdfRM|HH3^7v?)f z@6T`a$UgozXTM-;@4+gu412R_^WQXlxUyfj^p3<#vjh7t?^c}yvL{V0Mi}Iun16oY zaBJ>Af64Gk*?!JLFP=?PDz|T%?{=I|r`=vb{PK(qKDqWwN?!io9#L!m;dPow?!oB& z?uVuCg6#1%ewXs%T_e0cVi4cUv_9?me*SB#W_8J3+V9l+-*oe|`TLt!IFw{epTA$% zswnEH|NQ-n^4wQy^Ox)oZsp}(z+bZ8OGN*N(xc-2_nhogc{4KYRrWUTF3`)g-{f2U zWWMh8eGI<GmQ<%_+P8;4C_cFGk=^4PEmpHU&)X{~#+4qk++nZVp31pp+I;)%+Pl3o zr_Z-9JmwX2%zys=88%YJGW;d>7gCgMr}3BAv)z^QmVQ*czd!e~8*fJDe%6Y6{Dpd% z_I+782Aj8(*b91Fc2}oo?%&_iZK3zBalc2pcHjr~^Y$TH;cK?a@342tOO^)Nvo7;V zPUiIa`zL>BH$CP*-`?kiKes7=iM{ZWPcbw3OZGDx*s9AtDz<mdh~5eEPtbaknFV^8 z`x9LgQuVgC?6)+XSO&I-GtuwR!bkfKJnYvq)I4v`@$9PXI)xqfE0}E!woIG9e+%Ec zV32>zE1SUfJX5{jr_NupKYVVa;X-~;f7nO`WRG7lw;jko<x&=4djh(bGnMWyv0rY~ z@C_7h&yK7D`Nv*zPvbm|+3W3H90WeUshDBE`J8$I!~1>qE|*Wb&40gd|BTEV5})qv z+y8uqv&`BhlkJ<Y&0E2@bn<?mvMHQZvnKD)ROaxMeb}*osl4v1&6hj&>%83(?&uV1 z|Mkqq$#Wlc*zYa>vtxpNlD(Q9bH=&r$M^fpyLCXZY}0=2j7<q3dqRIb4xayhpS}L` z-D;oi?XzFCi&1OklF9Z8lb@#YE}d+Drs}6d)vU?(E6s1X%0BF{&&{9qbo1p7`<GK^ zXV}=>*k{)2QUJ22XU)6NyOW>T-AgwQ(iYude>d3h@yX{i>`(7(2w`}?Z~yOgQ#?Ta z3Ay~x5M)o>5$ScSmrUM2ZOWYh!KIV;%h$BLET1*mey-_0H;{ikt}W!-a=F7kKe==F zZlO^7*rd>+xeq$_ABs`fIl(?@|LO^jhj-0fZ||<9wDw~8412AATY^FMcw7lf0r}_G zl2~z2xE)^Ok-uum<o(-z@f{XgI@w;d%I{bCtjYUXe@SM5{Noe*@ii#i>}=z*6Xa6& zYpp0O2HErCnwtIH$xrry+H@d%@z0Cu9W&<G`U>6GG>Lj@+qOdG=AYK1_DdX+Stjl~ zvOi1yhR}me75i6C$ubbAUSOZ5Vyn9A$o&0NZyU=ry$sr~{P9}wge?{JOya$rnpwK` zm$ghpwZ4DYTXAf|5B^1!`%6356{~-HvbQ#6y2MlKPkWCBub-v&Qq->Zgspz+n<MtW ze)yk>jXS<STJNd<gH09okxqNF%2*fJFLaU;JpN?i{_I^0R!uL1?Ekc{n!M;rwSA~& zZ>d(6?tZI7CnpCz`nXqX_D@BX*A@FWthgh5MBve0ce5WyzsT6{%llBOaZDq@u17hv zD)r5g{m=X_m`~hy#6H86>sY<-wfzpZyZ4l_F4*rlX>H=JBlGR|KhR6NF1&GnaQQLO ziCZf6|8mGSTHW7euUBT?8PxR0w)?!&Qq|WL_7bPk>(}`_*~@a8h3T`NgI((q_dN}Z z684?CQpxwX_2~Yit2c_o#T~cTzms^p-uIgQYQwAk0@VxlcmDQIIR0dzy%z7;^y|VK z?N4q{=3H{6djA{gYRlFAUHemKo^;4MbKCZ1q1YjTMV0n7CTl0pko;#W^^CRdZ-}wo zxkYE~OVxt*{n(fko)~{^zx^p~(M*G~{jt^+aVrfQ_P6$=xVI~A+&_1ZOGV<}fc>F0 zRR%3*BKGeW_Bt$(tF^zBH!s@Q(|Er&kDx5q^Z&N`iS^&VE$Q07O@+I_I+fS1#EkRl zN~KAAH*V3Bz53_s-ivGaDiY(b**~3UdGyo7tNWiZwaigI-MOD*zfDNH;zoP6c{~1C zwrtq%ub_NJ@<8%_ttE4g`EBUj-*+zb-iPgW`;`o4ROCgA?@M?hqxgMEm%Z<69{&<{ z-hJnuPOQ(dezx~L+oluyI2r8(g;a!#tV;H~i^}B28I;*Sf67*-e7e*AkA8~gzn^~l z8+7iu1pW=MH=Z_siR6J~`=6G#z3$~`*{}85wf@6)JNr*v!aF6dUfZj%s{fU=xW@h! zpG`Lui_PC#7<r*Q?TFF7tt|q7!(;>P63*XLsI@Aw4_&cy>8FWT?K=<7*IH@VU_Vhj zalia8Kl?eGJwI5sY_Jztbu7C1OoY9`!uolB8#?V>8n$E@c^ccBEPr<YhK_*U+K0}r z0pc3=PKo-`kJt~u$Gv?nYD}$XJ-m0rOob&?{Y?8jx0Gqvi*2_r*r?1Da%87{&DoaE z-;z)4uV?T`=vh0_zDl!4^vSNt_5wFrghkgb+CMYqAM=&<dHYWZT(H&WskQ%^JV`Y} z^r@|9$XeC?b=mtT-0xl}|Ng^XUNOrJt5TfywHotJ@Z6eYH!ZP8xNH99{T?OHBc{B% zwBP&N9sTdgC+veezIW>LuG>GyDx;16$C~|`_xG5Ju3cpB!tlP(qP@-j*6wO01D@Lb zjT?m`9b-Q1RaZJPd3Rj9z04Gw%mvwhZJCz+?U-d_Y`5pQ#OjLdkbM&CYyWr7zicnO z@~UIVk)8W5@4R;<dP<x99k0D{`n>Dxzx39aKiM^Tzeh!kos?IieXhl%FIU#*+0XOX zQ7ZYwYQOit3g<_2KJGpB#71X#T>Ji?s~?3Kt^Kf9ebPx0ITu4am#vwVVi7_6c3o*p zv=`gHU)!MFY09fh_9~xGH-}AW+doy5v!Q40#Qj_E1c&keSYtm=TeDfjD{;S*MCd>B z_O|_(M0axidu(O@M37Z~_U~`Dwo7gQ?yt+XUtDqYp#{Sg+tW<*jUv9Q*e#u8bNH6J z)4pvh7ks+@?fCxX=TE+>$gbOeI9#D>=>sGC9h>>O+aee3uLyg${D@ch{!OZ_wh=FR z?42KHnybqi?Vq&%AOF|NXLhqpx4r9V{AxSTb>DoeuA2RePPZ^^W&UgHKINa+X~R0( z33oy+I?8{s?H2DAz47h1J^Rxsb7velWxvW``>Mo*oBQWEu+3<TTx5S|^IH|ZnkDuu z!g;PSd7Jki;ZFVFe5c<2tM2iRGZuaOpPp^Im8!~X_w2998LO@u`xDdoTMo2}?>lW) zw@hJ;>Au=}hJ!gQp>`n^N5n6;)Z2@G@M@{guCte5Thp4DaMM0FO>nxAQiQ#g{H~?v zyu$4d8qYL~$=ht-zrX92nXHk0efYY_GZuaJUsZm+-uRYjU(D9Zy4^>N?6)=+wFXrR z*u8(dBP`BM-LA$hS)_fT%f8i7-g;+S>i1{#PcNHs;M9K6cJ@!X4~_QQ9qU-CsT8sQ zcZ0teU(J&J)~~pi?R?F%-~ZvZea?64_rL8+ox7^`*}l3TOZBFUvFwXyK32K*h|&I2 z=`TAnyU)Y(5x>lVXbr`sd&@#*b$^@vZ0|&mSH{;qIPG7-!PW6v-^qR^$AR$WYYy3| zmX+J6Gs*5x4c0%DdP~}V(huKbY=`FB@o*(wd9~u}zHN!mS14PT?GrS365sr>(3b0( z?XgvT-|QB2+?_Z}J;nAYU-bLRPPe_Y^qOBjJ^OI)Lk-mi>E9RkXERHc7XP_uKaZ(N zarv4<`<etQM1nh4?`M{pa6Yknwf)msoiq0zoV)L%t^S+&*{%DpK0f@($-2yrxwY2p z>UrUPALl1~v1PT|KliT|>v_etZ}%aa1#3@kvfViM`sQY@7kfo=yg>HYn+LAi_*&m- z|B6*gi7rnr?qBA)MKri`wY@q|g>~93>HX%9N|#5bt=?~O&HCf(6<_Tx-t*Xg$*<L3 z>0xu4x^LxP?RZ^pwyf6uMGn`h7rqm-v*<piKCkug-ev37*@zn)-K*$*Cgs`(C;OxT zTabVDpM3$g$G+!!NSX$dti7{a?1jYc)%%Z^8i4GvUv2jQ<R5$IG(C_#`;WdYF+6@n z*lt=c!>-kR-}e1gmE%eM$F?tb0bdD2%vRg}PY-@6AN*unveHe_U`2qvnR(*bf-e{M zzs(SIum5z>zJimpvA9XoKI7B{G0FPX`+xY(@KmT-ZU5=7LZ)a^>;6CcO&=bKZ?$*b zZ-3l(y`tTo<97}xUt_ix6FDBFS0l2o`zc#NPebC~JzX9xDK~D}9+5fbXRsn*|Jtlm zEmFz>_OrCAQx>;z+l#(^!B*a+xqnkc@)_0b8ulh%n2t{odTN)m&+F}Db7p(SS380Q zQx4c|?kEa9wX?^z>fC~F$=8_ozue8l6|LB0>ri^zi@9u)?X!X%i8~ap*?OO9F@OB! z;{M9RlR*BlFYkUZWl<aVe!H5EKPT!}?=Sh#EhSmM+J4!6w#oHEPxm$W?gjbB{{JP; zlY%J+_RZdB^{#5G;=YuE|EuOEweBxcSFUa=5ZQO2x_b8a_w#Mn8Y<u2bLWC>cVrOA z9{a~f>cXXz1NQqgae~5azrf@_KTg!Iwl^%g7_78iW4~Y(lZQgh>izp#ICefZXWswK zR2dX*`&)`xPE@BJ$62eQZhz!#R4Q;SvL^GAjcopkJU_3W&&V`MC6MeJKS-S)-K zUl;d_c*<X$_~+t&zl95zwm-SJ-&9cg=YfvZ`wRbWlAheXdcVoab?S%GR`2&X`|Ql& z>{k1=naetr{9EmR#kc48O{m|ySnAcGiJ7hT2^<Wf@AC?6m&Vq%teGXeZ^i8=ic=3b z*r}d6GRNw()BaT_KJn%8Iqy$iAG*B#$wm7vroHE!S!MT2&OXp`OkZ~Y)F%?ghtgKt z-&>w2cy;49yZ@s05=#E9`!8jOu32p_Vi#=~wzY8LH@kNs+DD|j6!%F)ty#z<CA_au z-ptGTwS(RBlPf31e|EB$$YZ=d@y|v3@A`)|FA0CLbG#ksAjBqXpLX@6=H%|x_KN1t z8)SE!wlj&{6MZ<lbwBetovU{`m+q686J+#9S!Lgc+UC$B6Tj{IaD(B-^(|9uMeP@E zj4~45xA*S(zus>g>>}F~U+w&L(Vjzh63=8lXZyFo1%@{KpZ0w@qRM`-W3_$$yXA3{ z4P@;-#E!BR%{#r%`uN?W>YKjp6S$oGeRlm)yTuc~8EA-$*v&TH8aOetb-%jD^=98{ zMLWT7iX{`6^zF{u&U$qBzTdv%^KuP!zFypa;b!dgdCC6!mj>kWRWIeSUvYEx5z~g% z`xnGtpJvOUvA@7dVu9}Or*@2y^AG&kz_kDS%rCmv`3~C2UDCfPCh*BN$oB5nFG;QY zB~BT-d!GDn`*C}@d&oTZeJm?@FRfdfW>@!u+ga!9MSGQbXO1?0xoE$WPozZk(?xr= zAEvKN8&=y-eJ66jpnkRe4UO|%ldD(TubFXifo4*x{nsaJUPs2a+JAg2zK*Ll-}c`1 z$*aF4wc2|m6h&W{*17jcKTi`wrM_L?=5_qG!hZW^6xUC`^CrMPu{z#jS+c)9%gUG} z)lV1q2OenPa(}I9-#eLg%N`C5`^Vegr%kS2z5niI0m*+GnCzcz3pb65Z{089_Umnr z(<j>{`roJNX|UK|zp}aM@vQ%Q<^CtUkh-pCx9YvbEpJo*eUpkd9lh};VE?3BhFcoH zT-<N|_P0F4Ql9;$QY(KRf2Fzqgh$aVgZkC`d82|Zd3}Dm?`es~Yt5wA{V%@P1yu7L z+}F?a`23n!rhU`yoOjgJVA=29CSrfyeu6EiZv?`A$_B?~_ZQgis5|_6`>AVtSL}OQ zsGfa!fA<U3UHr<I?PV0+SH7HGXm7VTQ2hO#wfl7v8&<UYt+i)%NJ#Yk?qz>pPx(d7 zskZ$Gu6{09{zGNI^6A58u5RbH>z(|&FG8W+eqX|>8JW2ZcKne|2hMw(+#7bU?c<$< z6MKu7ra5rE*lzzvSi&~)+m8L-%FQ!g&Mw>^TeP&>&~K9cV|&fr+Om`Pdo5^9_Wka) zpPy^#jn-Yc_S>e-YhC_B#r`F?<O8c{c{{t$&TO+Pa`yXuib|hbAZS;ddEH$7ugboQ z-Kh@uYpv`g8kT(Ke6f9h$W*7>{K}X2`&@kW$v^$Xe)sS9HyHR$+P^!`YE!%4+WpmQ zrY0_1zi9ut<sX05oNBX|`rm1+BT#Fvp}1FF_^$H4+5QzVGb(cI`4-KL?FkUFQ~7>t zLb{Qv-Qt^e_vTEp-dCAlJy|{bvi+9hDIX)g?XVZ*O9a_tpI$mw;r*Vq_NjL#Khlz& zY=1B2R{XN{i|ixs8JV{1%H5xEJi$OmpmslpK`VpMj?a5FtA8(yP-x#@)U>xe__FH0 zjCD8uNzYT;_x5~!!JC92yPkgmvML<)_V<Gsr`x<fx!*=+7XOQ+oBI<Nihc_@9Bwbw zbKQN@jm7&#Lt3l+3O4U=X<pN$b+2y!rK%M%=dAno7ni?pKAR|D7b<UexS`6(J~&^V zzbN4N-shXQKi(6Wy;s4e;Ap(fHQNJOVc#@4>i2Io-l(zuVV(WHqH2jJ{7Uv*M*9^* z4u|je2@5=UWNEnlg}!NzKCbv}cfarCvZcvJ_FvjodR%V0YWMBxg7;4UH1_RaTG-rQ zWwig{y#=2hY?9sQ`edHX!h~bChY|&hj6c5ETlfCygN>I?>|djLeW}gsllI~)m)#!n zEA3BPrnT?^<Kq4D&a5t*ZY;L<aAcVLdd2U3%XV*3(Y{w_UtMzP$I+&%`})!gGz^YP z+6AdD4dtIuyMKzgw&GM)K0BtYrM#WxYIf0{hI?!81nzTIxOs8or4#lQi*Bc^d|0=? zfA%g=xY^5>o?h~Raj|_h_vQYhOT+hXxW`fLSFqW>rb<(Gda}{}jhT02&sq1`Kfg37 z=k>L}wo4-lHu6oVwKon*GTV~R2d|GNyI;~uOW0-mciV}ZjUS%d&QGd;_<zk&`>4Kj zVeZAp?8`IQxu5%A+s~SDAm{Cj1@>)+3#~#z7uxSgwY|4ebmM-7Bg*MLU#s@}aMfmZ zOzN_KbK1c8PbGuhn%&IBZ5Jx{_pEq--hLhTzTZ-^M}?a-><o7tv;C}Lw=aK%{GTn} zNA0<iUM>>;d&ItjzuEA)|22D=%-)>^ehcg;%lS^T5?WwybWiNbO3{t>Gct=NGySjF zze!|cUB{%Z{YvL_ABQa#+b70+FmriK<^EOqC*P<Ds@Pq9=;@*RO=922SG5agI~&-o z+&4#XtM}3U|L>f(;$D1g|JIxUKL3Q|{g37+eJ}7^us_4*n0-j-!u<;zCgnJ<vEBds zLv>xx*D8BS!v!xwg~jdH_1o<FwnKd1wY@P9m&a7vPl{k?zi>ltpLT2t>*rVE`yO+? zk~*lbZ#S>Xh@D~W(f#J<19uAlJ+l9H#d^u8gk*c>IUz-FXDrzNW?#T{E1?DZ3njw; zHLkI>cZv@ZVE$iWf2aD`!%4#8`x&&o#1@z<+i~B??`*qJX@BI^E+#(V54Hy3C*5m# zrR?&UCpLZccG;(`@h^qza>@SF-RtWPT)Ju>CRe)bq*uv)gJ;u=o1ggae|TQx*zK(w z>}T40f3Xtw-+#Aensv?JPWyR0Tqk#EY3!dHyDRu`v54I=1&bwL&S>m+p7^6Mc%{HT zx4;QCsrRq%Wni2V=<NG_?=OwXTX-***c;TWQfi!8YX4SewWX+4!~W2aD=N)T{O!MH z2a8#D1=!!;s2=~}X2kyLXH(|J?$op|d;i!wDA;&E|FubW`!>t&YZWTdU_Gm0@Ae@l zL+GD|-4^biZK6Lk_oZgft9`mD#Lm|4*UV_;Yx`wq7*`#*balVqGxjHfRt@&;GdUBd zeBQ8Mr2Wp(+gms6KhO00^T(SJ_V*U2ao7Cq+%L1~Y>Z>Dv3>j#_tN--pKOnrKDDWi z?b?5|;TZ>W-rv1hVYg1pos_bBT(McVa+d2pj+bupqLr`NZ@VkAu6AbW{>XcdEyuk| z?B~Y?g6y&1@9@UivMXT!1Fw8OTVa2DrS;u$K07t{|4#p4w@^#No@LSc(5x9BZFwK@ zWLL*_+3)69v+BztN%(l(pUCAet&6kwnr?8@*<QlF@8_(WU(z0*w(n(PeJ8MEv%P=O z^sd-Dm-b)oiP`$@(3$<QQ&;>{2-~vXzA@%iXYplwe)n%}@6%fMuf4|CTwSrmzO7O5 zg(<tlK9jaR3ddaL?bm!@r}65tnH}G@lwH>{()NioH`p-9bl5puU$_%w4^P*pvmkrc zJv@;KvL}8r-=;^0&g^GDtaOblY|H+?+diFZDZXs)y*e!iWKYSvZP_4uw9IAoZ}cnf z%j$E>0NFDo#CNY>fVf>5zqeabmAjqxEzQ!8uQK<gaz9=PvS+W9>Mf8xIR!q^AbU>S zke>MI(3$<-hUaX=!nW+MveMt$QGD4xW1dbt$exCWA=g3ntZb0%{r7=+U%5cePLMs% zJ}+Dv8(_GvyHoJyhNc9&9X7MKeY(}SZ)4&6!ytRS3rhBb>^W6&CJtoJPrk$6&kvp1 zFPJ>}zIND_{r+JG8(NDm+nc>k(*2Ouy1#o}>T!@gg*olI=dbJ7O*y>o-%*!&`?Y)% zR(n=l+}o~j{$Ih-gnbGz)gL}*HrXw)Kk+>6@#+1h`c6j#c5L21EqZ!w?43*YCnRqj zdVJ`Nz3}|QijrYl?4vAJXtft#-tWV9*AnEP<31a!K=zc$91ji<G}!mBv8f5<pIJ&l zG8cA;+ud(5Gpll)ZEL_>ZQ4?-xX*KsL><VU0(Yi)AbVC$-)9c8XVHgiMIR5Hv5$Pb zEmbjWi~W^5O)GnfFYlk|opuXkPgLypn;?5S51B_vItkecfBQ8DWY1y_OQjn}4EIf* zxZ%^14QY0F*SU4xNbK0hbf@JV$e!MccPxTCHt$~~)X@U6NBPaaWiJn%vF9s4B`6oR z#h&B$)`rI7%lq}RcC7;0b8^?pLm+$j5|aa7K9{r0GIEgv*~86XAv|&6y}gE)uC2d5 z(S6_Cz?4ItR#|o>Hs?-(?799p<2}ef9-RK=AbSMd^}XI4I%Ch&_{mi?Y>RyjgHk;x z+_=Lw901vKBKk%N$ev%^F)k-(y|HbT{l*BgM`54-^o1ej@cQUTs;_v{*J*oMTNZwd zI;gcz^hT9*OxfxEBE=f_c742Lzxc-dq|axz?7!FNH@h<EjQysKMe`?hUf%!XsNTQT z#ar#2IQOqw-M?i2vK1@76;E!n|Ka?Nf9@o+eZdVcmgJn9YabiCDd&Ed_`bqgu_N(E zJ@!4H(YwJ?DAR6|+lQ;sWvBN?SIpZ8vggFaBs-8je-(vSR0o~0e>F?VenRKv{p_}` z3CoJN+BZr)`VF#2Y|h&|AbY+w@13^qgO**5M7l-pxw-boUmEjV4pX&T?7V5Us&|^* z^vN5S@6_zr$8qk-O^`kNjlS#!*|T<)xCh9dx8IVO3xdwr?~G6|@9n(2-|}hPrHRE` z?T>wKp9HceDbV;9$eskIg^oI}mG-H!Jp2Um&&|Gl65>Y%?cVwpsOP#Q?2Gzfe(e%Z zquoop_b);A==r=@53=WKtfLsno-6CzFBS!zv1dvDDBRz9d4H3t)2IE#TkYAWsec05 zQy1%cpk#8JeY#j+_1Rytc00|Fwt(!JIpyiH7q%Mvoc7JTHFaKsoxsZ+hSXV&`$8^V z5sE20Z9gSC{M@dOm-aj7JW~05W{drbb4S~AgU;+XpA>R&Lg!`snlPu6tBSYof1;?i z6l9Ov;i~;0{~VdJZ_9>vTKh7t=Q4uAP3lEcap7|{JMqjFA`f#P*?u;9YpNyx#+FUg zWl?n5X?x}ft{#v*-}Z#Ef$WJb<Kn6dI<tT2pD3xRotN#iL#qnb7jNDF@#7VSHT_HW ztGB%h2l>aebIDF$W+OYv_o6Bwd*Y>{><(lJ+37x<WGNEhX{YB<lpm3mwQugW$3H># zgo&4e{UdfJH5Oz~iByGuNzj@75Bcmq^mSgguVqScTvWVuf84&<nIQiZC+@Qb*^?#q z|Cj1q<9+J__SAvGEkplrvBECReLAb3-4MT;wlAN#N%llxr`_MssXU-?`?x#n2FO2c zElM>Yd%mBZ@Vz|f%>He^12wxkFWax#7*?>nc<X+JM4eY4dyF|_BtZ69vmdQkZzo_E zw61k6$ev{f{o+rmO2F%*TN-nH%)h7Z6PMIccl*?C$EOqa<^P(~`zNbMoU&bcdjGya zuFdZ^oZg?vx!>`A*BN_$<%tLXww<xR%eqznW8WG3x4SP#FJql&FMFS{fq`$Hy>j%8 zAl9#|ZQTznS+k3Cp8ckSPECGJyZ1VZ27gRmA7-~r=+Ac%p~`&?g;CEz_VlVv-vhE| z?t9+9AbV`hio(x#ow5I{w?Xz@+Zp>`N<vG{^_{WrxodF&WRJnhY(J1aI~e(#p3A80 z`xIxPu$yz9z1N@W1sa)#c7p3}eL36{X17W_?`F{SihW}1PCNqHQ+(>(e~>+ijCWRp z?732Vko98M8GHRzd#*idJ7Yg3{ZiZgzBBfE#@Y@bd)_%mzX93P_VVWeWepR%&7FLz zAbU~+?;PIucZY5KKK)Hyondz6>g}y+Q!Do6nIsx8tUbNo_Lk^okUiN+AHeo(W?KCI zVAmOYflr3ZKDM2)w_onw0J3M+&xBhbdyYDLgYD^cl5^0DRM}U!TC@ukZXExbKG|m( z?n`Lj*kEwaabMNd=wqH!YwQ?;Rv-Jn=Cu6{Z;xlTD^J^BjQxG({f5)_y-&4|-s?KE zf1j-X(sylV_WR7aAoZZ{%znR<9KIlX_O&j|0{Q2(Iol_b->P;~mzWuX{L?02Cik|> z$gado$IkC~uARKdlND>rI`(aR6wJi1_O!hc-<j(md-_%tp9b0E7x?DH!LBp=?|)cQ z_@M2~ey5#+Ru}rt?0;V%dH`h4nFrslgX|Hz)2hF1+u^;6Tf+oE_ISN5jhr2L!8ZRx za@?y-r+xpb+?LP1P;2)q)`}fu&w(|2?||&7QDyE1*<)Zbv*dEunf<=!_ltdMJG1|3 zVB)uTeP{N+cD@(}vd3lprN1D1wrsJx=#ZsqcauYAIVju?-MD(Jufu3xkb`f(!2`#A z%9qwu^Yzx)ZJWDS3}jD5rTa9HJ*HhduY&9;zSVyFY}c9nE3c(@{0I4Gb>x-%eP{Mh qoAU7)$R0C~JDDJRc5VB{I<MPw-?<sa(?RxlJq~~6uYL|Z?hODSCVBDz diff --git a/test/examples/Si_bulk8-RESTART.wfn.bak-1 b/test/examples/Si_bulk8-RESTART.wfn.bak-1 deleted file mode 100644 index 51352fc1f4e1c3dd50ec6357065b7910a8139073..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 14132 zcmWe(U|`^2U|?WmU|`5#U|?WkU|?WjU|;}=D?s@mF<cm=#sR7aWClntIW))&kRw2P zL716=fdPbB8BhU84unB+AUxdUL17@kz`y{*f(#4{pfCpUL2@Tr&TjD+Ilo_~{?A26 zzwP^NguWjK(UuA?!Ssgh@ryzH)lVY9^o8(aGeP`?XY==1IBmC&cB}%^^4C3%g7~Ek z_Lo7lSFF?>5S_YC1x%lqSGMZ<qLcPtch`gInp5&8L40UzfTN3nfr0)sBLjm2GXn$v zB<XgZcNgrVvJC=+RTtV<HZ`d8WlP%``CW~vnI~n}qc^3aB1qaUYwZH-|Nn*dnKIu$ zIosK6U*PR0=fcy3_xV>lG)#Z_Y44K19MTQ7pZ0b<V$3-?t8DMWiT+o;OFr#Adf-Oh zmy~APhW6<jLH4vCb=xSSy3qc^L{X+l9+Q2ieSWreIGF5<P7Jl-m5{dE!)bl<k%_Gx z_s?TeApz=ke{Q~D-sz;cFX{Z%3-0SB?6_xbd+&Af(_W{g_D}WCUbcO~aYmu^jD+37 zE7#@csXg3l@^TT_9`1@$8$k9n?+E3OK54RV(!x_&o%f{e+<!DT=4P4fV?Mg3INrx) zpGhU_(U*r5?UXDT#HWRc?UM{!YkEKC)7}logMM%Aldzj}UC=(^*W<myLT4*<6C~^+ z)Uudm4>#?7{upA<pL-|4_MD8|Ss`yAZFlY4md`!&O!l3ucjJ!8HrW@El>a-`M{(a$ zy*<7c*v<DzOKf{49xuMHPXA`+lu8M^m#1!DIngI!_v+nMnb#FpY_A39PV71GX>Uj2 zffoUq4{d$;{O5qeZRYO^Mv#9lYkqNkx75V$YSM&_MOD)KB%Xa!=9(_G&x|#gRsX!8 zUH1R)vY(@r_XVW>&e*4{U>7JN(lguk<6efpy0i3ue%#BqW#4{l&Wn4uxwTxllO$nx zKPo><Y2Mwvm&3V_gY4m3@;pXVb)h{Iqxa{a78ASB8*>({-)LetuOyk-OVY%S!AdE2 zpSgmaocN=gK2er-O$Q?`Z(U-wFM65Y<2iB?c6x0Sb~-dl*sV}^*(<a`&Gxy<?&I6c zCF~l#J-%k&Q?Whw#~bXQiPrUC|A-t~-Tbv#YTt@Bud?SCB=<f4<Y0Q9$J8#tfkSZW z6?wbQ^Dn+%U1DjsVaoRx`?fpVIj&>pU#;<JuWQ#b*1hHucF&I~yuT!Naj#kK(Ti;R zKkgOMljgrZ>#psr9!R*|T~Q6TXQt0|jqXONeI<LcDh$(2?1H7zAD);lwXZ;6#w^uL z1-syt4QnJKwf3cK^Am6tRIuCo*1>B^lZ2hmQ=uM~pC9)=kWxJ{N93WcTE<ZZ+k+qX zw)Ag!5WmBBFDMUzuzQ~P<#X4M+rLf<=2t4)WWVvmr;h3u`}Ti0_;t^WSNrzQd3keR z)2n^^r?3Za>{>8+zeq(;^}~ge_vdFz-@LhS^8QWU{wsv;cG#CD<-BLO+hN~Uas0|8 z`-^r8^$Z#Dw>#_|lRrq^5Lj)uE_xTs&8zG8pKOjzj$1Zkzu*3nqt!3=*@p^;{F?D< zpZ$+zywc6D_U-?$A?APkg30!|Og(H57EZQ5aco2O@r9H3A34_cOz3XM{=FX!oEYzR z*eAB?cb-xY-7k6}Y<c|cj{Qp|wx>^ZOtM$F-*oNp^7Z?DUZlNQ-!x<Yrl;{oK=#Zn zlvSGfYTy2vCj#G_UhT6l&-hu?wqUaT*$or--d;F)zn^W-iwg@U+qcY%{Q&aM)+`+p z#=9N+H?J|?VkVnvzh~aE74f$__RAPET)ix?df%)hWAk|r*YAJP5?ouoe#ZV+Z|}YV zg<I^oqd!3Yc`q#iwnv=lw0`q~$@^QmgO1!?IN3g~a`!usJ)9=BFF^ib=dl9YbI`-4 z=#g^hei^N)8$scg-g*_}pZ(Lr#k=Fp&fD)OUg7#!Ylr<8e;<i6#q;;8N2&CbPMB|h zbH&HJg_ZN|_goI$`RqsW{<DnoQ&qT2?C0?WKiTxB*nVfIbXHkXrhT!)`DCepOnY7P zZR#A`O6+9|Opj%HXYOZD|L|?e`$l_jCWXQcY0vk+l;srfk-M~i+Y-jdAbT!dPTNyD zVgCLt5&tGEten4p<E)xsiJ!&$4?R4#!<V~cfBw^(iO2sG@3*?rTUefyx!+B?xK=tK zbANNFQ~mU*%lGl7)v;uIXYM~6dGIUul}-D&IQy#Hjn3PzUzK-wyZR1$SD8o0&lJzM zFLil87vvw7&D+l|sGPt5%Jk#^czzb!FFq?(W6ND)A9c&6WB#Aw{jwqkp&<V}cR6V) z6_B|<=gTjql<g(<r<(GgWO`@X?{)kH^3VPY6B>A=jn3OQ@s-6)Ro`Ks6g+zq$eyf< z&u>d6%->)7qbU~@Zo2kG-0y!B+Z+4{a1!S(*`LKze`3R*V*Bq~&zO}ZW$ynrb3-a9 z+{|7x?+@SBvVU38xxXNLSnuo4<hruS4wQ#L*fE)N6Thg$ewUj)iv3HO_V3(NRPxzz zy?s&P;pB`R>-Qh9yn4&_dw~7Z=w*KbZ#LUkWiH&%rP8ur$tG|8)NWOKexslA)0-;o zzj`zlUiXyVuW{rj{~hTj`(*;NErOLR_up9)!28s`#QsK<v3yu?(|+j>pJU!9<m|sS zeQnlf$MyS{XKjzl*s<Qe&{W>v?t6g!7qzY1{BJhz7iIB!+@aE9?^!QvT+yRyU)NUV zIlZZJKcm<0eb+pt_n$I;vU&5X#{F;8l;VSxEA9DDtmf8#RkHs^iqZ>VgGcu5F&1^9 zT4(HE6<qHG*%Q1XWgW;rlXPSZoxTU`fBjYaas18Z{nx8^ub!gPvOl{s&xL)W>i!z_ zMzDW!#jjVr^pv*$yO^(}S}DU`vdDHV$etZcr8-aTOZJ}$*LcqQ?BV`nm8a8}Z$G_% z`7Ne>pl}P6D60VZ$IWiVP0Q~A`!D~0`62pdvpw@dN$nXbE%tZ6K8)dLSKTi%>DOeC ze;i{%!jF1N+b`lOPLzF~vA^)QB{<xAR^4LPe^p}d`1VlL@#gdPO&S^Dm(q6Fdvxd5 zRqZ)w|EuNJxoDO{_NMb*wj2ph+Q0tm&fnc_)9l$d_0GQhaGE`L;M{myFZ2Be-4<4D z`qg89{P&XYcN=u~%Q$)WJbM*te{^rOyJ&ll{q`3|J9iW&?cZOa`&vKG&|dO#_)p<~ zQuZ36H<yC!>2TWoI-2Fs{x$WJ*$#y#+24(BnBCAe&E98&QN`(p)Alp-|K;84ZEkNC z9dvNhub%xABh_TDH|W|6Ex)h6@mc<Uuey+YkUjj%BBz^QN!%Z~;q{&i90m51UiC|? zxDsJ+*Jn4mYR|#_|NfTW0QpDHQS1J}@Fe>sN=CVjZPWIz)cjC)=HWDZ(J5l0hrG@0 zk7v722l?mI{~I5#HR#$mI=k*Z(4B7|r*qOxw7q9Ps|&O5j>07Sq|?DQNo*JP@0@tH z=E;Q}`|a-R?*!Q+yproM$Ul=6SfciaC+$z}i)!g=o3{Tw&sDFh52x+-Osy2R^fli< zQ{7<}$Um)nHBO&z(B0qpzE+CUG<5&_1(E6?dnP~q#BY8j(H@kCK-lV=Z}6>dv;9)@ z8MjN7D(>GkDW&+x;zj#YTV!v~JF;+pmDOVY8skm&xd(YK_^c`4@AW`_=c5T_`{nDq z8W+p9+h+(cZqS{QW550pyJ5Xs?fyL5El2b1O6+y^Yi^PFook=9wki3hRMY-zn-yB~ zFF)Du;!;>Du5e+0)U3nZM;0%#ufE{-W!{m6_S1K6umahWHCOp+{F-um-{@;*-zSvW zpZcr*biHi5eJNAw4818i`}bT}s}HiL)S}K}=Gs{M$Wxo=Oa9K?pHM2a^K3@?{;KH? ze2c@b*t@iy{QJ6Yzda|zv~Z9;lP}-+1M-i0=2xK_<4yY~?+n_exu)EHt?uOy*Cv$h z-{*ZldzEba{-fPb4}tu{`QTD#y<4sQ=3j?DeSH~a&$j*cEQ#N_`~PoZtG_ALWY2wR z#;1FsSM5K3@ih`_Jz#(1=C&!IaJ#o>)&Y=zO1#sALH4}Zxn+;ln)3b2Z7ZUWOenLT zezoQm$eu|$v%7Vt<m|sZ>2^2Bo^ZJtH#LnK_kS?Bat0J`iuR(%&t|0Cmw8BCZ#aKw z|6Jh-E6o+x?!VWxLH>d4*8OH;f<^xow%FS*o{-17^vM2ci*f>IU+mgn-F-|hb3&(m zB7689X^UC=AIY)ls~9xcD@Dht2hN(XKiEu6;(mCZ{RiFR!%e>%_TM^Ywfowfko_JP z;`pjxeA~C<MtSMMJrDQ&?`jDG*%P>`vHHKlmi<fbtm|Z5dc;0Q`Fib)i(UJ}cD`7e zKcRE~MaN&yL@Z|6=N?EiRWWGTUr;~IBXHIP`}Ds~%dXr?+`ou@|22?3b+@i7+<l*K zzx%)PO#Apt_8D$bGZKsU*ncostN1{6tNjJeHgAxB(xvvtfb7}(??Bgri(U4`TeyFf zP3W}mFf=%=Vlm5pG4o0WkbmNicp8B03GX>xW3248KU=7GLDTOB`|6Knb~omP*uPCk z6nB4jaQ`=_@SXWbR_||}mhTR-=cmv*FOYw#_eSJ{!cFD*zSKz<yX;?Hy5v$gp>zMa zyR!vlEoSY%sb8uL@()*>4mjNM-#bYCy5PS5XhhaIkUid^+In~1=kEvQArM~s(fNP2 z*f+cX^?kO6QV;F!MZHpresRz~^{RqG;LU^kZxy<2I4x<vpV#JnuG7A0_SWa3)@!Vu zwx4D53F}h_h4=sH-o3bdYmfbj$9zt#%$)Y!Wv!C2GxGOuu&hyCakgha^UE5Q>9I-s zFSxC1;<@X+UxPu!M{u&~{%NZ(J4C-YxIf2B&o1!hL3`7QN7+tF+VAgMcy6=fzG?fv z6a`vote$3{@chQw$9skMCv0y#UB0zvzm$tm4+k@+eQW968>y@E?YBjE`>r_KV_*4J zus6ak$-b3Y;#%(R^ZO^>`(pEQ;*R~k?CUi__PB@c&;<FXbBVR%4@rCbt2>ex*zKFP z|3_R}qRQ%N`&GLbx9>kFY~SKyv=-!_lKb5AO_@3O|DRYV&Bh<Pf5oLU@*sOSj+f1v z7Mo;mU}mE~>)LtyNxL`Id|tN0es58U4k+BNOKwsF`N!H}<B~a&_V$KTtnS$En`Xbf zH2uEH>S^|uG-GVu?G?7argVHR$Uhd|Va_+0IQM7wY@PbyRfzpv4K008xV<aeq8DkG zwEt(M%B?`*NA{B@<)oeYd)j_pyoW!3&-(q>dEdQ1lD~exRm?Q2C$F{kdn8WXpTg0y zzo0+1spMPpeqWW0ba#2i{lCJZgLOPA_j~PDd$4BN*L}OLT@Ji>Fx~zevunhgxXS&% z{F{AVZY|k=)8tQ4U^ByhhPuilF_T{IGjiK|3uKRBF>C+P{Pp%?r{8C)zR|LmtZrVO z!O^mRa{7MZyl>6++dKFDogu+!KkdcENje^t_VQ7!+Gm!1wR2n+bne~z#{GI?`#*u~ ziKvO*&v2_`zuVk*cW>Q~+rPHqQS#I0?)!_RjH~&3*4s-mJ=Z^)zka{zZ=-iAZ?x?H z#4b%q=V-D2vGu5X&bQ|M(tlNzW=SyG&;J))0`gBX8|$Le%f8yZzr9&7>Q|#ZtM7!% zZ{jNLpRHKf{&H)Hy^wSM{rSv~_S^iNY+m>D^nPadccvhFj+cjVg8Va2>(9p2H(LA4 z{L4K+_Wb=esj%={v;Eb}R~NhyWZbX#YR4ySk4pQMpW4-5FZ;T0s%o`<?d$aYvwL3d z1lc3?u;&iLtrB>BB%*V0nZ9~}eMgzL(gtQrd#9xh_x#(p?9Xd?6m7L<%l<oCcmyUb zIb!co_dYL8sB8a~$vboAg?8<4GdNS+Wj@P3r26aUUuzodSzUXRzRaAkzeB;~6lY4( z{sPu5b8?y*?5}^iuw%i(eEZYeS3TE`x@6C8e(6?A_8$A&OI%L*w{Ni*wR`SrwP%a{ z+sAKZK=$-Gnia<gb=iw$HPrNlcG-V)F<uC=$M)0So4?jH?9ZHiDHUW-8(aP|achtL z_dVy;=QK6!-@DCR=6P<Y{mLHye9`7Z`=cE{F+NzeW`Dp}L+OC_E%qME?<Rr#6I!aE zHfhO`{WX?PwbO;V?1Pk=UiXG}?LS|q-q~e7Yk!r4-#?Il1g0uvf0;SKe#tYj4a~wG z`y0M;p9k6V+kEe-`3v*+FI2GDKOykY{-RCE!iI5c_CJ2{{5dGx*jG8u0Qsla%N1-- zGjoedtWekf3;!3qoEzF@-?4r5WRN|6eUH_CuW8s{S8fZor{CyaaesN9{b#}Vi$URb zVvE*?C%K{f%TH+9OYFX~-{`@!LOb3A`_q@Uq`I*!vhP@$>=$OfX#cG<{O%X*HrYSQ zGui3;q|82M=WVWbrSko<$K%_1b=vKVJ)eed^vT`-Pfm$(M?|gt2@9pkr+pgjGtJG! zkN?fFXZEx0trpAJ-<7cXyZBr+`_^FZgGEh(_D?pfSqZWya7JlNnE4|6TmP*lU9j7< zzxJ{JIfEx<_Ir;mIo7IFZg1E%m78CueLvID^tg>ax%M9?tTNsaQM<qX=FBORZ<Op` zKj+>PkUgC5k9og%*Jyt}EcHmzl_&O7E=bM0tZ>2p$Mrm2H?~Fl&Ao2k4mDr2U%oxi z31p9!@GlX+CuRG2`%YR;S1R9M+JE*rw@&;1V@B%lHu&W3uPVM(zdfRM|HH3^7v?)f z@6T`a$UgozXTM-;@4+gu412R_^WQXlxUyfj^p3<#vjh7t?^c}yvL{V0Mi}Iun16oY zaBJ>Af64Gk*?!JLFP=?PDz|T%?{=I|r`=vb{PK(qKDqWwN?!io9#L!m;dPow?!oB& z?uVuCg6#1%ewXs%T_e0cVi4cUv_9?me*SB#W_8J3+V9l+-*oe|`TLt!IFw{epTA$% zswnEH|NQ-n^4wQy^Ox)oZsp}(z+bZ8OGN*N(xc-2_nhogc{4KYRrWUTF3`)g-{f2U zWWMh8eGI<GmQ<%_+P8;4C_cFGk=^4PEmpHU&)X{~#+4qk++nZVp31pp+I;)%+Pl3o zr_Z-9JmwX2%zys=88%YJGW;d>7gCgMr}3BAv)z^QmVQ*czd!e~8*fJDe%6Y6{Dpd% z_I+782Aj8(*b91Fc2}oo?%&_iZK3zBalc2pcHjr~^Y$TH;cK?a@342tOO^)Nvo7;V zPUiIa`zL>BH$CP*-`?kiKes7=iM{ZWPcbw3OZGDx*s9AtDz<mdh~5eEPtbaknFV^8 z`x9LgQuVgC?6)+XSO&I-GtuwR!bkfKJnYvq)I4v`@$9PXI)xqfE0}E!woIG9e+%Ec zV32>zE1SUfJX5{jr_NupKYVVa;X-~;f7nO`WRG7lw;jko<x&=4djh(bGnMWyv0rY~ z@C_7h&yK7D`Nv*zPvbm|+3W3H90WeUshDBE`J8$I!~1>qE|*Wb&40gd|BTEV5})qv z+y8uqv&`BhlkJ<Y&0E2@bn<?mvMHQZvnKD)ROaxMeb}*osl4v1&6hj&>%83(?&uV1 z|Mkqq$#Wlc*zYa>vtxpNlD(Q9bH=&r$M^fpyLCXZY}0=2j7<q3dqRIb4xayhpS}L` z-D;oi?XzFCi&1OklF9Z8lb@#YE}d+Drs}6d)vU?(E6s1X%0BF{&&{9qbo1p7`<GK^ zXV}=>*k{)2QUJ22XU)6NyOW>T-AgwQ(iYude>d3h@yX{i>`(7(2w`}?Z~yOgQ#?Ta z3Ay~x5M)o>5$ScSmrUM2ZOWYh!KIV;%h$BLET1*mey-_0H;{ikt}W!-a=F7kKe==F zZlO^7*rd>+xeq$_ABs`fIl(?@|LO^jhj-0fZ||<9wDw~8412AATY^FMcw7lf0r}_G zl2~z2xE)^Ok-uum<o(-z@f{XgI@w;d%I{bCtjYUXe@SM5{Noe*@ii#i>}=z*6Xa6& zYpp0O2HErCnwtIH$xrry+H@d%@z0Cu9W&<G`U>6GG>Lj@+qOdG=AYK1_DdX+Stjl~ zvOi1yhR}me75i6C$ubbAUSOZ5Vyn9A$o&0NZyU=ry$sr~{P9}wge?{JOya$rnpwK` zm$ghpwZ4DYTXAf|5B^1!`%6356{~-HvbQ#6y2MlKPkWCBub-v&Qq->Zgspz+n<MtW ze)yk>jXS<STJNd<gH09okxqNF%2*fJFLaU;JpN?i{_I^0R!uL1?Ekc{n!M;rwSA~& zZ>d(6?tZI7CnpCz`nXqX_D@BX*A@FWthgh5MBve0ce5WyzsT6{%llBOaZDq@u17hv zD)r5g{m=X_m`~hy#6H86>sY<-wfzpZyZ4l_F4*rlX>H=JBlGR|KhR6NF1&GnaQQLO ziCZf6|8mGSTHW7euUBT?8PxR0w)?!&Qq|WL_7bPk>(}`_*~@a8h3T`NgI((q_dN}Z z684?CQpxwX_2~Yit2c_o#T~cTzms^p-uIgQYQwAk0@VxlcmDQIIR0dzy%z7;^y|VK z?N4q{=3H{6djA{gYRlFAUHemKo^;4MbKCZ1q1YjTMV0n7CTl0pko;#W^^CRdZ-}wo zxkYE~OVxt*{n(fko)~{^zx^p~(M*G~{jt^+aVrfQ_P6$=xVI~A+&_1ZOGV<}fc>F0 zRR%3*BKGeW_Bt$(tF^zBH!s@Q(|Er&kDx5q^Z&N`iS^&VE$Q07O@+I_I+fS1#EkRl zN~KAAH*V3Bz53_s-ivGaDiY(b**~3UdGyo7tNWiZwaigI-MOD*zfDNH;zoP6c{~1C zwrtq%ub_NJ@<8%_ttE4g`EBUj-*+zb-iPgW`;`o4ROCgA?@M?hqxgMEm%Z<69{&<{ z-hJnuPOQ(dezx~L+oluyI2r8(g;a!#tV;H~i^}B28I;*Sf67*-e7e*AkA8~gzn^~l z8+7iu1pW=MH=Z_siR6J~`=6G#z3$~`*{}85wf@6)JNr*v!aF6dUfZj%s{fU=xW@h! zpG`Lui_PC#7<r*Q?TFF7tt|q7!(;>P63*XLsI@Aw4_&cy>8FWT?K=<7*IH@VU_Vhj zalia8Kl?eGJwI5sY_Jztbu7C1OoY9`!uolB8#?V>8n$E@c^ccBEPr<YhK_*U+K0}r z0pc3=PKo-`kJt~u$Gv?nYD}$XJ-m0rOob&?{Y?8jx0Gqvi*2_r*r?1Da%87{&DoaE z-;z)4uV?T`=vh0_zDl!4^vSNt_5wFrghkgb+CMYqAM=&<dHYWZT(H&WskQ%^JV`Y} z^r@|9$XeC?b=mtT-0xl}|Ng^XUNOrJt5TfywHotJ@Z6eYH!ZP8xNH99{T?OHBc{B% zwBP&N9sTdgC+veezIW>LuG>GyDx;16$C~|`_xG5Ju3cpB!tlP(qP@-j*6wO01D@Lb zjT?m`9b-Q1RaZJPd3Rj9z04Gw%mvwhZJCz+?U-d_Y`5pQ#OjLdkbM&CYyWr7zicnO z@~UIVk)8W5@4R;<dP<x99k0D{`n>Dxzx39aKiM^Tzeh!kos?IieXhl%FIU#*+0XOX zQ7ZYwYQOit3g<_2KJGpB#71X#T>Ji?s~?3Kt^Kf9ebPx0ITu4am#vwVVi7_6c3o*p zv=`gHU)!MFY09fh_9~xGH-}AW+doy5v!Q40#Qj_E1c&keSYtm=TeDfjD{;S*MCd>B z_O|_(M0axidu(O@M37Z~_U~`Dwo7gQ?yt+XUtDqYp#{Sg+tW<*jUv9Q*e#u8bNH6J z)4pvh7ks+@?fCxX=TE+>$gbOeI9#D>=>sGC9h>>O+aee3uLyg${D@ch{!OZ_wh=FR z?42KHnybqi?Vq&%AOF|NXLhqpx4r9V{AxSTb>DoeuA2RePPZ^^W&UgHKINa+X~R0( z33oy+I?8{s?H2DAz47h1J^Rxsb7velWxvW``>Mo*oBQWEu+3<TTx5S|^IH|ZnkDuu z!g;PSd7Jki;ZFVFe5c<2tM2iRGZuaOpPp^Im8!~X_w2998LO@u`xDdoTMo2}?>lW) zw@hJ;>Au=}hJ!gQp>`n^N5n6;)Z2@G@M@{guCte5Thp4DaMM0FO>nxAQiQ#g{H~?v zyu$4d8qYL~$=ht-zrX92nXHk0efYY_GZuaJUsZm+-uRYjU(D9Zy4^>N?6)=+wFXrR z*u8(dBP`BM-LA$hS)_fT%f8i7-g;+S>i1{#PcNHs;M9K6cJ@!X4~_QQ9qU-CsT8sQ zcZ0teU(J&J)~~pi?R?F%-~ZvZea?64_rL8+ox7^`*}l3TOZBFUvFwXyK32K*h|&I2 z=`TAnyU)Y(5x>lVXbr`sd&@#*b$^@vZ0|&mSH{;qIPG7-!PW6v-^qR^$AR$WYYy3| zmX+J6Gs*5x4c0%DdP~}V(huKbY=`FB@o*(wd9~u}zHN!mS14PT?GrS365sr>(3b0( z?XgvT-|QB2+?_Z}J;nAYU-bLRPPe_Y^qOBjJ^OI)Lk-mi>E9RkXERHc7XP_uKaZ(N zarv4<`<etQM1nh4?`M{pa6Yknwf)msoiq0zoV)L%t^S+&*{%DpK0f@($-2yrxwY2p z>UrUPALl1~v1PT|KliT|>v_etZ}%aa1#3@kvfViM`sQY@7kfo=yg>HYn+LAi_*&m- z|B6*gi7rnr?qBA)MKri`wY@q|g>~93>HX%9N|#5bt=?~O&HCf(6<_Tx-t*Xg$*<L3 z>0xu4x^LxP?RZ^pwyf6uMGn`h7rqm-v*<piKCkug-ev37*@zn)-K*$*Cgs`(C;OxT zTabVDpM3$g$G+!!NSX$dti7{a?1jYc)%%Z^8i4GvUv2jQ<R5$IG(C_#`;WdYF+6@n z*lt=c!>-kR-}e1gmE%eM$F?tb0bdD2%vRg}PY-@6AN*unveHe_U`2qvnR(*bf-e{M zzs(SIum5z>zJimpvA9XoKI7B{G0FPX`+xY(@KmT-ZU5=7LZ)a^>;6CcO&=bKZ?$*b zZ-3l(y`tTo<97}xUt_ix6FDBFS0l2o`zc#NPebC~JzX9xDK~D}9+5fbXRsn*|Jtlm zEmFz>_OrCAQx>;z+l#(^!B*a+xqnkc@)_0b8ulh%n2t{odTN)m&+F}Db7p(SS380Q zQx4c|?kEa9wX?^z>fC~F$=8_ozue8l6|LB0>ri^zi@9u)?X!X%i8~ap*?OO9F@OB! z;{M9RlR*BlFYkUZWl<aVe!H5EKPT!}?=Sh#EhSmM+J4!6w#oHEPxm$W?gjbB{{JP; zlY%J+_RZdB^{#5G;=YuE|EuOEweBxcSFUa=5ZQO2x_b8a_w#Mn8Y<u2bLWC>cVrOA z9{a~f>cXXz1NQqgae~5azrf@_KTg!Iwl^%g7_78iW4~Y(lZQgh>izp#ICefZXWswK zR2dX*`&)`xPE@BJ$62eQZhz!#R4Q;SvL^GAjcopkJU_3W&&V`MC6MeJKS-S)-K zUl;d_c*<X$_~+t&zl95zwm-SJ-&9cg=YfvZ`wRbWlAheXdcVoab?S%GR`2&X`|Ql& z>{k1=naetr{9EmR#kc48O{m|ySnAcGiJ7hT2^<Wf@AC?6m&Vq%teGXeZ^i8=ic=3b z*r}d6GRNw()BaT_KJn%8Iqy$iAG*B#$wm7vroHE!S!MT2&OXp`OkZ~Y)F%?ghtgKt z-&>w2cy;49yZ@s05=#E9`!8jOu32p_Vi#=~wzY8LH@kNs+DD|j6!%F)ty#z<CA_au z-ptGTwS(RBlPf31e|EB$$YZ=d@y|v3@A`)|FA0CLbG#ksAjBqXpLX@6=H%|x_KN1t z8)SE!wlj&{6MZ<lbwBetovU{`m+q686J+#9S!Lgc+UC$B6Tj{IaD(B-^(|9uMeP@E zj4~45xA*S(zus>g>>}F~U+w&L(Vjzh63=8lXZyFo1%@{KpZ0w@qRM`-W3_$$yXA3{ z4P@;-#E!BR%{#r%`uN?W>YKjp6S$oGeRlm)yTuc~8EA-$*v&TH8aOetb-%jD^=98{ zMLWT7iX{`6^zF{u&U$qBzTdv%^KuP!zFypa;b!dgdCC6!mj>kWRWIeSUvYEx5z~g% z`xnGtpJvOUvA@7dVu9}Or*@2y^AG&kz_kDS%rCmv`3~C2UDCfPCh*BN$oB5nFG;QY zB~BT-d!GDn`*C}@d&oTZeJm?@FRfdfW>@!u+ga!9MSGQbXO1?0xoE$WPozZk(?xr= zAEvKN8&=y-eJ66jpnkRe4UO|%ldD(TubFXifo4*x{nsaJUPs2a+JAg2zK*Ll-}c`1 z$*aF4wc2|m6h&W{*17jcKTi`wrM_L?=5_qG!hZW^6xUC`^CrMPu{z#jS+c)9%gUG} z)lV1q2OenPa(}I9-#eLg%N`C5`^Vegr%kS2z5niI0m*+GnCzcz3pb65Z{089_Umnr z(<j>{`roJNX|UK|zp}aM@vQ%Q<^CtUkh-pCx9YvbEpJo*eUpkd9lh};VE?3BhFcoH zT-<N|_P0F4Ql9;$QY(KRf2Fzqgh$aVgZkC`d82|Zd3}Dm?`es~Yt5wA{V%@P1yu7L z+}F?a`23n!rhU`yoOjgJVA=29CSrfyeu6EiZv?`A$_B?~_ZQgis5|_6`>AVtSL}OQ zsGfa!fA<U3UHr<I?PV0+SH7HGXm7VTQ2hO#wfl7v8&<UYt+i)%NJ#Yk?qz>pPx(d7 zskZ$Gu6{09{zGNI^6A58u5RbH>z(|&FG8W+eqX|>8JW2ZcKne|2hMw(+#7bU?c<$< z6MKu7ra5rE*lzzvSi&~)+m8L-%FQ!g&Mw>^TeP&>&~K9cV|&fr+Om`Pdo5^9_Wka) zpPy^#jn-Yc_S>e-YhC_B#r`F?<O8c{c{{t$&TO+Pa`yXuib|hbAZS;ddEH$7ugboQ z-Kh@uYpv`g8kT(Ke6f9h$W*7>{K}X2`&@kW$v^$Xe)sS9HyHR$+P^!`YE!%4+WpmQ zrY0_1zi9ut<sX05oNBX|`rm1+BT#Fvp}1FF_^$H4+5QzVGb(cI`4-KL?FkUFQ~7>t zLb{Qv-Qt^e_vTEp-dCAlJy|{bvi+9hDIX)g?XVZ*O9a_tpI$mw;r*Vq_NjL#Khlz& zY=1B2R{XN{i|ixs8JV{1%H5xEJi$OmpmslpK`VpMj?a5FtA8(yP-x#@)U>xe__FH0 zjCD8uNzYT;_x5~!!JC92yPkgmvML<)_V<Gsr`x<fx!*=+7XOQ+oBI<Nihc_@9Bwbw zbKQN@jm7&#Lt3l+3O4U=X<pN$b+2y!rK%M%=dAno7ni?pKAR|D7b<UexS`6(J~&^V zzbN4N-shXQKi(6Wy;s4e;Ap(fHQNJOVc#@4>i2Io-l(zuVV(WHqH2jJ{7Uv*M*9^* z4u|je2@5=UWNEnlg}!NzKCbv}cfarCvZcvJ_FvjodR%V0YWMBxg7;4UH1_RaTG-rQ zWwig{y#=2hY?9sQ`edHX!h~bChY|&hj6c5ETlfCygN>I?>|djLeW}gsllI~)m)#!n zEA3BPrnT?^<Kq4D&a5t*ZY;L<aAcVLdd2U3%XV*3(Y{w_UtMzP$I+&%`})!gGz^YP z+6AdD4dtIuyMKzgw&GM)K0BtYrM#WxYIf0{hI?!81nzTIxOs8or4#lQi*Bc^d|0=? zfA%g=xY^5>o?h~Raj|_h_vQYhOT+hXxW`fLSFqW>rb<(Gda}{}jhT02&sq1`Kfg37 z=k>L}wo4-lHu6oVwKon*GTV~R2d|GNyI;~uOW0-mciV}ZjUS%d&QGd;_<zk&`>4Kj zVeZAp?8`IQxu5%A+s~SDAm{Cj1@>)+3#~#z7uxSgwY|4ebmM-7Bg*MLU#s@}aMfmZ zOzN_KbK1c8PbGuhn%&IBZ5Jx{_pEq--hLhTzTZ-^M}?a-><o7tv;C}Lw=aK%{GTn} zNA0<iUM>>;d&ItjzuEA)|22D=%-)>^ehcg;%lS^T5?WwybWiNbO3{t>Gct=NGySjF zze!|cUB{%Z{YvL_ABQa#+b70+FmriK<^EOqC*P<Ds@Pq9=;@*RO=922SG5agI~&-o z+&4#XtM}3U|L>f(;$D1g|JIxUKL3Q|{g37+eJ}7^us_4*n0-j-!u<;zCgnJ<vEBds zLv>xx*D8BS!v!xwg~jdH_1o<FwnKd1wY@P9m&a7vPl{k?zi>ltpLT2t>*rVE`yO+? zk~*lbZ#S>Xh@D~W(f#J<19uAlJ+l9H#d^u8gk*c>IUz-FXDrzNW?#T{E1?DZ3njw; zHLkI>cZv@ZVE$iWf2aD`!%4#8`x&&o#1@z<+i~B??`*qJX@BI^E+#(V54Hy3C*5m# zrR?&UCpLZccG;(`@h^qza>@SF-RtWPT)Ju>CRe)bq*uv)gJ;u=o1ggae|TQx*zK(w z>}T40f3Xtw-+#Aensv?JPWyR0Tqk#EY3!dHyDRu`v54I=1&bwL&S>m+p7^6Mc%{HT zx4;QCsrRq%Wni2V=<NG_?=OwXTX-***c;TWQfi!8YX4SewWX+4!~W2aD=N)T{O!MH z2a8#D1=!!;s2=~}X2kyLXH(|J?$op|d;i!wDA;&E|FubW`!>t&YZWTdU_Gm0@Ae@l zL+GD|-4^biZK6Lk_oZgft9`mD#Lm|4*UV_;Yx`wq7*`#*balVqGxjHfRt@&;GdUBd zeBQ8Mr2Wp(+gms6KhO00^T(SJ_V*U2ao7Cq+%L1~Y>Z>Dv3>j#_tN--pKOnrKDDWi z?b?5|;TZ>W-rv1hVYg1pos_bBT(McVa+d2pj+bupqLr`NZ@VkAu6AbW{>XcdEyuk| z?B~Y?g6y&1@9@UivMXT!1Fw8OTVa2DrS;u$K07t{|4#p4w@^#No@LSc(5x9BZFwK@ zWLL*_+3)69v+BztN%(l(pUCAet&6kwnr?8@*<QlF@8_(WU(z0*w(n(PeJ8MEv%P=O z^sd-Dm-b)oiP`$@(3$<QQ&;>{2-~vXzA@%iXYplwe)n%}@6%fMuf4|CTwSrmzO7O5 zg(<tlK9jaR3ddaL?bm!@r}65tnH}G@lwH>{()NioH`p-9bl5puU$_%w4^P*pvmkrc zJv@;KvL}8r-=;^0&g^GDtaOblY|H+?+diFZDZXs)y*e!iWKYSvZP_4uw9IAoZ}cnf z%j$E>0NFDo#CNY>fVf>5zqeabmAjqxEzQ!8uQK<gaz9=PvS+W9>Mf8xIR!q^AbU>S zke>MI(3$<-hUaX=!nW+MveMt$QGD4xW1dbt$exCWA=g3ntZb0%{r7=+U%5cePLMs% zJ}+Dv8(_GvyHoJyhNc9&9X7MKeY(}SZ)4&6!ytRS3rhBb>^W6&CJtoJPrk$6&kvp1 zFPJ>}zIND_{r+JG8(NDm+nc>k(*2Ouy1#o}>T!@gg*olI=dbJ7O*y>o-%*!&`?Y)% zR(n=l+}o~j{$Ih-gnbGz)gL}*HrXw)Kk+>6@#+1h`c6j#c5L21EqZ!w?43*YCnRqj zdVJ`Nz3}|QijrYl?4vAJXtft#-tWV9*AnEP<31a!K=zc$91ji<G}!mBv8f5<pIJ&l zG8cA;+ud(5Gpll)ZEL_>ZQ4?-xX*KsL><VU0(Yi)AbVC$-)9c8XVHgiMIR5Hv5$Pb zEmbjWi~W^5O)GnfFYlk|opuXkPgLypn;?5S51B_vItkecfBQ8DWY1y_OQjn}4EIf* zxZ%^14QY0F*SU4xNbK0hbf@JV$e!MccPxTCHt$~~)X@U6NBPaaWiJn%vF9s4B`6oR z#h&B$)`rI7%lq}RcC7;0b8^?pLm+$j5|aa7K9{r0GIEgv*~86XAv|&6y}gE)uC2d5 z(S6_Cz?4ItR#|o>Hs?-(?799p<2}ef9-RK=AbSMd^}XI4I%Ch&_{mi?Y>RyjgHk;x z+_=Lw901vKBKk%N$ev%^F)k-(y|HbT{l*BgM`54-^o1ej@cQUTs;_v{*J*oMTNZwd zI;gcz^hT9*OxfxEBE=f_c742Lzxc-dq|axz?7!FNH@h<EjQysKMe`?hUf%!XsNTQT z#ar#2IQOqw-M?i2vK1@76;E!n|Ka?Nf9@o+eZdVcmgJn9YabiCDd&Ed_`bqgu_N(E zJ@!4H(YwJ?DAR6|+lQ;sWvBN?SIpZ8vggFaBs-8je-(vSR0o~0e>F?VenRKv{p_}` z3CoJN+BZr)`VF#2Y|h&|AbY+w@13^qgO**5M7l-pxw-boUmEjV4pX&T?7V5Us&|^* z^vN5S@6_zr$8qk-O^`kNjlS#!*|T<)xCh9dx8IVO3xdwr?~G6|@9n(2-|}hPrHRE` z?T>wKp9HceDbV;9$eskIg^oI}mG-H!Jp2Um&&|Gl65>Y%?cVwpsOP#Q?2Gzfe(e%Z zquoop_b);A==r=@53=WKtfLsno-6CzFBS!zv1dvDDBRz9d4H3t)2IE#TkYAWsec05 zQy1%cpk#8JeY#j+_1Rytc00|Fwt(!JIpyiH7q%Mvoc7JTHFaKsoxsZ+hSXV&`$8^V z5sE20Z9gSC{M@dOm-aj7JW~05W{drbb4S~AgU;+XpA>R&Lg!`snlPu6tBSYof1;?i z6l9Ov;i~;0{~VdJZ_9>vTKh7t=Q4uAP3lEcap7|{JMqjFA`f#P*?u;9YpNyx#+FUg zWl?n5X?x}ft{#v*-}Z#Ef$WJb<Kn6dI<tT2pD3xRotN#iL#qnb7jNDF@#7VSHT_HW ztGB%h2l>aebIDF$W+OYv_o6Bwd*Y>{><(lJ+37x<WGNEhX{YB<lpm3mwQugW$3H># zgo&4e{UdfJH5Oz~iByGuNzj@75Bcmq^mSgguVqScTvWVuf84&<nIQiZC+@Qb*^?#q z|Cj1q<9+J__SAvGEkplrvBECReLAb3-4MT;wlAN#N%llxr`_MssXU-?`?x#n2FO2c zElM>Yd%mBZ@Vz|f%>He^12wxkFWax#7*?>nc<X+JM4eY4dyF|_BtZ69vmdQkZzo_E zw61k6$ev{f{o+rmO2F%*TN-nH%)h7Z6PMIccl*?C$EOqa<^P(~`zNbMoU&bcdjGya zuFdZ^oZg?vx!>`A*BN_$<%tLXww<xR%eqznW8WG3x4SP#FJql&FMFS{fq`$Hy>j%8 zAl9#|ZQTznS+k3Cp8ckSPECGJyZ1VZ27gRmA7-~r=+Ac%p~`&?g;CEz_VlVv-vhE| z?t9+9AbV`hio(x#ow5I{w?Xz@+Zp>`N<vG{^_{WrxodF&WRJnhY(J1aI~e(#p3A80 z`xIxPu$yz9z1N@W1sa)#c7p3}eL36{X17W_?`F{SihW}1PCNqHQ+(>(e~>+ijCWRp z?732Vko98M8GHRzd#*idJ7Yg3{ZiZgzBBfE#@Y@bd)_%mzX93P_VVWeWepR%&7FLz zAbU~+?;PIucZY5KKK)Hyondz6>g}y+Q!Do6nIsx8tUbNo_Lk^okUiN+AHeo(W?KCI zVAmOYflr3ZKDM2)w_onw0J3M+&xBhbdyYDLgYD^cl5^0DRM}U!TC@ukZXExbKG|m( z?n`Lj*kEwaabMNd=wqH!YwQ?;Rv-Jn=Cu6{Z;xlTD^J^BjQxG({f5)_y-&4|-s?KE zf1j-X(sylV_WR7aAoZZ{%znR<9KIlX_O&j|0{Q2(Iol_b->P;~mzWuX{L?02Cik|> z$gado$IkC~uARKdlND>rI`(aR6wJi1_O!hc-<j(md-_%tp9b0E7x?DH!LBp=?|)cQ z_@M2~ey5#+Ru}rt?0;V%dH`h4nFrslgX|Hz)2hF1+u^;6Tf+oE_ISN5jhr2L!8ZRx za@?y-r+xpb+?LP1P;2)q)`}fu&w(|2?||&7QDyE1*<)Zbv*dEunf<=!_ltdMJG1|3 zVB)uTeP{N+cD@(}vd3lprN1D1wrsJx=#ZsqcauYAIVju?-MD(Jufu3xkb`f(!2`#A z%9qwu^Yzx)ZJWDS3}jD5rTa9HJ*HhduY&9;zSVyFY}c9nE3c(@{0I4Gb>x-%eP{Mh qoAU7)$R0C~JDDJRc5VB{I<MPw-?<sa(?RxlJq~~6uYL|Z?hODSCVBDz diff --git a/test/examples/Si_bulk_energy_force.inp b/test/examples/Si_bulk_energy_force.inp deleted file mode 100644 index 316166b..0000000 --- a/test/examples/Si_bulk_energy_force.inp +++ /dev/null @@ -1,63 +0,0 @@ -&GLOBAL - PROJECT Si_bulk8 - RUN_TYPE ENERGY_FORCE - PRINT_LEVEL MEDIUM -&END GLOBAL -&FORCE_EVAL - METHOD Quickstep - &SUBSYS - &KIND Si - ELEMENT Si - BASIS_SET DZVP-GTH-PADE - POTENTIAL GTH-PADE-q4 - &END KIND - &CELL - A 5.430697500 0.000000000 0.000000000 - B 0.000000000 5.430697500 0.000000000 - C 0.000000000 0.000000000 5.430697500 - &END CELL - &COORD - Si 0.000000000 0.000000000 0.000000000 - Si 0.000000000 2.715348700 2.715348700 - Si 2.715348700 2.715348700 0.000000000 - Si 2.715348700 0.000000000 2.715348700 - Si 4.073023100 1.357674400 4.073023100 - Si 1.357674400 1.357674400 1.357674400 - Si 1.357674400 4.073023100 4.073023100 - Si 4.073023100 4.073023100 1.357674400 - &END COORD - &END SUBSYS - &DFT - BASIS_SET_FILE_NAME BASIS_SET - POTENTIAL_FILE_NAME GTH_POTENTIALS - &QS - EPS_DEFAULT 1.0E-10 - &END QS - &MGRID - NGRIDS 4 - CUTOFF 300 - REL_CUTOFF 60 - &END MGRID - &XC - &XC_FUNCTIONAL PADE - &END XC_FUNCTIONAL - &END XC - &SCF - SCF_GUESS ATOMIC - EPS_SCF 1.0E-7 - MAX_SCF 300 - &DIAGONALIZATION ON - ALGORITHM STANDARD - &END DIAGONALIZATION - &MIXING T - METHOD BROYDEN_MIXING - ALPHA 0.4 - NBROYDEN 8 - &END MIXING - &END SCF - &END DFT - &PRINT - &FORCES ON - &END FORCES - &END PRINT -&END FORCE_EVAL diff --git a/test/examples/Si_bulk_energy_force.out b/test/examples/Si_bulk_energy_force.out deleted file mode 100644 index 70c39df..0000000 --- a/test/examples/Si_bulk_energy_force.out +++ /dev/null @@ -1,760 +0,0 @@ - DBCSR| Multiplication driver SMM - DBCSR| Multrec recursion limit 512 - DBCSR| Multiplication stack size 1000 - DBCSR| Multiplication size stacks 3 - DBCSR| Use subcommunicators T - DBCSR| Use MPI combined types F - DBCSR| Use MPI memory allocation T - DBCSR| Use Communication thread T - DBCSR| Communication thread load 87 - - - **** **** ****** ** PROGRAM STARTED AT 2016-02-04 10:13:38.137 - ***** ** *** *** ** PROGRAM STARTED ON lauri-Lenovo-Z50-70 - ** **** ****** PROGRAM STARTED BY lauri - ***** ** ** ** ** PROGRAM PROCESS ID 30677 - **** ** ******* ** PROGRAM STARTED IN /home/lauri/Dropbox/nomad-dev/nomad-l - ab-base/parsers/cp2k/test/examples - - CP2K| version string: CP2K version 2.6.2 - CP2K| source code revision number: svn:15893 - CP2K| is freely available from http://www.cp2k.org/ - CP2K| Program compiled at ke 4.11.2015 08.48.42 +0200 - CP2K| Program compiled on lauri-Lenovo-Z50-70 - CP2K| Program compiled for Linux-x86-64-gfortran_basic - CP2K| Input file name Si_bulk_energy_force.inp - - GLOBAL| Force Environment number 1 - GLOBAL| Basis set file name BASIS_SET - GLOBAL| Geminal file name BASIS_GEMINAL - GLOBAL| Potential file name GTH_POTENTIALS - GLOBAL| MM Potential file name MM_POTENTIAL - GLOBAL| Coordinate file name __STD_INPUT__ - GLOBAL| Method name CP2K - GLOBAL| Project name Si_bulk8 - GLOBAL| Preferred FFT library FFTW3 - GLOBAL| Preferred diagonalization lib. SL - GLOBAL| Run type ENERGY_FORCE - GLOBAL| All-to-all communication in single precision F - GLOBAL| FFTs using library dependent lengths F - GLOBAL| Global print level MEDIUM - GLOBAL| Total number of message passing processes 1 - GLOBAL| Number of threads for this process 1 - GLOBAL| This output is from process 0 - - MEMORY| system memory details [Kb] - MEMORY| rank 0 min max average - MEMORY| MemTotal 8070384 8070384 8070384 8070384 - MEMORY| MemFree 1463328 1463328 1463328 1463328 - MEMORY| Buffers 870024 870024 870024 870024 - MEMORY| Cached 3864964 3864964 3864964 3864964 - MEMORY| Slab 674124 674124 674124 674124 - MEMORY| SReclaimable 635132 635132 635132 635132 - MEMORY| MemLikelyFree 6833448 6833448 6833448 6833448 - - - *** Fundamental physical constants (SI units) *** - - *** Literature: B. J. Mohr and B. N. Taylor, - *** CODATA recommended values of the fundamental physical - *** constants: 2006, Web Version 5.1 - *** http://physics.nist.gov/constants - - Speed of light in vacuum [m/s] 2.99792458000000E+08 - Magnetic constant or permeability of vacuum [N/A**2] 1.25663706143592E-06 - Electric constant or permittivity of vacuum [F/m] 8.85418781762039E-12 - Planck constant (h) [J*s] 6.62606896000000E-34 - Planck constant (h-bar) [J*s] 1.05457162825177E-34 - Elementary charge [C] 1.60217648700000E-19 - Electron mass [kg] 9.10938215000000E-31 - Electron g factor [ ] -2.00231930436220E+00 - Proton mass [kg] 1.67262163700000E-27 - Fine-structure constant 7.29735253760000E-03 - Rydberg constant [1/m] 1.09737315685270E+07 - Avogadro constant [1/mol] 6.02214179000000E+23 - Boltzmann constant [J/K] 1.38065040000000E-23 - Atomic mass unit [kg] 1.66053878200000E-27 - Bohr radius [m] 5.29177208590000E-11 - - *** Conversion factors *** - - [u] -> [a.u.] 1.82288848426455E+03 - [Angstrom] -> [Bohr] = [a.u.] 1.88972613288564E+00 - [a.u.] = [Bohr] -> [Angstrom] 5.29177208590000E-01 - [a.u.] -> [s] 2.41888432650478E-17 - [a.u.] -> [fs] 2.41888432650478E-02 - [a.u.] -> [J] 4.35974393937059E-18 - [a.u.] -> [N] 8.23872205491840E-08 - [a.u.] -> [K] 3.15774647902944E+05 - [a.u.] -> [kJ/mol] 2.62549961709828E+03 - [a.u.] -> [kcal/mol] 6.27509468713739E+02 - [a.u.] -> [Pa] 2.94210107994716E+13 - [a.u.] -> [bar] 2.94210107994716E+08 - [a.u.] -> [atm] 2.90362800883016E+08 - [a.u.] -> [eV] 2.72113838565563E+01 - [a.u.] -> [Hz] 6.57968392072181E+15 - [a.u.] -> [1/cm] (wave numbers) 2.19474631370540E+05 - [a.u./Bohr**2] -> [1/cm] 5.14048714338585E+03 - - - CELL_TOP| Volume [angstrom^3]: 160.165 - CELL_TOP| Vector a [angstrom 5.431 0.000 0.000 |a| = 5.431 - CELL_TOP| Vector b [angstrom 0.000 5.431 0.000 |b| = 5.431 - CELL_TOP| Vector c [angstrom 0.000 0.000 5.431 |c| = 5.431 - CELL_TOP| Angle (b,c), alpha [degree]: 90.000 - CELL_TOP| Angle (a,c), beta [degree]: 90.000 - CELL_TOP| Angle (a,b), gamma [degree]: 90.000 - CELL_TOP| Numerically orthorhombic: YES - - GENERATE| Preliminary Number of Bonds generated: 0 - GENERATE| Achieved consistency in connectivity generation. - - CELL| Volume [angstrom^3]: 160.165 - CELL| Vector a [angstrom]: 5.431 0.000 0.000 |a| = 5.431 - CELL| Vector b [angstrom]: 0.000 5.431 0.000 |b| = 5.431 - CELL| Vector c [angstrom]: 0.000 0.000 5.431 |c| = 5.431 - CELL| Angle (b,c), alpha [degree]: 90.000 - CELL| Angle (a,c), beta [degree]: 90.000 - CELL| Angle (a,b), gamma [degree]: 90.000 - CELL| Numerically orthorhombic: YES - - CELL_REF| Volume [angstrom^3]: 160.165 - CELL_REF| Vector a [angstrom 5.431 0.000 0.000 |a| = 5.431 - CELL_REF| Vector b [angstrom 0.000 5.431 0.000 |b| = 5.431 - CELL_REF| Vector c [angstrom 0.000 0.000 5.431 |c| = 5.431 - CELL_REF| Angle (b,c), alpha [degree]: 90.000 - CELL_REF| Angle (a,c), beta [degree]: 90.000 - CELL_REF| Angle (a,b), gamma [degree]: 90.000 - CELL_REF| Numerically orthorhombic: YES - - ******************************************************************************* - ******************************************************************************* - ** ** - ** ##### ## ## ** - ** ## ## ## ## ## ** - ** ## ## ## ###### ** - ** ## ## ## ## ## ##### ## ## #### ## ##### ##### ** - ** ## ## ## ## ## ## ## ## ## ## ## ## ## ## ** - ** ## ## ## ## ## ## ## #### ### ## ###### ###### ** - ** ## ### ## ## ## ## ## ## ## ## ## ## ** - ** ####### ##### ## ##### ## ## #### ## ##### ## ** - ** ## ## ** - ** ** - ** ... make the atoms dance ** - ** ** - ** Copyright (C) by CP2K Developers Group (2000 - 2014) ** - ** ** - ******************************************************************************* - - DFT| Spin restricted Kohn-Sham (RKS) calculation RKS - DFT| Multiplicity 1 - DFT| Number of spin states 1 - DFT| Charge 0 - DFT| Self-interaction correction (SIC) NO - DFT| Cutoffs: density 1.000000E-10 - DFT| gradient 1.000000E-10 - DFT| tau 1.000000E-10 - DFT| cutoff_smoothing_range 0.000000E+00 - DFT| XC density smoothing NONE - DFT| XC derivatives PW - FUNCTIONAL| ROUTINE=NEW - FUNCTIONAL| PADE: - FUNCTIONAL| S. Goedecker, M. Teter and J. Hutter, Phys. Rev. B 54, 1703 (1996) - - QS| Method: GPW - QS| Density plane wave grid type NON-SPHERICAL FULLSPACE - QS| Number of grid levels: 4 - QS| Density cutoff [a.u.]: 150.0 - QS| Multi grid cutoff [a.u.]: 1) grid level 150.0 - QS| 2) grid level 50.0 - QS| 3) grid level 16.7 - QS| 4) grid level 5.6 - QS| Grid level progression factor: 3.0 - QS| Relative density cutoff [a.u.]: 30.0 - QS| Consistent realspace mapping and integration - QS| Interaction thresholds: eps_pgf_orb: 1.0E-05 - QS| eps_filter_matrix: 0.0E+00 - QS| eps_core_charge: 1.0E-12 - QS| eps_rho_gspace: 1.0E-10 - QS| eps_rho_rspace: 1.0E-10 - QS| eps_gvg_rspace: 1.0E-05 - QS| eps_ppl: 1.0E-02 - QS| eps_ppnl: 1.0E-07 - - - ATOMIC KIND INFORMATION - - 1. Atomic kind: Si Number of atoms: 8 - - Orbital Basis Set DZVP-GTH-PADE - - Number of orbital shell sets: 2 - Number of orbital shells: 5 - Number of primitive Cartesian functions: 5 - Number of Cartesian basis functions: 14 - Number of spherical basis functions: 13 - Norm type: 2 - - Normalised Cartesian orbitals: - - Set Shell Orbital Exponent Coefficient - - 1 1 3s 1.203242 0.269412 - 0.468841 -0.102290 - 0.167986 -0.147195 - 0.057562 -0.015996 - - 1 2 4s 1.203242 0.000000 - 0.468841 0.000000 - 0.167986 0.000000 - 0.057562 0.083755 - - 1 3 4px 1.203242 0.085242 - 0.468841 -0.143473 - 0.167986 -0.083408 - 0.057562 -0.014565 - 1 3 4py 1.203242 0.085242 - 0.468841 -0.143473 - 0.167986 -0.083408 - 0.057562 -0.014565 - 1 3 4pz 1.203242 0.085242 - 0.468841 -0.143473 - 0.167986 -0.083408 - 0.057562 -0.014565 - - 1 4 5px 1.203242 0.000000 - 0.468841 0.000000 - 0.167986 0.000000 - 0.057562 0.040189 - 1 4 5py 1.203242 0.000000 - 0.468841 0.000000 - 0.167986 0.000000 - 0.057562 0.040189 - 1 4 5pz 1.203242 0.000000 - 0.468841 0.000000 - 0.167986 0.000000 - 0.057562 0.040189 - - 2 1 3dx2 0.450000 0.406941 - 2 1 3dxy 0.450000 0.704842 - 2 1 3dxz 0.450000 0.704842 - 2 1 3dy2 0.450000 0.406941 - 2 1 3dyz 0.450000 0.704842 - 2 1 3dz2 0.450000 0.406941 - - Potential information for GTH-PADE-q4 - - Description: Goedecker-Teter-Hutter pseudopotential - Goedecker et al., PRB 54, 1703 (1996) - Hartwigsen et al., PRB 58, 3641 (1998) - Krack, TCA 114, 145 (2005) - - Gaussian exponent of the core charge distribution: 2.582645 - Electronic configuration (s p d ...): 2 2 - - Parameters of the local part of the GTH pseudopotential: - - rloc C1 C2 C3 C4 - 0.440000 -7.336103 - - Parameters of the non-local part of the GTH pseudopotential: - - l r(l) h(i,j,l) - - 0 0.422738 5.906928 -1.261894 - -1.261894 3.258196 - 1 0.484278 2.727013 - - - MOLECULE KIND INFORMATION - - - All atoms are their own molecule, skipping detailed information - - - TOTAL NUMBERS AND MAXIMUM NUMBERS - - Total number of - Atomic kinds: 1 - - Atoms: 8 - - Shell sets: 16 - - Shells: 40 - - Primitive Cartesian functions: 40 - - Cartesian basis functions: 112 - - Spherical basis functions: 104 - - Maximum angular momentum of- Orbital basis functions: 2 - - Local part of the GTH pseudopotential: 0 - - Non-local part of the GTH pseudopotential: 2 - - - MODULE QUICKSTEP: ATOMIC COORDINATES IN angstrom - - Atom Kind Element X Y Z Z(eff) Mass - - 1 1 Si 14 0.000000 0.000000 0.000000 4.00 28.0855 - 2 1 Si 14 0.000000 2.715349 2.715349 4.00 28.0855 - 3 1 Si 14 2.715349 2.715349 0.000000 4.00 28.0855 - 4 1 Si 14 2.715349 0.000000 2.715349 4.00 28.0855 - 5 1 Si 14 4.073023 1.357674 4.073023 4.00 28.0855 - 6 1 Si 14 1.357674 1.357674 1.357674 4.00 28.0855 - 7 1 Si 14 1.357674 4.073023 4.073023 4.00 28.0855 - 8 1 Si 14 4.073023 4.073023 1.357674 4.00 28.0855 - - - - - SCF PARAMETERS Density guess: ATOMIC - -------------------------------------------------------- - max_scf: 300 - max_scf_history: 0 - max_diis: 4 - -------------------------------------------------------- - eps_scf: 1.00E-07 - eps_scf_history: 0.00E+00 - eps_diis: 1.00E-01 - eps_eigval: 1.00E-05 - -------------------------------------------------------- - level_shift [a.u.]: 0.00 - -------------------------------------------------------- - Mixing method: BROYDEN_MIXING - charge density mixing in g-space - -------------------------------------------------------- - No outer SCF - - PW_GRID| Information for grid number 1 - PW_GRID| Cutoff [a.u.] 150.0 - PW_GRID| spherical cutoff: NO - PW_GRID| Bounds 1 -30 29 Points: 60 - PW_GRID| Bounds 2 -30 29 Points: 60 - PW_GRID| Bounds 3 -30 29 Points: 60 - PW_GRID| Volume element (a.u.^3) 0.5004E-02 Volume (a.u.^3) 1080.8451 - PW_GRID| Grid span FULLSPACE - - PW_GRID| Information for grid number 2 - PW_GRID| Cutoff [a.u.] 50.0 - PW_GRID| spherical cutoff: NO - PW_GRID| Bounds 1 -18 17 Points: 36 - PW_GRID| Bounds 2 -18 17 Points: 36 - PW_GRID| Bounds 3 -18 17 Points: 36 - PW_GRID| Volume element (a.u.^3) 0.2317E-01 Volume (a.u.^3) 1080.8451 - PW_GRID| Grid span FULLSPACE - - PW_GRID| Information for grid number 3 - PW_GRID| Cutoff [a.u.] 16.7 - PW_GRID| spherical cutoff: NO - PW_GRID| Bounds 1 -10 9 Points: 20 - PW_GRID| Bounds 2 -10 9 Points: 20 - PW_GRID| Bounds 3 -10 9 Points: 20 - PW_GRID| Volume element (a.u.^3) 0.1351 Volume (a.u.^3) 1080.8451 - PW_GRID| Grid span FULLSPACE - - PW_GRID| Information for grid number 4 - PW_GRID| Cutoff [a.u.] 5.6 - PW_GRID| spherical cutoff: NO - PW_GRID| Bounds 1 -6 5 Points: 12 - PW_GRID| Bounds 2 -6 5 Points: 12 - PW_GRID| Bounds 3 -6 5 Points: 12 - PW_GRID| Volume element (a.u.^3) 0.6255 Volume (a.u.^3) 1080.8451 - PW_GRID| Grid span FULLSPACE - - POISSON| Solver PERIODIC - POISSON| Periodicity XYZ - - RS_GRID| Information for grid number 1 - RS_GRID| Bounds 1 -30 29 Points: 60 - RS_GRID| Bounds 2 -30 29 Points: 60 - RS_GRID| Bounds 3 -30 29 Points: 60 - - RS_GRID| Information for grid number 2 - RS_GRID| Bounds 1 -18 17 Points: 36 - RS_GRID| Bounds 2 -18 17 Points: 36 - RS_GRID| Bounds 3 -18 17 Points: 36 - - RS_GRID| Information for grid number 3 - RS_GRID| Bounds 1 -10 9 Points: 20 - RS_GRID| Bounds 2 -10 9 Points: 20 - RS_GRID| Bounds 3 -10 9 Points: 20 - - RS_GRID| Information for grid number 4 - RS_GRID| Bounds 1 -6 5 Points: 12 - RS_GRID| Bounds 2 -6 5 Points: 12 - RS_GRID| Bounds 3 -6 5 Points: 12 - - DISTRIBUTION OF THE PARTICLES (ROWS) - Process row Number of particles Number of matrix rows - 0 8 -1 - Sum 8 -1 - - DISTRIBUTION OF THE PARTICLES (COLUMNS) - Process col Number of particles Number of matrix columns - 0 8 -1 - Sum 8 -1 - - DISTRIBUTION OF THE NEIGHBOR LISTS - Total number of particle pairs: 2952 - Total number of matrix elements: 498888 - Average number of particle pairs: 2952 - Maximum number of particle pairs: 2952 - Average number of matrix element: 498888 - Maximum number of matrix elements: 498888 - - - DISTRIBUTION OF THE OVERLAP MATRIX - Number of non-zero blocks: 36 - Percentage non-zero blocks: 100.00 - Average number of blocks per CPU: 36 - Maximum number of blocks per CPU: 36 - Average number of matrix elements per CPU: 6094 - Maximum number of matrix elements per CPU: 6094 - - Number of electrons: 32 - Number of occupied orbitals: 16 - Number of molecular orbitals: 16 - - Number of orbital functions: 104 - Number of independent orbital functions: 104 - - Extrapolation method: initial_guess - - Atomic guess: The first density matrix is obtained in terms of atomic orbitals - and electronic configurations assigned to each atomic kind - - Guess for atomic kind: Si - - Electronic structure - Total number of core electrons 10.00 - Total number of valence electrons 4.00 - Total number of electrons 14.00 - Multiplicity not specified - S [ 2.00 2.00] 2.00 - P [ 6.00] 2.00 - - - ******************************************************************************* - Iteration Convergence Energy [au] - ******************************************************************************* - 1 0.191310 -3.618313869735 - 2 0.731569E-01 -3.691159009622 - 3 0.405574E-02 -3.699900512584 - 4 0.328704E-02 -3.699908407293 - 5 0.320845E-02 -3.699909118998 - 6 0.316809E-02 -3.699909477757 - 7 0.331859E-05 -3.699923449535 - 8 0.110258E-06 -3.699923449550 - - Energy components [Hartree] Total Energy :: -3.699923449550 - Band Energy :: -1.012729790251 - Kinetic Energy :: 1.397012768229 - Potential Energy :: -5.096936217779 - Virial (-V/T) :: 3.648453567279 - Core Energy :: -5.703543362687 - XC Energy :: -0.980691562795 - Coulomb Energy :: 2.984311475932 - Total Pseudopotential Energy :: -7.145739758818 - Local Pseudopotential Energy :: -7.987908627736 - Nonlocal Pseudopotential Energy :: 0.842168868918 - Confinement :: 0.451836279031 - - Orbital energies State L Occupation Energy[a.u.] Energy[eV] - - 1 0 2.000 -0.378230 -10.292155 - - 1 1 2.000 -0.128135 -3.486734 - - Re-scaling the density matrix to get the right number of electrons - # Electrons Trace(P) Scaling factor - 32 31.129 1.028 - - - SCF WAVEFUNCTION OPTIMIZATION - - Step Update method Time Convergence Total energy Change - ------------------------------------------------------------------------------ - - Trace(PS): 32.0000000000 - Electronic density on regular grids: -31.9999999982 0.0000000018 - Core density on regular grids: 31.9999999939 -0.0000000061 - Total charge density on r-space grids: -0.0000000043 - Total charge density g-space grids: -0.0000000043 - - 1 NoMix/Diag. 0.40E+00 0.8 0.75558724 -32.2320848878 -3.22E+01 - - Trace(PS): 32.0000000000 - Electronic density on regular grids: -31.9999999956 0.0000000044 - Core density on regular grids: 31.9999999939 -0.0000000061 - Total charge density on r-space grids: -0.0000000017 - Total charge density g-space grids: -0.0000000017 - - 2 Broy./Diag. 0.40E+00 1.4 0.05667976 -31.1418135481 1.09E+00 - - Trace(PS): 32.0000000000 - Electronic density on regular grids: -31.9999999907 0.0000000093 - Core density on regular grids: 31.9999999939 -0.0000000061 - Total charge density on r-space grids: 0.0000000033 - Total charge density g-space grids: 0.0000000033 - - 3 Broy./Diag. 0.40E+00 1.3 0.09691469 -31.1974003416 -5.56E-02 - - Trace(PS): 32.0000000000 - Electronic density on regular grids: -31.9999999897 0.0000000103 - Core density on regular grids: 31.9999999939 -0.0000000061 - Total charge density on r-space grids: 0.0000000042 - Total charge density g-space grids: 0.0000000042 - - 4 Broy./Diag. 0.40E+00 1.4 0.00245608 -31.3378474040 -1.40E-01 - - Trace(PS): 32.0000000000 - Electronic density on regular grids: -31.9999999888 0.0000000112 - Core density on regular grids: 31.9999999939 -0.0000000061 - Total charge density on r-space grids: 0.0000000051 - Total charge density g-space grids: 0.0000000051 - - 5 Broy./Diag. 0.40E+00 1.4 0.00235460 -31.3009654398 3.69E-02 - - Trace(PS): 32.0000000000 - Electronic density on regular grids: -31.9999999888 0.0000000112 - Core density on regular grids: 31.9999999939 -0.0000000061 - Total charge density on r-space grids: 0.0000000051 - Total charge density g-space grids: 0.0000000051 - - 6 Broy./Diag. 0.40E+00 1.4 0.00007565 -31.2972158934 3.75E-03 - - Trace(PS): 32.0000000000 - Electronic density on regular grids: -31.9999999889 0.0000000111 - Core density on regular grids: 31.9999999939 -0.0000000061 - Total charge density on r-space grids: 0.0000000050 - Total charge density g-space grids: 0.0000000050 - - 7 Broy./Diag. 0.40E+00 1.4 0.00009004 -31.2977293749 -5.13E-04 - - Trace(PS): 32.0000000000 - Electronic density on regular grids: -31.9999999889 0.0000000111 - Core density on regular grids: 31.9999999939 -0.0000000061 - Total charge density on r-space grids: 0.0000000051 - Total charge density g-space grids: 0.0000000051 - - 8 Broy./Diag. 0.40E+00 1.4 0.00000186 -31.2978454163 -1.16E-04 - - Trace(PS): 32.0000000000 - Electronic density on regular grids: -31.9999999889 0.0000000111 - Core density on regular grids: 31.9999999939 -0.0000000061 - Total charge density on r-space grids: 0.0000000051 - Total charge density g-space grids: 0.0000000051 - - 9 Broy./Diag. 0.40E+00 1.4 0.00000252 -31.2978835492 -3.81E-05 - - Trace(PS): 32.0000000000 - Electronic density on regular grids: -31.9999999889 0.0000000111 - Core density on regular grids: 31.9999999939 -0.0000000061 - Total charge density on r-space grids: 0.0000000051 - Total charge density g-space grids: 0.0000000051 - - 10 Broy./Diag. 0.40E+00 1.4 5.6405E-09 -31.2978852054 -1.66E-06 - - *** SCF run converged in 10 steps *** - - - Electronic density on regular grids: -31.9999999889 0.0000000111 - Core density on regular grids: 31.9999999939 -0.0000000061 - Total charge density on r-space grids: 0.0000000051 - Total charge density g-space grids: 0.0000000051 - - Overlap energy of the core charge distribution: 0.00000000005320 - Self energy of the core charge distribution: -82.06393942512820 - Core Hamiltonian energy: 18.06858429706012 - Hartree energy: 42.41172824581675 - Exchange-correlation energy: -9.71425832315954 - - Total energy: -31.29788520535767 - - - MULLIKEN POPULATION ANALYSIS - - # Atom Element Kind Atomic population Net charge - 1 Si 1 4.000000 -0.000000 - 2 Si 1 4.000000 0.000000 - 3 Si 1 4.000000 0.000000 - 4 Si 1 4.000000 0.000000 - 5 Si 1 4.000000 0.000000 - 6 Si 1 4.000000 -0.000000 - 7 Si 1 4.000000 0.000000 - 8 Si 1 4.000000 0.000000 - # Total charge 32.000000 -0.000000 - - - !-----------------------------------------------------------------------------! - Hirschfeld Charges - - #Atom Element Kind Ref Charge Population Net charge - 1 Si 1 4.000 4.000 -0.000 - 2 Si 1 4.000 4.000 -0.000 - 3 Si 1 4.000 4.000 -0.000 - 4 Si 1 4.000 4.000 -0.000 - 5 Si 1 4.000 4.000 -0.000 - 6 Si 1 4.000 4.000 0.000 - 7 Si 1 4.000 4.000 -0.000 - 8 Si 1 4.000 4.000 -0.000 - - Total Charge -0.000 - !-----------------------------------------------------------------------------! - - Trace(PS): 32.0000000000 - Electronic density on regular grids: -31.9999999889 0.0000000111 - Core density on regular grids: 31.9999999939 -0.0000000061 - Total charge density on r-space grids: 0.0000000051 - Total charge density g-space grids: 0.0000000051 - - - ENERGY| Total FORCE_EVAL ( QS ) energy (a.u.): -31.297885372811063 - - - ATOMIC FORCES in [a.u.] - - # Atom Kind Element X Y Z - 1 1 Si 0.00000000 0.00000000 0.00000000 - 2 1 Si 0.00000000 0.00000001 0.00000001 - 3 1 Si 0.00000001 0.00000001 0.00000000 - 4 1 Si 0.00000001 0.00000000 0.00000001 - 5 1 Si -0.00000001 -0.00000001 -0.00000001 - 6 1 Si -0.00000001 -0.00000001 -0.00000001 - 7 1 Si -0.00000001 -0.00000001 -0.00000001 - 8 1 Si -0.00000001 -0.00000001 -0.00000001 - SUM OF ATOMIC FORCES -0.00000000 -0.00000000 -0.00000000 0.00000000 - - ------------------------------------------------------------------------------- - - - - - DBCSR STATISTICS - - - - - ------------------------------------------------------------------------------- - COUNTER CPU ACC ACC% - number of processed stacks 11 0 0.0 - matmuls inhomo. stacks 0 0 0.0 - matmuls total 396 0 0.0 - flops 13 x 13 x 16 2141568 0 0.0 - flops total 2141568 0 0.0 - marketing flops 3807232 - ------------------------------------------------------------------------------- - - ------------------------------------------------------------------------------- - ---- MULTIGRID INFO ---- - ------------------------------------------------------------------------------- - count for grid 1: 360 cutoff [a.u.] 150.00 - count for grid 2: 4516 cutoff [a.u.] 50.00 - count for grid 3: 7804 cutoff [a.u.] 16.67 - count for grid 4: 3048 cutoff [a.u.] 5.56 - total gridlevel count : 15728 - - ------------------------------------------------------------------------------- - - - - - MESSAGE PASSING PERFORMANCE - - - - - ------------------------------------------------------------------------------- - - ROUTINE CALLS TOT TIME [s] AVE VOLUME [Bytes] PERFORMANCE [MB/s] - MP_Group 5 0.000 - MP_Bcast 12 0.000 5. 0.21 - MP_Allreduce 171 0.000 42. 37.27 - MP_Sync 4 0.000 - MP_Alltoall 183 0.000 15743. 10992.20 - MP_Wait 264 0.000 - MP_ISend 88 0.001 13472. 1474.84 - MP_IRecv 88 0.000 13472. 12454.02 - MP_Memory 264 0.000 - ------------------------------------------------------------------------------- - - - ------------------------------------------------------------------------------- - - - - - R E F E R E N C E S - - - - - ------------------------------------------------------------------------------- - - CP2K version 2.6.2, the CP2K developers group (2015). - CP2K is freely available from http://www.cp2k.org/ . - - Borstnik, U; VandeVondele, J; Weber, V; Hutter, J. - PARALLEL COMPUTING, 40 (5-6), 47-58 (2014). - Sparse matrix multiplication: The distributed block-compressed sparse - row library. - http://dx.doi.org/10.1016/j.parco.2014.03.012 - - - Hutter, J; Iannuzzi, M; Schiffmann, F; VandeVondele, J. - WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE, 4 (1), 15-25 (2014). - CP2K: atomistic simulations of condensed matter systems. - http://dx.doi.org/10.1002/wcms.1159 - - - Krack, M. - THEORETICAL CHEMISTRY ACCOUNTS, 114 (1-3), 145-152 (2005). - Pseudopotentials for H to Kr optimized for gradient-corrected - exchange-correlation functionals. - http://dx.doi.org/10.1007/s00214-005-0655-y - - - VandeVondele, J; Krack, M; Mohamed, F; Parrinello, M; Chassaing, T; - Hutter, J. COMPUTER PHYSICS COMMUNICATIONS, 167 (2), 103-128 (2005). - QUICKSTEP: Fast and accurate density functional calculations using a - mixed Gaussian and plane waves approach. - http://dx.doi.org/10.1016/j.cpc.2004.12.014 - - - Frigo, M; Johnson, SG. - PROCEEDINGS OF THE IEEE, 93 (2), 216-231 (2005). - The design and implementation of FFTW3. - http://dx.doi.org/10.1109/JPROC.2004.840301 - - - Hartwigsen, C; Goedecker, S; Hutter, J. - PHYSICAL REVIEW B, 58 (7), 3641-3662 (1998). - Relativistic separable dual-space Gaussian pseudopotentials from H to Rn. - http://dx.doi.org/10.1103/PhysRevB.58.3641 - - - Lippert, G; Hutter, J; Parrinello, M. - MOLECULAR PHYSICS, 92 (3), 477-487 (1997). - A hybrid Gaussian and plane wave density functional scheme. - http://dx.doi.org/10.1080/002689797170220 - - - Goedecker, S; Teter, M; Hutter, J. - PHYSICAL REVIEW B, 54 (3), 1703-1710 (1996). - Separable dual-space Gaussian pseudopotentials. - http://dx.doi.org/10.1103/PhysRevB.54.1703 - - - ------------------------------------------------------------------------------- - - - - - T I M I N G - - - - - ------------------------------------------------------------------------------- - SUBROUTINE CALLS ASD SELF TIME TOTAL TIME - MAXIMUM AVERAGE MAXIMUM AVERAGE MAXIMUM - CP2K 1 1.0 0.003 0.003 16.366 16.366 - qs_forces 1 2.0 0.000 0.000 16.261 16.261 - qs_energies_scf 1 3.0 0.000 0.000 14.562 14.562 - scf_env_do_scf 1 4.0 0.000 0.000 13.643 13.643 - scf_env_do_scf_inner_loop 10 5.0 0.001 0.001 13.643 13.643 - rebuild_ks_matrix 11 6.7 0.000 0.000 8.995 8.995 - qs_ks_build_kohn_sham_matrix 11 7.7 0.001 0.001 8.995 8.995 - sum_up_and_integrate 11 8.7 0.004 0.004 8.711 8.711 - integrate_v_rspace 11 9.7 8.583 8.583 8.707 8.707 - qs_ks_update_qs_env 10 6.0 0.000 0.000 7.943 7.943 - qs_rho_update_rho 11 6.0 0.000 0.000 5.934 5.934 - calculate_rho_elec 11 7.0 5.683 5.683 5.934 5.934 - qs_ks_update_qs_env_forces 1 3.0 0.000 0.000 1.052 1.052 - build_core_hamiltonian_matrix_ 1 3.0 0.000 0.000 0.646 0.646 - init_scf_run 1 4.0 0.000 0.000 0.592 0.592 - scf_env_initial_rho_setup 1 5.0 0.000 0.000 0.592 0.592 - fft_wrap_pw1pw2 120 9.8 0.001 0.001 0.547 0.547 - fft_wrap_pw1pw2_150 54 10.1 0.027 0.027 0.500 0.500 - build_core_ppnl_forces 1 4.0 0.348 0.348 0.348 0.348 - fft3d_s 121 11.7 0.339 0.339 0.348 0.348 - ------------------------------------------------------------------------------- - - **** **** ****** ** PROGRAM ENDED AT 2016-02-04 10:13:54.586 - ***** ** *** *** ** PROGRAM RAN ON lauri-Lenovo-Z50-70 - ** **** ****** PROGRAM RAN BY lauri - ***** ** ** ** ** PROGRAM PROCESS ID 30677 - **** ** ******* ** PROGRAM STOPPED IN /home/lauri/Dropbox/nomad-dev/nomad-l - ab-base/parsers/cp2k/test/examples -- GitLab