From be13b968df5367d878812987feba6cffd8c21c4e Mon Sep 17 00:00:00 2001 From: Panadestein Date: Tue, 8 Dec 2020 18:56:36 +0100 Subject: [PATCH] Fix nuclei formula --- Makefile | 2 +- electrons.irp.f | 25 ++++++++++++--- jast_coeffs.txt | 70 ++++++++++++++++++++--------------------- jastrow | Bin 48224 -> 39776 bytes jastrow_provider.irp.f | 2 +- nuclei.irp.f | 5 ++- orders.irp.f | 4 +-- random.irp.f | 21 ------------- 8 files changed, 62 insertions(+), 67 deletions(-) delete mode 100644 random.irp.f diff --git a/Makefile b/Makefile index 24b8b55..9ae4c71 100644 --- a/Makefile +++ b/Makefile @@ -2,7 +2,7 @@ IRPF90 = irpf90 #-a -d FC = gfortran FCFLAGS= -O2 -ffree-line-length-none -I . NINJA = ninja -AR = ar +ARCHIVE= ar crs RANLIB = ranlib SRC= diff --git a/electrons.irp.f b/electrons.irp.f index b83b836..38186bc 100644 --- a/electrons.irp.f +++ b/electrons.irp.f @@ -6,6 +6,14 @@ BEGIN_PROVIDER [ integer, nelec ] nelec = 10 END_PROVIDER +BEGIN_PROVIDER [ integer, nelec_up ] + implicit none + BEGIN_DOC + ! Number of alpha and beta electrons + END_DOC + nelec_up = 5 +END_PROVIDER + BEGIN_PROVIDER [ double precision, elec_coord, (nelec, 3) ] implicit none @@ -22,6 +30,7 @@ BEGIN_PROVIDER [ double precision, elec_coord, (nelec, 3) ] end do close(fu) + END_PROVIDER BEGIN_PROVIDER [ double precision, elec_dist, (nelec, nelec) ] @@ -48,22 +57,30 @@ BEGIN_PROVIDER [double precision, factor_ee] ! Electron-electron contribution to Jastrow factor END_DOC integer :: i, j, p - double precision :: pow_ser, x + double precision :: pow_ser, x, b_one factor_ee = 0.0d0 - pow_ser = 0.0d0 do j = 1, nelec do i = 1, nelec - x = rescale_ee(i, j) * rescale_ee(i, j) + x = rescale_ee(i, j) + pow_ser = 0.0d0 do p = 2, nbord x = x * rescale_ee(i, j) pow_ser = pow_ser + bord_vect(p + 1) * x end do - factor_ee = factor_ee + bord_vect(1) * rescale_ee(i, j) & + + if (i <= nelec_up .or. j >= nelec_up) then + b_one = bord_vect(1) * 0.5d0 + else + b_one = bord_vect(1) + end if + + factor_ee = factor_ee + b_one * rescale_ee(i, j) & / (1.0d0 + bord_vect(2) * rescale_ee(i, j)) + pow_ser end do end do factor_ee = 0.5d0 * factor_ee + END_PROVIDER diff --git a/jast_coeffs.txt b/jast_coeffs.txt index a064ab3..5b23c66 100644 --- a/jast_coeffs.txt +++ b/jast_coeffs.txt @@ -1,35 +1,35 @@ --0.00008444 --0.00003393 -0.00009230 --0.00006224 --0.00007552 -0.00009380 -0.00001181 -0.00009928 --0.00007969 -0.00001493 -0.00007319 --0.00004341 --0.00002337 -0.00007838 --0.00006033 --0.00002512 -0.00004281 -0.00004631 --0.00008111 --0.00001758 -0.00002855 -0.00008134 --0.00003731 -0.00002581 -0.00007723 --0.00007482 --0.00004745 --0.00006915 --0.00008698 --0.00001448 --0.00005086 -0.00001145 --0.00003374 -0.00009208 -0.00007611 +0.00000000 +0.00000000 +-0.380512 +-0.157996 +-0.031558 +0.021512 +0.5000000 +0.153660 +0.0672262 +0.021570 +0.0073096 +0.002866 +0.571702 +-0.5142530 +-0.513043 +0.009486 +-0.004205 +0.4263258 +0.0828815 +0.0051186 +-0.0029978 +-0.0052704 +-0.000075 +-0.0830165 +0.0145434 +0.0514351 +0.000925 +-0.0040991 +0.0043276 +-0.00165447 +0.002614 +-0.001477 +-0.0011370 +-0.04010475 +0.00610671 diff --git a/jastrow b/jastrow index 51f88277cddc49658e22f1f43ef7eb789d236655..a86a7df7803acc5fe6f5a3485907e153ea93009d 100755 GIT binary patch literal 39776 zcmeHwdwf*I+5g!jM}x50fW+cub)`XrkQgq3M4=?GfsGA{2#SgXMcUM5+v|f|}+M2+Gk`*Fup6bcuw5 z{b&pxPcH*1IhsegLMc}$$;)TTNB$b7`&=fn-QcsV z2p$EM+{KWioj>-;mW)$nc|G}MkF;MwS-f^`UG1DvBhRg?$*-%eZ(f|gxOi0lsF4MW z8VZJsdQ(4nCtWyQl*IOyssz11^2CvT>NS^#?>n?8GcE7Uhwi>(&(3kr?%PON^@HS( zewrcFF9r9w9`H|rHx+*Y@U9ekhW4OmY7h8e5BNXzfTt~~%KcXl{9}9I-voTBdNuce zzZLoif|mTG_fYOhJ>d61L8@{Y&Z+RB9`N_|fN$-gTn+I!Z*D`Azp1L8Q8-fbR!p1h zt?@Pa=G89p`13N@0(UNr_QH&z4I0})O#2CtD5{?ufdv{QR(zj@uIoDCU29k zs?K|DQ>`CrOloyQeRY)|3Y%){=UMQJ{N%6cj`7yOm_&3#qpv;@BT7l;G!3bqU)5Aq zjjDCC2BlWlH7rWd@AG>b8vV5m^@~!JVMT$ls%f5;%Pg^~sj2EZdU_#UJWBH|Ze*_7 z$k8!%4b|F019-jFi>thIYwN4(YOklu$thjFu&TCRo7?2`X_G1|#!v7LFBo1>q)nQ5 z(bQ>E%P#beM&o0&NEagne`km>gJYeM<8SrP7y)er}>mxjN}>GYT)>g z&ZMHx`}Wr3scnOWz5r=8E*kP_ro5Ub)cD@jTd39YFu?26PhgvLnmaU#;7Q9?)2vUay-$yY6-f6+xE%=odyu*Tjz=Bt8O7bcTzORLU zwFPe-r#db86D|CkEO^dk>glrJa}*NL)4DIy6V!ckU7Dx8)1MaW^0aQqT(8B7*7;mo zY|sY$cAQa)6!Y^b8Ql}ZGiVL-jM8YA$Un|Jn=0BV@((f3CW@{S`9CwC!Td^*|0DBk zs%VGE-^x6jDB2?Ozh$0H6Kxdvo0w+^MrVus^~^K$qEke^iFr0zv{dA;W}YD!Efo1V z%(LmDE|I^2c{W)zTjZxQpT)c;@)t1Ari&gp3dA7F&+R~7JoH=S8xTC9qAri7y+BJN z*MKo{CY&;?GxAmuf$*l|3~ku{ion)cp1}K_)(;O%nN}InFMyOMWPg1U!PrMwX3obo zF6}GF8E%0ZhBle*3F(6YJtc=dfsgzrqp4m351Yz>F`!k(+vvz6j`sQdy3 z!@5rarh5h-g}i6&(R5GXfM;!FoX57&v-PNd02D5j3j0uDPx@NMJ_tQp^SPeZv6n#r zDq&9_PjKuN$VJ|PCn6qHVWYkjd^VcXLc_jX50!3%{i^6(a2WZrY11$9EZw&Uz6$D> zL8QI2`FMD(4G2%*gUZmLU(hknQ3->>6WHwu49yL@bF&%Nm4TYv%*Y)-2>%y?bIkyQ zDU!}M3G9wMLt0?BU{sO}{Ky*M}}y|C#x zlyr&3D5C z&@hMmQ-k*JncD0x4(2pLzI~hj^q{^HLTENeWH)>s)p=tf7@I~-rCAQC-Dazm#_7-c9MC_<4OI&3*CWLy{c^o(n?`b>%HPw!LH$K!J;9t60PbMU z>+mO9q5}g~WzacT(8gYYoTWt;09zT%|FDXzMMZ*j&=%%wuT8eSQyHpUmx%@|3yqIm z6msv&th9aVu^p}q)g8^O2+iD{?FrqmG5exW?x~*8$X@aOS>C?wMlbx4omUxnsWK2b zIKh#f`v!lXw#V>mHUe`U{n{TcoB1LVeGVY8-SGCzm*LP({A|O|X801Gp8>yK3|!yw zXo_|iu08|Z zeD+eX1Bw#_Pr#kqh%x@_@n8*n4)hJ2qQY&t8vru!;|a9owh)H&x4{HTF;0qvcR6&u0S8nIP|#$E+iuJp_TbJVSEMLTcFN%D^Fr)4aly=>jB5 zMG;8kK?fRQy$uD!aCfc;dSYHGaOW0R1dsQGGP>lz=gxISHhjY&@AJq^mRK?~H`CEx zg@{(>33H%!c01GpPNDk9ZVD~kK>sI%y=8th;u(!J08GK0+i5{?>_mX}ZI1SPX=FRb zwEqGE4Co3CNU;)*wygj#dOQ%Q{0M+y+dS3IOPOb-9G-P0>zhvn%V>io&q8TEobK-* z$%29HyD-LWo#lc~Odzr9hSOCTxx@NH$O^A5-9Q^-gP zdJw0Dk!wJ~66i(h#az$TT;zS?AoWMU=*^h9#QNlgLdu0-gzIC^u*KmB_f>|cA{S9a zxM&bipqvLF+9u#U-@v#OCaf35lf#jeE;^Nbf&6`_o1@*qMT&IMTfabI$u4RTwXnMA z;je^?-X?d*Zoe|Mi*|s^>Y}faXJ!9EY2~6uD387Y6FMyY-Pi}{FfXhR!EaH$U`~G| zF{I}ZUy7au z79uHIZfYDm8c7x$`R$i*MOZ%<`Ccdz@jd{_=t9u2K2I7C+Cjp2aOet-2hB}lEMWQ6 zAJ(gpqSYc$9Br#NQv67 z-`k6H<;MYFqFW->P-evXx0j=wu$~2qFfibBAPLM@AP7uYe-lZ1*U|oSq21{9$B#h? z$sr9KinWIH8()Z!E+KcwK6M5|+66ArPt9JR27|02J&*<>q~D^n3hA*Bi|)cNH9kfj z1{t-z2RPbX3`0fzT8S5Z`DCMh1*$dy=eY)^zruv|Z()>CzX_Hud-nM$>NkP~CDiZm z0a3rp$sMvcU6#6j17VQ0e&cB%>h~R`RsBW^<-3hhZ6zoe)y|||jB4D&M*2vp4WLE2 zFA~uf@84nlaa4uji*DxD;_;pbTyX5|GZ7;(M_W7EHCQwPtnG&!ZS$D3e~V~u-|T3s zU@m6?Y1_C=tda`D`Xxw-fC!FlLlRTcX2@`gYI6gRFutAWA^Ob#m@LiUZ5{~PEe7o- z0!7g?k#n>?#HzjbxeV`I;P^OJ+=|~p{X=+^fi1?Q>Ii({S_mGgcswzezCcP8j*%ju6no#n?~9lfRQm)27FSO*E);pQz}?)Ecs<^-$4I z!DX%J3DV?slvWje97?18*ie|?-9fuP&83yIGd}`?JE(sIE6O>qa|?_0FIoi?PgAgH zCmiv9?p<(ZR@qF|o*m@IeWq%+YE2rpunU7qF4M6k)iqw{S0GY^x zpT=VDU=1-MQmNVX87LHM0d%94`$vZF72!JsHV@lu9vLl%Fytb+q)K_gPbHsQltd1a zGgS02!ptIHDdj&jAbZLHQ=aB2BO;F*^1Tz~`$+lwQ^|iJx=#7Q*7t3JF9vVj8ZFs*L+|pCJJVLS zhE`%S!?D4R&|;grWOws(vNLwS0v;U*1H%8hJuQLfPLsp$r#p8L`oDt zHV3Kp7^%HY>JL*SwJ3&E`yo=lM#?b#6GgS(auKK)a!yA02Xl6jcblW_2zobG*nX%B zI@bXymJ9jE5v$(}ES3w`0u0+0T3{R&a9{nRSTQUy$c+T%is5<`<7j&y-ikc>iCi)C z0|y_+vX}5Xs9#D`uwr-?IZyRFmDPJNKC>QTLmYS$yJwu=z2d>%9GmFSm{%$6A7L&l zmU!;+d~?LJ2EONW2De1_c6v&-I__BKc&gL!)HYj}qy05jqCF;#;ye!b;fMqF7%Po- zRy3^t5rB@xbn@Q@rbS`G`t34_S@I7g2AoR-$e#k>X#1yd(byxXO3=9!zlG!W33uvt zat_gt0u8q%ZbM#lw4a4S;dC#hs-4r6@b3`o@OmTjL-ZGI<*q(Nel zHs)jL@)&S@)FSPJ$%e~0{`I_Uc@t9BWlJ+F$z{v-kCf{k5{1nqA*?@*jC9}{B!Rh4 zVZ!TlVF-bGfZO4xu!7NB5pM6&{PK8|JUxDKe^ zkHUTnoH*KQ-;XHajL5{ZLi)45&vbBb8-={2~K$io%5T82ri@ zbhHnlS(v1F-g4EsaL;osum{>TdEWAOB*o4F69HES#gSPk3g<1ilRIQTZ~>Q6-OpQ= zfy=rxd4L9DPMJ+<>cJdPYOK@51XUP0jvS$)E2$StseULp(j=wM0gNtDJCjK7Bgh`X zFLS+n9k(vKfZ@oKn`*G4s;kF4q(U-1JgAq!WftyfgHrw0?2-~&cOU#VZ!=yB#n?;jaE`2_vuh5 zPY$_N$_Hl1?R`)9Vqm`T#eaEH`yvMhS$$DN1L2GPlvciYjJ{YVB8hRb?6dE{RGohv z{M)|Yv0^u~Pt#C5OW>ekIF8;#cZBu6sDc>l@2UvA!TsSkN3bufbUgZ^aKolj$Gsal z#{Y0ZXv~Xrm#5?n$K4p%JMq3$`T2~Bl1|5c2GY^Khl8hev_BKJ8cw|f9Sem%0wA6G z4+C?n!i4p|t7IVm0TKhwTL3Uv%>i(<_hV&b=WZE-l9E0Bhi;x;x<`0gaEI(CPD<_R zhrwm_^vlT8)2)=Fb)k-G+xjK+EBwUp#rXOqw>WQx%T9 z@vf-oI2shPUtONMqG>SATG4^hAUcM_ zZdmu3OwPGv!rEJ5!g?W+qPV|`2>_eAh0R}OP&)3x9%&(dd0(b2qnp%ukb3ljk+N>yE0U>jvu`AVOWKx!%riPex9yxMFV&O ze!M~p7Cpzbcq@z$z#HK*8+LQKcsC{93Uh0hXUtpW?fY*W=LsS|Iq-31Xv_dPKEo31 zq{CyUabBMsy1}WzRq#(nZg=-g&x3n9c3IuCL%JuE?#YIGPHSER_v9tGrz2?J zIl-u}v&VQ@%|;F8ECoQkoMqJ4`6(g+)(ma3O9^%f&CT2e$;zYbYuycy0qq z1y2jp#}1wj#?7?Uag)j~W6y?P7QZFF{=mAZo14s(BQoX>}Vcfs^X*&A}FC8e=eBB3Je=; z9>4E2Y>?84$Yx)wQ;{@|&gFOhVd$KTgO4JV@* zYYojK_)N2grwb$H5q!QfE8r|Z5;Z?jVZwTzN(P*RNenpq36TFm8D`d#F=c{(FVfOe z`;Y`?tHNO2rjh~YCnN@(T>|9)i2w_nQmkygIaBMlG@CMA12W+V{;Zv1#qtTcL-yja z7zR@u!A}8~^<=scd9FS3D2)}1ID)@aDBmsgN4D$$N2utZ)QhEAf6yZ5NvVUNMGN6} zSe`|b1I{x5VDxyw?z9ofdZ6`REk)Y8#Dp)1vq2-9E3E%p<*az%{>}3@_);9i{{^N7$Ckr1 zoVRgso_Z#Gc%ifv14Idu7z4hbS<85fG^I6u7zwNo^*(lyxb8~RmyG_;nJ>h2cehwK`e^IP7jufZvGb6xd?K79K-sm&EDXsd<-9ot>93I;Wj!@BNTA|+0Ii=JN z(4qs-2C#f9oDZ;#Kv<6Vb4X*LI9Z71ZxRgM_H7}cg^TfiE}s0IOqw>`uPPjIp-@y* zXbsu>^-$3Zz-6uI)vN|8I*`(;qJM=_F*qcQ2f?Cy;C}XQc8s8XG;^E9Xy8Uw!umiU zjHu9uAQ{Y=P7d_=9AMGoXR)T$&I!PZ9)G$(>kUK?(PslB^mxwj3k=#!0@dTy4Daod zy*&+7KF6Hl_n((I|A|Ou=gx#0qjNtG76jt!07f;O+YQXq3KQ11sARzT8i@htY60>C z0FJhP5De;9AuT$0AAH>t)-OO-C`{xWe`;QW@vfb)F;^8W=89b$CZXW~rHk{Z_k zq%dLqDU}R3|4m}Rxl(}qUlGVId-f?>_b$5&Etfc-U%Ev^@^8o;vgZ}C%MJwB8Z3n{ z3VikgUO{V-7I3X$O1@HGi992D2BlRbzYRJaC|;=_eH9#`qIXXgIAcD)T1pKBjE*e2@ zghA-C_nr)cLiTYZQ&;p3a9Jz*Z&m{pZJ@NO=tofcALjGm*sI{m`1dr|!Z5jBZ-WJR zHzzH{Y`y@=p#64onDh7Xl#b8eg`?qwYG*TeL~kxd?a^b)h#sO32T176oWI8!v;qRv zoAoD9Z8?93UXeYzYm^K1VUNBQzk_-Qc>eSJ{Sc(Avly2hV*VZky+#-RHCWI*9RRX3 z++tv^QJAnERLMa8auNg1W&q4S69F7;ywe)gPeEFA@#ubb@#B#d3KP0G=kF2>tk(HE z+6xV(=5Nq@;aI`?JvzCQnDRnuF*Vs`$O$jJKs2}yY1 zW`zms|E-dN{JTjEINJeG-%A0aeT~jP%VctT$%OAv6(+1-u9AWLD@Y7DCjdZ0f97Q8 zKj$PZvD0wA7E?gNV*c`%M4xz^+#!4KVX6B>E4ZwE;xEXvPh3f9)hEtCmCep?BS)y{ zWa?G(^+G8%2r!ywbpE9XHB{gx0y$q_h%^R*+8m1JFBS~l=IgOQ3m1p=bK*(%X<|*A z#NtG1;QDBj&~>%UgX=GGu4Rhr1%s#2A9gC#9QKAk1MXp23M2JgNyMabuCw1XB%Am z;w273hm%}KhszRe2G?;Q@qw=0ab2sFu3cf_iB6dZ*L6@Jbxl>eo;0|wka=)@8Q1lS z()BxoYrE3bijR-uTQi6FKGM98gfFKS2Xmf99dKI4D{5=s#cRo;dFP_^u)bQBgvHv5 zU1=PzJ|-cU^GpFYbi-o7vCo487hd)Xv~!;T`F5z|Au@*Qn_vYmAx#c7XX5~Cs3+9t zMd62ZeH?AaU87Qb{XB%}kQ;Y@g7$gP%`<3MXd-5{n&0BG&nBSf8!R_vhAuP%|4Ed%M-jIcM2~? z<$RJ2ROp6c9If-_35c_K1~-Eo#h%uPz|ZI5xHuVhhi2HwG=RdWx;uA1Z%1X%Epz1k zz$5Y>Jb{c{;$2jTnr7tAEqCO-*m^^5o{d}wDJcrdWXg50+;RB^682G?zh)h9tT=+Y zls%2cKm&Bj^EK*usVsEe(P4)~`?>?0%L3mG{-XU*;3c=CeCyWlFyXG4h&oc;=GJvK z@%aWed)&zy-|}>^kWI~B1va(ru{qlLIWd{v;;(J)z}CTQw?^G1T}x(HRQJ!FT-}8e zuN@xSTNi}f{cM%CZDqKaWZPB|a=)nI7T8W>ABkIC>o7y}6n6)19IXf3jh|9^_X{_W zF2_}`8GD*{8Xp$DfI_v_jW*W&Wp~LIhq!K0*1FbK&f8u&oI_eO;Ja2hP403g}^YzL3ckMnxaXE@rcNk_YtRCyehAbDD?n*aJ%)zTd|8 ziPQ6$$3o~4k=Jh@%6zB#J~0b{Z+xFv!<6_w@d}t@e4p3_Ww@D`8`MuHj^82{;S>+w zEFnw3!D@7L(f(WyU@WPdnPbR?q!rY?lrz2E5i*A5c_O#tO%f3sF$}CyA1!K zReQLnSn}F`^(t@Mg}{E%-;2N5YZ%AiUuzRxsGZSf4=$xAY}u8PnNo!OHEU1~f*{p3|ejw;!?oA~>*1Lk-G-y>ZA_xBYi+up3e zM~EF?!cC>t^G7+qVYOvL`ycDRVl`crbYF3QA2`Xlub4%0%KM5t!34^Wysx?v_Z1%@JjQ*+dD2DzEpjfIf8>2d@jfWWcMht6`$+&8J>&NkCj)T;RqNxBV(6Rq z6%SHmi2eZpwu{xyvzS+VL*u^UCkD$K7K(9San|#4ySNl2KD@73g|rxCjr)ohBFh3% zk%Y?*#{7went5OGk--rB(f1WkhhBC3!c(Mm&%s&Zi1;w{U_iC2g7JO;NnpMNL6{QO zGZitc>qw^4z<~2JN(S;@05B#{wRQR$d&Y$G25y}^Ldbg5c{*&saYBsTA$wykcNAxU zEAc9V2VB+@!1=TsI|>)2)sCVQ;|Q<3|qM!0APP@sGYARH2K?lmx73KP~piYMPkQo876@?oC}5WlZj zis52gCJKA51x3mu8M^3VT9B~2Jo1ci(c9z>*(d)zwTlM9Qmc!~X*pc<1ErOVrb2o0 zeZ{bT#o5LvZjAq7eJb*Fgq**`-x5&wd+kqxDBg#1iJo&da_BuXh(zyk!8M4_um?x+& zgGq%v-;o|i8yNc=*QkcQo8^b~#{gu*K5k&{P#Ao7rjh~YKS&HX{~$pAJOD@Aeh3Ej z5lD+({2g9!gmo9PLSe$!u(*H8{`(czh}F)0gGB#bkA~`gUW(VcVf_SX5fy;-EB2$h zD7znU06zyn6fHkaAax=`_;wYcwiRj<0t z7Lj)gO9irL_83Yqp zNzwgYsA2g=e-ur1D)osz_yvXsN85MwU0Byt6#`BhlE8cnvZ!KMKcbQW=l4j$TOSFK z{}e!UilO28IMWNHhV}au2H!=iWWc$G#DMcJ0_0=rG#m}z<1bnyIvPE=v{m^JP7w8v zl<htRQ#DzTk}1eyj$U)sL;n(~svMiVVdTAkgCkVrpkDQ|ce0ey0i*Na zWn+J!m!B!u|F=OyaO{?Wh`gBp`e4z_pCcu%985(Dsxalbk)kJq_PJcr3lrc5tDn9)zQn^pKO5E2wQ4$jO3Af<%ckh$&b$`{nu$4*Z8XAg%RT#;ut>$!- zBZhLR`Ky!pYibw5K3TX~mbFGrzD3nlbw0C7J@{9$d40Ylf6;4(!&273l&epYi~IT5 za=v=6&uEPQ$X9PRnl;g)G9$$E!8jwoC!t+g3v`zHc++YcaYhWPwlySE+7V|?tmlZ- z4JFu9L?7|=rtQ@9^LZLI3Tgl2u~)8KJ{J3)5L=P~8Q4Fqi^Zk^P6Om$ zZ*y;q#U23ck0aC#fM4P8_Av7}OzDGRcqE_;a2nu9zz+dEfVnvMpGi7kBcL0w4R9LZ zoq%8AgOaxaGj_s0%-kaZO96B7%~S{Av^Qe0J%INEeg(MUP1N&v=m(qu_!VFs;P!W7 zu}6rGfVtnpKY;y>Fwtz+Pt|OT zvu&sLIX<()cD(o_W~YG`c`OziB6$c>^FtqKYrqAL-q&WfXyZ;f z|C|xI1A!ItZonzfi)xuqAyPB(vkdX`H)6+UWiL-F%NlTVdRdmMHDgj%OGersvRuHF zWo3`gqU-?tbl|5ybmaqT?yT(jf`Y$0jO-)$i9p}yq{a1lsPAIx+tOX%9Jh`#|H%l{pV>F8T2%b%3hnUQw? z30ZkC35{TwQPG=i&>#JFHdg*WAWfg8ErHWeS?VnW&5gdh6|iT$@N$CQ@xAG%I>m61LV-cL|CkveE=C**Hi8;iX}jkNVT*gB!NXoI&w^I#l# zOZqc0{yH+!Uox6N1t4wAh8-|R|5F*W)7BoV4C>MZpUU9;vv^x9b}|dV z!*V$nwE>3%aq+MZa9(Nw?sc?F4{eG*aA|TMn8Y@v{07LoFn?vEPCd*2I$3_alFz{W zb^!C)zm8eHCt04upD5oA`9{ox$L_-mlI2DDr$I*>#aGgv z8qQmF=Dd~m(FtnCx~TWWtW0zR(XJlI?E55fjF<`B;Uruma2b(U?0uGvhvR-5aC?CJ z|1s`&L6`f}V~+b>sm6WAW*&!?$FS)%< z|N9}|aUd3Z0^-8|{}uPY8>rhrFZ?1F^O25;Nw#;nIWLQO(8GCoa#m?ms6_~aPC{LQTnpV9`GfND?CBIVUNDz;zZ`I&?f ziBR}FnNW6jNp+oLc4B*F_(RgaO`;px z|Ns3~{|{WU#y`l4{Y}EBB;+S2e0ECsv4n>tJWe*yNfMqV;V20wN_eS+b0xe+!kZ;D zO8N1>)&EZ<-z|FMuj&$YC^}!{iC2Y^tnl*9nkl64=F2#fVk8vZd}W1KPu5pTspw*V zoK&tUZxselclksE6y1C`lt?z^b_b;`?YVGX4|scC8~tJf~|}+RBu8ho;8W zL_OJBsT^+;@qNwzUojEiPun!ZiZRmZ+KJ}>8)>Dc<86={pA+#p_EW5U4ee4W+l7AQlPmEaiFX)iZMei|OME|xFYQ79r9I$%QqSz9 zau@Z$ALs#pm(;T-$)3MUd}k8=Nr_)6@oGJ=UgBGl@I1Rt<+n%=`0snbXJK4RrRVG( z@Fl>rU0kwXC_hh^5#ZCGl!KsEXYz^mHFjwqm??r5HE1fxj>M-LhtvP$8Ft+SkDML;otg(9mpl zV=SPQ4^^3jeSxQ*m425%d&F-|8z%UnM~)}LY|zRif9E2DUzIzr2mLqqfZu_7(N6XM z5OX}@bH6NiwjBDzyBFYB%W~(-`l|Kg|46(`+EXaB;O_%T{K_Vmz$iQU{{V|}CmM|6 zT?mqAO1#TphtG|t5O|fAC;`Sv{+96uR=mSO@*?2r&&s5B@k{;#a)Z!fGvrzY-bh=2 zpD+B-eZKp15By=Ne^s`jMZ9xCa;GZSW#CHY42_ z^xR|s?FZl`N@5IfGUBi==1sMJQ1!|a@tMFXGUo~Y0a~8yFN(if@MAumYREhxMX!&ssK8rIXqt9_lq zJAm(}G8Q;Y;1QoL11W)4D)FUB?e$BEZ}b@a?;Ak7LgHsv7zj1aTrKgPGX8y%|9Xip zbr}c=w1C7H78rQ({R5JB0WZdxBtNf_{9T0xzgih@miPk_uhwbXC4Nu2!7tv4BN>tU z`y~1QsKnnX%Z1y;!*9~qF4{07FTQ0!vQXknxnRcgxdF5Z5+A83)fgdjMO6FO~-k&C_O--Mk@2apNSUKzq0== z;8VqU1`gzCe}@bN@tqoM?GHTdoSkIP7|Aa#CxOo;>vf^z?@aPrt;DaA{#SmxN$8Qn z7Vub_(O&bfQJM|nN~xzbY2DQ+^kDwYG5Ex{YDm5+@h&-@qx$06CGj1Ad{i9rVj=7M zM3O%>V5cMoJ*uQQH532q z_*z~qsJhPUtH;#^{D<*X3w_?2=7kHdgNP;P1tV*w8&G^=yw=pET`wfbe z1dAY!JFQnLaJ$1hAJr`wF>2%}EuJeHRbt7F9$BPW1TE@tIly|OA+BK5utLp3w1~(9 zPV>ruMW(RWQs}7B2__T`v(So%jn=GsEb6Jw-_Tq=-xvQH-&g15F|V)AEHn{g5$Ew} z;-~y^`-)0NTWlCLx>$VVALq8H@p>mtEt_2KEx*w1^`ea@T{zua?vY53d#dK0RC&?( zvP$no6DM9$KFvF=YOG_tJR*awE9gdhUg19CUp6SpS3-n!a3I2&KoECTpg z{PkY5MDM)04Nd;0s(SAtA0v6*B9uDuqN&rSmR;x_ty#!^Zv!sfBv9jL*laFdaG*=f zwlZsJ+&sc=}bkK+gw){uhzmU-nfEU3DG2?#c#*EMGN`HNcX2)ZMLwc1@@+6jw}}?2QlW(gpgyf?PMr_5i_pq7we{ZSMZOx0 z%|14Z#rEz^9k&_pc;bFhl~h)X8YMM^h`}H0a>H&l86{fnt*P=?8R2QRXFLQ=lo$_; zo}|n(txp+}nC9wxip%PT`syk_3T>*bpO+ApUT^i{DvlFXb+~Ys=syv%qU&{UUNQnPJm&@FhCaFxSMzi^5eihH&>m zIo<}A!jzmEK2$!{m==HM*v&}ZIE>d(IY;#gn1iS>XQ7%!4c_@x^)-z8_*jH-TCRl> znpEjBMonWXF@{DqXuM8>BC^smmWhM2ndvFR--(?*W$Vc>woYP}$Y$wA*3_(G0rOUh z5QB$OHj61`>M`2@4as8SdC|5akm>If0vre|PBlmH>iJboxF?ITJjHWT|tC-b0ecH$#CIE^=ySvR`)Oxy-m{!7G1Z{Uo{8N z-(a71Uf;4}Avo zHyM;`uoJ{3{Dh1bv?gC&6$NCf5&Mz?F)FesiA;rUnrc1-|)mXFR_K zN>t7eH5SH(gi51?s)e=Hu!$j6fYDPc!1jD0w%qCn;C}=8J(~?bIOYoBDf!hhpWj}FzuE!63s^r!8LJD?* z#wXc-s{U=DQGPs`@Tl*K6l{?draWJbsPdKk3gE~>`}oYqxtJCH%*dL*S(fz)AY+wR z-zh1m&f%0ywoE8^Tavu`o=HJvuaZ}J1@B3cSLb{R-s(}322|xL>|c`P)%Q^fN_QxB z;cGL4#Dh{^wV!$qs^Dfey5&hO|1r?$KPBJiGJ{#cXJz{-70Ld8AxU0+m!+U9sRGJx zO5bKAaI`66lQ{KVg=-;m@pv(E@x0Qr0|dS*P}Nl5dFA~``S&2J`UcBT--8uOd29bC zU6rf+cqedlkt$z(hg;e}ZcXbHzvI7@yn>&CWR=$%T?TWb0r94Y;A{M_h{8BZq?G*NB)l@m3RUb_S!c`gJLS1O=NUJCRJ~RI lXLu;vt@5rmBNtU^`lI+^dV!O^LbuA#LYeX}ko^Vi)U3%{p~ICiq6>vD7l%dJ=C)~j+FW-5Pa zsQpQC{MoAV?XU=Wd^994pBW$e%8B3IAd)@cL(1_Q>T;K(9QFKnAMN&;s=nU*GEeER zp{m|7s=8{?xUr+E%d@MiY8soeoASnGj~koQP@8jsuzL^`VV{~jZI0vT-+sJsX>ev! z*5zHtMwKr8=(B{`bE!9Vpgi*MA-}3l(3z*}WB$MMOajxXaTv>cp`D1-e@2LG6}rivZhjRWmQdS zb=3;`8&RW7`FmMuRgI&vzM{ebqV^yING=J>h>a6jWhlYa7^bdkEG=JF>RUp>XcEdK zogxG#@Aadq(ADBJIEZE`;Yq$8k$OW4MmP zOTN70Bx-Q73V?ck((}neH3*+n^N*&dDf+q&RDws-b$ahFB%O1Em~Um$(Pg3QGRz8F z=z1z7{g8#;p9J+Bv(VMlD{=;RNCp@MPeh;04vCYlV?8OuI{4e+Af5AoenxjljC8h< ze)2jbM!I$0E40v=ryoy;#7I+15$k6@x?K-Cr7}HCb*WDGt zGxRN{Ipucm66xJcbE@s`6zQEzbBgVLQlvLC&8fA!L!>t_Z7|&~(oZqXDYm;=q#t9N z(@J-pNI$?dmk`|xMfx75xpe5BA<}m+%_+CLP^5p%G?x(Fqec2=raAR?XNvR+ra9$y zr-^hu(_C6~J4E_grd>=QIu63nMND(b?e0RFz6<{WbWbq%&;o~JBQ*`*1wwcroHC{} zeC>Dw;ZfH)4*Gk{KZ^aYE%5l?_q6OkG-Gy2(5RdS*9TJ$P9-Q<$9mHzU`u$U9mz56 zq#EBN$Uo*wM^oh>k4<%Qq|0?qGir2fQ4su2!JqpT^z;lr?(rY;Z2I&P&!%IE9_KdC zYsYF1t3XO?I6Czu~nEb*6TC42lsvchxG5E}-Q z8qB?n?EZt{zkG`}0X!g~U}j-1T1269*?Rti8`zal?gN+>K3f%fmM|doDff(!3N~hV zR&S-^S6$=2TKhXZ0}MP=f!sH4cR1RP`o0%P*-vbTZ%iP4 zCWN=`_T>eP8IXw1eTf<~%}|vC>6aIy&!qfb6li=Egv=@IL}!X}joU%kh!tiTjh>GNNX<&WPC`4L&iq{d?vQN?OMAGt>~VMfl&>=Rk{xDcVFQN%y}XiQ*BId^%!*zF1Wwxzj)mqtpQoh88`djvzU zdUf)cPB>4FFDe^AMsb~^L3HYXEQZ%&&UEZRJ{#CfPVIQvC} zZH#Z|F|g7#Z1Z8){LL*9=c>#Se<#N$)5sg_Y3WP`z8B5=k|%$2sD^ zuw8oo>K6~Yu6`XqFTr)4_}Pu09dI-qzUi>*a?l1&IED$J@lVBQD{7M;Zfs{PqGPwf zuWV=cT=y03IU))@GEYRI)t`hh`^Lncy&zyp3KnGzVL%whp6S69a~nz!UkJQi;H_Dm za635%d;G0g9{8MzOpiZA#$*Jo(W}`+5Cwp7Rg^VvfIwQa)&nNtC)AoXgRO^$`I(vS zW@q@kW06R}nEzu8l|XgYUvXOC0!+`5*%N#`s}o2`uwbvB;OVSpBoWg}{10Vq z2MV#f09wEdl)TMAi~NsgZ3P~F|BFaO^iPlftuZ@1BaFX74944Al;%GQL9~b*o|Y|6 zG~w@-iRaLCm;3P!FVkxCW{O3T^?2=@A~}iix-F*R0og9t?8mJpqg+ zd)gwd)*4z8Fm9rh5k?)rnBC~VPET3JZKSz8oAb9co&^@4WT_w^19Do9CHe-1&!QC} z!+{JwZExRH9WdVd5y!l%m90+`sg7x)t1Xd^T)mYhNHf{lTSa-Q)mBCCm1PWo%(k~z zJ?Zhk>j_?bmFJ5+o=qPocJBmZl%E;M7Y&p^xgnz!wI{P?aC`j5Y9vuMT>3er`i(1@ zq;jn=K46SuGGx>N2)Fn%u9K82G$mv#jiO6{q89XnjLXUA&pwqPWSj%QN+G%JWG){j zV2%tKFW(~Nc7Kk^3-My(Xci8=h*>)F99twt^PS)pn){8HfQEm2M5z2IaH`xkm5%Xe z|AJ|(G=f9#n9ZT|cvc!jL8WI!1_naL)5M|+mXL892*SXSu`-I@9!1vz6_$mJSwN|V zxw@#OqFk9h4~0sCDIQeoUVk>J!Ib%Exeb`9!T#ySQGnq;e1_EyF;F-zP8Vhz43A=I zuJSMt(&~`k{W5;RpjR+Dq(Pfz#TfLV$e=+V5VIQe0tm{WUq{i-DEd55Wzfw)V+^Vk z^}?V}pehWy0o9rYt!c3v^wHs1gC?OA4O_s{+MvZkCUk|^pGbxcF4BUz7eRb@pei&R zuzN3hM|c!4+xlAW-3sQwry(PN-@-(W)dwglU@Qd<3zuqw;0&iHf2Z#%biWyzb@(Qv zrXVFcSI8KpGeSlZP;?*;*MUHjD6r*AXXA(stYl1QDCtK;LPe&^Hj4bC{}i?SC-Oywexp^g%HGcR(@DL@iW&0qx47CK19QN*6Ioc(HpYo=k3zxN=!S8RxaKQ^ zRquFV|+S1-7=h`qgQSA7zVHjKs^{Ti%>5{;-xb< z5;I8+rhKOv4ZtyEhL)=j#*V~4A8jj%zY6Gg>O@Z9tG^qgR-#cEV4A+%9O|cGT5+CsQe}pj)@?aJ62~ zRC*pHx9#RmAV<{^GBSY)2L}q=Kr!ocJc`AMA9ReEMI?_fW&xe309YmB5->!aDRH9Ct*j_!O(oh7 zyXfVkFd9~g#*j4u>x## zTRof+b=wQEZp&>4m)&i3peeWAbZ)H6b{r5cn*a{6$oaE~9@#fD%|(v6x|t^^Tnt=6 z?73#t+YdkRR0;?PH!YWe9x~2DN_abDq(#x8QS|VQw)Md8AV!SFcc2VLTv7e$8xrEATK*NHk2T3phY!B%4eE9}OcY%`_`X@-{7 zfTm0_P`Jf{*!7ZO!)dH1nA-$B!(&yU62R_R3@tN(Rkm|IQxr(4fD0$$TxelYz*r7_ zC!!Zu1PVR{iixp0ko)vXU`Krofs}jL$9MQv1k#g96Y=6WiZdpKj6=Zajev0&7^Ex^ z=RoFuAbbFOB~&sJc20$_bbfDm7QCH)gAbfl}2)S-iWn&P8XCXY0k|}C#yaj}b zBhhv%(96U3B7}sD)v8)qnRz-|*>1+iyAvtmAK;5wc&q^#8?haXdb^Nm2rD)#ae?%w zfk_)5Gjs3LxyD-NinEheGxvIx>%qQ#KuNGVD~;RAf3PjE>CVGezoD4zl5sge37f1a zYX%PkJfuFJRR{NR7yBejKn9LnTeI>6@#?I>kj3pjIv$P?IP-F&%MJw}a$8xwZlG$E_4fQ3SS2* z``ynz7SJZ^U@x*>0(Xq8x6r}ELzS$lab+FqMb=GrSvi#TS1BXx1Qh<1Zp!>8rN9dS zFs(dmQ25|S&>Y)cHoCdpZA?X2X75m2-EcnSnEf?iJhoi7aR)HjKkfhqM~U+>;<68{ zG4oeg@`=m*Ix~NfnUC|)sYry3tMHqTyX9-AYc0Ayur1DsV|Gs*vtvx;1&%JR5B-OV zf|HyFhUxvdfAjF2ul+rL`^r;-E5|$Cf23ALAa2foW7T;$rSorU`I|F;TjQHW{@vts z=kL66&w-Cj`-=SA+%21(6Z7{rJ~Xj)hp(T%U{8}+T(m+GhD`cnNXik=0ZiEsj5P&} zhnU;xOOBq;V^D?+zhc-DGOksn=c!@z{b1%C(D7GZ9enY&x(=3ae0aIU6sJC6u7X}57gs@|v*#wDoP9Qb9~yGX!IAa1 zav_9N?5zIT2OPzQv<0UvaMf}GZrrM+tJ29h!srLEK~$q05}f7C-|TAR*M;!&Tv@n87aoOP2S+wB zJ;JC0=&{oD$-D{z)5!;Y=1OxA&BcD9k+ljZtCi+U;OV*2G<_^Q=+Au@r06h<=>MFH6#}7e+u36A$gSth3GCHs z@n|luFu|B>nnx*>3!0BnFSlSQVIBhg{yKD_f#5K=o`36#-GaG=wgIEWqyTuujtJ?cZbF5;Eq)ELzyT26iwWV=Zz9c5;D~dmG@) z9j>-LUn7;?$dp(bEd?fn<~6`rHk|wwMOCU~_k$>5#$&d}^Ei777#Gxw7GBGjg1wLa zBq(I>nEN4FI9!f?`+EM%eWGJN1un6k?}rMo1U}5P|3kFzE`c4iu6eV!f~H#w{qaVT z9J;m;WkSY(BE<)2Gf>q9eiB6&0i{M}Gu%aPHthg38k+XKLQB)`t+h98Kd5Xyu&fj& z91O1lG-du2GP|9KGi(?7@hLbqZoe*rYlVaSMmf;%$M1=L-2@!{I-ek9ECApG^fI8z zsLUue8>o=Nem#kGi6h??mqa}kb)mIiivMJ>gbC!FbEqOw1hT=tIjvKpdFz1L5P_kWxNg8AWf8q6>h=^r4s7ICSH7 zC<8ZM4MCH$jvxh*ix6CC>@_#DBx5uD4P0XF|izR2}u;gxl+6KVP30bY#{0XyC>mtb2K zyn0(&(Cx%FaEK?Ew1)zwI*)K$n2E{`ARb^vgWX$E*09R<0NxC@b*?J%0Ld%~jvP#6 z*tFoq>boLnge@apf6t$Y4+UOKE^Y}WSGwjJt+f5szvBB6w`<_!W5xOJHoo5SqLURR zUY(Wp{pzeCjSph?oM(Bji4{9m8aG2ztfon0<4tu|52B5?xf>vQ4<@n(xC+ z>-I`FO@f>DHSQ5^`XuV6$L(&i9MHB0QjT6L`_FmsQT+ashn^iMxB_az@Qe#29D1%$ z+uDKlLK1q`4vc=?hn@$?>_bh+5FJeQAh=;jRxjPC2)xMWF@Ww&;`~wdr}-!7Pf1xl zyHoCX_^J6#4xC}zVTIw#qx`G!<-~m1|5PX2XzyfePvAZFxAapv{&tK0wlAOk?OUj> z`ddGaJmUxm=#l5*%wJ_ugToh5g9XO>(o$nLAirPvRB*f2o+8@gU-rkdPz9UiULhuJ zXaErXvK+wHA^!-Hr9(a-lIW1=zQ@tRqBStMMW6go65xJxSF>~Sy;JY~Q~Euod%gfK zhfjVZ>h^EP+1)pC$otBDHk(^9d{8@ z*a`lI{k<&Xd!UMK<16}I>Gv22LPAdDJt%={KO`m6DC!0sEsYR-q@=f2@) zs26B5R5uqA3y=_t zlPV~K#mQr+(tPFcOc@p@rwVrU%AssGR%Bx&g^Q9cNE5ydXv*{onR0P*s<83khC!qS zb4Npbc&#e*2w?XBY4lV0Zb2IT&kJLW-XSt*^qT_9=p?WjeJ=>g=rvLF{wR7a&={lV z2nDes^HJ?!$}_0aG`e$<-RPRvVvT+UX&U_zplS3KklEc2eIvXNWXoj9eVXmO&v+S* zjT2CCp66=iTf9KQVSF2c^E<9bHbT3Q@h%9$O26?DP@IjBL-XhUohh8?-NVu;vmu}H zWhGLGFWD6UWf?ocBjU?%1lgEI@(AM+fSCBgia)1R>j0K2mDm7-D;6}xja10_$03JYuVV;OEWrA!xPVxxjet?(=!xhXrnfEgS9bxnersQQf zAY_a|9-kAeV9o}M=+$xsa~7rbSiy);5fe<_eo6HB*B4S>1krv8K=57{X(uW~pZF)U zLdF}K5;6{G)SvMQGyNHR1jxRNAY?oUKve=pBQU%!8XhN@T+51Lf(he7_}BkrJjq8n zR)1idiYJI4?o{;QdsMs#Uzm>6)kw*9#rF`m$61S<3^^93^h46_6nhBRxl=e~3kI49 z0S9FWc#dhY7P)d{%(T2)WN`3w2`pm=!E?>fz8`+zj)y=HzQlJYQS^~0+5nVVnmv3! z3-X}kT_}SPumDxM*PlyjFs1uSj2c)E_D?rr6XXcLxdS19nEL$TB4NV8@N|~eA>c!n zR$mOzrQ`!=m{vrP;n(oYop(ROt}`B1v2(u zvE#pB>Q(W_zxn-?4_d1Ij2rr3{A&v^{(%NkR*(ZLjCxobNZ*gfVKkVH6rw@KH2`H9 zSAs`GgEm1nen9dFV-!G4G~jw;yv!Rxpre6tJH)E>#(3C`PPOrBn1(01-dM(}|K)n) zL6oykh%KOrQNm5m*<KA!2<;BNCQI$$*5r#_F0+|GL~r6pIyaF zf5uDz%#Noqh40J&SS7w2K@%WEOo|gC*7OJw?T7>U4C@^#4j?`J6twZpDjgny557Xi z&mfe#9r%@Ak8HsqPwe+ot?)H^6;c5TO}Oh^aEJi%njA}u zn2uRPTXbD#?pU)M+A`jpbW;+ z#9@}Pv<55m*Z^|Zi?L(rDAF8DsVuF>(lE%>Ylxqq!dPRL3JQ&>B4#xv69gD@t~6#G zQp$q!qG(aRD#+)N$LS4k0Cdr`J z5P3rp4Ua$wJum=y-E+ByI2D-uehTB7zQ=&}of$Gdg}Qvgzhn`dQK>YzZUU;fmT0aY zO0I=bt}kJM;@SyR$eJU$-co7Ex)6<~xRQYiu2GU}kV=Cq80D(bTqjE|od(yVb4^)I zK!vOiwg~Ols5H2~8|9j*xt^C?2~nV{wUXtK$TZ3HP>Gy*89F=WhB`tr@f50HVtejNk`3Z^4t;+Xxs zpNF&SZ)Y(p`wDgeoaE+#jlRZxxfYOc09jHJT$WZWzZxq^6EZ?kUse+jqj7|qw2eY$ z$eMtOlm6t&iLPN=S>UUdmAE}eX;`Wwm)Nkw6J7IPWbWq_!e0xHxz-**tGWM(sRixO zNdY^@>~IIKJ1)-8iu|t*f7SohYhQ6n_aaER{a>ST*YdjyO4`9g2}ksP1qE-%G@ZX| z)e--(mM*8OmG{!g`VN2XDe~_d9(wIyQU2bW78aK!OfB16?0RCGCpc*Elwi@B&JyQa zr~$u=gGnJ=C*8xVR7DUVu7~o@=)5gBV=p|0UvJ>&UHtIS{Z(Xi@gnIP6N5#0PWPXv z7VUkhWs|cge}}8>0?cV{B;33=n$8KX1t<0JDSi_jmhAaH9O`QOJ?E#EB)YowX9X{e2r*c62`l3TF0wcz@OXD;!9QJFUOq{y=`e zn*Q(%e)F}P->;_paVGH{zQsiW;~$q3@ATCK3hqG{M}%$+R5Mdn)A%!xCP6XOc9X2F@eYu-;N?i5 zV5q3HF#y6H_^>;W>l1Z1wt}!4JDWhtWm1K<8_^;lR-@el>A!_J(4>v)5D@1gV&k5* zRH&^M#Bl4^T-JOAX>91IegWe+Xb3nF06rYAiQqwS=>Ta~urA{EZw?gfLltym@fU9u z1?M?E`u=lbx|;$5DHp&P^tzWZiVz^OK*V+3p&%A-dOs1zGuf-p24B$kMp|(0@xv2b zX)AXm6Zm$dO$!$77+mapZAx(Vj!gMYqHa%IpT2!c3@jAl?@sCBGKz1ctRE7WGd7qhI;{=(i%@ zwRS7gaKu%@JLy9J@OdR}HW!>MT$AxW_~0uzsY7795Jcr|XIyUn2TAJc@^aw3PVC2J z=6T6V55x_`C;ZAfB(E^SQA zF^3?$0UPiBas_(;V(&tHzP#%!;{VS5<>_>>_5L#U@#wMt;{I~ka~uQTMgRMk_m`W% z)cgHqWX8R}d@d>w!Cc>8UJg!7nJXwQ-u>mvnLi?ZJOQ--zj%N7iw%gmSbUxX*xUVO zbVA$xWhYDFYiG7z^lRsM_m?MwQ}x5x{c!qk&{*i==Ob9e^a%MCE<@C>S+CIsJ~?!5 z0472b-vO%o%jcjtwzp{jYI|FjF_)5~_m>Ar7XCR{y~~yNm#_LGtOyw$VCT~&?k}%J zim~|%P~z+p|IPc$uTP`A|N8yqEY!{o33tUfcmFOe!Tsf5;^iQcyjmu<#jS8d!1#bk z+;ahwaEm`gPb3mKnHwe9E%gvHo;eG7%Kqnn0m-THcu--8|9so>!G z))>I>7tdlxM@*e>%@t-G4A-)>-qHO4`TG8H4Jwod-FZfgL2V*~2DK2g8ngfeWzc0& zbYT>|2Bh`r=UY0!h4B!z-7oLREnB{&FoaOSI>_8#3OyOb%7Ko{1grDY+%_2RhbBZ()TTNWTxI5fGjM z;OL(N4~2{$DUl&#IWXZgj{fyLt_c{MSY(9p6o5JU=aVK!KTo`$l`M}E=rb;T;(g^a z@I}bD1bKW;aN<1@EMi>B6Yt@a*5kzchg2%&IIH_29{{)`B- zLdJeg!FS{u^=Ev=Ox%ALAp0Q#y#E5ADgmPf7^SvAR-919eMju6w_NXI6Y6T)$hCG^ z#z=~Id)4nPLs`FD&sX)R&lJWX*~8?u*RX1T_Q}i?qx!b*#0&!dT_5x|7?C5ejME9M zW8rYIITrRs(IZiG6VR9e&;;Yf02l%9VnknwD$P&(eu{;3?6=!Pe~%G83T4!YZbO;_ z;C?`J04x_W5jY%rjJ>C4d3r6sW{t`6F&G%o}0h5l6Xj?)FH7x~s%q zqtFpf{1T&P2=*|_|LY-Jyf~CmG-SMjA~bPhIcxPBJCRm_Hnu?zTo?qEI@coIUp zSE_o@4y&0V+To=1n0B~ZWUw6`5m<(}5@v0OYe5%%GGxq+qSa9}A80c4GksPoltpBH z5L%%fcEeX@JGd^jx5NFvi*1K*k)~_YSz1TdI>_v9fI;C@kYfT2hA;l!RKd83@O{KL zA&|R&60oDb;y}u^Xmsq#W(3mbA^ypa~F8;{P zS(qQ(R}}l-aL+0Bf9{^`_J1?S6U-g}@~o2K-_nlo^e18crf(B{=Ztwj$|4v;k`&)}kSl<@dMM#;#Dx~|KJ`ePHpw9z+9_aHxp9lIp(C2|Z z5A=DU&jWoP=<`6I2mb%{0RNq#;+Zoh<&XCIDwfray0#Smm{jczIaT#_IUqV_FR94% z)%r@SGpER$%u4*@Pqp>hBSkpfiMV8n2ORw2bB+Y@WP9uCYp<`uKapitVpX7~v8=kH z%FJ3^QM;_dSHC>R*W}~B_Ehc|pXug5`Bch(R0_&eRacac&qOJ2S#527xg)ay|MJuL z%waVkjl6JpOxlvFsIDj*NuJ?S(YR85%kghxK{EaUDrHXm+AFHbRgQn~OR*igL?T_6 zSWSa^74;3KX>>J|uBfbR5Z1+2VG#;>udgWcnTCuWUqj;X_^Ook;4iOQhFX-Ak;$}*Zhn8iW+al3EI(= z1HKyFk}~EswV(mt!i))Q2q&rFP9I_ib&r1CMLl{gK zLa{aXMIw8Fue(1I8H@Azy?~ua@BMis!asAi{uhzR4AMLB$JYVp{W20+2YBp(NaPsc z;0Gg-0eCZ%_Ye+E0X>gIBKH9fdkk^`=lwboIRd!uw~@%`B*+6S0z8J()K;gQ7uRKSe1b%$yTLd@^a2{Y24r1>HJoYN=Ap8sD zAs*%7Q`&iiZ$%neU(x|z z2lO0>LhT)U*ky$)5gr>E>w1PxZZ#NH)t9rKT-Og~50z z?{WMrMBh1+^$K|&Q{Ls2cf=~s-9Hm-UHB;b@VO$rl#h#oRpaAJSBB>aq6%MshJZ}2SW*H z?(Pq9TeiLY=T}}8|2`eRO+SV{t>=)7A$AF9N$B4~_Zg{aYZI?d3)NLr~?@ERQ+!3;G@LC-JaThjY;`uOol`YnY3b ze)0ROij|$t#52(IdX4FFgJ(shY0vDwul%~i5N1ef5JVem2(#7^u2o`k>SiQ zO-*Y_NI0CDT)c(-W+C!Ec_SJd$nOK~7-)Y1?pP@Je}`FQ3;R?%^8WPZe_x#Mg1nXZ zn)bgM=c)S;Oud;{d(aPj*mp*QHUzW@C+s_`>}@CdUM9-dq5Mc(UA_~z@Z}UScO+b@ zTCfP@JHfLKJQd`b#CTK8STK>XU?S%_%$2wF^T@fbuqIyHiShIf zw)CHE>GxE6fD|JS&48z<^f9kQ^qJba`(+xpE5(B&0`cI2lz5(0`8@aGqx134S=Nd) z&y?`!Gldxn!I6`AhN*n@dc`q_XUcf=84AyT_~52O%%h6oz*~}4F~{A?zCnr`$2sCT zT|vGV;-l%wDy{YI6uAxuUIU3o*Z*}y+5@6sJ+9JxAH^r4(poS4V|Bv*a+QzU9U`ss zWvT!9-|JO-*=ruh11vm~6ud&gN(HZ1@D>H{Qt%N4pH*<1f^R7Jk%C7R>^D%>^Bo1x zS8$wyG8gAr9THRD4tRUKa#xHxpD*O)o7C8H>ixIk#(A3j7Qq{R^KE7cMCbSR3OQbW zg+)YyIp#=tyCsSIcqy zQp(l!@~_gx<6qMYUF+GbV@9hxEK$9CwjriRkB{|?|AY(;Ail($FZw~ zoSyOII>CQxk9c*nqTA!!PlTMFacmM6UGio}-*YYJF^|yz>R!!l=+Gpi`fY zha^Gk^91y#oI2I8qZEIgl0V~INf7UPfwn3+^ORkBo%*1ncReWi3ng&8py&$~UDx{_ z=(>+7RslX0bi}Kfk}BTG0_D4E>bYyUq@N>!V*muklQRnROvPgbv_r_5DyiZfGtg*o~dvv_)N%&u4lCT1krD;oDv@dI`tX+ zuuNPi+=vpA@;{7$y8H&H$=KsZten|OK_w(C8&y*r!e*X&ev^a6#SHa|QEI`uhp zwv;2jj|2KL`4die>{b0l=Hh!2MPJvUDH49D?$ z=UWwj`*_KJh6Ij#75zy?*Xz1R6n)0!l2Lq{19XFu|BlW7yA<8?sN~n}dr-(xm07@P zIKZKvbxKdTQ9NTnr{6l2KH?iJn9H=euo+tx_`bPRD|Yj~6NW5?lK&Q1n8ZKUXMvlTFS8ir%H-fVOKR=t3Vge~IsDq0_5c zzO8*fQ}h`&ea^y0j`}y-{5DR}7ux#IbkO6)$r}@7yL6~}wX3gFax!iD-!9}}9deoE z6yLA`{gtA3a>Iy6d|M0nY%lfhQT#Jb6YLJ{hmS$0T|S#Xk1P5in_cInNPnh{kaFNI z@#KI`IpsFLc?7>=vVd!poOTswbR4J#o%5Q#-!>@udu{UXRrDfTJ3cGqTPjrO9>w2e zn=g(l`VyNT&Q6u~cBufP=Y=tX-c$Y+ioZ@Z{yxRx1zq$XrN8csz<@|>f zAigyQ`lzCJU~eqJ&La5m7{b~Rd!9SuOg?+;q@Zn z^_4DmRMk}ZK(Jqys;%Rtqy};5R#xw8@HOHuuar4(>!)HdE(cZAdwt8iW!1Gc6%7zp zUh7?4UAw5X+FOnbRSn+K#wJHu?XtS+3SUJzj_u-dcq>siFC>*N_g2*4Vw9t@zI0iI zx4d!LvgIgZNqMPd469Y&WwqsvxVpnTUW(lFURK-;OFPbsYs67- zRfBg?V^y`!v8cLSk-)7_u{EzwSBMK#-o=%*^*&hR<*$ZH^jF-ttM=YdU*)SPESx-P z`pnrg-P62zkW^i@2nXy9jiNs~q9w~}Yj8Esi%W!6HH&R!bZtFquko7dM$fz1P`oWg zmzb?&S}AWK;^tTlvSpVs^P+d=dSumzn{Z}r6-{-}v7~syrQQp2F31!*_Av z4tkjWfmXbWh#_2Kcg*5twKXsw&4gB`^G$g*w8ZX%rKNRsrD*c&qb{%(LK)G3`fWy3 zr;)RK_ zrcy?T(rR4Fjdh~v)?zgFY-RRB@$xR}N{k#4WzC+@P=U757H9{@9b%?=6;^TfW)Lp` z%`VfUQ<*iGbqYsHElpaqDqKiq_w5;;qyD#4rpxu{PjRAK^i5?CCB^VH#cJp1@c?zu zW6b%eyrI^+1U^NN*PJps%7`O^s6e9{N0Z9#9uYCNw?t72R?|?k1YMd?BEFTqC>}v@jYABtk_gVHz~syHI1f4*AU%b zL|g5AD}(TcE%fPTmpo-ly=6<5Vjilhjt;@GkYO7FX2*)U zL+cjJGVN2I?;%CZl@NXcA4TP#U_R`rUu+0ByLnF;QBSH4OF!$SM7-wbpA4`nb)w~l zmCxot^OikykKcIK@LJ*Xs{Mj(&}j901OrYkVmYlNs(0*#7mUuo;G(=-*~fGju~k!A zzqlc$4tBoSiMMAnvdMee5$#)1KZvd@hjhGV7wb{}>4u(b1l_&llGWj@;1m|O?}!-m z?@S3^wx!j6HH>5IKT((A=(O}LTNBHbpqXzPU`|+;s?8};ufR2lui zLqF;AgMst80zZlP(dG5`LK^Dt1If!r9skG#e)Dw`zs1w#^>;-Y<|t0Bzvk0$Df0Pk zpHA!Vk2KVu`A~oAsr6T49V>v)r@FlUPD#T~#i;Ao<+c5-s{91Sr@v>?P=DXZ@;uVD zmtPASeM0@jqu-w^sK4XX+4lDLBg0xgQ$Ypu$fAhPnI$7*)Hk#$J*S`*V*7Ewh2MzW2w_1d~{%39F^>}PS8{Fqt(>kORNjZ0(O3#*6?Mp@tuz@ufH36^4~1K7iD$dW{vv$vh}Kb zrgE(AU%FoH$G@Qr{io~K-=XcjjneRU0VFVg>+%|YgdA&mM~7l(9jd&3ucpguy>WtR9IR`utp_ZJ+;#rV;7>-d4!bEfX1p!%iXs~lzLvOKz8P2%ynwq4T*zw@+A N@t<^Zh diff --git a/jastrow_provider.irp.f b/jastrow_provider.irp.f index 9652d85..70a0acd 100644 --- a/jastrow_provider.irp.f +++ b/jastrow_provider.irp.f @@ -8,7 +8,7 @@ BEGIN_PROVIDER [ double precision, jastrow_full ] print *, factor_ee print *, factor_en print *, factor_een - + jastrow_full = dexp(factor_ee + factor_en + factor_een) END_PROVIDER diff --git a/nuclei.irp.f b/nuclei.irp.f index d11f19a..7a89720 100644 --- a/nuclei.irp.f +++ b/nuclei.irp.f @@ -62,11 +62,11 @@ BEGIN_PROVIDER [double precision, factor_en] double precision :: pow_ser, x factor_en = 0.0d0 - pow_ser = 0.0d0 do j = 1 , nnuc do i = 1, nelec - x = rescale_en(i, j) * rescale_en(i, j) + x = rescale_en(i, j) + pow_ser = 0.0d0 do p = 2, naord x = x * rescale_en(i, j) pow_ser = pow_ser + aord_vect(p + 1, typenuc_arr(j)) * x @@ -76,5 +76,4 @@ BEGIN_PROVIDER [double precision, factor_en] end do end do - factor_en = 0.5d0 * factor_en END_PROVIDER diff --git a/orders.irp.f b/orders.irp.f index df09f13..187f1cb 100644 --- a/orders.irp.f +++ b/orders.irp.f @@ -47,8 +47,8 @@ BEGIN_PROVIDER [integer, dim_cord_vect] END_PROVIDER -BEGIN_PROVIDER [double precision, aord_vect, (0:naord, typenuc)] -&BEGIN_PROVIDER [double precision, bord_vect, (0:nbord)] +BEGIN_PROVIDER [double precision, aord_vect, (naord + 1, typenuc)] +&BEGIN_PROVIDER [double precision, bord_vect, (nbord + 1)] &BEGIN_PROVIDER [double precision, cord_vect, (dim_cord_vect, typenuc)] implicit none BEGIN_DOC diff --git a/random.irp.f b/random.irp.f deleted file mode 100644 index 487982a..0000000 --- a/random.irp.f +++ /dev/null @@ -1,21 +0,0 @@ -BEGIN_PROVIDER [ integer, seed_size ] - implicit none - BEGIN_DOC - ! Size of the random seed - END_DOC - call random_seed(size=seed_size) - -END_PROVIDER - -BEGIN_PROVIDER [ integer, seed, (seed_size) ] - implicit none - BEGIN_DOC - ! Random seed - END_DOC - integer :: i - do i=1,seed_size - seed(i) = i - enddo - call random_seed(put=seed) - -END_PROVIDER