From 274b648e4c2c01c7888cda1b14a3f971f0a5e59e Mon Sep 17 00:00:00 2001 From: Panadestein Date: Mon, 7 Dec 2020 10:55:37 +0100 Subject: [PATCH] Added input files --- Makefile | 2 +- OUT_GF0 | 5 - OUT_REF | 5 - el_nuc_el.irp.f | 8 +- elec_coord.txt | 10 + electrons.irp.f | 18 +- geometry.txt | 2 + jast_coeffs.txt | 190 +++++ jastrow | Bin 0 -> 43912 bytes jastrow_provider.irp.f | 12 - nuclei.irp.f | 37 +- orders.irp.f | 76 +- orders_inp | 1810 ---------------------------------------- random.irp.f | 2 - rescale.irp.f | 2 +- 15 files changed, 261 insertions(+), 1918 deletions(-) delete mode 100644 OUT_GF0 delete mode 100644 OUT_REF create mode 100644 elec_coord.txt create mode 100644 geometry.txt create mode 100644 jast_coeffs.txt create mode 100755 jastrow delete mode 100644 orders_inp diff --git a/Makefile b/Makefile index 3988e26..24b8b55 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 crs +AR = ar RANLIB = ranlib SRC= diff --git a/OUT_GF0 b/OUT_GF0 deleted file mode 100644 index 8d63a23..0000000 --- a/OUT_GF0 +++ /dev/null @@ -1,5 +0,0 @@ - 0.17703281415941130 - 0.17332587478375222 - 0.51043161486617861 - The total Jastrow factor - 2.3650290467901534 diff --git a/OUT_REF b/OUT_REF deleted file mode 100644 index c932497..0000000 --- a/OUT_REF +++ /dev/null @@ -1,5 +0,0 @@ - 0.278125485678554 - 0.308815718675770 - 0.533040790952379 - The total Jastrow factor - 3.06479902202983 diff --git a/el_nuc_el.irp.f b/el_nuc_el.irp.f index 39e3c58..84cffd4 100644 --- a/el_nuc_el.irp.f +++ b/el_nuc_el.irp.f @@ -3,12 +3,11 @@ BEGIN_PROVIDER [double precision, factor_een] BEGIN_DOC ! Electron-electron nucleus contribution to Jastrow factor END_DOC - integer :: i, j, alpha, p, k, l, lmax, cidx + integer :: i, j, alpha, p, k, l, lmax double precision :: x, y, z, t, c_inv, u, a, b, a2, b2, c, t0 PROVIDE cord_vect factor_een = 0.0d0 - cidx = 0 do alpha = 1, nnuc do j = 1, nelec @@ -45,8 +44,9 @@ BEGIN_PROVIDER [double precision, factor_een] z = b endif do l = iand(p - k, 1), lmax, 2 -! cidx = 1 + l + (ncord + 1) * k + (ncord + 1) * (ncord + 1) * (p - 1) + & -! (ncord + 1) * (ncord + 1) * ncord * (alpha - 1) + ! This can be used in case of a flatten cord_vect + ! cidx = 1 + l + (ncord + 1) * k + (ncord + 1) * (ncord + 1) * (p - 1) + & + ! (ncord + 1) * (ncord + 1) * ncord * (alpha - 1) factor_een = factor_een + cord_vect(l, k, p, alpha) * (y + z) * t t = t * c_inv y = y * a2 diff --git a/elec_coord.txt b/elec_coord.txt new file mode 100644 index 0000000..0d33424 --- /dev/null +++ b/elec_coord.txt @@ -0,0 +1,10 @@ + 0.47107089936082125 0.36934601008972345 0.13813129163530469 + 0.11734479025488309 0.36584274529092520 0.92672173271006275 + 0.35754797210216693 0.24703841579308339 0.17042939824115788 + 0.31813480208419698 0.81855922566341688 0.93785308019180225 + 0.69675380000569165 0.48537277058577866 2.8246759780116260E-002 + 0.42646034508932040 5.1852406712776689E-002 0.93402091856761427 + 0.85418292218111624 0.51831660694819415 0.26738545106704503 + 0.34322893096420737 0.80132196716714288 0.33887619652274625 + 0.69754618540855695 0.39691148679374866 0.12716214781484358 + 0.57090264858892203 0.70538842967387250 0.31300897404191264 diff --git a/electrons.irp.f b/electrons.irp.f index 8915894..0a7bed7 100644 --- a/electrons.irp.f +++ b/electrons.irp.f @@ -3,7 +3,7 @@ BEGIN_PROVIDER [ integer, nelec ] BEGIN_DOC ! Number of electrons END_DOC - nelec = 100 + nelec = 10 END_PROVIDER @@ -12,14 +12,14 @@ BEGIN_PROVIDER [ double precision, elec_coord, (nelec, 3) ] BEGIN_DOC ! Electron coordinates END_DOC - integer :: i,j - PROVIDE seed - do j = 1 , 3 - do i = 1, nelec - call random_number(elec_coord(i, j)) - enddo - enddo - FREE seed + character(len=*), parameter :: FILE_NAME = "elec_coord.txt" + integer :: fu, rc + + open(action='read', file=FILE_NAME, iostat=rc, newunit=fu) + + read(fu, *) elec_coord + + close(fu) END_PROVIDER diff --git a/geometry.txt b/geometry.txt new file mode 100644 index 0000000..9c8041c --- /dev/null +++ b/geometry.txt @@ -0,0 +1,2 @@ +0.000000 0.000000 0.000000 +0.000000 0.000000 2.059801 diff --git a/jast_coeffs.txt b/jast_coeffs.txt new file mode 100644 index 0000000..593d63b --- /dev/null +++ b/jast_coeffs.txt @@ -0,0 +1,190 @@ + 4.7107089936082125E-004 + 1.1734479025488309E-004 + 3.5754797210216693E-004 + 3.1813480208419699E-004 + 6.9675380000569162E-004 + 4.7107089936082121E-008 + 1.1734479025488308E-008 + 3.5754797210216693E-008 + 3.1813480208419699E-008 + 6.9675380000569162E-008 + 4.7107089936082132E-006 + 1.1734479025488310E-006 + 3.5754797210216696E-006 + 3.1813480208419699E-006 + 6.9675380000569168E-006 + 4.2646034508932041E-006 + 8.5418292218111638E-006 + 3.4322893096420741E-006 + 6.9754618540855701E-006 + 5.7090264858892204E-006 + 3.6934601008972349E-006 + 3.6584274529092522E-006 + 2.4703841579308343E-006 + 8.1855922566341692E-006 + 4.8537277058577867E-006 + 5.1852406712776692E-007 + 5.1831660694819417E-006 + 8.0132196716714300E-006 + 3.9691148679374872E-006 + 7.0538842967387255E-006 + 1.3813129163530470E-006 + 9.2672173271006290E-006 + 1.7042939824115789E-006 + 9.3785308019180231E-006 + 2.8246759780116260E-007 + 9.3402091856761429E-006 + 2.6738545106704506E-006 + 3.3887619652274627E-006 + 1.2716214781484360E-006 + 3.1300897404191267E-006 + 2.5932160037276120E-006 + 8.9781735169153588E-006 + 4.4307945771984676E-006 + 3.0357326911046815E-006 + 4.7057217460270767E-006 + 9.1888685615387477E-006 + 1.3674140578661999E-006 + 2.0898016957168522E-006 + 8.3133097492563700E-006 + 8.1930838386368430E-006 + 5.0761544134895837E-006 + 9.6316690786857319E-006 + 7.3265410239880617E-006 + 8.0282246083863469E-006 + 1.7840082668643953E-006 + 8.7402085917097983E-007 + 5.0021478344850869E-006 + 6.9080963589020322E-006 + 3.6461563934505928E-006 + 2.8079367396118341E-006 + 6.4561942747179404E-006 + 2.0143276334222562E-006 + 7.2282413910517550E-006 + 4.2803612481228009E-006 + 9.8791262118408863E-006 + 1.0651141027935296E-008 + 2.1266859661958340E-006 + 4.7254708764414140E-006 + 6.7960239791661172E-006 + 2.5361974009558710E-006 + 8.3022917885344452E-006 + 3.9304102414302272E-007 + 4.6287756683245075E-006 + 4.8818348901487870E-006 + 5.4217534412207837E-006 + 9.2515623629987977E-006 + 5.6776909044188313E-008 + 8.5143283852193907E-006 + 8.4916543326313185E-006 + 2.6803664649533946E-006 + 3.8619825277037558E-007 + 4.1645240139124852E-007 + 4.7076790048846678E-006 + 3.7400607062400029E-006 + 6.5903674693304497E-007 + 7.0511416830091020E-006 + 7.1970472254075054E-006 + 8.3779145190230263E-006 + 4.4720525973186892E-006 + 5.3527028821450744E-006 + 3.2469611963900162E-006 + 2.7474685177806738E-006 + 2.8766396048376057E-006 + 3.0832590779442896E-007 + 1.1970794212078929E-006 + 9.8875410686779822E-006 + 3.8833096533945047E-006 + 6.5493531824651745E-006 + 4.8669577823889157E-006 + 3.0536887907279188E-006 + 1.5453212737038092E-006 + 4.2010386366569512E-006 + 3.4608479854625175E-006 + 5.5863755844856395E-006 + 3.1626318255858814E-006 + 8.9535689383160287E-006 + 2.7295097461499133E-006 + 6.5761403856116881E-006 + 9.9777710552222946E-006 + 5.7483205746228453E-006 + 7.2938932476699272E-006 + 1.8975608710106486E-006 + 9.0217837038189077E-006 + 9.2351027820514674E-006 + 2.6491154807706522E-006 + 1.5921008779667090E-006 + 8.6057966599747186E-006 + 3.7446440900687352E-006 + 3.4272640538904732E-006 + 2.0533004904915601E-006 + 9.2668744161565381E-006 + 5.7347639398840740E-006 + 2.1537537756204565E-006 + 4.9515420151455622E-006 + 2.5237493063925278E-006 + 8.7662988099822665E-006 + 6.5819974241574931E-007 + 9.2929324306096999E-006 + 3.4693265321361314E-006 + 4.5904275891192394E-007 + 4.3467474285524594E-006 + 5.6775174141423421E-007 + 4.2394412605675084E-006 + 6.9012994065703607E-008 + 6.2566316036101626E-006 + 3.7163201120581581E-007 + 5.4155288561946149E-006 + 1.8786819174573623E-006 + 8.4967241033663055E-006 + 2.6961279691572182E-006 + 1.7416078313976980E-006 + 2.1226880367527847E-006 + 5.6985152337693010E-006 + 4.5338244070194161E-007 + 8.1224288329971305E-006 + 1.2583441000010966E-007 + 4.0920562284569457E-006 + 7.6329426691450179E-006 + 1.7160963129884001E-006 + 6.6076648110477140E-006 + 1.0731735278990562E-006 + 6.1073303601923391E-006 + 1.8860699966918250E-006 + 7.1975926015998006E-007 + 4.3303800633873823E-006 + 2.7581618400452237E-006 + 6.0028583379676688E-006 + 7.7027018498727923E-006 + 3.6999785309263525E-006 + 7.4314767294185678E-006 + 1.4088729129277634E-006 + 6.8721711355692943E-006 + 5.6081301435941245E-006 + 1.8073198064823138E-006 + 3.2030269920557565E-006 + 4.6370278037429165E-006 + 3.2322542793190805E-006 + 8.0124486124572907E-007 + 4.1995613708513501E-006 + 3.4756871649772262E-006 + 8.0509458007326399E-006 + 6.5522790234765574E-006 + 2.5617201839219564E-007 + 9.6199052717720368E-006 + 5.1882924569606104E-006 + 4.4263082100511957E-006 + 7.7019916315718217E-006 + 7.9953275316002850E-006 + 2.4508656586831879E-007 + 1.9660391811454640E-006 + 1.1556395663871900E-007 + 7.2925513202755867E-007 + 8.6334898352408556E-006 + 9.1764148916426310E-006 + 1.9877147326789627E-006 + 6.6911730452935485E-006 + 4.8461832878537443E-006 + 9.5998598209927511E-007 + 6.5297643919457324E-006 + 3.3382674828396801E-006 diff --git a/jastrow b/jastrow new file mode 100755 index 0000000000000000000000000000000000000000..4279f94c5ec27e1837b73f5693e89a43fd756116 GIT binary patch literal 43912 zcmeHw3wTu3_3xP^M$ z;~-kK)OuU$tJEj8ts+znA`ajS@kM=s59~y&P((yUbAP}6ICJJ?g8l#R{qA?~{WyHd zT6^ua*IsMwwb$8upL4=>lZq!NxLk(55{N`b2I-lq_N0%MHP8^1G+{>CvV zCkdLYFOUGW@~TRRg~^JZ2*`G;RlO7TS9}_34GGz;g<%%3Lt zG}QH0fAeL%;ee48N!QTcj=ouajjUfCiiJaz5o46Tx*grYc579;wW^(lS;}7;YJXB3 zf1X$MPFRRKUK)~@*Yw@~a^fq?M6nyZT{&Ju-R^v}qn`iirQJS5_1BYMW-I+QRNWh= z)l@GSJ?gZY^4yy0y5_~Xi^q)49X)DzQ~mIf!tQ=hgfVsQ)C-KO2Q*#&{^qA=&%SYE z?rnjthZdbaZ9ny<4wOe8UgTHZ2|DX^f2{vcuE}6JJ`Q7f5Amm2Yp`cfqpII4@8~wS^`2myPoWUKCcHo%^v8( zd!X0%KtH($x~B)a;VYU|;w!IctXNpx4KUH!{=LATVLmE@|QOH zeLl$=pK)PjeWSmzw9dz0O_ddmzQ&5u8s8O-)qaSvYRl^D%1ZqZ*jQb+&_-|alfS$> z#a9kxV#)Ol6?L%`(TkJQQlzY^w6U}drghT>y_VI~H^s=W@cZf;{MGe!O>z3LlfYQo zxX{jJ^;p{2SbC*VvA6*agzw5O^OaRy=Bq5Nt}*D9n)))MRsx@|Y;mctvbwIcrur)S z8(E`F`Mb8Xy3VL9;e-#ONB$P?UbBl{6 zO!SQ$KFZM8NC>j^=ve9*}XKjv+&ijDA@@#aDk-o;*==S?wfsM|-_2sqE z^(~%t({1!Vni+7mjgEn~uK6~)+oEFpZ1gl6y}?FTOE^)x*haVK*A^T7CpJ0DY;@jR z^|iuAPuE1il{WebHu_QYs0WUE;HU?Vdf@*~4}9r9;@V0mO6S}r} zTc1zfC=~9>{|Mx+L7(8adq4p&maACTxvLA;pm$j2Ew^)a&gr6j z8_T@qb{2^8O)PT>(U~X8*R#x9Z)cV$U&S(Sxt(dE+{iMQ7M+GDU&=CXy`B4h0AbJq zmeW|?g))5?c@lJQD1YBn!+4gOM(zY5(ict{(Gj`oOakFi_W*en)U6WSWhr(fW09lxCi5#GRVe~LFaei8cZO78$EQU)DgFz*JV>ynL?!`%bgzy~>% zQJnEYd7TTvtMv2Cx?rYv>E~GpkYM^kv@Z7slLqC8KH9P~KpxT2zTMq+pP=k64h@<| zr>t*5nL!aQ%ueHa#2eUK94OCD_6GWAM=ruZtji!Zlz$%C1A8N{e~&Q%+#{i2W@9c| zM5BQmy};ge3?;Pt7^X!|QH`D;3`%{FY(OP{*DaAD0j|)fstF$ zXoUyEIp|pE05dJ6mBDdLg=3W&d#T(Itz0`~Y0Zvnl@={MEG;@wDrs3XR5zM#Yt#p> zrTN9G(I*mGjku2&ZaCT2=*H`1qo1lqj{)v?rUj<=$5rgO?(XM$j2wj+#L9?%D-n=2wfDsOxKj27#=}AaiQ5ekJO@+he zrvSVrwZH3LUJDM-MVJ^5!*7-Dy`Ok4^ak^HibBpC#X<8hq4meNFFDe^EMNyf{>wXLZi$Z-) zFLoW;?{TO3pGzjFc`mIaRQOy*v1?yZXx1jAZOm`zabjzi>)M>+0EFdi!q1C)dR49f z9iHh(-tV6CB7Qbtc(-}SKG*zDFl_X;ZeYdx_;~>pJNLWiZ3NtfpO3s_Uuk}DEgf6b zx*@69^=Z-APIvnymsRrb644r+nZFl9WKc;)Jol8_Wy>FAn+NPWFcSoKPIhJRM?-gDD3{c75y} zn$iocysjN^<@?D~Lc@Mu;`&EPsAdxeXhN5FXzqEGjbRHd38!ID|H@)$cJie5?brG^ z3OV~Hy3?|EO$qgTu{f}gl>QOqfBG9g{q}iVFiOD9B0o3BEyCqHIJ*JvC3!l?ZIl&g%U)y9$UD`%zA-%iFw0q*{cDQSlO0c0A60KPs{SOxFXqyy1EL3H59ei8@M}r zF{QO-w*Vj`5TybwKnnv8XZwLi`l88*Eix;jQ7FITN>DEYm5J_GUw3s;b(%x9kTk{3 z0stavvin7rf`f^Q9*8FOR3s`oNjh?4F%bZ4<+xy(ujm;;6Wt3zY zL_h}Av>rR+#Gv(duYJ%onD7AOE+r9>cCQ=>CM0a~4?1Py!L78dg z6ROo-jads7UZRwDuU+E}yzdQ-yV(2no8AroOzhkWp;3NTpnfeoPoUheS%`MY>>1qN zfayUIZ6iOUrUCOL7O7ktj1QW7uMlnWKoxEYWR4(1*gQc~!sgH@ngbNQApA1>k}r_E zT?*IT|@xt(7j`D4C4o5#lCO{nZz&k>2y{3ay?Y6gKu>b?>x zF9c4N+e_$}K(3Ev+^Iu@W=vw%csM%^qM*_uQGo!?^b=d(f)2P!fjQGa5C(?LaZ&XA zC^{9Wuq?#}g;!IjN5P^dVRl8kOWYfm9HlyYLTVXWiXPI2IiXHNr+TyX%a z1!ACZF1(F0;a)&XW&`RwU&1dKw0sc(4O)tvVKeA{Q9*<5C1yA1S`d^$mqpQ*D0&T0 zWzYnmF$PtNeqqqVR4J5l22+(b$m@3+^uU*~1|6deI>m0#LLn2nLhJ>SVI9}$$P8ak zgjZgJ5PG2#uyY4iAd$VmY;&9IQ@QsUpAMgf&GYbEn8>+$4@Cvdp`hW$Ka3zW-Q^v- z)qgRDH%+q+-hfgsltgfa&Al*2x#@YJ2p~?^*MKq?dRxD7H4j}!h7ld%q#KEZlYT84 zv|`V-0tXg-nZ0G^9cYJSvHT(qb!c3>AO|+?Igb7k;}bIqPe6(qg(rbVzTJmrfy@_yBMse5 z5H|k|phn|&QM3uD7~inj0#rz0ntESsWjLK5FR)GL=cxLz&6sZa~ZKQ9`E7K(9-Nbt$BU^4puf9cW5IPCjrngV+N8u)cgs6TAavL$w1Vd zBJ11)YE?4UyaOWLZLdOnGeeVsV+z76suix;lbCJlzk!Ur{;=?+u zC1{>0`gbpHgfpUUdj<_{Zp)tzE~ne_KvQnJZeXm-0-p<)C4fULa$XkMBXjHJFXZ3-&Lb9RV1*a}Q=CF}1)sh%#qmDqHgyh?#8->N_{%7Zy0p=rGop z^8|&)6cMuoDd4(I-;Nt56HZ?0_K4nDZ`k8nflI zSYrk&V?1oFjrpUHDO1a0$*}Gxq=oXwKz!r?gwO+hP}ez=spWWJl`FWODGa6zf(s|% zTxEV?(7Y0=PeK$|1jl~?6gS42V7?zl-~^^Am~t1Q9h_GM(>J4$%ooo9llfu|FnS|s zKF-<>|2%PyV(kYad$3UHoZSNXk)88r5jq?(kKuRhoVU_OH!PFY`E3R%BVvNrsa^8q%h z9sNs`f_V;r9ep`=^pAmqsrkq|`ftEx#cR;K6n)4so)1jM#|&VMa-J^{m-Sw&{!Cjv zaao^d)ep4lu@_KOc|Q-E36Ao&)v{N0zrVtOh2qzht1ys z6Eib^6|nY$ti98!U8ZYOZep#N4=q-0qpnRqj9nDZ+m~Qbg@t@Al*O9hm#7zOf{-Aa z%Sj$;UI)-^O>mQ}yPCjS6PQ)jnqbJM7^JZ2Mm?|Sn&3lr-EB?ayt#b|Rh&1s*Fod$ zM63xu1r~dQTfLhBZT>+RfrE$V$q_c6Mp-2GK;|>7!a7A$!sc3y26CTbWgzn&0Nm-y zSPGjB0HVv__|7U&$})!&a1Zk*>7!m`9BO_Ca<@Iqr|1{nb<38Y!*J2E9?^$YRUXMh z&7T2`*d991Ww~)f^C_qvEUQl<&<@_`98W^w3hS5e^o5eL0bae;;Fx{?5%5t%nrr@-bC_hVnf+D4R;A{OsgFi2JKelN~Xxt`tSUeVMxy=RcyI&ydcg17(x$o|9Khc3FP1Yi= zXl#f3_mb#t?+dMAdvpQ!-N$0E4nSS8fZIo(V1KPF!a($;YGLyrIK=RU&EC<?n7Tgm36gp$way435 zhcX{8Kpgi4hqV>#3+?P%!OMlVdm{lrOXgA`Q?6S&B*VH?(n9%XL40JKYBUS5vlE9h zYW<5fz#q6vTKt_nXNdd+=epPrR!HMJ8{tFlr|HkdXcjc@DHG$g{8E&nYh&k@;;#_k zOdq#{CLA_4>|7hye#j9pE5RYw#+#rHmXi;#9J4lFA}Fvk6Ym+uE+52`!htuiDBKavkneTRke6aYL^I@cl!TE?d&rQVmIS^k`uhEnz^;`Aq4lbDr~+1;36e#c19~d1Bwj!Ishsn z04~{zJjjqjZrJ>%Rgp^`q{i0-$Xx({hV=k6xbvNFv7EA24xCxA;(^0!IKZbd=M>w(71NgiT}RQdV) zA~df>XI5ygy~G)spS~9xnv>Crp*fGO^%n1ULZ(cWcSweHr;--RUj*@yzN*n+z|L}d zGmx7}U&JJ$FQF<#!(IO-X7kga)gZtzSLT9_^EZ@eFwlphXh#$c0gW*l_Z4aM*Dwe3 zaWw>4MsJ+!H2UA~#u|M-TG8lgw$?`D<|@_2ZD*BaSa$+xq5LX{kMviKMgex-BAy`S z&PKGz*f>UDG`63>G6xZ?Z-Ibg@PlCJ4+3m%AI4Jd1M^r?2O42>J&L^Gg4IBkH=b5h z#KxmQQMpl7Vs`vhQnn}xDl37;xak2(|?!MURrmgS9x@ibn(M=QB zTD$3wklDHOd+dkL0y#R}Rs>T{MYQ3JZ}|^k@#x*d)3jwKHveNV<*L9D9RI?6@zg)% z9>C&Y{$J>+(5Ow|GeQ@2;W@)`U<)ka!(Trh0|)c}jS-bEBH&SZG1&3Mxft8b(A*c& zLLL{kwf()Jq&F#`#C3q%+$?lfFOo8U_JoU{&oTMGanWF?&Z5- zhv%Y@w=1;R<^6gGUM1MDD^YC4?9NPo7wk|vvPMIuffmhMQG~{HNs>5>>TB%aoSI#y7}o?x~Qo21>vTWp>Wg6 z9nwvc;ijF9_Y&jrqw#p2Xr7k&qHtK+p!I-= z%iE=|?E%>dk}V*&iy{IN;X53-7!6Er5tHeX0QVzYt>ENyG&zss#i#UpOn5#DFGs$4 zGwSw#;f5CL_Q{8e$9~t`;SAbeJKe6^RSRbt; zHF`cM&r;q5G45H4yX_tL8|R>4z}X42hYNDTrVj))npHr9=7r?IesBT#Q)*bk6P7}h zu%myU@s9HnP{oe^4na15AbF^{7a-=mgxmRVW!)Zv=yIv;85k4!;SJmi!e%S#cyX~R z?k-|=Djt*scJl2DoJW8q87z1S>^m4=OKEst!`gI6p%0T8X1jvod5v*SUP?cbM7#LpRyDZMOc@;Z~d0Zmz z2iW(0gZqEy8nwUUc(2XoK)eHC)QNNE1K=F_R{%S&0IPML7g>Lp`@|x}yr)gYMdaiD zVC^qi3#BVB#s)<$^kZB(2yJby9EhURmCkMAd{6bW4-Evr9nY#(rx4^t(@B@E$fFMR7Y(5f2pNpcm0Hv1J5&`KLonhHA2g+dC zaTmI@o;ut=6U&a6M<_G5W7%PQgz^`Zx$Jlm(2|LJ7G%mLX@anE@49}Zh4N2@_{g=Y z(Vc*ueWcMJUBost`kz0KG5Q-(L8Fg_?g0Bpup9j%2+HWSQS{{~dKb_bqc6}#ms4## zH$s<|(Q9Tnjh^>ftkKV)Orzfdw2ZzGGCPk!Y(!oK*>6W+x1^K79s}mB;1WqEM<$&|P=+~-j7M1&Hh-llSf6V& zkozZA1~P90zyg0JOJQ?90K3Hdzrc`mPK}dv{@N|+tcZMvJw@335N>AQd$#C&(TX#N z&aE0>KLS;BM6_3D~_U+Gd+n|`!z2%@O_c{~N)2|5Ejs}NFME9dlB%(zu$7~ae z1%;Dw3NdRI!lNz_Fnb5X8aVc+D2Y+QbMPp7Vif%pVM1-J`Ewr%a{mKmFegVqkToX{ zPj$}8tuM#U$sxi7?0?u=?|<%xO#LkPPv|h#nCk?E##}?pZpgW0^^uqngK0i?iMm-R()79=x4b#naHX+ zs!;)8=TL6?CIgfA$Jl2>`6k98PWp-FkNaQ{p7wA(GXkoI&0Dmw0dpnLVETjPz!RTe z!ZFA&>k*d7Fsu_9rc98{_edUUz6=nPVYq_XA?vmf=nS*bEmtt(;1EPnDe8C~VFfc4 zEF!Aq3g$dY>$ZXssUl_t^Wuvl;*-eEBzlZYCRr$pkq%`3j8$Rt7)=SACu=m2*`Jkx zOjCf|7p4o=*8!+X(7YcQ9&!(klT0q{mQ0uzA`ktW`J@`{*!`XtR6aqz_*p7q_%4+% zA~XA;G#4c~u6V!W);MdCu}TxQtQ(D@(<#mrux_hx$d}W=3Wv1sl__8k%i^{?9V;Q* zS|ouo#eq?pz%sufr!|M*xi>%%z6_fiqUc*u^j@Gbv*r_+E>ghDPzEVrD!TNn9Y|^@ zq@z%r8v>_+cIQ)To^QS`$o`Y6yCqu&(@Vga@Tn!xBe z=+ZL!lOm_l+c(D=JyPjEg{`&G??I+qi%bL=`>)vfKYsjq@#nwwKG=QGQtf9(Kp)J1 zZ2{&#&|u1gAhTS1{QzaeHfW;zHkCPNxiiIOsAW)!W8qWM5$j5$Rc^DNZ~rBp$X zl^7nHjn~%HVnomME^3Kn2%jl52%3gR2u3D6Vxt1=l3WwN;hD^>YlG z;yR?c21u?XRPh4WHBqj4iW>Q5qtHv2!S(jJmaNr4c|~5ATx(SsT<1r*uGd_TNUo7l zu3)ssxti-n$yF28t3B$8dD;`Dl8f=c3q9`iT6%2;$}4h~1N#c*K913wG~yuMyyAocB_#(H!OTju;8%P+Vk`hVZNWk0 zP_W~TttyheApcg72h9SxBy?S@cMc9vllVdAoX@Gtbh6?XSSu&FbDk&1x2;REhq>6Qf0)%7M2Cg~|p2DT2~>IrQ3y!I_` z@|Qq6PhcCqN^~zT1y12x#LGA8KEsYS*o4e-6G%!`F)VM08HUKKyzMKm2x-hwiVU zVi!LwearKV@EA;ev2}wBAC4A|-QsR9#2w9pa*=C~XY4=R?frR&Yo&XI=e%dhA#k@p z%DWRP((oq>{Con%!EpI87+`6Ef86aCg9RKpXu1Sp1=r2@x|U?;d0lvu5d1l3A!OX` zJCG$OwH@>?fO><9)5OkwP#dEWUj#ys?MYMAzjVqw{I$LFLv*vk6FAiQMHG~2D8)CgtGHiK@8|BGAm7jR;z;Qq9{GOmG(5ALh^g7-zc85JA7hO7 zb1w;|d`oPLzdV?p14bUN2F+mdK-6X{|{~4rIkb>hIgw*B- zL3jjjqX+Y!CE;N6y&&9;?MpD_PSHp6-5^|s(FvwMUBIrI?*QRaB<-MCE_!TU9(m|H zu2;@S5sWa0$;NscfX~;aq6!IR8UQcO!bGs*!Ph=?J2TYK-MK1wVJc@l60lN(9}IplvUEp*Wj`zN^5 zmTXBT7`36;wRLJJX@ln}p=D^+YZ-Fh;B)Db90xSEe-8~Owh46+3y9yeG{iDLecpz! z-Rk@v`U6yiQGwXDzK8x$d=JeC_PRcpH1@^jI~i36Q`9HWlc2u%@)@T?gT}C+=1tH% z2_^YHX&SIz&;sv_0;A%5KJ)`Zg-(1PWjNv?RKb@u0Ql$`UowxM0~8-Hk7R{wKRmY& zE$|SoI5bOs-HeZ!$Fov>s=Ni&{CqfokC$gBNAmzaS;kk&+56-t%iCD!&u3@!)&Gv~ zm5cE0^2&?xp|U+*)cz)v8RSU*SHD;GqfwmimEY(MSv`HPjE!IH_sWogfc=rTuKnwK z{9YNm=-4-<)%VJ!XcF_*^%?|kWId={sC~SSOV#N2$~cg)z1__3m0tq~HyKBMyZL0u zvw1XjzbNK0LS!_0j(2Rd-L1Y?zN0{g-i^S-oZ+jj>U-tW&>Y*pEC99rE6ZF6KHK-o z*^(ulK{=M=j^qve|xNF3@^avDj z$P&5yDUOWUi$zR0^}X_!=mh#_f$Cx#oG+09=S!M`^Q9=70W{GNd@i6ki+xTCcol%1 zLbCmPtgY5~%Ha0x zWhjIl&nf8Eqd(`FGZs4@EV{t1KIaKPi5<@fN#R`hg)jr>KWwd!Le56L{$BZ76r@4r z@QKZ!J46Ky`W3PDEoh_+sslk8ba51Ih@!PXl|dtc#uzjTa)d#Lp(+e2LARDc%f>hj zI`~AaK^LMG4Xa~oZBV|DDYtO`^;-$~$&{b965Lf8a=HDUV&@KD(NwGtUN1HgpSIQ!S41s372 zvB^;LIRI<+SCb~c3*?FL3zFqY0(}OgPlSiSH(U>(j@J=RgvWt}qo58y$5L9i6XC0= zRLpVDQ-2eaVECn`7z&K6kaF9_^&Aprzs`Q9@vG$~4b4b+$ku!!Tp8U5W~#RM3F zF0Gf^79x~m-wprjFPH$Qp^ZAxy#TE^0czP=Pk<|gOgRDWLxn!hebHoJdm~w@(HKCi zi=Yyd<&AJ+bj1+aBc9)VgPcNxI;+LnyTFJf@@3NgSklXNUewnzi-ye)(1a#FTQ1Zv z-$7X=+SquihAl{{_pqJ{tbx$ZC8{5c!>z0kuaG_n=E_9Cn`W9EaM+W5?lTlMDDnZw zF$qTOPw@9mM#J}A{t3bS17`p`=r0PUoQXllK4f|@y$bnHenU48m|U?=2S&^A{S<3E z{G-Lzfajkl;)xB`twK=qlOW8-yPd(5&yZvBq-qTaH{u9?&mYCUHnRRp?k|p0;=1iG z`1iy-7Zr!b?->q9Tv!x%%X2|d;49B8PvGzc-cas)AkQow{5|(skpX{(^&7rV^bZ*E zu3G072lgXG3j_c3boJZmZQY1J=``ZY+2q{x2g!2C{7yr(1(TwSn z$L9I`6}1hgU0RAiCS8BU@ao2f;UF5bsw%Sl_5RYDtSPc4s}g_wxxP``Ux;Y!a?29O zpFg|?NDxo9uc5JiQ8oTVw$+JUfx70hnu=ur|yu!yTu2lb( z_}kc!j6Xmvxq9HMs3BK5{@}M_+jovccAaB44eC`iHd&_8)lj;kva(577gvP^XyjW| zQRa^tLek*)3YS^jafOvy)s~k?8m8=|D|AKPnROI4OhZ|+80}^?Ry38B)>K$_X%ThO z$Nb%cx~ioQf0H||zE*wRkJMMx`6{Hxk1Pj#On%+37EL-`Tp!1b8!f=3sEZCADKg%~ z$i-TWpo^=8XrzOrXAT+0u9Fkb;@@=ODtXz%U0s`juYRPfD;w*L9e`_4-tlNx*FM0t zk9BnwkiMp?%MUmkX>TRqA!PRmUSchl@kC)W;9|h{09OL;16&K3g!|HFzzo10fWrVIfCYdV z+i_3^nDcs9*Dk=>fCmWQz@8eb^_6dSb=?ja0bC83{B~E@A;28KK3I2*0n7qi_HI|# ze86P9sEFRVz-55?|LMR}d~vz1nqjyWr@2nP@G}>2bsZq()%0Rd36qsr0}nrqB*q+g z__C{UNov|H2_>l+%M&k7%?c$=OU=32^rYt9a7;<+nAT$xQ(e7M^NLb)CZuLfNX>vz z(A?7-;x@kU(#x;BD*nDbpFJldcJ)023v%%nW&439bf2D@wmjjNsTrZf>8V*aCwWtI zZZHc|^IDHdcri5x?ukkkQrUhrWFA2KNz@VDh7!PZGktkOk|#B%)ts8zactrsXh#8X zs#Xo`SbK=Sv!8_d^BU$w-1g_DwwMVU9nGhRG2t9tjX5@qeZz+dH^GS$dJDOnr+>$e zXyOs&-s+HB*qgz$1Z}rr&SG;dF50;Qw3k3T8}qP7Jrn=vv~x0LWZ`EU+he}J0bJ<1 z0J=`;olvF*oP8Vs&$4H_y3PRXp%3`-lmv&;9`-@KN5R%*9ntY6%?sM?4%%GMRyk-5 zpgjzlwwZ0)KwIsg@y|irXQMGrA6D`P^+5sSdkbS=c_KnFlr+JLQxD_x+|(6j;zNnO zQpX^65g+z&ox;8;?=a-Gywuf|!f{1R#+Nsl@(NQM5_ZKC=CQ?nF~*GVeEl5ac(IT9 zpv?vC012Xx^QiOWSe>mt)*Y!2%3A|@eO~G68u=gToM^g^AdtT2A51t5c|1ks#XkCh zHs;l?t`n_3j;Qm_BXl9+r5y5-w?UrhLB`Nx(E5RWJ;vJKJL>%l3V6kD9-utLuBjMGy0$#Y|l0Y(Gi2FGl+;?1hdzcOP=L z7yb7~mp8NhuCA`H>38Ir#QiRvrp35CjCup#S>LSmiyyBl{I<6dmwuRRQ=KBcK?tqBSHQvTN`u|_$5eu*T6r#(2bCjQTlwWp~KTzd9Qk*<&1CAG0#g*feP@l!= z|0nijRWF}57{*DetpC4Qo)_`b_4w~N>P4AnQMmNk#B_!5sESvPs^?n;y!!FX85hr> zc=6nbSGok&Z`Djb5CK&+Ggd15`YBEvkYI?FNYrQ)ZA07jVOZSgQLc$&p z1uNdB5apFh5MIO+Wvv&Eenh!k)#GiIn0h7t|2WY9ue`J4(c-Wd$50Ahs9>dnixj*; z!8;XvK*1*zd_lpt6#P`dg9;wgSN42@f~P7tT0yPX>J<`G?+*C7J#tr!x}Fc@<&)Id zcIpYX7Heq!=%aC~N+7!4dKB$=gdDG4oZl%i-Rd3BpS#uD8$*V!k9l}(Q!%Dq7<}%; ziFo}1|Hu8{tPy#`M-I=kamPm!4VST6PZ%2u|Hta*-p9y_OHVNh z9P!9@iAJiS^JA^K2}=pAm09 zCn8(!xW&#X@kHYUD_`1c6OEr5tJHjprDqtsZnaai9DR3))-mlS8UOBzQY6Ovj0=|O z_{mZX&?GK@W)E~e?{*n(qhLM|T>5^4_vhlu;rV(3?rN(X<2YH-w>ju@6kXqMQ4|-G zwhMP5?FV&xF|HQ$?&ER~=@yTSqUeZVaZOP4HHyynxOtWIAfJCeB3{2*j-ld5_sg%)Jv$uwYdu#fdSB^&qc`$z zJiqbXgLw3w=Ea%dcNu++PimE58;B>z3w|a~XCAv)$myO>E*JdAcgt7TE4nkk-67<3 z&tuyF;7?BDf#)>j@JgzkH>sm!wHT{uPa|t{2g~l z`fLe|&jr6K*uV+6&rwd6S|6)rtZ}NK+ZYu(1$64OOFRJN7|E*=`cqDW8rahmo$o!6 zzT#(+Af6Ee<+~4*vr^fm*Qxg^dX|!3pyWKG=<|Oi3A*18gnY$h178SwiZOkLB#LLb zK>3WD`lJn(bp0H+4+QEtXcJ?c26~pApm5qjIlD?EK|GrVIu&%n81lfv{8bvo8^yQ#uN!FO(4+?%{@%fUh z^?wTVczW_3eD*tA#lPm~dwOEr=Sev#yo_TJx8yH(#MLQ^zEAm0$LC}*Ph(URh54YT z#hJggptIjR<)3M)X^WzFsQD|N$D{67ik@|jWEAgo0R5w&tHw64sR#N2&DlR+1LaE!}2ioWXCQbK_ebb*qyZIC3qq3C?yfpRh&@xNTr8ytRqRM8`< zU+w2DTK)x6uz3Fg==+MkYNDizcUFKNP;%Bf`aKx}sgL2vhkTEL^cEEeI#142^c@aA zG=eV1@f6vwc<&8pEBOO+i;oHt?n5U!&r`Le(YWKtR~#$cM)(`aXQ_z)QS$1vn6N`fahp z&l45@Dtzt2>qH5R8H(Pa=z2YNiJ~{0FBx_GT(0E5=kWiHioQ(6p`On}==HnNNzS=Rq=`dW(7l%JbD0-HnYyGE!PQR@=L3SwKe+4>E%W;h3^@_g2 z(eHzrU+K9_>Ay+Q=R4x$1JI9GJa%vm)<5+BF4eD9hm9QVn(xRDS)fz@UFS(o@opN> zF^X<*Lx)Se8w)r|&;=>xSGnS!{u4>oepm!L`}I5gA5!#v4!a&#avIe5!c1{Jt@!VA z`0Z`Mub6D$b0w!`rX=e8^Bw5Ce>vkbAw}w$;TW$Wir(g+PXaxj|H~EsTE~4kr07)+ zKRhDjU|p&biO&BWJ;>j!_#0H`JC)5}DES?Z`0SIa;#v7w$L%jb7jb*46d>Le20B{P z^CdmSC1H`G7xb0oTnUV+N`98|kCtDn=*yHoIxiWe3#xPL>ireN%M70n1)slkp;28| z?FYg6F;jg5KN@Nhhhk-o{w9AjK7K1R@KsL5LVUba(dhHn`pRnR>nfTcti0a0u%>=N zX^pQOA3QbrN}CrOW%acUH5LAfavZD0zP0t`&G@K>-(V?n_wT0SR@fTi5Y+xH6~5{6Rl(NbBMl2ZJ&NWTHh!BF zd>?L)X5NX@x%VS(LB3$+w`ny2fuT zs3}(@aO+cM&8y25;sYn&!pizaKdkZbS5qb86<@^F_^xQI_E!`XoI82ij9D{0Q+;C~ zsit}X4#1n5MLZhOmSy#I_;k*PkNm3Z7CPGK-n#W(7gaZU-pqmG*GzPYHAez9Iw8?!CA8wVf|$Z} zPRA^)t*?Xm7$%H5oo~sjqa{usTvpo9P>La66jjRJ2yMgw>Zc75VGR<$?yIS9id8`S zkDng0Kj##-oTMu%n`>&Keb<)q+ee#EhUY_w7;I@%ZE1B~Y_;w&MxfAI+gedUP2}Xa z`qyr;jeR#B{e7er%)X*oCBA4nomDcC$#hoff|?4$S2PWNFR!liH8)k1BP~|YGkR== z9a^Do0mz3gQ734Nw2DG~M~sO;`cM4}-JrL=3B+pcD=+nz%7C%rEgBRS$$CB$jY@3| zl#h*ydVEz1({0n$iWzHSmsOQE;zL_Z@p!YlY;h?kc4-a1b&YkRh(_@lX!jvy$caZp zQCA}1MRK%;y{Q6Y<8*M`nPHh%VHf8d2Ju`i8XrAInboIoq}0-)MZ3axo{ZJ*SvBf^ zTW7jmw|I(^l%fv+yM?)!tCm>p93A{n2N7f4M#`J&eO2%&B3^UKq$U%E%v2nusK%nk zZg)%$VIqr-mS}S4tti@!Fd=qk#~nEpzVSttf{&u|k8m4<2<@_N zv!X%;BK@cD?D2+`U)b9_JJQ|X&gXDU^kemIzSd9aZ4;rw?^DZn$7JZ)9DeN0Ri@6A zzEP)-%o{uMj4}B#cJ$EOvopFWZCuzC(+7hp_U_qzcsP{stz&LPVnXY0Z-#ByO~)msKw zHzHft)~!!$o?;Zx&hOb)@mYO0_cM^?OCWT*ZMLzf=dX_Jxba#Z9u7)wZ z>B?Gv=>kB1ql8r&tFFUMyulbQ_NWW%nul}QU!GfCu1KBDQX9~! zahXqZhV+ojL{G-)&tdwczUf`(a|Q9)<>LeTN&1M0Q@ z`aKv8wf=aVA};IvLY1ci@@z`C*YC?%P*vdZgs@$=&rt0z1WdHG*YDM6sNbVg?A`pw z<88jPru*0L-)J~oHPrfRJ`FEJJ>P-T<&|_ME)5%GQT%Z>>Z$ctWaBCz^r>#I-{;Y= zLow?9b$e}pn`%En@#*(|G}P|_vi%A8akgI$8sm@t6_RpcZ`n@C#ahw-(IrU$KTK2E~Z_G!5*L?BY z;P>_R%h~^GRNC9?_X9NK-!idZ&i-V*G zTn8eq90{!7y1j;4{{9zDsx9}JKu7da-{iWM!xEEFS_WJ#=RV}iOevYLjXuWm$ zD3%hd$Jp`OPr6rf@J}GyFGtb%pPNeT#Q*>R literal 0 HcmV?d00001 diff --git a/jastrow_provider.irp.f b/jastrow_provider.irp.f index ab1e3e8..81a8635 100644 --- a/jastrow_provider.irp.f +++ b/jastrow_provider.irp.f @@ -3,18 +3,6 @@ BEGIN_PROVIDER [ double precision, jastrow_full ] BEGIN_DOC ! Complete jastrow factor END_DOC - ! integer :: i, j, k, l - - ! do l = 1, nnuc - ! do k = 1, ncord - ! do j = 0, ncord - ! do i = 0, ncord - ! write(*, *) cord_vect_0(i, j, k, l) - ! end do - ! end do - ! end do - ! end do - print *, factor_ee print *, factor_en diff --git a/nuclei.irp.f b/nuclei.irp.f index 0b431e8..090132a 100644 --- a/nuclei.irp.f +++ b/nuclei.irp.f @@ -3,7 +3,18 @@ BEGIN_PROVIDER [ integer, nnuc ] BEGIN_DOC ! Number of nuclei END_DOC - nnuc = 10 + nnuc = 2 +END_PROVIDER + + +BEGIN_PROVIDER [ integer, typenuc ] +&BEGIN_PROVIDER [integer, typenuc_arr, (nnuc)] + implicit none + BEGIN_DOC + ! Number of nuclei + END_DOC + typenuc = 1 + typenuc_arr = (/1, 1/) END_PROVIDER @@ -12,12 +23,14 @@ BEGIN_PROVIDER [ double precision, nuc_coord, (nnuc, 3) ] BEGIN_DOC ! Nuclei coordinates END_DOC - integer :: i, j - do j = 1 , 3 - do i = 1, nnuc - nuc_coord(i, j) = dble(i) - enddo - enddo + character(len=*), parameter :: FILE_NAME = "geometry.txt" + integer :: fu, rc + + open(action='read', file=FILE_NAME, iostat=rc, newunit=fu) + + read(fu, *) nuc_coord + + close(fu) END_PROVIDER @@ -43,21 +56,21 @@ BEGIN_PROVIDER [double precision, factor_en] BEGIN_DOC ! Electron-nuclei contribution to Jastrow factor END_DOC - integer :: i, j, p + integer :: i, j, p, q double precision :: pow_ser, x factor_en = 0.0d0 pow_ser = 0.0d0 do j = 1 , nnuc - do i = 1, nnuc + do i = 1, nelec x = rescale_en(i, j) do p = 2, naord x = x * rescale_en(i, j) - pow_ser = pow_ser + aord_vect(p) * x + pow_ser = pow_ser + aord_vect(p, typenuc_arr(j)) * x end do - factor_en = factor_en + aord_vect(1) * rescale_en(i, j) & - / (1 + aord_vect(2) * rescale_en(i, j)) + pow_ser + factor_en = factor_en + aord_vect(1, typenuc_arr(j)) * rescale_en(i, j) & + / (1 + aord_vect(2, typenuc_arr(j)) * rescale_en(i, j)) + pow_ser end do end do diff --git a/orders.irp.f b/orders.irp.f index b367df3..de52116 100644 --- a/orders.irp.f +++ b/orders.irp.f @@ -22,63 +22,25 @@ BEGIN_PROVIDER [integer, ncord] ncord = 5 END_PROVIDER -BEGIN_PROVIDER [double precision, aord_vect, (naord)] +BEGIN_PROVIDER [double precision, aord_vect, (naord, typenuc)] &BEGIN_PROVIDER [double precision, bord_vect, (nbord)] -&BEGIN_PROVIDER [double precision, cord_vect, (0:ncord , 0:ncord , ncord , nnuc)] -implicit none -BEGIN_DOC -! Read Jastow coefficients from file (NEEDS OPTIMIZATION!) -END_DOC -PROVIDE naord -PROVIDE nbord -PROVIDE ncord -character(len=*), parameter :: FILE_NAME = "orders_inp" -integer :: i, fu, rc, maxord -double precision, dimension((ncord + 1) * (ncord + 1) * ncord * nnuc) :: allord_vect - -maxord = (ncord + 1) * (ncord + 1) * ncord * nnuc + naord + nbord - -open(action='read', file=FILE_NAME, iostat=rc, newunit=fu) - -read(fu, *) aord_vect -read(fu, *) bord_vect -read(fu, *) cord_vect - -close(fu) +&BEGIN_PROVIDER [double precision, cord_vect, (0:ncord , 0:ncord , ncord , typenuc)] + implicit none + BEGIN_DOC + ! Read Jastow coefficients from file + END_DOC + PROVIDE naord + PROVIDE nbord + PROVIDE ncord + character(len=*), parameter :: FILE_NAME = "jast_coeffs.txt" + integer :: i, fu, rc + + open(action='read', file=FILE_NAME, iostat=rc, newunit=fu) + + read(fu, *) aord_vect + read(fu, *) bord_vect + read(fu, *) cord_vect + + close(fu) END_PROVIDER - -! BEGIN_PROVIDER [double precision, aord_vect, (naord)] -! implicit none -! BEGIN_DOC -! ! Vector of the `a' coefficients -! END_DOC -! integer :: i -! PROVIDE seed -! call random_number(aord_vect) -! aord_vect = aord_vect*.1d-2 -! FREE seed -! END_PROVIDER -! -! BEGIN_PROVIDER [double precision, bord_vect, (nbord)] -! implicit none -! BEGIN_DOC -! ! Vector of the `b' coefficients -! END_DOC -! integer :: i -! PROVIDE seed -! call random_number(bord_vect) -! bord_vect = bord_vect*.1d-6 -! FREE seed -! END_PROVIDER -! -! BEGIN_PROVIDER [double precision, cord_vect, (0:ncord,0:ncord,ncord,nnuc)] -! implicit none -! BEGIN_DOC -! ! Vector of the `c' coefficients -! END_DOC -! PROVIDE seed -! call random_number(cord_vect) -! cord_vect_0 = cord_vect_0 * .1d-4 -! FREE seed -! END_PROVIDER diff --git a/orders_inp b/orders_inp deleted file mode 100644 index 4437df9..0000000 --- a/orders_inp +++ /dev/null @@ -1,1810 +0,0 @@ - 4.7107089936082125E-004 - 1.1734479025488309E-004 - 3.5754797210216693E-004 - 3.1813480208419699E-004 - 6.9675380000569162E-004 - 4.7107089936082121E-008 - 1.1734479025488308E-008 - 3.5754797210216693E-008 - 3.1813480208419699E-008 - 6.9675380000569162E-008 - 4.7107089936082132E-006 - 1.1734479025488310E-006 - 3.5754797210216696E-006 - 3.1813480208419699E-006 - 6.9675380000569168E-006 - 4.2646034508932041E-006 - 8.5418292218111638E-006 - 3.4322893096420741E-006 - 6.9754618540855701E-006 - 5.7090264858892204E-006 - 3.6934601008972349E-006 - 3.6584274529092522E-006 - 2.4703841579308343E-006 - 8.1855922566341692E-006 - 4.8537277058577867E-006 - 5.1852406712776692E-007 - 5.1831660694819417E-006 - 8.0132196716714300E-006 - 3.9691148679374872E-006 - 7.0538842967387255E-006 - 1.3813129163530470E-006 - 9.2672173271006290E-006 - 1.7042939824115789E-006 - 9.3785308019180231E-006 - 2.8246759780116260E-007 - 9.3402091856761429E-006 - 2.6738545106704506E-006 - 3.3887619652274627E-006 - 1.2716214781484360E-006 - 3.1300897404191267E-006 - 2.5932160037276120E-006 - 8.9781735169153588E-006 - 4.4307945771984676E-006 - 3.0357326911046815E-006 - 4.7057217460270767E-006 - 9.1888685615387477E-006 - 1.3674140578661999E-006 - 2.0898016957168522E-006 - 8.3133097492563700E-006 - 8.1930838386368430E-006 - 5.0761544134895837E-006 - 9.6316690786857319E-006 - 7.3265410239880617E-006 - 8.0282246083863469E-006 - 1.7840082668643953E-006 - 8.7402085917097983E-007 - 5.0021478344850869E-006 - 6.9080963589020322E-006 - 3.6461563934505928E-006 - 2.8079367396118341E-006 - 6.4561942747179404E-006 - 2.0143276334222562E-006 - 7.2282413910517550E-006 - 4.2803612481228009E-006 - 9.8791262118408863E-006 - 1.0651141027935296E-008 - 2.1266859661958340E-006 - 4.7254708764414140E-006 - 6.7960239791661172E-006 - 2.5361974009558710E-006 - 8.3022917885344452E-006 - 3.9304102414302272E-007 - 4.6287756683245075E-006 - 4.8818348901487870E-006 - 5.4217534412207837E-006 - 9.2515623629987977E-006 - 5.6776909044188313E-008 - 8.5143283852193907E-006 - 8.4916543326313185E-006 - 2.6803664649533946E-006 - 3.8619825277037558E-007 - 4.1645240139124852E-007 - 4.7076790048846678E-006 - 3.7400607062400029E-006 - 6.5903674693304497E-007 - 7.0511416830091020E-006 - 7.1970472254075054E-006 - 8.3779145190230263E-006 - 4.4720525973186892E-006 - 5.3527028821450744E-006 - 3.2469611963900162E-006 - 2.7474685177806738E-006 - 2.8766396048376057E-006 - 3.0832590779442896E-007 - 1.1970794212078929E-006 - 9.8875410686779822E-006 - 3.8833096533945047E-006 - 6.5493531824651745E-006 - 4.8669577823889157E-006 - 3.0536887907279188E-006 - 1.5453212737038092E-006 - 4.2010386366569512E-006 - 3.4608479854625175E-006 - 5.5863755844856395E-006 - 3.1626318255858814E-006 - 8.9535689383160287E-006 - 2.7295097461499133E-006 - 6.5761403856116881E-006 - 9.9777710552222946E-006 - 5.7483205746228453E-006 - 7.2938932476699272E-006 - 1.8975608710106486E-006 - 9.0217837038189077E-006 - 9.2351027820514674E-006 - 2.6491154807706522E-006 - 1.5921008779667090E-006 - 8.6057966599747186E-006 - 3.7446440900687352E-006 - 3.4272640538904732E-006 - 2.0533004904915601E-006 - 9.2668744161565381E-006 - 5.7347639398840740E-006 - 2.1537537756204565E-006 - 4.9515420151455622E-006 - 2.5237493063925278E-006 - 8.7662988099822665E-006 - 6.5819974241574931E-007 - 9.2929324306096999E-006 - 3.4693265321361314E-006 - 4.5904275891192394E-007 - 4.3467474285524594E-006 - 5.6775174141423421E-007 - 4.2394412605675084E-006 - 6.9012994065703607E-008 - 6.2566316036101626E-006 - 3.7163201120581581E-007 - 5.4155288561946149E-006 - 1.8786819174573623E-006 - 8.4967241033663055E-006 - 2.6961279691572182E-006 - 1.7416078313976980E-006 - 2.1226880367527847E-006 - 5.6985152337693010E-006 - 4.5338244070194161E-007 - 8.1224288329971305E-006 - 1.2583441000010966E-007 - 4.0920562284569457E-006 - 7.6329426691450179E-006 - 1.7160963129884001E-006 - 6.6076648110477140E-006 - 1.0731735278990562E-006 - 6.1073303601923391E-006 - 1.8860699966918250E-006 - 7.1975926015998006E-007 - 4.3303800633873823E-006 - 2.7581618400452237E-006 - 6.0028583379676688E-006 - 7.7027018498727923E-006 - 3.6999785309263525E-006 - 7.4314767294185678E-006 - 1.4088729129277634E-006 - 6.8721711355692943E-006 - 5.6081301435941245E-006 - 1.8073198064823138E-006 - 3.2030269920557565E-006 - 4.6370278037429165E-006 - 3.2322542793190805E-006 - 8.0124486124572907E-007 - 4.1995613708513501E-006 - 3.4756871649772262E-006 - 8.0509458007326399E-006 - 6.5522790234765574E-006 - 2.5617201839219564E-007 - 9.6199052717720368E-006 - 5.1882924569606104E-006 - 4.4263082100511957E-006 - 7.7019916315718217E-006 - 7.9953275316002850E-006 - 2.4508656586831879E-007 - 1.9660391811454640E-006 - 1.1556395663871900E-007 - 7.2925513202755867E-007 - 8.6334898352408556E-006 - 9.1764148916426310E-006 - 1.9877147326789627E-006 - 6.6911730452935485E-006 - 4.8461832878537443E-006 - 9.5998598209927511E-007 - 6.5297643919457324E-006 - 3.3382674828396801E-006 - 7.8892927193568410E-006 - 4.0545438529473023E-006 - 4.2315175175221882E-006 - 8.7779323931576882E-006 - 5.6698502492352151E-006 - 7.4488720279602900E-006 - 4.9746787237257920E-006 - 2.1470924036132623E-006 - 8.8211617238327347E-006 - 5.2716457263083544E-006 - 6.8881507160486730E-006 - 3.0561586993413626E-006 - 1.5321103895247068E-006 - 9.9056678698256946E-006 - 7.9463311803198588E-006 - 5.1831686319353133E-006 - 8.2256337617059838E-006 - 1.9240358430525620E-006 - 4.3989125583153533E-006 - 8.4753375583306775E-006 - 7.1936192205316367E-006 - 3.3918338802167370E-006 - 1.2905767459500207E-006 - 7.5281621970489112E-006 - 3.7785650698919820E-006 - 2.2122923833237741E-006 - 7.7614354308813529E-006 - 4.1614475853883898E-006 - 8.5036600804519213E-006 - 2.4987303115094565E-007 - 6.8330576000932780E-006 - 1.2311783444282410E-006 - 4.3933911060159918E-006 - 9.7024813994546786E-006 - 4.4590121260581168E-006 - 1.6335850567541500E-006 - 4.4622481978343948E-006 - 3.8405737088863379E-006 - 5.7049415063669120E-007 - 1.1042471688324885E-006 - 3.5515926390478426E-006 - 2.2348927953418466E-006 - 1.8358616934339679E-006 - 5.9649244369617832E-007 - 6.1201026056940337E-007 - 7.4853729998439442E-006 - 5.0600912174293041E-006 - 2.7255819526353399E-006 - 1.0487437382193289E-006 - 1.4215965228411642E-006 - 9.7884335122164506E-006 - 2.6820607301194102E-006 - 9.3440711101772464E-006 - 7.4365300915620062E-006 - 5.8988807480775675E-006 - 9.2888697326872487E-006 - 4.7916106213613335E-006 - 2.6373735460162985E-006 - 3.6164373569987609E-006 - 1.7238452556617725E-006 - 1.6864056345231927E-006 - 2.1336502239386412E-006 - 3.8421864748963126E-006 - 5.6779742564259207E-006 - 1.3007671573054048E-006 - 2.4877906781527903E-006 - 9.4647435418263192E-006 - 6.5466261265971966E-006 - 9.4705960839120536E-006 - 7.4943583639177683E-006 - 4.6729538330085988E-006 - 6.8424762992843149E-006 - 2.4405158477784641E-006 - 6.4346097450801584E-006 - 9.1974391422440142E-006 - 2.8942506600420397E-006 - 2.6212237639989923E-006 - 7.2356378940722717E-006 - 4.8474759421637001E-006 - 7.4066347457883273E-006 - 9.5711273457649725E-006 - 9.8751337648048810E-006 - 5.4034072577735666E-006 - 7.8488630011203323E-006 - 9.6970968664011973E-006 - 8.8668458846582551E-006 - 1.6753682155565765E-006 - 2.3492418917994343E-006 - 1.6072421320941479E-006 - 9.9730332532416325E-006 - 3.4302013913563894E-006 - 9.4293990164209319E-006 - 8.5786998873232962E-006 - 8.0167629631948003E-006 - 4.1528134932671270E-008 - 5.8536686255071749E-006 - 6.5280026244230991E-006 - 9.7940668217482297E-006 - 3.4568838540586625E-006 - 8.8830849834772234E-006 - 4.0617230137762616E-006 - 8.9092811238552202E-006 - 9.7211480331855960E-006 - 6.1260399811390646E-006 - 6.4364166314542788E-006 - 8.3577791604662015E-006 - 9.5735890038188523E-006 - 6.9549554817249781E-008 - 2.0273457901530702E-006 - 3.9316803918448329E-006 - 3.4719181209234741E-006 - 3.0598491070621661E-006 - 2.8997338259331874E-006 - 9.5609434851333089E-006 - 4.1902842408987451E-006 - 2.5876591577206554E-006 - 2.1056421237099854E-006 - 6.3031497968762332E-006 - 7.2530039318275781E-006 - 3.0824909535908942E-006 - 3.8174191834999297E-006 - 5.5798118841099312E-006 - 8.5685502329479496E-006 - 4.3548244933106222E-006 - 2.9476121287621352E-006 - 8.0064796094441557E-006 - 7.8513946164704174E-006 - 8.0179169392501073E-006 - 7.4121643549080100E-006 - 6.9677723073585737E-006 - 6.4730088550863103E-006 - 1.4643807550873478E-006 - 5.2555310942560499E-007 - 6.9875455169783247E-006 - 2.9220313292350654E-006 - 8.1920611535312833E-007 - 5.1092435844601988E-006 - 1.4117169620739435E-006 - 7.0062454798629108E-006 - 2.6114634387061346E-006 - 6.2502145309488072E-006 - 8.2414956266183235E-007 - 1.6989634112870924E-006 - 6.9163231634036845E-006 - 7.3625087778373973E-006 - 7.1854188114099136E-006 - 2.1876749780264074E-006 - 5.3164852566265220E-006 - 4.2869517755284544E-006 - 4.6706183902426171E-006 - 8.8179293356220365E-006 - 5.9855486604178359E-006 - 9.6213647458493329E-006 - 7.6901164105689705E-007 - 6.9828784852540924E-006 - 8.4007551004112864E-006 - 7.9193852119960885E-006 - 7.1871976549804277E-006 - 3.1443314579486983E-006 - 2.7049139924710863E-006 - 2.4551632352758481E-006 - 4.4319108212989790E-006 - 8.2150584255019515E-006 - 7.3221458354131411E-007 - 9.7581691743915488E-006 - 8.3088933725356665E-006 - 4.1243971205565112E-006 - 4.8727533373002145E-006 - 4.0433561671611718E-006 - 6.4542388633113676E-006 - 1.4131968384385075E-006 - 4.5597263389556965E-006 - 3.2495487461874640E-006 - 6.3579329131857670E-006 - 9.0657649454017264E-006 - 1.5902503444944728E-006 - 1.1665741417061117E-007 - 1.3034697366008397E-006 - 1.0111475074454601E-006 - 4.4430436524002764E-006 - 8.0518792720348742E-006 - 1.7205121282508120E-006 - 9.8271985692119369E-006 - 8.6041883859727464E-006 - 9.8368140288376928E-007 - 8.0618628191384117E-006 - 1.4519390717557713E-006 - 6.5193698821044862E-006 - 6.1482313738382069E-006 - 1.1018120846688985E-006 - 8.1392522621615715E-006 - 2.1141354802198844E-006 - 2.5321742921602476E-006 - 1.8185441265292093E-006 - 7.0286836935569023E-006 - 6.4334736478863998E-006 - 2.8441079813898254E-006 - 5.4476384023351865E-006 - 4.4339428751014498E-006 - 3.1054834325132120E-007 - 3.8210964878851242E-006 - 5.9037515812345194E-006 - 2.3174686286700132E-006 - 1.8614261606145445E-006 - 9.3108372673720453E-006 - 9.3167693370933180E-006 - 2.2029222128113159E-006 - 4.2923268075947221E-006 - 6.5901041091479434E-007 - 4.9343920043145524E-006 - 7.7509955783571225E-006 - 5.8188027633112042E-006 - 3.4954640642238212E-006 - 2.0327320845174714E-006 - 3.7563342659519164E-006 - 4.2293741814419894E-006 - 9.7682955605932446E-007 - 5.3653287409997299E-006 - 4.1191296123816615E-006 - 8.8650234749779739E-007 - 6.3855388662400045E-006 - 2.2404236069495443E-006 - 4.1567242413025124E-007 - 7.1912007354129528E-006 - 4.6439723519071521E-006 - 4.4162525823315682E-006 - 1.0015315334899677E-006 - 7.6846587049407034E-006 - 7.2848589540208949E-006 - 2.0740471362235805E-006 - 7.9341917190575158E-007 - 1.4981874761044423E-006 - 1.4086030732149160E-006 - 6.3295593385373713E-006 - 7.1034220767948757E-006 - 2.2280528393357215E-006 - 5.4441045531281290E-006 - 9.1285786521128165E-006 - 7.6508233893637148E-006 - 1.3720058190188535E-006 - 7.7728102853102702E-006 - 7.5526858374867220E-006 - 6.8227749057264477E-007 - 4.9003762735521456E-006 - 6.2217941385036517E-006 - 9.2119971979416092E-006 - 8.2007164574454961E-006 - 2.3058997093825641E-007 - 6.1093717461052296E-006 - 7.9220619093335316E-006 - 7.9302861734307610E-006 - 7.9404693130767433E-006 - 4.5041748361358435E-006 - 1.4862400280293355E-006 - 9.0123873070502122E-006 - 4.0160707570629397E-006 - 1.1304766436378045E-006 - 3.8571412542580865E-006 - 7.3922177194562794E-006 - 7.3415807350903177E-006 - 7.1791956480965729E-006 - 5.8612335174691540E-006 - 7.4495611328806270E-006 - 8.8913751041091356E-006 - 8.3627778369869377E-006 - 3.0207312579973046E-006 - 4.5758541313740212E-006 - 4.0280255763489096E-006 - 3.9479653013089425E-007 - 2.9424347284904064E-006 - 1.0138898055719559E-006 - 5.7030611614821766E-006 - 7.2407458450613533E-006 - 6.6811109421482286E-006 - 2.0963536820871545E-006 - 2.1882603457290708E-006 - 1.5356667940231651E-006 - 9.1853464815623413E-006 - 1.4646224738664472E-006 - 3.5659476383331726E-006 - 9.9297289599346218E-006 - 7.5050908858094537E-006 - 1.7036697674509629E-006 - 5.4828749299534386E-006 - 4.0954357667561397E-006 - 6.2502304506433063E-006 - 4.0561906404871678E-007 - 1.7746772528998924E-007 - 3.9532755375658396E-006 - 2.7708729285534798E-006 - 3.8674007788210922E-006 - 4.4515561343786654E-006 - 2.7832168731980837E-006 - 4.8213597427338634E-006 - 1.1231865627351179E-006 - 8.4647312247123575E-006 - 6.9303175309768698E-006 - 2.3198166070162732E-006 - 3.9278379171194834E-006 - 8.4747625176460778E-006 - 5.5026908888332863E-006 - 2.5233481940069526E-006 - 2.4629510762943308E-006 - 8.6287648595744044E-006 - 6.5579340270033100E-007 - 6.9281244186382047E-006 - 9.8786708264740045E-006 - 8.1076948739910911E-006 - 3.1279349601389121E-006 - 6.4749258884428165E-006 - 1.4206970407339416E-006 - 6.7516151386178279E-006 - 5.5188452946307513E-006 - 5.2639487108979502E-006 - 3.6705951748748136E-006 - 4.2892754262151677E-006 - 4.1934860793496178E-007 - 8.2605024344381491E-006 - 9.6333038927667375E-006 - 7.0157369943310041E-006 - 4.3857381631390606E-006 - 1.7500434418612511E-006 - 9.1116052713451256E-006 - 5.5184817650060753E-006 - 9.7849312617351406E-006 - 4.1147190802448335E-006 - 6.4837034405345659E-006 - 3.8199481156537847E-006 - 1.6549433732671505E-006 - 1.7244421507131926E-006 - 1.9677276756987339E-006 - 1.3896070611666534E-006 - 1.7092126136549492E-006 - 8.2569203929582356E-006 - 1.7972223946802846E-006 - 5.9810911256113826E-006 - 1.8535754671079109E-006 - 8.1308111061466828E-006 - 4.6016702067046518E-006 - 9.2083252140628064E-006 - 1.3688267976915836E-006 - 6.6641188188803389E-006 - 2.2073156156343934E-006 - 9.7751378494659615E-006 - 4.7797661330075549E-006 - 1.9949437430821551E-006 - 6.1204453752567364E-006 - 6.8018672566254541E-006 - 4.3702897074440973E-007 - 8.4314508495841939E-006 - 3.2469779067809158E-006 - 1.3800343794834258E-006 - 2.3376741489602663E-006 - 3.8565933932378338E-006 - 8.7870400304783869E-006 - 4.7825482879772912E-006 - 7.2647923189680271E-006 - 1.4793931150084018E-006 - 2.8849568217425250E-006 - 4.0717399346220929E-006 - 3.9214194200621866E-006 - 6.9825407754900465E-006 - 9.7749432034137996E-006 - 6.2651537853722025E-006 - 2.8742739745690028E-006 - 8.4319478937106408E-006 - 5.7139521710019677E-006 - 3.7092505367835207E-006 - 2.1921470869513549E-006 - 5.6702780395655133E-006 - 5.5662996267948240E-006 - 2.1693964822983306E-006 - 6.9707927589932982E-006 - 9.1017685202695225E-006 - 5.0149124958578777E-006 - 2.9224793513132397E-006 - 2.0977285577082974E-006 - 9.0928016027554417E-006 - 4.6326465340089556E-006 - 7.2716075051055678E-006 - 1.2626623010888029E-006 - 6.6449383590433073E-007 - 8.3285210433376944E-006 - 5.1110338579068963E-006 - 1.3879023603262965E-006 - 1.3696070852155885E-006 - 3.0548085163653083E-006 - 1.5279291063412227E-006 - 8.8626809636765222E-006 - 5.1025781464889339E-006 - 4.7696282493556347E-006 - 3.0814972922931708E-006 - 2.8221053369720731E-006 - 6.3526787932251454E-006 - 1.4412221988209563E-006 - 7.8291796684353121E-006 - 3.1978327471581262E-006 - 1.4465854832854486E-006 - 7.8826815612032306E-006 - 6.4105181152943131E-006 - 8.7894185731033918E-006 - 7.5914500315749954E-006 - 5.8576114272097047E-006 - 3.5201188527664532E-006 - 6.9136632893321433E-006 - 1.7051588150521658E-006 - 8.9218048226862820E-006 - 5.9199714907871322E-006 - 1.5949535778504344E-006 - 6.4533361772453261E-006 - 7.4654228341951049E-006 - 6.3119143526271627E-006 - 9.5410567978127187E-006 - 2.5761378938737979E-007 - 7.9319386741114218E-006 - 9.2451871848932844E-006 - 7.3432798150616621E-006 - 2.2719344293257516E-006 - 8.7441905652593814E-006 - 3.8269184080218420E-006 - 2.6126610019214671E-006 - 2.3041828077598795E-006 - 9.5226503335866445E-006 - 3.9835786659719596E-006 - 8.1948240279883798E-006 - 1.6190771372773894E-006 - 6.0115501868825309E-006 - 6.5098609795497889E-006 - 7.1894751289390824E-008 - 8.6676120648817015E-006 - 8.0356617068274680E-006 - 7.2803510882423364E-006 - 5.8085518216496768E-007 - 9.2454375990990717E-007 - 5.0773576231745913E-006 - 8.8186685457599672E-007 - 4.7560206227399793E-006 - 3.5091571969479065E-006 - 9.9953170045124347E-006 - 9.3802020430725106E-006 - 7.7421618244391093E-006 - 8.9996420942488742E-006 - 3.1463253304570141E-006 - 3.2687591639620619E-006 - 5.1717798740693074E-006 - 6.4958342357222357E-006 - 9.5186252820914489E-006 - 7.6518095003100796E-006 - 3.2812681864244242E-007 - 9.1562890889687178E-006 - 1.0409405655086858E-006 - 8.0620348877191582E-006 - 1.6782210532972575E-007 - 3.2571291016199368E-006 - 2.3736398566310537E-006 - 1.7016508496239891E-006 - 2.3788453676848900E-006 - 6.1617122533147287E-006 - 7.1484531810520677E-006 - 9.5497968190890952E-006 - 8.2340408908746413E-006 - 5.8704101296757954E-006 - 3.0673987141744399E-006 - 9.3837407263787899E-006 - 7.2295545278004454E-006 - 2.9070689091327020E-006 - 9.6073295570431147E-006 - 2.8602787183258172E-006 - 3.9450734004638715E-006 - 7.2630651128304849E-006 - 5.7290243087668742E-006 - 8.7004472608817813E-006 - 6.7896014988095546E-006 - 9.3612949075366715E-006 - 6.5830712337494538E-006 - 7.7175105697898607E-006 - 4.2489019964479874E-006 - 8.2532265919616471E-006 - 2.8038391163108558E-006 - 4.6679871676240708E-006 - 9.0353295696701374E-006 - 8.2258827833684375E-006 - 4.7323995047753101E-006 - 7.7694070407929929E-006 - 5.8800097387830701E-006 - 9.7763354774013993E-007 - 9.3590576764904533E-006 - 3.8050773452382841E-007 - 3.6288606868416708E-006 - 5.7665411023096473E-006 - 7.2936819781210018E-006 - 7.5728515196123091E-006 - 9.4073164054273133E-007 - 3.2918829668373264E-006 - 2.8807781501643682E-006 - 2.3991267616282710E-006 - 2.1011887257152587E-006 - 4.8791580977329560E-006 - 7.3178893639086346E-006 - 9.3185350920641632E-006 - 6.1813531942559649E-006 - 5.1232563881662119E-006 - 5.1807841203848634E-006 - 5.7139457032909507E-006 - 2.7218008967915102E-006 - 8.4943876582210064E-006 - 4.4609164920857630E-006 - 2.3021677607560600E-006 - 4.0668976702136067E-006 - 6.8815479592920528E-006 - 5.8080046867981451E-006 - 2.2170266764376936E-006 - 2.3356418926271030E-006 - 1.8011678022604884E-006 - 9.2100937238280369E-006 - 4.1152943125626077E-006 - 7.4497587201785084E-006 - 3.1813647145382630E-006 - 9.4722110615235561E-006 - 6.6386567064107567E-006 - 5.6801030565481600E-006 - 6.0212964458227239E-006 - 1.1972487134288313E-006 - 9.2307385628462531E-006 - 9.0716788490252633E-006 - 5.2593673365604233E-006 - 9.3203434061147718E-006 - 5.9363305760960832E-007 - 3.8533607338406785E-006 - 2.8126977989834692E-006 - 8.2414758692229325E-006 - 2.7334699321976242E-006 - 2.8058911457102956E-006 - 7.9727789587638969E-006 - 1.6959666483085092E-006 - 4.2649145238949884E-006 - 4.2317227749571857E-006 - 5.0721716921555251E-006 - 1.6703056744866219E-006 - 3.3222465510053716E-006 - 7.4150782823261925E-006 - 4.9777703477543600E-006 - 1.1324488856773841E-007 - 7.3734931204165550E-006 - 6.7146374802920735E-006 - 4.5274489315057734E-007 - 2.0927916702039198E-007 - 4.9481395800683938E-006 - 5.5432647885081415E-006 - 1.1479965505291358E-006 - 9.5176829741217779E-007 - 3.5176711526672190E-006 - 4.3972188554581350E-006 - 3.1299014912162051E-006 - 2.9980591890211972E-006 - 7.3516784769216876E-006 - 9.0548875247421690E-006 - 4.2908830526258216E-006 - 7.2889483073860401E-006 - 6.4553073295146899E-006 - 2.4698222954374739E-007 - 8.0447821265135427E-006 - 9.9957528908258199E-006 - 9.7875338493318533E-006 - 4.2954684342167687E-006 - 8.7944091592383499E-006 - 5.3915503615560094E-006 - 4.1195470246451095E-007 - 9.1116702812983235E-007 - 9.7916822375035039E-006 - 5.2689459387834462E-006 - 2.5820890208911009E-006 - 2.5652924144170977E-006 - 4.1330712253845415E-006 - 2.8010198437469162E-006 - 1.2334688713102161E-006 - 1.3673216695937508E-006 - 4.7039302864866610E-006 - 3.9295170992148809E-006 - 9.0355529469744078E-006 - 5.3349499988383515E-006 - 6.9713787631783004E-006 - 5.4096545313859213E-006 - 6.9957488343259771E-007 - 1.0692315314445301E-006 - 6.1491945482560568E-006 - 3.6833248487244210E-006 - 3.8429670710544430E-006 - 1.7797384199507771E-006 - 1.0727027019526682E-006 - 1.7571777123010169E-007 - 5.6160677511581427E-006 - 8.5900360682434126E-006 - 3.3195927064466061E-006 - 2.1518441717076622E-007 - 4.9485007898054093E-006 - 6.9060919658316942E-006 - 7.1556682038077009E-006 - 8.1382556329416546E-006 - 1.2296450155594618E-006 - 6.1320606241899747E-006 - 9.0656903683155766E-006 - 6.9059410698403219E-006 - 9.5522458033691963E-006 - 7.5523545819834385E-006 - 8.4361944144224205E-006 - 4.1830151811710349E-006 - 2.0301530037840844E-006 - 5.4974740557041274E-006 - 7.9824777762066499E-006 - 2.9864196608947650E-006 - 5.4079306219120050E-006 - 7.5135698219664917E-006 - 1.8887600493820079E-006 - 6.3169325997806076E-006 - 8.4688570465885979E-006 - 8.5029078619960082E-006 - 3.3771052836482699E-006 - 8.5142633382212922E-006 - 8.0353847284950897E-006 - 8.8266542396942208E-007 - 9.7217005981712610E-006 - 1.0970499359459641E-006 - 5.7935288853447345E-007 - 9.1066789688770869E-006 - 6.7847012962133719E-006 - 5.0347323374549110E-006 - 2.5437864364323083E-006 - 5.3258023085156840E-006 - 8.3822906283691646E-006 - 5.9458320372612766E-006 - 7.0322476744309002E-006 - 4.6858139485818891E-006 - 7.3808599108874224E-006 - 1.4170854674197466E-006 - 8.7367866408382490E-006 - 5.9007462955448421E-006 - 3.1925302261879177E-007 - 2.1233633762369765E-006 - 6.4131473584731955E-006 - 2.8565127889766462E-006 - 1.4142656136639843E-006 - 4.6524149382595210E-006 - 6.7748715139604748E-006 - 6.2762670399471801E-006 - 6.5449681844295607E-006 - 2.5210057620847694E-006 - 5.0457993714846989E-006 - 6.2076383473705955E-006 - 5.6675644776084978E-006 - 3.4756773268587984E-006 - 8.7150289515203705E-006 - 2.6488324372265896E-007 - 3.3142676114723273E-006 - 5.3092362751209871E-006 - 8.6174091869786144E-006 - 8.6170920806065731E-006 - 9.4068476564332628E-006 - 4.8587286174710516E-006 - 9.1579112007577565E-006 - 6.8272942157492783E-006 - 1.0968888222089414E-006 - 1.6558222738749553E-006 - 8.4890120975622906E-006 - 3.3344374623189633E-006 - 1.7749536724856121E-006 - 1.6183411352739386E-006 - 7.5647029474941182E-008 - 4.5964289049907681E-006 - 9.7700328745461568E-006 - 8.3196969788377080E-006 - 2.8036763959393811E-006 - 7.2771024434559390E-006 - 9.3610024976622848E-006 - 8.2968680751815544E-006 - 7.1614603714475151E-007 - 8.9907450021787470E-007 - 6.2509546104340225E-006 - 4.6650462304913168E-006 - 7.1533038027811327E-006 - 5.9567736089789061E-006 - 9.7567540418407530E-006 - 4.6303547201364085E-006 - 2.1865371883583364E-006 - 3.8783071327148050E-006 - 5.9024413881056744E-006 - 9.1331648373527147E-009 - 9.8170147634459921E-006 - 4.2713511210451353E-006 - 6.1421389593279521E-006 - 9.4283547615098621E-006 - 7.6419366019089262E-006 - 5.7508122822917448E-006 - 8.4037892596697396E-006 - 8.2752308251127734E-006 - 7.4526283022058129E-006 - 5.5953950117570252E-006 - 6.5570654968112582E-006 - 8.0363549364410482E-006 - 6.5427164034340539E-006 - 8.8376733613674839E-006 - 6.4971476405215343E-006 - 4.0958307971730307E-006 - 5.4499819902860655E-006 - 5.7086916806622353E-006 - 8.5640123991328676E-006 - 6.1999247929754599E-006 - 4.3385279521170030E-006 - 6.2429200584044782E-006 - 6.1641335204622932E-006 - 3.9125348211047062E-006 - 5.4160280011500201E-007 - 3.4509850667053258E-007 - 6.1067326889398401E-006 - 6.3127579244976075E-006 - 7.0448872713331128E-006 - 3.9693658236536543E-006 - 5.4096593496557126E-006 - 7.6655385203150316E-006 - 4.0053297588244707E-006 - 8.1164220186892628E-007 - 5.8581318932685470E-006 - 5.7821363417527113E-006 - 1.0404711345236395E-006 - 9.5214305728543610E-006 - 4.1097610262129050E-006 - 9.5741677670723699E-006 - 3.8970912371438810E-006 - 3.7370379299487567E-006 - 4.1148673934621776E-006 - 5.7491572956776794E-006 - 4.9363542980599851E-006 - 2.0833641858562070E-006 - 4.3860870396094981E-006 - 4.1216329901786842E-006 - 7.9188372976502176E-006 - 3.1880823591282083E-006 - 7.5039102997899273E-006 - 9.6308626124475831E-006 - 3.1455537514544332E-007 - 4.2987565716056611E-006 - 3.0081580961684265E-006 - 4.5902134135375609E-006 - 9.7054737029855315E-006 - 4.9820310740211942E-006 - 5.8365573058147006E-006 - 3.7943229693631934E-006 - 5.3809618360924535E-006 - 6.1876738966907268E-007 - 7.4871851772067616E-006 - 5.5647256715266392E-006 - 2.5326181670632476E-006 - 6.2320874994514606E-006 - 9.8105210231680188E-007 - 7.1707173461081185E-006 - 2.3091580664577772E-006 - 3.0173915430344412E-006 - 6.8831583808246685E-006 - 2.9698824721568452E-006 - 8.7285346322726131E-006 - 8.3558422656395131E-007 - 1.7142672672676652E-006 - 2.6157511707064798E-006 - 1.6058365500638529E-006 - 4.2017274901290306E-006 - 3.3960264694119892E-006 - 1.4577508452807142E-007 - 4.8977024081625744E-006 - 6.5192036290515484E-006 - 7.0768532269719766E-006 - 9.6202613053151672E-006 - 4.8623500590195152E-006 - 1.5029688670805898E-007 - 9.5426827870006070E-006 - 7.9866969247833483E-006 - 5.8760150633193183E-007 - 5.8360723228560438E-006 - 8.3583251614909151E-006 - 9.7812771343452563E-007 - 9.3060265972534443E-006 - 8.2923686847665385E-006 - 3.7198212604682728E-006 - 7.2899053267077302E-006 - 8.9956851221067346E-007 - 3.9630182224978828E-007 - 4.6098317204160789E-006 - 2.2268463921917383E-006 - 4.1855226799259737E-006 - 9.9215935024748208E-006 - 7.6421283412483165E-006 - 6.1372582684420931E-006 - 2.6293523668669907E-006 - 1.9102447248883791E-006 - 4.9686962798366756E-006 - 8.5349871440823083E-006 - 9.2663917426525807E-007 - 8.7627418758392254E-006 - 2.9596385845705502E-006 - 7.8239807591789071E-006 - 9.4030789287705184E-006 - 6.8463930899688202E-006 - 6.1276853324604636E-006 - 2.8351141592305933E-006 - 3.5223007241232609E-006 - 2.0038358215186095E-006 - 5.7148718813902273E-007 - 9.8629942006308705E-006 - 1.7290023537211398E-006 - 3.4783844841073200E-006 - 4.3453750985071297E-006 - 4.7703369792323694E-006 - 9.8584950083110824E-006 - 8.4275855816088520E-006 - 7.0800334533836688E-006 - 2.7618768501060899E-006 - 6.1763550487575899E-006 - 3.6402824354513657E-006 - 6.0955205338064797E-006 - 7.9699424802248899E-007 - 2.2912659604282660E-006 - 7.3429440291633079E-006 - 9.9919851099089768E-006 - 9.5271696258345175E-006 - 6.0339590101256983E-006 - 3.7257120498249598E-006 - 1.9245141156775925E-006 - 7.2792156990022743E-006 - 4.4178798435237157E-006 - 4.4134503831591011E-007 - 6.0411133579401546E-006 - 6.6931935494002001E-006 - 4.9381158183378499E-006 - 8.7634154191481138E-006 - 9.5889555823886713E-007 - 6.6205666208233290E-006 - 2.0481664749723285E-006 - 3.8512063667030722E-006 - 9.4886697195597840E-007 - 2.3582803351473328E-006 - 9.9542833807752867E-006 - 8.6210497933333799E-007 - 2.1272056919086480E-006 - 3.5145660903996325E-006 - 1.0439834581823537E-006 - 1.6042870879019822E-006 - 9.6664202128684810E-006 - 2.5906972627440753E-006 - 3.3744422640900877E-006 - 6.0304878062134279E-006 - 4.7051731960480675E-006 - 7.1342898689939107E-006 - 5.8515515345103955E-006 - 7.3211837677334217E-006 - 5.7665848906819564E-006 - 9.0671223127506870E-007 - 6.8945100576994215E-006 - 1.9012833318032520E-006 - 6.3998738799253972E-006 - 8.4655231200517026E-006 - 6.1489991383975100E-006 - 7.0459816304512793E-006 - 9.6567931566437236E-006 - 5.6149855941763019E-006 - 7.9313150855101690E-006 - 2.5651935332197176E-006 - 1.2703115275606792E-006 - 5.1469668735580579E-006 - 3.6738447705898550E-006 - 6.8880055087411098E-006 - 5.8722805842964899E-006 - 1.4259943015099943E-006 - 5.7327092092966149E-006 - 4.0521053921986493E-006 - 3.0803362623225218E-006 - 3.1520666886477800E-006 - 7.5766806904470469E-006 - 6.8594900104155570E-006 - 7.1519793539062772E-006 - 3.3035087988884263E-006 - 8.4408849539486463E-006 - 2.2490612554058443E-006 - 8.0776897099750277E-006 - 1.9325320303639117E-006 - 7.5459287002044123E-006 - 3.0381414541275187E-007 - 4.8891196058571679E-006 - 3.9412958858803773E-006 - 8.9572063321824569E-006 - 5.6989883556243965E-006 - 5.6702544140060158E-007 - 2.1500452786956363E-006 - 9.0580822334166944E-006 - 5.5987658990316576E-006 - 6.2197232543354757E-006 - 5.6344516913760363E-006 - 8.6350347705369435E-007 - 6.0817223543351973E-006 - 7.8523759216874492E-006 - 5.6759309724299059E-006 - 3.5987208579869303E-006 - 6.7089307367273260E-007 - 5.9461196924942289E-006 - 2.8443842636596329E-006 - 6.2974372473781606E-006 - 1.0384046015580185E-006 - 2.9370504176287316E-006 - 1.5165889270241629E-006 - 3.8508346276775569E-006 - 7.2988958618404793E-006 - 8.5599276115660436E-007 - 3.6487265442875931E-006 - 9.0420180650845572E-007 - 4.4034978702604427E-007 - 5.0719889405213742E-006 - 5.6162035604853760E-006 - 7.5815389514699446E-006 - 2.7272643433954281E-006 - 8.6225826643504735E-006 - 5.2694535292707350E-006 - 4.5169978681277906E-006 - 9.7290798278468093E-007 - 5.3297249655747790E-007 - 3.2467298364985310E-006 - 8.0875780633561384E-006 - 3.4720840598493598E-006 - 8.1644167324332534E-006 - 1.5115391057566187E-006 - 7.0755238430751696E-006 - 8.1028296917690098E-007 - 4.0610181471893182E-006 - 5.5876003745966800E-006 - 7.0333832743440889E-006 - 1.3710455955214475E-007 - 2.8607543576234185E-006 - 7.4595255055165935E-006 - 3.0952856246177189E-006 - 7.0703471222053847E-006 - 3.7877709134118357E-006 - 6.6524051486223910E-006 - 8.7273038009468250E-006 - 5.0819220542005674E-006 - 8.4282767132542081E-006 - 1.9755270938695194E-006 - 4.8993835653672672E-006 - 1.4347479655538742E-006 - 9.6204834813192731E-006 - 7.4823375084742042E-006 - 5.0081381568443685E-006 - 3.9968306770740483E-006 - 1.7598053240707513E-006 - 9.5821337985292856E-006 - 5.0295506062042474E-006 - 5.1551678188221930E-007 - 5.8839913858864626E-006 - 1.0696363937074871E-006 - 3.8547553258217392E-006 - 6.7199028484472783E-006 - 4.9502559209942321E-006 - 5.3215583404135360E-006 - 6.9514475280658124E-006 - 8.7089277933322860E-006 - 6.0232494163330092E-006 - 7.4240800253452562E-006 - 2.9978153326255978E-006 - 3.5800898834011922E-006 - 3.5238809658236884E-006 - 2.8782688610685826E-006 - 4.5298950025128470E-006 - 6.0268971683547160E-006 - 5.5449080284348023E-006 - 1.4467532371413028E-006 - 4.0509268590054733E-006 - 5.7698511503747609E-006 - 4.7846203399041154E-006 - 1.7949434486131189E-006 - 9.7826536528555668E-007 - 1.7181269826929669E-006 - 5.8937029099159029E-006 - 2.9461422175346343E-006 - 8.5730804571858347E-006 - 8.1090275681754034E-006 - 8.2032105883622874E-006 - 8.5030912756547757E-006 - 1.8796998403626998E-006 - 6.8523796509765145E-007 - 8.8772546093851193E-006 - 4.9741972173673956E-006 - 3.1962262033093596E-006 - 6.1486118655297065E-006 - 1.7319836262930490E-006 - 4.1423448344554521E-006 - 5.4337346634496299E-006 - 8.4880281980610252E-006 - 4.2587893265138960E-006 - 1.9408936157187851E-006 - 5.3789292685536302E-006 - 7.3336518392311825E-006 - 8.9265484261212724E-006 - 2.4720039212752645E-006 - 1.3577213319486061E-006 - 6.9687325646473185E-007 - 5.8908525106763648E-006 - 7.4864419193612421E-007 - 8.2759215543230726E-006 - 5.3896270751065882E-006 - 8.1724742380708433E-006 - 9.9614341173527291E-006 - 5.2409203436121667E-006 - 1.3225416121458633E-006 - 7.0140959847017602E-006 - 1.4577516510457102E-006 - 1.1446078321404830E-006 - 3.0495592155207821E-006 - 9.9831890097100853E-006 - 6.9169260695793047E-006 - 7.6903545238989234E-006 - 7.0170216281443378E-006 - 2.8985706717687524E-006 - 4.8034855029188520E-006 - 9.6648449942650016E-006 - 3.2612278919133923E-006 - 7.5126150279838172E-006 - 1.4629981249208470E-006 - 4.7707641024284597E-007 - 4.8607566368304429E-006 - 3.2940601969705756E-006 - 8.2554287543526119E-006 - 6.0900458089054641E-006 - 4.6945735434236276E-006 - 2.5225347965076374E-006 - 9.7253524218168246E-006 - 3.3437404234577854E-006 - 3.0160561321606308E-006 - 9.0181085783816702E-006 - 8.0527334323919788E-007 - 1.2954655266924065E-006 - 1.9130312090834655E-006 - 1.1670833747792409E-006 - 3.6820963102448237E-006 - 5.8932060940729338E-006 - 4.7772644163706486E-006 - 2.5155934275576730E-006 - 3.3019204162939577E-006 - 4.0092128915163453E-006 - 9.3650421790753715E-006 - 5.6438112979934247E-006 - 9.5343057080408921E-007 - 1.6639324140519476E-006 - 5.7071787452969039E-007 - 4.9955903491799264E-006 - 1.9765826216985983E-006 - 8.7444261911886600E-006 - 9.2021824672544231E-006 - 3.9747385593103783E-006 - 8.5217915542331144E-006 - 8.3389370128126052E-006 - 1.4584911720492556E-006 - 7.8464448307670910E-006 - 5.6584938679990406E-006 - 2.4750709327596567E-006 - 7.1567968520532774E-006 - 1.2429211160735112E-006 - 7.7740705398940642E-006 - 5.8840124189581635E-006 - 4.4742871172840594E-007 - 3.0896226853820411E-006 - 7.0669456310584403E-006 - 7.3937120368012246E-006 - 2.0696119067652297E-006 - 7.5510000201288920E-006 - 3.9681929329582481E-006 - 7.4538733916146347E-006 - 5.7676424566192678E-006 - 3.3992605845233283E-006 - 5.8672253112392711E-006 - 4.8192937607780887E-006 - 6.2237359627088192E-006 - 8.2359300054629633E-006 - 8.6696415490972341E-006 - 3.2381201507407202E-006 - 7.9238558153627571E-006 - 9.8238791298868292E-006 - 4.6591969827824099E-006 - 6.6944669835806421E-006 - 9.1777089370095545E-006 - 5.5712166692827865E-006 - 3.3514838170167142E-006 - 1.4631555266139575E-006 - 9.0450506770011889E-006 - 2.4792227612834797E-006 - 7.4120543764922916E-006 - 9.1493707890689093E-006 - 9.8895430853857882E-006 - 1.5577237374578359E-006 - 1.2274032277505953E-006 - 4.6482869558107683E-006 - 7.1184623814928465E-006 - 3.3568788125388508E-006 - 8.3648402639242171E-006 - 2.7281803194386058E-006 - 3.2540458090654881E-006 - 6.8009297824167057E-006 - 1.8932093694583420E-006 - 7.1333672269985939E-006 - 4.1584457282247430E-006 - 5.5746955950587387E-006 - 7.3021663355727320E-006 - 6.2511438121767405E-007 - 1.8532399367838650E-006 - 3.9970486634575545E-006 - 6.2491817353279497E-006 - 9.6514973001220733E-006 - 4.0742471038408913E-006 - 2.8818195322459129E-006 - 7.2168032803314018E-006 - 7.4904169526870140E-006 - 5.5844780867528539E-006 - 9.1586651752828112E-006 - 3.7943051251608220E-006 - 5.6429499275122196E-007 - 2.9601939848836633E-007 - 4.3060330583904519E-006 - 4.7535014549838464E-006 - 5.5812655593995746E-006 - 7.8910807196809711E-006 - 9.4185790685844805E-006 - 4.3836094097243254E-006 - 1.0911914350935448E-006 - 4.7136642737248236E-006 - 9.5945755851524721E-006 - 1.8575116080475198E-006 - 3.0620068517583046E-006 - 2.4796711592762058E-006 - 1.0758270269525850E-006 - 4.4902803294147532E-007 - 9.0824287934223376E-007 - 4.0352533060261911E-006 - 2.5406763923111787E-006 - 9.6964735351798146E-006 - 9.2569320753151952E-006 - 3.9415812369842292E-006 - 4.7925457379795602E-006 - 3.6150047309438765E-006 - 4.0773330600999338E-006 - 7.1093726696095892E-006 - 9.2370218090873762E-006 - 1.8861523316405206E-006 - 7.3342488988522102E-006 - 3.5078919855484673E-006 - 9.4007634983548373E-008 - 1.0333977512471838E-006 - 1.6353325329283109E-006 - 5.4283745120826766E-006 - 1.0154482275442568E-006 - 5.8835257162739222E-006 - 4.2332700175356990E-006 - 8.6019595122777519E-006 - 4.5506366336532804E-006 - 2.9714297389755843E-006 - 2.1502422810165478E-006 - 2.0728182068696136E-006 - 1.1444840386016642E-006 - 7.2559968004449014E-006 - 6.2712972258245993E-006 - 2.1613195276838128E-006 - 8.8466784701495974E-006 - 3.1791166094859882E-007 - 8.1818931630121933E-006 - 9.8551400269488082E-006 - 3.2759613232626888E-007 - 9.9101418920251138E-006 - 9.6632215682464522E-006 - 7.1096768867832557E-006 - 1.6319978996714624E-006 - 2.6565654589378310E-006 - 5.1493418000553985E-006 - 2.0258860693598549E-006 - 7.9322787656428101E-006 - 4.8347291428090491E-006 - 1.0295202892894185E-006 - 5.3553231134387094E-006 - 8.9740483253379236E-006 - 3.1925313990879601E-006 - 9.1703412981822279E-006 - 3.7101295191265974E-006 - 3.5408071776176299E-006 - 9.7883187505898987E-007 - 9.4885812055681732E-006 - 6.0739543275338245E-006 - 9.2624677979732186E-006 - 6.3244161100350910E-006 - 8.3412552384762406E-006 - 8.0877767202348239E-006 - 5.2692994939630725E-006 - 4.1550946252980960E-006 - 2.6900386848119465E-006 - 9.1551349948460507E-006 - 4.0392029296010209E-006 - 8.1607050311477494E-006 - 9.0181633140077834E-006 - 8.6030439901574983E-006 - 3.3841619419229919E-006 - 9.3071195257930848E-007 - 3.7965385810158417E-006 - 2.5406809874012339E-008 - 5.2867644456207680E-006 - 9.0207401787000462E-006 - 8.7883720017929202E-006 - 3.9440212691982502E-007 - 4.6701018635106900E-008 - 4.4890055181903669E-006 - 5.8766169543128356E-006 - 5.4556423322769804E-006 - 5.2458256553284822E-006 - 8.6933761664259413E-006 - 2.4130964267977862E-006 - 1.0947284560192905E-006 - 7.8766916184996048E-006 - 8.5748306256951490E-006 - 1.4682515380584060E-006 - 7.2487387259592140E-008 - 7.4882552332083742E-007 - 6.1714665452229469E-006 - 9.3006616731552119E-006 - 8.9270495802793058E-006 - 3.1960319996013169E-006 - 4.8894879633129507E-006 - 5.5824963153002501E-006 - 5.5913337321623069E-007 - 4.9546038625796278E-006 - 6.3197973300598959E-007 - 7.3907888881071276E-006 - 8.7037299630839746E-006 - 4.0074327688552037E-006 - 2.2265372976455758E-006 - 3.7555631515893265E-006 - 5.9332343284169683E-006 - 6.0457799818582416E-006 - 7.1021977999043412E-006 - 6.3716238544415850E-006 - 1.0176822174006085E-006 - 7.6456382883083031E-006 - 5.9872000559753592E-006 - 7.8564057375470510E-006 - 2.4338428657479219E-006 - 5.8529874532399126E-006 - 9.5948560537506643E-006 - 9.0572712420806724E-006 - 1.0829897093770914E-007 - 5.5116198521580592E-006 - 3.6210517581316730E-006 - 1.8599722748506788E-006 - 6.8782839574117245E-007 - 7.7854216303197701E-006 - 9.4516937709005118E-006 - 3.0624183259898190E-006 - 9.9462811926706167E-006 - 2.5602722212233366E-007 - 9.4085919185321083E-006 - 7.4911433440907662E-006 - 8.3833651569746252E-006 - 5.0190426079542951E-006 - 9.9970216161645241E-006 - 7.9939099628971216E-006 - 1.3427980176971723E-006 - 9.8378069894889988E-006 - 5.4984188172990328E-006 - 2.9747261717269548E-006 - 3.0709948661250476E-006 - 3.0248008209207090E-006 - 5.6471175687948700E-006 - 7.9114628901370832E-006 - 1.7934727473158641E-006 - 1.9161195700971413E-006 - 4.1644061681877370E-006 - 6.9168024336410822E-006 - 7.1439469405878069E-006 - 9.4541176203181022E-006 - 2.2974721639932261E-006 - 5.4495125688409503E-006 - 4.5890577367630605E-006 - 4.4685846936281420E-007 - 6.9166255557825403E-006 - 2.5035067309196371E-006 - 3.3002075937968692E-006 - 9.2498050881214956E-006 - 7.7399674087531046E-006 - 2.8828937387694912E-006 - 2.5144501869083261E-006 - 4.1740631609301859E-006 - 2.7395997020684085E-006 - 4.9980719853880232E-006 - 2.0121192512764553E-006 - 2.2158934179387103E-006 - 2.6792635777367194E-006 - 4.9619535646235058E-006 - 9.5268241364740554E-007 - 6.1638895996524073E-006 - 8.0649751593661656E-006 - 9.2243318818358697E-006 - 5.5530451315918977E-006 - 5.9747927886237293E-006 - 8.0825785573839435E-006 - 5.3418020576429560E-006 - 5.6702335475885263E-006 - 2.7937569409673437E-006 - 8.2012127791223465E-006 - 4.6624464652695842E-006 - 9.1876148131641382E-006 - 5.6779390460022174E-006 - 1.9274398800914764E-006 - 3.9726887546478156E-006 - 2.5324145561984924E-006 - 7.1756751062174222E-006 - 6.4056437817905338E-006 - 4.3096421797891128E-006 - 1.3597009710929222E-006 - 4.1233683811999724E-006 - 3.7336659871990912E-006 - 7.9145040673133181E-006 - 1.6593440410178164E-006 - 5.7324949207856442E-006 - 5.7803368645898698E-006 - 9.8167043770714308E-006 - 2.4610331235207895E-006 - 7.7799344131020319E-006 - 4.9176175416996993E-006 - 6.0396704550420571E-008 - 4.2610979643224836E-006 - 4.1503290808773418E-006 - 3.6574129934713729E-006 - 8.8074628314760261E-006 - 2.4705698248979726E-006 - 8.2455543330362858E-006 - 2.0914509105203874E-007 - 2.2225119942607688E-006 - 4.1453634661663076E-007 - 9.9861842115538247E-006 - 6.7053400876108541E-006 - 3.8213738718154769E-006 - 2.7981729854686488E-006 - 3.5990066813869041E-006 - 4.0795625177334189E-006 - 4.7434095601959959E-006 - 3.5282570273191453E-006 - 7.0876641408921111E-006 - 1.0378466695751566E-006 - 4.0485033009831240E-006 - 3.7780449413287389E-006 - 1.4111475789490415E-006 - 8.9999600859721153E-006 - 6.0977782312550581E-007 - 6.9509153384643040E-007 - 2.8489730181244689E-006 - 5.4387475955551480E-006 - 4.9677616475719162E-006 - 7.8273064452031905E-006 - 4.6273361476156817E-006 - 9.7122903738002639E-006 - 6.5217381174586445E-006 - 4.4888374422779936E-007 - 1.0910739158238471E-006 - 2.1057075394456172E-006 - 6.5643463593145814E-006 - 7.2016191067057193E-006 - 1.7728377625398463E-006 - 4.9268552213974970E-007 - 9.0618738404710574E-006 - 7.6848572004055937E-006 - 2.2919021347331548E-006 - 7.3778546173441945E-006 - 6.7647633689106342E-006 - 1.2704702456574381E-006 - 7.5501719105537208E-006 - 9.5572034066428019E-007 - 8.4464252103265027E-007 - 5.6667124282726219E-006 - 5.7183929358800271E-006 - 7.2084913443054789E-007 - 7.9473913428260109E-006 - 5.1135381493593474E-006 - 9.3325401865813770E-006 - 7.3947884364900228E-006 - 7.3631342853859199E-007 - 6.4124688451152164E-007 - 8.8863577216174942E-006 - 7.8575225434531159E-006 - 8.1742729443385683E-006 - 6.4266124263806517E-007 - 8.0408409367312753E-006 - 4.2656476100603646E-006 - 9.9549107906656499E-006 - 2.1410296360847514E-006 - 4.0166360966619520E-007 - 7.7826693303183340E-006 - 8.4360531407711755E-006 - 8.1350732473109122E-006 - 6.3787644838239969E-006 - 7.9850239839249185E-006 - 2.7033303750663610E-006 - 2.6677996085532975E-006 - 8.2048174169855968E-006 - 7.7539174052747859E-006 - 1.7798894049733539E-006 - 9.4668136173175826E-006 - 2.6238094507658871E-006 - 2.0518739858332794E-006 - 6.9461091802814006E-006 - 9.5333579198823631E-006 - 3.9263223322018963E-006 - 2.2910307404480525E-007 - 5.5484321495958908E-006 - 2.4372771731056897E-007 - 2.7543548709753841E-007 - 9.2449316433872075E-006 - 2.3452199576700516E-006 - 3.3040145886289897E-006 - 7.0539316008862114E-006 - 9.2132390656499274E-006 - 2.0143932271291012E-006 - 3.9513985442912665E-006 - 3.3045800015298579E-006 - 9.4092585714261156E-006 - 5.0243991088406974E-006 - 8.5888372762979150E-006 - 3.6108936208795130E-006 - 9.1037786028275809E-006 - 8.6797952085047011E-006 - 6.2049909421052644E-006 - 9.7002347103223856E-006 - 9.3616287729749765E-006 - 5.4632616316872674E-006 - 9.7008468457909254E-006 - 8.9021081352941680E-006 - 5.8273048461723662E-006 - 5.5313508198790166E-006 - 5.7563492408580047E-006 - 3.6174827299334468E-006 - 5.1064091974378086E-006 - 1.1267257830786603E-006 - 5.1643695295098274E-006 - 6.4832433101455971E-006 - 7.7029629047827180E-007 - 8.4598828948286507E-006 - 3.0837412160781767E-006 - 8.6171850564530151E-007 - 1.9992770409162696E-006 - 4.3932009745401115E-006 - 6.4838135271593626E-006 - 2.2032487376624467E-006 - 9.9450758186931110E-006 - 1.2041241667241732E-006 - 3.7782452470753305E-006 - 1.9352611760538553E-006 - 7.8907234650457283E-006 - 1.9747873164957797E-006 - 6.8633941136977632E-006 - 7.7315656600696432E-006 - 8.3597614538040670E-006 - 1.8213023873064772E-006 - 1.8646352845056225E-006 - 3.3313078966323597E-006 - 2.7630090587072942E-006 - 1.9266812306558678E-006 - 5.0100545255013933E-006 - 9.4721689948853668E-006 - 5.9077872669413202E-006 - 7.7978379170890019E-007 - 1.2043191154853194E-007 - 8.7703413295809198E-006 - 6.8330098747428874E-006 - 7.0525759982688250E-006 - 7.4790206057674339E-006 - 5.2188051897429588E-006 - 3.6772010878625565E-006 - 3.5801025118413559E-006 - 9.5604213539523732E-006 - 6.1348256015774166E-006 - 6.1673723025140551E-006 - 4.4084698068380673E-006 - 7.2213096355415744E-006 - 8.7092894300476089E-006 - 9.0051372318745467E-006 - 5.5248241486719571E-006 - 8.2588104956707316E-006 - 9.4433152971805664E-007 - 6.2976097442807434E-006 - 2.0439584649719025E-006 - 5.8943426454419624E-006 - 7.8566150565667521E-006 - 4.9053650908610732E-006 - 6.5511252347633644E-006 - 1.2881129397886217E-007 - 6.6514277601692620E-006 - 5.0005193320231461E-006 - 2.9339512350674471E-006 - 5.9989143090770830E-006 - 8.2832136595370302E-006 - 5.4053788587545252E-006 - 7.7147325446757429E-006 - 6.1756961342632593E-006 - 6.7568497638126710E-006 - 1.7701472666555441E-006 - 4.4030931202520514E-006 - 9.2793075523899770E-006 - 7.7722437458390214E-006 - 9.1172925020563260E-006 - 1.3268259528519678E-006 - 1.3064550207204318E-006 - 7.3321470808282959E-006 - 8.5108739619454192E-006 - 3.7904504976632106E-006 - 1.0070167731075497E-008 - 4.9221158415753887E-006 - 1.9306551169136023E-006 - 1.7050057736800907E-006 - 9.3598811409357856E-006 - 5.7099048980995500E-006 - 8.5919945079172475E-006 - 8.2795343421958611E-007 - 3.1350956993806859E-007 - 9.8944000195447665E-006 - 7.5687885678191747E-007 - 4.4781579655541074E-006 - 1.2643371281506199E-006 - 4.8563626054561844E-006 - 6.8630129195451531E-006 - 7.9515447601045482E-006 - 1.5228433929167842E-007 - 6.6246050612815873E-006 - 8.2598521334630774E-006 - 6.4095162850159817E-006 - 7.1382807850096833E-006 - 2.9169832946684750E-006 - 5.0238115594255275E-006 - 4.3720188209703066E-006 - 5.3048596983153291E-006 - 2.4497256639095902E-006 - 5.9940925770866358E-006 - 2.9146260033898897E-006 - 6.8146645907927464E-006 - 2.7117863065010019E-006 - 3.6413906913983954E-006 - 3.3053175530287309E-007 - 4.6585902559337038E-007 - 9.3081948381043902E-006 - 8.6286091040796142E-006 - 3.0992569391390281E-006 - 3.6906536969217975E-007 - 2.6351493129068417E-006 - 5.6814049674629842E-006 - 8.6254666377660895E-006 - 4.5033490611198515E-006 - 3.7610397272784818E-008 - 3.1583839236883362E-006 - 9.0801496119848945E-006 - 2.5805079016684910E-006 - 8.7628702297840890E-006 - 1.4702574858379682E-006 - 2.7189065699615733E-006 - 7.1083951901805901E-006 - 6.5226232438446720E-007 - 6.4035856830137988E-006 - 6.6707349938719337E-006 - 1.8445756623320311E-006 - 7.4486505293180022E-007 - 1.4560485613848363E-006 - 9.0354708775863402E-006 - 6.2953177275228397E-006 - 2.2861907314901620E-006 - 4.7759449135379774E-007 - 5.9095339503708306E-006 - 8.4618946517459773E-006 - 1.0232756542813026E-006 - 5.8291153152581235E-006 - 6.6312394927168286E-006 - 1.0091428633736810E-006 - 1.3250374534665278E-006 - 9.9877327613500726E-006 - 2.3112895980115336E-006 - 6.3756228088726990E-006 - 3.3146243010303657E-006 - 8.9615733975543922E-006 - 1.5870596245116565E-006 - 2.8577174207016502E-006 - 2.1075481052671432E-006 - 2.1407084328738748E-006 - 6.0007050777105443E-006 - 3.9108150471029204E-006 - 6.0474715378470453E-006 - 3.4416868154076479E-006 - 1.4152604036617401E-006 - 7.2523903636304901E-007 - 1.6019122223381689E-006 - 1.3376383808919945E-006 - 7.5664621385558501E-006 - 7.3454579920553751E-006 - 4.6350608163851104E-006 - 5.2084064382030503E-007 - 9.7686845263668865E-006 - 8.4478264623560142E-006 - 3.7344678659735799E-006 - 7.1968607502661646E-006 - 5.2542451607961033E-006 - 8.9194326763856855E-006 - 2.3920159034567912E-006 - 9.2592621135113264E-006 - 1.4892778056850622E-006 - 8.1688319608476567E-006 - 5.3353594155786688E-006 - 7.2322185067889536E-006 - 4.2065534289878781E-006 - 7.3145267067486493E-006 - 9.7792026350530453E-006 - 1.1448952539633851E-006 - 7.1680930138367762E-006 - 7.7430321297221974E-006 - 4.3989152813334601E-006 - 3.7887254154745199E-006 - 1.4090649449227667E-007 diff --git a/random.irp.f b/random.irp.f index d5e928e..487982a 100644 --- a/random.irp.f +++ b/random.irp.f @@ -19,5 +19,3 @@ BEGIN_PROVIDER [ integer, seed, (seed_size) ] call random_seed(put=seed) END_PROVIDER - - diff --git a/rescale.irp.f b/rescale.irp.f index 69b7185..ee91ec0 100644 --- a/rescale.irp.f +++ b/rescale.irp.f @@ -3,7 +3,7 @@ BEGIN_PROVIDER [ double precision, kappa ] BEGIN_DOC ! Constant in rescaling END_DOC - kappa = 1.5d0 + kappa = 0.6d0 END_PROVIDER BEGIN_PROVIDER [ double precision, kappa_inv ]