From 1c6991b3237b119af464d09506c730be431e5423 Mon Sep 17 00:00:00 2001 From: Anthony Scemama Date: Tue, 23 Oct 2018 01:13:32 +0200 Subject: [PATCH] Doc: Benchmarks --- docs/source/_static/cn3_energy.png | Bin 0 -> 8812 bytes docs/source/_static/speedup_davidson.png | Bin 0 -> 8082 bytes docs/source/_static/speedup_pt2.png | Bin 0 -> 8335 bytes docs/source/benchmarks.rst | 162 ++++++++++++++++++- docs/source/perf/davidson.dat | 24 +++ docs/source/perf/davidson.plt | 194 +++++++++++++++++++++++ docs/source/perf/energy.dat | 26 +++ docs/source/perf/energy.plt | 173 ++++++++++++++++++++ docs/source/perf/pt2.dat | 20 +++ docs/source/perf/pt2.plt | 176 ++++++++++++++++++++ 10 files changed, 767 insertions(+), 8 deletions(-) create mode 100644 docs/source/_static/cn3_energy.png create mode 100644 docs/source/_static/speedup_davidson.png create mode 100644 docs/source/_static/speedup_pt2.png create mode 100644 docs/source/perf/davidson.dat create mode 100644 docs/source/perf/davidson.plt create mode 100644 docs/source/perf/energy.dat create mode 100644 docs/source/perf/energy.plt create mode 100644 docs/source/perf/pt2.dat create mode 100644 docs/source/perf/pt2.plt diff --git a/docs/source/_static/cn3_energy.png b/docs/source/_static/cn3_energy.png new file mode 100644 index 0000000000000000000000000000000000000000..38450961ab4f9d40cd6b5448a06991540cbaf89c GIT binary patch literal 8812 zcmaia2UJr**We@sNa#&Kia=|Z&X)gfc0U(`^w?N0E0VxzXLRWw^2uhZhNmY%@W8sq|(yH0&yP@BS%ZaNO zFG70F&dzppbP$Qe^78WZ^z_isP-kaneSLi?DJc#P4#*Ce$sQ63aLOQiA|Af;})cFaem~L(b3z z%1g`YJ^*x|PkykJb99^l@B(dh6;r>ArB~G9Z{IL!to+&K*6eO^sptxYU&8S)m|~R7 zj1KC%^dIHl9&m_2Z%wG8oxX?GQe)m&BX~vLA~DA1jJUjFG|;xv0LdTXfPm8}m{z+e zJn@;?{|PWlIVsn1aF9vh9<6n)(l8MY=1uaLA;(T^UBTUK8JBh6u)Q3o(JUPjwy4IE zVEXYB;QH!BS?cb5ydRuSQ5Ktf$!CRHQ`Ek}ok&Y8n)BRix$67`2n^2);`GmruXB{9 z=btM0W80|ZJO`&;V zTHpl>5@SC5Z~-@A81aV*?;e&G5@pk@u#+lIqbb54;}0wRs6bkg>|Uw~xlTA!PgU{P zml4sV;Ofxi^^6IPQ@TgtWo1sG&2EOk#b1gDRZkj4SqYhB?N7 z)>y4DL*c1Tu?~VNUdSGh4rxE2>%cVm$0MUUwc_ z?0WYe_5P$1$9)rzsDW|SMO%c-ag>^vbLZpSv66vjZvAw{)C=pea#;%cmwl(mI6s23 z4h<*#F*F;!fXq)fT;SYZoAsR$8jICruWS3wh6diJnkRdyon;0~eSE5q-_0M&ev_*# z?*G=v_|~%?J3h9w@3c7T32Bm@+qCAU&hC1*%RIs8J0(6v_Yezk_)77#jgIV1NmbKwiBf(%KS5?bD69dS1^?#cZ&j}!)6wmQw_+qX3gZ)G&NUxoLG`I^kotd~0G z)+D~k8h(}Rd9!S<>NE;d{wdLNK?N?Es{`w?c%*V~`q3M1;_SWZ@I#|Z&lPC zdT}!wA|^$}I*IO=FUY<#X<6CunEtb)Lg7~Cw{PFT{Frvr)wiw zZfJrx;fDA^SU%&obS ze#V{sM?IGEl$?OQ8|k89Yz@PT3H~Wu`R5u2TLf681wyO7G2bYJr8g+A+K&?^7)C|Z zuSL(4Z#-2rM3O)WZEMt2&I5GznPla+o+%1^hibflU%}pA4J%)1?ITEG)X5BAyNu6s zg`?_;XUG}rDSp2%9lG! zd3MRNas$?kIKzvK?Bv8z&wl-o30>DGPZeMVWu3?7`GEb}L$}=;*Dkwt{SXOzg4xxe z^pfkJPQiAIrGt}>h*E|jCGeqF&0^qi_}DYo7(BYteI4&Y?&;0tH{Q151k9*^zc9#k z(91YlC8ej@FD+YnllcirItx3t$?XLbLMA+4G1RJ@9gy4ccywy4JHfxq!Ab;tH4oxU zrBhg#+4=j`kf$m;BPPV=LFiC5^d#&CGY+fXHwXzOo0OZ+y;elQ66e;a#}t+{in9_X z5MV|Vyp^y^g23p*R3UKAJhXG&<{pHJOTf&BsSR?|WNBR21DSEIi!PnjG48czV6Ei8 z#!JgDKB}lT7XbyOSRsez*65o~>@J?-@E$Y5abQG#;nNB0w=95}>NV{`+~vOU>j-?^ z?TkB^6gc)WwKEzX@epiM(rHh9NV$SIYH5RxSu+0hV?OEH8XYFCz_T7HHn%ZvgnznjC20cSjd~%Z5G|z)yH#id z1s55H8MISlk(|*od-W(xklH!~zo;8bh0q^fp`9rF&9^tY-RZA7kgr4z7bYSG6a+%Y z8srUILtg#J#Rxl9VTIh(D*m!36YQ>C`i&ZXzRxlbg)DJ84D?Y>^le$<>ycc#KQKst z)x+b4&7vqZbYo*`QRfHBYT|>AeSMZ~GP%@oh$s8z zPtL@Mc;-}gjG$p`87PF=EZ3nC`1gh2jc@Z_D5bn%Y!RptPxDr7jkHzY2`5`%Gtp$n z>A)DCjoC12u@=XJM~pryys$!OO~W4A44EUaT#R9{YGjji)NFJ1A0Q}o;+jl8GDJ~*!U4i zeBgkf90X7+uIOeXlwKr z5yjE_tB$N`gqa6>!R9r2QzJI6bt&+{z~KN+Pbun~z&!^u`Hjd{`FpRZKUY@5E?m$; zh_nt+bx68!&Hg<9?ikY@%4LYQqoag=;OWVn za+u}Cv~T%8ksW63EDn$q2v4#%ABvFWM^_CHIlL(JiR@pq0|%)?SZ$<)*jSr_C5EmV zSx2h#T8MygPNHC>yls0XhQ*GmH$0fdT>+oeSWHNgfrNR(g*m5k|K`1^cHjvA!pmq$ zgF8>bYZjSq+s*{RU#}&vJ9A7ml{Ha%j3p#-{|AEyMX^D*P%<9r3dMEt`7zdu8)WXnj$z*n_< zU_ncq28d{%qKQCo&Cg|)_3*g3)->E9G}k`cv@YV5ITeoYS#PxG-Z{wEBD8KRz>do+f1y>L>K?DBlz(tN^=h;0wnD3xWz-C_*WHlxD9f}5M7 zr)wcX|9Eh17CQxJ3}JYwl0_*ci(8|gJ_ORD5@9lYkRI8=(+zy=qxq0=SNqz<EqkheD60ZeJ`fB1S9-K!S!;O2xJ%m0YS3T(;q#SKQ_)}_4z<9 zq~p_Nh3;hNP}O>*?A*M!O*(jau*)}#a{1_U=c|)jK~Z>2SAhzJm#`!?&M=`QW`vs#5EPLFZY5MHxda18c)omD zD5oToNGYsioIz~)cc}z!+VL*!S53Tl+MlNnJn#IRR`=lZp&SbzREm|2<4q)zwYtf% z;i$9pEAIh*>W}17vjz)B$&3fSs-o^>_hvKN?;RP6AOnEhgTsrHEcpsZv9zc?M)sy~ zUqy?EsOzgl%8{L|Dwu!T;8q$QoOkMQU#6S9Tveg0+{LPORxAyrK45WZu-Hp|eC%Cm z7n@2Rj@jg9B5-(dI`F7yiKVp-N2^({zBz}WuAI!CQ`0%PE}{O-!KIHzxh4H=qbh8s z*@*HsTaFYenpVoe`^HCxBLgp_`UuA(MS%&cLNPz>R{WJ!#eCFg~j~wUS_1I|xn>=z& zHUqGP)T5kiVz9vH)=@4bTYQrd4?_5stvWOuQ~>QcPnnq4Q^( zp>4E&FBZ_SSDA-D=PjW0VobPq`2v!g*z2iTrRv(Kmw~~We@d8&@xS}^rUJJ39B$$F zO%H_zrJmiVqc{e>%|h5g_RFhvg$D|Je`lL&=c5%eJdH;v zO%s72fkW;5c`7`Umg8H+!4Ujhct&doPmX7M;e!{*I`Ao(y}(W+5=*0^%F{PcJB);}d$ypGajB($^GdANaKavgkcGj%2 zSU-&-*F4vwy(;F{i;DDYv~UonAu~Yyx+|RLc_oC$KzEB%T@dCc#2|rFy%OyYd-OK^ zA8~rBz{;-#W+rU}FZ`WUe#~T#X9zJ?v^~7p(Lc)9>$S@Jg4c=eZJX(*VZUN;dl7nt z(uqv4pQ#xO&3?eA{I0%tf^bhQIsIFvI?0N{j-=v#M>6g)BXK=xT<`@7(;w{JEqoi2mZ|IVvw+Fj-e5u&%j^_=F9x37UrAl}-E~mm z$T1~Y-f;a7dG*H6VVag5<1#SLv} zYQwz|++9~>c-L4e*%Ydz z*&D5|gz31?nU{=>wB1yHP@#*nDTICR{BE?s81_NuG$ZY2kG$nbSiMX~I);~$O|bacH%TLd zZF<79Z7;>{6wFylHvVE3;Ly6w%2K;}-U1|SxuaM5rDKOh6{!)|w#RJHUQgL@7E;!) zwh{}ndA>^IyuK!LTio(zUrl|USd=_t8fo%Zlm_`KpFSDon!*=WW1#V6$4lFa?mnO7 zPU9Dvcq)8Fs))Xt7AS9DSie{+X@LfhT{!k-Gsx$jPIS`3z05d(M4vRxhC&w4`(~#` zkIkaMa$)E^hL>YsO4EL#Pd^Ec7(+zFICtT4|7IgHXC_A2wg*L%;?2a$8T zL+}pSNB2WMeYJk=Y7H*v`2A2W#1YcC?W_Tv=)CZG{vKAY_5i^j%S{M;CAz1@R6Y#K z-$pN)?uRcgrT-_g(NuK3KE~$bnczQgEzwjSx5m=u?-ZOG@G1>_vIyEwQ~I%PL4&&% zH2qSzO=Ww)1+0Ood}5 zCrWc9l67qvp?L6I^`zvj&Z_13Elyj8e^JC{>L%B<8*?&6N}NzIw4~ksotAL!$5pHd zXjr99u1taoC;1d1{D{ zno!$Qe+iWsB7*G@SIe^tTLF#1@Qh0x1k;ryi8R4VVqdgNIi z8)`>Wm#=T{1bUU9y<~vEUoebaZ>+?8Y%fSyFd3J=>T(j*TzeWV#~Q2cpVJu+i6{6a z*rwou>RqogP|?Pm6Cwv++K4<9Qx5(0PS;i$^gknz!ES22zG9N?*iP?3ayQ?Du4Rv6N6{vJ( zX51{w%oSROYeV`$)Z?ES)$e*Ccn>8@YCL6jZzaXa^^6}jL$G-Zc|v|}K14fiv(MF12Xx3YF@T4~-5HvcmAk^wS>*pFqU zbKxy;{us7bv+Xc}sKoSlQZXtm>;(v~ZyXm_<;#aU3WG2C*4u|^6?wDYZ_U(E1v|xj z2hZH}2UnG$^eVDRH_W=!X-I>Nv4MQWO3epQnjdW4$zgy??ks0EE4lh=Ofe|o4#u~& zukF1olPXq)i{;baE~e-9`0E5Kb8&jXjk|;NuKR{+KxTxpnpVN=DFUlQxF+$if#tcD_@HnG^3| zygm5?9fqgbrQsD;9ql4U1bOURzsL6MQ66zXsp9@)5fBI;L%I#}O>>+KKVk~7Q4 zsNr{{m8>6s+0Gp)mfT<(qy8!&oau+#yi>MLINBPM<@V79zn?3b8z{(AGT`)1U0k>SjZ{Tr$Hq0YQ@D#RcSmsBcDq$< zi?yT105O$M=ddm@(?juK3e_Ud2TWrA&?w%Bnw$IyWlZcdrolhZaLFiSy)Km1)715L znpzokLiIaF9AERya%A=ZKkz=eW1FM-5h@B8dpy_H(v=5IrJRl9g1CRFIM$(iM;dvO z^`}|Gl|Eo^LAdH2{!q0+YD9nYWx3p^ihI;u-XC7kJO|$)rvB;!2T+<w%AxM5Y3wxU^J^FlcirBI3Kv8UT+cKQpl5O&yY@uGEp5CK&nux?3P zSB-4|&&3Eg`vaD9;5 zl)T-P;|hGKp(5*%3XT+Ew;LVnMN%k;N%d(LBXDGi^SZ^K-HH5Cu_F4{d+|$fXvb>u3qr5 z8l(9ooPo)A!Z;`Fwvh<<=(-RPOJB6DbGgU>2=?diMXp1o%_Uk)BW;xK)LExlRyW_> zY-e|&f0mgnZ@sPZh7O>VN=3q&o%<}Q89j^Q15-eJaeqKX2QykfEfkS_$5;UP zW`5tf!|BB(q!_M#@vKf};&u>YqWVXpM+$sup`1Es?zl8zM;c2#<02g&i8S?#5C2>H z56j50x#r1wuNEuj;&v;G+FeY~UI(P&LpUk}vr>~r#NXbktPgqpL&&X<0x5>To*s$$ zQhVoTiLVs)qbK-PPcN>pud@U) z-#vAQ`CgnimHdjumHP$EQlVq-WJ}nUjq~J+{3x3-{t_UOAIZv+YacA6NgdhN==J|K)O2#!emi22YTi zgh9uSu(>9FG|I(yX{o6~7C$?_UIiPoeuq}&S@h(&!5DbzF7rNnPRmFP`}Vg8&TW3W z$1bkap6xFDlHyI$lvUIVH13%A(yB|e{bclbf9Vq2SdB%q&-(K>W6!wq+!$K1d+sLl zaO!hK`2ea6CB>DcVR^HL?l$v4CxsLmbTUo?M@8F>0gnV4KXzCwCDA4G8F{J}kxmUj z6>g2yAsOrV&gIW|q#D3q9-0vUjNTUgCA4`#J)x5R!lEy2ra|Z1*AZ%r(SUfflo0se zc}vj$W|MqI_q)F_8>vfj2GWdhvu^r_WD($)K0-_ zTm!2P1aZ%K#IqLb`o4J;qJeMc<3h)(*Y}5YTe4hs8zklC#AT+H=9mux9<#=C{ONt3 zsn099VjoOk@5=uBdocfS0%nN2prz}#JJhw$f~?Uv;WfPO2!^{A%W8Lq7H=)Cab(5@ zW^K5f`6^R0w?QUjoYHcJ;(#fzxL!$4DCW(}iS3qS?uw(`D6hoSj-VnwS^*D6T-|S( zmNYoF+1bqP^Y=KZx3!is!yfmQxM`Kn3|BDWZnOmzcM`6t{pehwuF@|R{O~7%&xuD@esO-*!wFu$HyN~kZp>L~`Ow}kpb)l)IIn{oR6@m3@xFzn z?B#s#aRo6hT(l(C25YxMTFb)bIU9NC1%bMgcL6?2bX<8)Q+bs#u^L7x!o38W)Zi@5 zN~oL*F8b6<4KxzI0IUnG3f$}YuEXA|M?Nff}my0b(wNY&KQq4+6$dr(aR(f1(o`#BF4ARX(SBA&}{5pEst; yr{8-Hl5Y5Y^xYTvf)lIJ@5z$2RQ|vFW=nZ!bVg;mAf^;j188d)sJ~amhW`&Gx8V^0 literal 0 HcmV?d00001 diff --git a/docs/source/_static/speedup_davidson.png b/docs/source/_static/speedup_davidson.png new file mode 100644 index 0000000000000000000000000000000000000000..29c9f7948a88274c8cf95b97f31e2c84b126ab0a GIT binary patch literal 8082 zcmdUUXH-*7*Y+gzjtT-IMFFXTC<3Ao1QDrHq=+;@ilBfHIwVIdC@5We6GTNyP>`w+ zM4B||ia;X0Na#Hz@4@@I-yh#v-?P5I-%7GFXJ*e_dtZCcK2x3=8EUgJA7KUnu<2a6 zWC8#Z0sui|q=9RY95Ms&gPM{4H7zQYI+X>?%Vt5Q4+@1M5D0d5c6D`iLPA1WSy{EU zwb@kR5UP3yzNDU7pIzRuL|y8r$fQ#7R4PQptZ z;)1`x%v~s910degf1%fFqHp<@KbIFd~(@WVTfF zV=9Ievvm42Y}e%EWP5u%9*-|AElo*D!D6w_&dvq~24~KkIdbF(%!4}dJCzFZq@=La zKG+&{XdmO0J%op-08*#og)iRz-O>W^>Hv^Pco+|^-~R*I+1YS+7(%s3>KI%?piK-- zf%?=BXZV51>&kT>0NB6o{}3;v+4%twSkSq2@mfIoLIzt{t_l~K1SD$q40|}=I&pSZbV6V zv4AWpO$i}erx-W=#BZJv!P_W?X1Iq~T0G7_6y;hW+j0O*Ou#LdXc(FH?TcW6R=lG> z4LmihUkS3c{r|px1}Ni@sK%Hh5qv;!Sf!<;-b=dUoBB%A7a#9ZcXA zQ9=z8zbqSrdF0itqLIn3F7Y8^ELEP!mtEwTP6y>&;}UnNPdfWhk5)8;uPVigZ-PjTr~;XV~UnujV5o+WPq7{ptrUU2KT`LOP*Dk37=VvPne zvhLQaeE^6jFL9&PKwRguK2tA%xl%73I~Aq~UF_xn8uuf!1BDi?|K{C;+SHzUnAb_{`#I3R@xOkrhipkv)C0qnX z36;N_8Pt1HVs)8Eo8bOH-7t#%`9=;rDiJA+K8|Vi=LB5Gino>#F=*%ZjIikC+z~bq z6Q664TRNj%!(Ry)7oK>o2+)G-k>@cX04Cu~o_%w*$v=o?3f}DhK|D#R^5+5*67pMl z$P%evNwz|gLlLvD&WtxoDu>QGlO!vp(u#w$ZR9f)c#@Da0(-9hng>C69&=hU?AZsj zVuXu2K;}&{kLO4ID>_FpPeDk((+%^f8TrG7Co=fty!*zrD}578yncn_tqw4PalY9Z z?Q-Axs2=dv#LMP`wS2}*bOfgA7w4Y%&xqNEK+i~ddl{yqy9T;QqG<-*AP&Z%+T7CY z?na`*igtU9Q|BSv>n^-OJ9qEo`lo_e?!A_eBzNuGq=#=y?K2+d7A0dWm<($lc*!My z_&N^r7+5tu>`H^uD2kk51SY;^e0wcrvZ*4nVAMeA3hY~9tr2?IbKlxjcP}$+9dNa| zr43?Ia~o!oKKphHz$rqu+IXKz$yq;`a3fpNwt_eQsx9mlsi^dvtWU5{ve6@__G5qx zeFJu{gxqHPWx=i1d3uoiE7oEyevr)?aJ~LQ_9yJez?88(xoY1lTuH$ZU*3Nr3}#;O zrvdHbZzTXxSQapw@?s9s0bG)U=N8R-gpnbsJGLvG?#EF_e~v~~9>3pDK0HLW9;H{gV(1vsxVZt_f-Yym8F`f%}<024TIc;@_J9uPQ{ znM+Kn^A$acz@wv-aUyJ>ibVvagusg;LyO0d!Ue)BG zy^6p;`SLWpt`f8eZ$7bV;*ksxn^RpBBD?z3^W)#x(&`s|=f2yB9@3a&Gf!O)d-E}N ztSH&gl6N~aw+uAhkk7b7GG4PuoS5;`EXh?Ua58UGGB$4$?O$CL(b8O^fqqBD6(tvk zW0nsWzRt@pq4!=K)Qt}#n4#9(Cy-AV+HxgFoOnOlo@(Ikr5@ejV{s4o8A2vcQVXlF zaK}-DVI~2Z;;xWGDiQXVei^6I74+-bG{FIuXk;mp+7pEshqcjt`RHu+)x4)r*DU>& zGVWf<$xM!wn6L(c-?k+c#_3M4p86m_KC^}07&F*m+R}mvdpePMm$EN{P#4P0*kS+; za8SfNjKPkxfNB=97$*;4nh-=f*-D@_;n4ZsKO;;*;Jvi1LA&e(K>-PQAhA!s(&5{y zrjt>xh+eQrqpORq4e3@b(g9XjqSwT}0MS9FgZ2lyNy^W$9 zD=+5o7YVxaE8$#6NfN?za$IAv2cY6*4)FN=7hinnF{-UFsFYsve^dzqLP^U+4oT%< z<1JpJEnoGlzqUTQj4$X4lRNA|wap`Pg)C>pxYLnkO}BWLhRz>LniBjVEh|jNZ?@>Z z_|W?2>geB*nbnvi3m$yVt!x*YAV@b(o^JO=Vw2!kYEs@q&c9p+BDdC`+jAr(7SF|m zty!>w^PzMlbMK~dmu{nz&OIate>^4RCAs8ZRoVWeR#c43;#_s#lk!}MQP~=sKVieb zYP-n>o4_0-RouIpqdB}?`N^_!FHH-PbaQT{zx)mP>`qI_;X-;$V8Wb={kn;ERE53o zT!Q_UNXCR~MO0q``GdK(bE*CMQR)?JW7!RVZfM<9t71I3e(|!Sj@xs9-g#Ry^Plg{ZX9WxV3*$L~vH$b6-rwGi4;qhTg2ksJ2j$@_1 zEKytiw(KSCd$%Po>OQI1kTdbB%oXLy8LuqkD%Fgbt?P?&zOPqRaQkimFVda{GRv#(LL^7j}b8- zwSjKbvE2-dh4Q}zvszyDy?X9BqF(6wHqtWb1?K2%?ut)HQev$xlNveJ&>N1B4qyCK zUU#oc!_DfW>bNcK9oSx#3x__OoUr>A$3c{P4NTx{@xU}p?jPOZmE%i+F@b;LOZ*7l z=q9*g9fu&(Vb>bz@G(mi{x3p*8O$;D;)IU<$G9BC+5paEIryoZgtY^Ll@!E+_Sm?zx_zCK7zW)>=T9^md`63 zs&ROsLwWo(QeeC&ah7j%9mM#NVh2@JBrWqzx;=wG+Y>Y;VVk#p(Vk;2Kk!4v1e7n7 zMlxk*#XvtLv}yHg9dC}ht==%2M@fx_T^!X+JhMpc$HV8$4-p zHoTO7LyYIB2t?}FOiz9^RDhZ_kj~NvKp3ia3`sN&AdLTu7(i`U7=&K_XSmRpjKc2A z44mZmu$%qA)JTO5wY8mSF-G&>--q5*vwp8lxNi`X`$G|V1BMsi9zj5yG|TiQMtbsJ z&mTp9t^}2Qqmhm~|7u%-xB##Gg8%L=D)Zt2<2_DC@L)Sw#`a}9tcptjYT8pt>_T-Y z<%;;Zp?wgPTDm6r*J4aFhAgrl#~H3*P^#LK9^#m=lNw{Nm#6T2BWN;-s7t~9kALFX ziKfBwOa5L?Zo=}?^mu*c#1~5pA6%-|6pWtSUSz-jcNQLEBNjTuyy+jme@{qfnZ`&Y zK#*E!&#=`3C+s$P`>rfa5VEEiR5GpoZ@ey%fA9Rcp?vc4-c~TJsBMp({KCLfdauM? zZv3sMn=VE^l#jt_*{NA`^k1sXuAVl90 ze^mgpD;s!=EgMVP?!HU`P3_O@{2MkEQK^r&fa07ODPp!ae3f|h=beKAym-H$}-B+Xdco?n&g%e4_mc7<5@8MIRIY7kAB^0dnpo9Q%@BpNC`lT{$X?l0^8>) z>s^;C0zV%A!n-stf=)Vxc_)Xu_=WLDf{i{1Um_(!;*<>r)XsWK_xpj2E;MUpfWe%r z9TVc{VWz0h4nfPC!PQ;k_sGEpQ9HrAn^?rcN<|81XO+WtH08Vte;=P$Tq&_QsCJzY z+YHGvZ7J`~UF8=GUgCM~am7cVq%gb=o)EZPMstwFKwJHRq5+o{i|UwaN_N zn>^2+xYw|lTHp2Ja+!b_qD0`DmbLs|R~D*3Gu?-67MJWdKLdE1C#_t3O4xhkyydPh zBbHB$W-%P$Al9bnyn`w2INKV2h0O_!X zz=;vhJ^HoxoGc((>LRLqdw2U~CcH5-^kPfLY^Q*RhXVAH9m4tqC~^2?WCnciRGpoC;dWe3L~z7Y5KMzgVAPCC>a|q_ zitqGB2gc8sog>fZXB7vYEVsP8YQ8>2TyF4g_W4=vR64q?8%yscjX*vAKp*Grx$efp z-B>pFfplbD`R4Va3Yq6qMc>M-R6)(3nf}tBxtVj@KdVFhfR`yhP*ceWsXJMqq3bjv zD;4oPt@`(UH=6999~4q6YJJ5Z=bgJkKOYM;A)kkGg}xL1nb_uXD@|Tjiu9}4&6n>! zWj#;G`Hx-S=v|wrvw{^7s+g1Cfb2H^hxF>5)A(tma6A{;LqCZZtUnHVC+B?k{j>G~ z@05)7(jTfyk;vz%1J2Hd?#z{Lj&+*(!sd-@Tmg$oLB@3w{l`)>9HP&A+i5*H1(}n3 z3Cqemy*PFTBEsIL)8_Qo_Z|~Ek4gr)u!0brizZ*H7^@Zd9U^7r z&ds#0M~!w~bBUT1z6oTr$WKgpo<&J%)&CIl2{3UyFYx6X##Q-i*6m2ATjg`gPKQwG@c(>hAxO>;PJq;+DwKg5$9#7nC6Qn6Vo-LJck|x3rP*F~^gWCLASu)TDQ_HyL zlQ3KI`z!V1gFA}cjm~B#gv$#~HnrRlXmVcg#&RR<3Pi50*Yg_PFGh?+K_2ftV<5&# z6&iM8Gy`PitgMu8Ea!IVtz$w?(;b*la|0Yx%8uon#>;avr+;YLZ)4W;Molq-(`$ zwKNXlIVjb25Q`k}y3zHy-dgpR{S9C=HOUh;RNL*R4eAM@mk!`hQYJdI!pq*$C^w>j z#5G#P@EAv&R1lDrLxk1gCj2QeTX)9wa280{WhGlol_qofiIKWymhT(s!+!tprsB;(8B(JT zvLul0YCCOVa_tp-&LYL5B^8aj;JRkT>1N1oebgSgzh^O?~FD@p16t89?9Oy14m7K^$v#X3;Xg32L|(YThQ(ht+!DKwvW| zQJt+drAcv?7bm2mspp*R=6!VQcGJm}iuY3!*lTR6MphsR_ePc*a{?5P6aFggFXjGj zaj3od!-)~B;ZKXy+%yVy)YGtb=jmB-YN^LcF(V@2E8*lVgNDT>MrHC;f6-sTiq~st zQT|!+C_~+~jCkr&cT;#Lg18Ye+bFz8S>D_%EzLQ%mxa4!@FF$Wpe7UTxmZwcWxBKO zf$>?XQ5;PV31@s8s=3JQ(buLB(MvT_b(xyD5l`x_sa3tjs<`d^Xp&Aal=k&}03zhIIPCfmgfFC3ibkED=5{>rE0qD@+-5 zyM@W4E(A4j_!oIZ0~VMLaCFR!_kYDip8xN0Px(3J#7{TH^|3zWd zTqVul^tpA9R=7I{wZCr$PwzG5Z+DumM9!jkj1flO>nKkrnRJ|zThr9u+)cNay9s0s z>b*mBz-Ft)l)2MJ8LOQ0kx()cs9T$9`RBwDkl9>!9+6z?MM0{&*ZEaB<$uxU7fJ-s z6}A3a@oVwY zp0puoRcvdbty4zWVoo8TP}Cxtai)M4(s&WbYH1Gwd?x9?pQK!Lc)_Fo0vgd(7N%n{ z8!aU+8;dd^@XFqMwQ{EBfW`X$x$gBxtGE*ia*#63IrGO~w*8Uc#lh$qIv0s_TwTUZ z6N;G6N-1*Y0F}fwh(qxzn+hIS`MI7HLq~$J z-~mpi9~YkoRh4^DSnZe6Lj>WaG=q}_OcU<~@vO${C7$5ZZRXlqxS4p?XT=BbnlP`w z2DY#7>_1qpEpuq9TYXa$N6ACKq+jziVkEE4*y^(EPcv04T@w^zZUy||`S61O0)jQ= z)d7X~56VP_iuC3rSn&(}r+p{sLBNYxR_&XhJaRD5vE}-yD)!bf_$>bi;KjUJQ!EpJ zw@v)un5(%7P-=!aCvb9|VT17TWlY@{Fl;R^5Ju(>5WnT#h-O@|_!?Z#IU)C7E;+iv z+R2hzZo}JrioK@Ror%@WUFAo%N?ahBAhCR6mW>gKk@ZFyl|)Yvln*0{4FfcMs{&{ zdvLMfJ9Ks!Ols!06)=w*XiC)n!A$ux?$(Rd5z&j^$26%?sXn8ap>4-25?y;rQ7!c4 zM#-kCQUsm-ciNMtTwVQM)}Ke4n=TL+Zr<`+FS#CSt#&nfwiaR zyP5orNH3p=xs*O?tlHz=x%K~Gz4Lq78x?NgURQlx7PC0a36`2zWl~rL%UxeLD4vRE0#$ILNcpGkRcT*WqaA_IhJwl7-!~x{W_9yANXMy^_Dr z-E=j_Z6dwieVpb;^xN;82C6#_D)(lQU~Mi^ix#~Y+0RhLuD&Qwyf)tHl6-b%VY0n% zNaO-N+`N2#eLhl>E_%`QIRd5D;_IlssI$`x{TM7+FO&bcapBpYGbp7r`BMh;9*H85 zZD{w9XXsH_2WEAiwCKN_suOf}UN>dnZqDC!bXy4IO}SqY-l3qt&o(T>(`Pj>tvFQ_ zEmBu&nekBhM+hdyC)d?q{Xpt8}J z}c#WOeU@KwV($q$0m#;A7!1}P=y zPtTPD_Ciar!T;`e9#y-SO^dcK2{(VfgQm7e_S7j|RE_I9T>E(|`nl)l#=ggO|Nhjw z^xCF6YHyiowicK>#wPILJ(hdm^I@sV4ldR-EBNkhwcLY*2NX|fr(2z_&hqQ2 z%TFxM>bmTTtUUvnvi?GN4c)~A4W73Bqu?fS?Cj(BH_<;lj@-HM>lfzFttxZzpS&~? zag?~-o5Xvxef}%N?XSEM^|$7T+GLKrboT??#~cnx79CYxz90d zqY!s?>gK1%${b!V_^`h96O5a@7V0#*+<;dzmbYDyx9{~JNnEXtwWY#VZ@Z{%Uq*>&kjxbw*v0|B|8L dTiSP2w`S*U6Yj%b;M-uJqh)xhSi?T-e*lNXa9#ia literal 0 HcmV?d00001 diff --git a/docs/source/_static/speedup_pt2.png b/docs/source/_static/speedup_pt2.png new file mode 100644 index 0000000000000000000000000000000000000000..e4bec94e94d63547b12c6d425dc055d572356ef4 GIT binary patch literal 8335 zcmch7cT`jDvu+ZCgbspq2pts=1rel%j!N%JwSY=5QlxHyPgIa1Md_$0y@=AJB!JS3 zNC{Gd^d1Q<0g}7%``x?NUF)3l=UHj{-TR%HXJ($a=edcI9wQwO9Rvblyr!>h27$l; z2n32lP=R+~%u=IZLD|H>LWe@3SbbT6!Tg27PoyeWB3lsps#Kmh;*aH0TCn3GEsJHQUUA(n13 zi6*G2yr~q(&}QT&c0&sJCLm;sb2<1YowQ#I3gc1R{2bLfdETm#{lTy>NITeYcQ1F)cva!X) z%|yJcEa=q2!a{FvF9w78{P}ZwdU`}ego}%dp`oFioE#4i4;Ta#?gxbe`EdSx1Z5cX z1_hkhxab66015;^Q80oU4`(|&As7?{vbl)??Sc1Cejs^yd7wCG0@l-41ztfPnHiZw zP>W+OV1d?K|E3=V!u0Lr15K4=;)Os2Hm+%FSOjHlWYb5!Kg+JuaqxVt@c15Y*`Z?# zoCO}mipZszma|`ZbUjk&{g^>cc$N+h*Lecuon>KRr-ntoL!nIA z5r`KssD|nP|M0sx=C<9pWMaL_{W)3)^ng}T{k{-xC+lE<3(P*t*s23ZtH< zL06cDoHipSIFA=U#h0_Pd*ccyr#7zXFQ}(IEa7Qw(AC4JHW?4-WQAK9!_HGa8Q6}n zOv@|5bjus<$qVPS>M~a0?j^*Dad`w+nSQw|${!Y40{}Wgd|6468(&6_JT`dXaaQm! zjGcvq-~)Pu4SQ<~_qj16+rF96at+J_ogXy4I{#}znjQShc2~qi1Rp%#+yIgzSS~U~ zNBQNre*NBuc7n!)XR`Oy7$MBOZ5tNt^VgJp!Z2`| zAg0FLNdnB7q0kOrOcF_ZXQDah1CHRYIXyF>e)NUXra3cPoxc|SfI*c6NJMM*0kYpq zuBu~S*idMf58nJ!{el1tj?xyD5x4pRiCjhI_&O<7?1fP! zp*)_favuu?o{!Q9nx&cIdZsYm5n3A^g1KaJJ&!e!i>2~`{KhOF=}Ivh3roAX+9|jc z-7Ie1wq0Sz1I*Ll;(#^dcl(KX`l3@XjHd6j8Oh>)lmCIEwy~y10965 z)?IArTV|N9WTYC7Gklw^547gd=*fWtWcaqza< zcJ79b(2iNreXMrVU5)|*mNQHpk(J32qrH!<-zz~iw3Wz??OO51f=;H_X_kx$4uKC!;U`${4Eh%{+pr*PQhZzL*LvUMv`6e-IX|`7gRe!i zOA$=J2@z_9a4Jk|7yE4_qS?PIhx0tAwvZ#p(ta{R(4Jj+ZH($=^$?5j4j8RSVxROn z$JKK=Me6;FAZ_#&_&xKT3v+bPYbZ$~ov$_~yj@8M9{#obNbQX}_b9Cf60a;|Dwo%b z29`OkUMk#7i2P)v|5jGWdUBk_gg@FTJJHu*pk~_Wv9Hc=JVB600W^Ve=(}Yd9m68jPmF>L@deV-bs6FGSArnPdYV zA8<3>`HfwCupAUNZZZrro?MZq1bvdncr562TQll&L3p$u*lpCLMz(d&*ZrrSK^R6! zhL~fP4JZ^yA6g-7xCRQlcHJcZy3&XAPQKb@PhW9D&i1IqS0@U+uD>T~U|JP*4YklU zzNkCU*E;=?&;2zcd-&cVKiD5OW$>{kmh}Vvnwk%<5b0zQol771md>n2!++F;--O@u zOSv2WK$+8YsMyi=s$dgfj>bkg^MrP1X{*yUa8zXMz+P0-3FXQ3NjPfwLTkUo`fE0V znUrj1UWt$`WveQ$;H1!pp(OVTiK++G+T(z5U%pe)lTSv^LefhnC#>j2IlXOpK~5?9 zCYR?5*8^3sQnuf%Q#)$pAe@_w z*t|9EDiE6q_K+^D`7a+0|_Exq)t zd9T8oF1)lMq9YI6?vXxFK-4vA;l436<^tb_Dg^yN zmA6o7k-u(TtF@KLkMvio3r&bTOuU`V;xpMWSTMv zGR5ySPACcWuhb?2w$s)8)E%Ym{|qlda|y~+0}8eh(Lwta-W3C)yQwSz?7c(yuue5^$%rIVN0{THl}R9Lzip5I7{s< z=>uIq&0)s%7IXE_YEE(Q6J5-_Ul;In$ohDWiAmz#v6|f0Ur>TmPyUp7spce>ZvBJ+ z%tn`1XQYj<+X~@~qorzGPs=Gre5xp7=Sx~OulzLb^h*tg9mA#oRg@S45VK%KzXz2x zXnOSnVY7br4ZEn-VcA&oZVbEK>UD$=?6vI(@v34$=9hom>Gq_#-}>?^oja(TbS10h z`K`N|U}i0suOqmj>Y9#=YEDz>)~ul{KbakC;_a^LbO(@mK6#blFEg+T6FE7unl{`n zQtQ1PXLnaU=-ZKwfujT9aX~_3?|u2{VB>9$@oG-WMb@Wfd7i#}mZR>uZD@1#(sh5z z_-C*0K^Vv#(vpCD;NjVMAZWxwLHO$rPahD^$d6|IiLYEA?Y-K)rsrFTWUXXM44kIi>n z*2z=t&GPGl_XN@gRsuFEoq2+=GUHOj5BhVtVy+J3f&1*eO1@iWQRZD_`O;9H^zqjZ zJHIT(+*ogD0Aly%PG#md&dV}V!Y>fij_NKDWLm3eyCx!sjmKmeZQ)rePv^1C23l#G zGcQSDN^iqZ9V@>ME+HD+RjT^#g2(_1vTD%J)AZs$G8@y~>0H=ml#NY5@6X5#;Z}1a z(to7&S_W96sxJca-HT{(4qL=0d#M;gLxg#tE^c zTXjxRINz+YXvUvAJuLOfGL;K!BC~dO?j<9tSEunTH8$iMsKm0tl_K58(URU+sk6iJf+|Ydi#!vo^kv)NnQ9#RU;Q`%k<*MKdA` z!f~}CVZHl&nzQV^*i5xQTlu{T6mXdE74Q!C-Mv%#)8rPWHY~iPKhzk#O!d@cDm5|a z5o{l}d@Lpy1;|xI{Zd(@S=h z%1%y8jLnAcY_3(iB`8Kqr`SB|1xo`vaX;y~=75fYwdl!B3HNLknVJwihhY+6dvw?^ z$Tr&>zQX zLQX2z5zF3?)>g-lAO&a-M2yG_W2}O8=>QYixoPd0x2HB;4Z*yj?tt?xkFPVr+D{T% z&2~{Dvlt(v<|JOshNs~>>iI0`3H<4Sej11;o5 zwh`LpzSp)FKMmVH&fvr!K@?lfywtT66@;lL{v0J{OWoPu8b*Dck)5@V0+5fYo-IH` z@5sN^L50DBg|;g;_7onGpLsJwyGVNf$`BJ>nc~~lnC|DU{A(@4n<6EL`AY;ja0J{e zWCPM9y#cx+>d>wQjj9T^s$Ad-wdexa@WyoM1RUFx@VjMx6>RV9A%FVO0wYIO(N)|Y zh|ZUHfVFu<6kNX`xt_mrXf)Bm!n z*VU$7{vBY1G_!|ZdU0@|w{I)n=ysXkX47i2O4$YhCnqh{TTlyS-urHd&1NQAC$8qiiD40x@@q3 zn1+ZokbK^{bTxO0iNE2!*T~hJv_;|?N$mj}A7KJUkZ^}%!O0_VTCc81;+#{}xRSwA z>F^&ThcWB(xiV8|t)xU(GmeMuSWj`>L{tCwe{q=Ip!eL58y?UaCKkWMCF;z;}C!@}9Iob{q^BXf!ht?kZ381ms zq)78JwpL{zwB;cU1)&gO+n{R#F5*RPB(z0wF}5Tz(6|1p@0@8WS9j>SuEIl`gh3J> zS~olawFQud=Qd_)nALg{2xXyPG9<6#`ex1VjmlAk^Djb?g<>7 zi{ES|zRjFO&9#Ap4~NV%@Q0D*K@r@fE?vAi(Wv$mLjP|n_oyd7-X<^r(V0jKG29hK zXE|D~VR1obK$@s`650sC#kajxJLhM;((&c;7xMV1mKnjU?;^O?L9Ldsnn$NhP$v~L z^(o`7nIp*opGmpnh!l13CkeMXGO+1G`P)H8!len4&lj}sOXrhIjnHZ-B(%On3!10I&{A6 zw_nPFPw9Uau|cP}JEYWeSTkH|)YO<1d+e6Df2P=Dq{0%)*c8Sb_^sP)z5II-GQoB! zsu`HF-J<)`2DWfoBN-0Vx7NDx)`sd&$-jdiBGR}2?Xyk?{N^i*h2!yhcKAf|&o$lE-qVvpK5tX%X zf;XDVA6~;~!91-|qZMz74v9!dZM$L}3>mBPnjV>4PLGic!-#6Ym*~YbMxwYSt~UGNuGS(Zdw{ z={~Gpxq(bf@1?HRLlulkJP?JAGoi)}TjW9#n-s~+du0T!@1?WhU3mpHmgpLP^>u7k z{e2lPy%ws^FQ`rY*-LrIK?3;i35raxy@{5J;0m{KmLH3bJ#nO)Uws-XFU4)o^KR|& z*IJ^d3!eO5UV^?(GG^Mnw^g}*&dRgEkv9DsFQ0H;fIuSJ+Q5unG0v8e3{DCvsNTfR zUVicGx9zy3=yU3<+|PwVPT4FW5O)RYrMoI?y@tY6nHJJGWv5u4AkOb@*@?+;2I$aj z^bf0Qp!iho)pyE!1#6ol&_>C6_%KpN@;<2r9fC5$c#(yuYjtiyH`bybMGp9j!6IIi zvBCPgs*UvZk=|jAyIL?K&2fVs$4%}&~RqJPd;w9$v#FTWT{CIndg^r{lT`f~4z0 zRJ>W^+j!4zoelG$wy{%9oA9_xu~jG9~@im@yh#RT7`ME9MgnbFhQervlWvD{e5J32Blf4>M9K zH!%ySg9nOkojtb4EOp#q`ebb42mi+X;AOh%JpTtY{R7WhzL_n0`F*I&af2n@g^SlC zI88hX9DMnb9SNp*e^~Qhm0PgljlI_==*0z3dvo{etw=mg1Mn5>1@CvvwXUOHi{L^& z;P&2%1BfB z+Edu!sAZC_D+8w5E*(vYP!rr->u9pAY9FN zV60c50kyRUgghc}B@J}clP?PA8kYTDMY2L?bNlgp78C2UD}J3#mz*5Wn4JI%1Jwl{ z-PR^|FI9U-Jp(qKJrg7KXYA;+mriJ>EWmFUinHCE!pMfs9YfF<7QU8w45e3(tV)~dAD<}UQG>Ib>_lGDt6?dBi> z6=C$(_fF~*5ZBcZFGIs!Up`+irm-2ihbQ<|b?1BLa`{BIg`9r57Mi$ghx?KwW1WgJU^_?$pncQM9=oW{7T3=;Y znM0qTIet&Y5Q(4DN?E2NG5vgfd}qKn*74u?5by&(6Ef#mvVzG=0+H3 zp!_dFL!Yh5{-NIC%dZgcAxPC8ocWs&(uF3)@S$nP`Nt zI%?|Q{L|sqnn$Nryq@2~Us$@Cs+7CspzE<>bq~!y*pzpEZkl1yKEbmuZzID^?g;bTd7hH+@VPC0UC}%G&gMC6vT9{RZSK(R zngZt(y59A`Be6NvbfC$HlPU8VV8~g?$rUsA#k9w2wwbjdb-fH-BQ%jvZ#>55TvrRU zp3%bW(?1LNz)8sRX%a0O) zj`}8eXb9AqL7COWvu38*jeYdWILnZ}$8CchCQdhK&s0-Shpz6(x*%^D<|EE=yMgyW z6}ElOC{8$Q>As2G!f0}HQ5@g`_4}G~%+)!j|Bohdk&OgM3$f33WE^>tjiF2XGed>^ zR_irkd(dC)!Z*8VZT1+|>E6(sFSh#3tSB?em>STK^ZE21 zEkCNsp-Ew@^d9?DszNtSUl?e&OSk`WEKcg^(aUiby_UXDMRGE~$Eg!;egsI5QQ>>q zv8e5%N%GG^aqT3w%1*|Y1fZuFJ=ZJ|OD=DJjPm(?Y6%K`#`Hg(NOx0&6bJMQILyF* Q7(gJ`bd0piHJu{=4@>P7=>Px# literal 0 HcmV?d00001 diff --git a/docs/source/benchmarks.rst b/docs/source/benchmarks.rst index 9df2aeb1..fa077cc6 100644 --- a/docs/source/benchmarks.rst +++ b/docs/source/benchmarks.rst @@ -1,12 +1,158 @@ -Parallel efficiency -=================== +Benchmarks +========== + +The determinant selection, MR-PT2 and diagonalization are parallelized with +distributed parallelism. Benchmarks for the [NH2-CH-NH2]+ molecule in the +aug-cc-pVDZ basis set are presented with up to 50 nodes (1800 cores) on +CALMIP's `Olympe`_ supercomputer, and 200 nodes (9600 cores) on GENCI's +`Irene`_ supercomputer. This represents an active space of 18 electrons +in 111 MOs. + +- Nodes of Olympe have two Skylake sockets, 2x18 cores @ 2.3GHz. +- Nodes of Irene have two Skylake sockets, 2x24 cores @ 2.7GHz. + +Convergence of the energy +------------------------- + +.. figure:: _static/cn3_energy.png + :alt: Convergence of the energy. + + Convergence of the variational energy, with and without the PT2 correction. + Both energies converge to the (frozen core) FCI energy. + The plot is displayed for the ground state and for the 1st excited state. + + +Variational energy +^^^^^^^^^^^^^^^^^^ + +================ ================ ================ =============== + Number of dets Ground state Excited state Excitation (eV) +================ ================ ================ =============== + 7 -149.489 186 -149.207 354 7.67 + 123 -149.536 265 -149.261 860 7.47 + 3 083 -149.685 606 -149.404 450 7.65 + 29 409 -149.826 151 -149.547 275 7.59 + 168 595 -149.900 352 -149.626 058 7.46 + 1 322 537 -149.946 655 -149.675 032 7.39 + 8 495 334 -149.972 032 -149.704 145 7.29 + 9 356 952 -149.973 375 -149.706 822 7.25 + 42 779 636 -149.987 370 -149.721 470 7.24 + 186 978 487 -149.998 582 -149.733 039 7.23 +================ ================ ================ =============== + + +Variational energy + PT2 correction +^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ + +================ ================ ================ =============== + Number of dets Ground state Excited state Excitation (eV) +================ ================ ================ =============== + 7 -150.161 107 -149.904 883 6.97 + 123 -150.116 958 -149.849 465 7.28 + 3 083 -150.043 5(2) -149.780 8(2) 7.15 + 29 409 -150.022 2(2) -149.758 3(2) 7.18 + 168 595 -150.019 9(1) -149.754 5(1) 7.22 + 1 322 537 -150.017 89(7) -149.752 55(7) 7.22 + 8 495 334 -150.015 97(4) -149.750 87(5) 7.21 + 9 356 952 -150.015 89(3) -149.750 66(3) 7.22 + 42 959 496 -150.016 75(2) -149.751 88(2) 7.21 + 186 978 487 -150.017 51(2) -149.752 90(2) 7.20 +================ ================ ================ =============== + + +Davidson's diagonalization +-------------------------- + +We present the parallel speedup curve, and the wall-clock time in seconds +required to compute one iteration for two wave functions measured on Olympe +and Irene. + +.. figure:: _static/speedup_davidson.png + :alt: Parallel speedup of Davidson's diagonalization. + + Parallel speedup of Davidson's diagonalization measured on Olympe and Irene. + +Olympe +^^^^^^ + +======================= ====================== ======================= +Number of 36-core Nodes 9 356 952 determinants 42 959 496 determinants +======================= ====================== ======================= + 1 775.55 11 198.70 + 5 169.88 2 288.58 + 10 93.22 1 213.95 + 20 56.86 626.41 + 30 43.76 445.65 + 40 36.18 350.25 + 50 33.67 295.25 +======================= ====================== ======================= + + +Irene +^^^^^ + +======================= ====================== ======================= +Number of 48-core Nodes 9 356 952 determinants 42 959 496 determinants +======================= ====================== ======================= + 1 572.98 9 154.30 * + 10 72.55 922.07 + 25 38.88 412.34 + 50 27.95 241.35 + 75 27.54 183.63 + 100 27.86 165.68 + 150 28.14 134.05 + 200 27.77 134.64 +======================= ====================== ======================= + + +PT2 correction +-------------- + +We present the parallel speedup curve, and the wall-clock time in seconds +required to compute the PT2 correction for two wave functions measured on +Olympe and Irene. + + +.. figure:: _static/speedup_pt2.png + :alt: Parallel speedup of the PT2 computation of the ground state. + + Parallel speedup of the PT2 computation of the ground state measured + on Olympe and Irene. + + +Olympe +^^^^^^ + +======================= ====================== ======================= +Number of 36-core Nodes Ground state (9.3M) Excited state (9.3M) +======================= ====================== ======================= + 1 7 883.74 9 829.19 + 5 1 629.06 2 022.36 + 10 832.89 1 029.91 + 20 440.76 537.37 + 30 303.31 378.69 + 40 246.12 296.31 + 50 201.84 241.55 +======================= ====================== ======================= + + +Irene +^^^^^ + +======================= ====================== ======================= ====================== ======================= +Number of 48-core Nodes Ground state (9.3M) Excited state (9.3M) Ground state (42.9M) Excited state (42.9M) +======================= ====================== ======================= ====================== ======================= + 1 4 935.81 6 152.29 24 586.62 37 440.59 + 10 525.95 652.23 2 458.66 3 086.19 + 25 237.47 286.06 1 041.69 1 295.43 + 50 144.39 174.12 588.35 724.25 + 75 109.13 129.17 446.74 537.59 + 100 100.75 103.43 367.21 450.32 + 150 82.04 91.77 298.63 358.25 + 200 75.62 85.25 268.96 312.23 +======================= ====================== ======================= ====================== ======================= -The determinant selection, MR-PT2 and diagonalization are parallelized with distributed parallelism. -Benchmarks for the [NH2-CH-NH2]+ molecule are presented with up to 200 nodes (9600 cores) on the -`Irene` supercomputer. .. _Irene: http://www-hpc.cea.fr/en/complexe/tgcc-Irene.htm - - -.. include:: work.rst +.. _Olympe: https://www.calmip.univ-toulouse.fr/spip.php?article582&lang=fr diff --git a/docs/source/perf/davidson.dat b/docs/source/perf/davidson.dat new file mode 100644 index 00000000..dae57ad4 --- /dev/null +++ b/docs/source/perf/davidson.dat @@ -0,0 +1,24 @@ +#Olympe +#Number of 36-core Nodes 9 356 952 determinants 42 959 496 determinants + 1 775.55 11198.70 + 5 169.88 2288.58 + 10 93.22 1213.95 + 20 56.86 626.41 + 30 43.76 445.65 + 40 36.18 350.25 + 50 33.67 295.25 + + +#Irene +#Number of 48-core Nodes 9 356 952 determinants 42 959 496 determinants + 1 572.98 9154.30 + 10 72.55 922.07 + 25 38.88 412.34 + 50 27.95 241.35 + 75 27.54 183.63 + 100 27.86 165.68 + 150 28.14 134.05 + 200 27.77 134.64 + + + diff --git a/docs/source/perf/davidson.plt b/docs/source/perf/davidson.plt new file mode 100644 index 00000000..8af90013 --- /dev/null +++ b/docs/source/perf/davidson.plt @@ -0,0 +1,194 @@ +#!/usr/bin/gnuplot -persist +# +# +# G N U P L O T +# Version 5.2 patchlevel 4 last modified 2018-06-01 +# +# Copyright (C) 1986-1993, 1998, 2004, 2007-2018 +# Thomas Williams, Colin Kelley and many others +# +# gnuplot home: http://www.gnuplot.info +# faq, bugs, etc: type "help FAQ" +# immediate help: type "help" (plot window: hit 'h') +set terminal png +set output '../_static/speedup_davidson.png' +unset clip points +set clip one +unset clip two +#set errorbars front 1.000000 +set border 31 front lt black linewidth 1.000 dashtype solid +set zdata +set ydata +set xdata +set y2data +set x2data +set boxwidth +set style fill empty border +set style rectangle back fc bgnd fillstyle solid 1.00 border lt -1 +set style circle radius graph 0.02 +set style ellipse size graph 0.05, 0.03 angle 0 units xy +set dummy x, y +set format x "% h" +set format y "% h" +set format x2 "% h" +set format y2 "% h" +set format z "% h" +set format cb "% h" +set format r "% h" +#set ttics format "% h" +set timefmt "%d/%m/%y,%H:%M" +set angles radians +set tics back +set grid nopolar +set grid xtics nomxtics ytics nomytics noztics nomztics nortics nomrtics \ + nox2tics nomx2tics noy2tics nomy2tics nocbtics nomcbtics +set grid layerdefault lt 0 linecolor 0 linewidth 0.500 dashtype solid, lt 0 linecolor 0 linewidth 0.500 dashtype solid +unset raxis +#set theta counterclockwise right +set style parallel front lt black linewidth 2.000 dashtype solid +set key title "" center +#set key inside right bottom vertical Right noreverse enhanced autotitle nobox +set key noinvert samplen 4 spacing 1 width 0 height 0 +set key maxcolumns 0 maxrows 0 +set key noopaque +unset label +unset arrow +set style increment default +unset style line +unset style arrow +#set style histogram clustered gap 2 title textcolor lt -1 +unset object +#set style textbox transparent margins 1.0, 1.0 border lt -1 linewidth 1.0 +set offsets 0, 0, 0, 0 +set pointsize 1 +set pointintervalbox 1 +set encoding default +unset polar +unset parametric +unset decimalsign +#unset micro +#unset minussign +set view 60, 30, 1, 1 +#set view azimuth 0 +#set rgbmax 255 +set samples 100, 100 +set isosamples 10, 10 +set surface +unset contour +set cntrlabel format '%8.3g' font '' start 5 interval 20 +set mapping cartesian +set datafile separator whitespace +unset hidden3d +set cntrparam order 4 +set cntrparam linear +#set cntrparam levels auto 5 unsorted +#set cntrparam firstlinetype 0 +set cntrparam points 5 +set size ratio 0 1,1 +set origin 0,0 +set style data points +set style function lines +unset xzeroaxis +unset yzeroaxis +unset zzeroaxis +unset x2zeroaxis +unset y2zeroaxis +set xyplane relative 0.5 +set tics scale 1, 0.5, 1, 1, 1 +set mxtics default +set mytics default +set mztics default +set mx2tics default +set my2tics default +set mcbtics default +#set mrtics default +#set nomttics +set xtics border in scale 1,0.5 mirror norotate autojustify +#set xtics norangelimit autofreq +set ytics border in scale 1,0.5 mirror norotate autojustify +#set ytics norangelimit autofreq +set ztics border in scale 1,0.5 nomirror norotate autojustify +#set ztics norangelimit autofreq +unset x2tics +unset y2tics +set cbtics border in scale 1,0.5 mirror norotate autojustify +set cbtics norangelimit autofreq +set rtics axis in scale 1,0.5 nomirror norotate autojustify +set rtics norangelimit autofreq +#unset ttics +#unset paxis 1 tics +#unset paxis 2 tics +#unset paxis 3 tics +#unset paxis 4 tics +#unset paxis 5 tics +#unset paxis 6 tics +#unset paxis 7 tics +set title "" +set title font "" norotate +set timestamp bottom +set timestamp "" +set timestamp font "" norotate +set trange [ * : * ] noreverse nowriteback +set urange [ * : * ] noreverse nowriteback +set vrange [ * : * ] noreverse nowriteback +set xlabel "Number of 36-core nodes" +set xlabel font "" textcolor lt -1 norotate +set x2label "" +set x2label font "" textcolor lt -1 norotate +set xrange [ * : * ] noreverse writeback +set x2range [ * : * ] noreverse writeback +set ylabel "Speedup" +set ylabel font "" textcolor lt -1 rotate +set y2label "" +set y2label font "" textcolor lt -1 rotate +set yrange [ * : * ] noreverse writeback +set y2range [ * : * ] noreverse writeback +set zlabel "" +set zlabel font "" textcolor lt -1 norotate +set zrange [ * : * ] noreverse writeback +set cblabel "" +set cblabel font "" textcolor lt -1 rotate +set cbrange [ * : * ] noreverse writeback +#set rlabel "" +#set rlabel font "" textcolor lt -1 norotate +set rrange [ * : * ] noreverse writeback +set paxis 1 range [ * : * ] noreverse nowriteback +set paxis 2 range [ * : * ] noreverse nowriteback +set paxis 3 range [ * : * ] noreverse nowriteback +set paxis 4 range [ * : * ] noreverse nowriteback +set paxis 5 range [ * : * ] noreverse nowriteback +set paxis 6 range [ * : * ] noreverse nowriteback +set paxis 7 range [ * : * ] noreverse nowriteback +unset logscale +#unset jitter +set zero 1e-08 +set lmargin -1 +set bmargin -1 +set rmargin -1 +set tmargin -1 +#set locale "en_US.UTF-8" +set pm3d explicit at s +set pm3d scansautomatic +set pm3d interpolate 1,1 flush begin noftriangles noborder corners2color mean +#set pm3d nolighting +set palette positive nops_allcF maxcolors 0 gamma 1.5 color model RGB +set palette rgbformulae 7, 5, 15 +set colorbox default +#set colorbox vertical origin screen 0.9, 0.2 size screen 0.05, 0.6 front noinvert bdefault +set style boxplot candles range 1.50 outliers pt 7 separation 1 labels auto unsorted +set loadpath +set fontpath +set psdir +set fit brief errorvariables nocovariancevariables errorscaling prescale nowrap v5 +GNUTERM = "qt" +x = 0.0 +set key bottom +set yrange [:80] +## Last datafile plotted: "data_davidson" +plot \ + 'davidson.dat' index 0 u 1:(11198.698/$3) title '42 959 496 dets, Olympe' w lp, \ + 'davidson.dat' index 0 u 1:(775.5464/$2) title '9 356 952 dets, Olympe' w lp, \ + 'davidson.dat' index 1 u 1:(9154.30/$3) title '42 959 496 dets, Irene' w lp, \ + 'davidson.dat' index 1 u 1:(572.98/$2) title '9 356 952 dets, Irene' w lp +# x title "Ideal", \ +# EOF diff --git a/docs/source/perf/energy.dat b/docs/source/perf/energy.dat new file mode 100644 index 00000000..2e1790ca --- /dev/null +++ b/docs/source/perf/energy.dat @@ -0,0 +1,26 @@ +# Num of dets Ground state Excited state + 7 -149.489186 -149.207354 + 123 -149.536265 -149.261860 + 3083 -149.685606 -149.404450 + 29409 -149.826151 -149.547275 + 168595 -149.900352 -149.626058 + 1322537 -149.946655 -149.675032 + 8495334 -149.972032 -149.704145 + 9356952 -149.973375 -149.706822 + 42779636 -149.987370 -149.721470 + 186978487 -149.998582 -149.733039 + + +# Num of dets Ground state Excited state + 7 -150.161107 -149.904883 + 123 -150.116958 -149.849465 + 3083 -150.0435 -149.7808 + 29409 -150.0222 -149.7583 + 168595 -150.0199 -149.7545 + 1322537 -150.01789 -149.75255 + 8495334 -150.01597 -149.75087 + 9356952 -150.01589 -149.75066 + 42959496 -150.01675 -149.75188 + 186978487 -150.01751 -149.75290 + + diff --git a/docs/source/perf/energy.plt b/docs/source/perf/energy.plt new file mode 100644 index 00000000..a5379a63 --- /dev/null +++ b/docs/source/perf/energy.plt @@ -0,0 +1,173 @@ +#!/usr/bin/gnuplot -persist +# +# +# G N U P L O T +# Version 5.0 patchlevel 3 last modified 2016-02-21 +# +# Copyright (C) 1986-1993, 1998, 2004, 2007-2016 +# Thomas Williams, Colin Kelley and many others +# +# gnuplot home: http://www.gnuplot.info +# faq, bugs, etc: type "help FAQ" +# immediate help: type "help" (plot window: hit 'h') +set terminal png +set output '../_static/cn3_energy.png' +unset clip points +set clip one +unset clip two +set bar 1.000000 front +set border 31 front lt black linewidth 1.000 dashtype solid +set zdata +set ydata +set xdata +set y2data +set x2data +set boxwidth +set style fill empty border +set style rectangle back fc bgnd fillstyle solid 1.00 border lt -1 +set style circle radius graph 0.02, first 0.00000, 0.00000 +set style ellipse size graph 0.05, 0.03, first 0.00000 angle 0 units xy +set dummy x, y +set format x "% h" +set format y "% h" +set format x2 "% h" +set format y2 "% h" +set format z "% h" +set format cb "% h" +set format r "% h" +set timefmt "%d/%m/%y,%H:%M" +set angles radians +set tics back +set grid nopolar +set grid xtics nomxtics ytics nomytics noztics nomztics \ + nox2tics nomx2tics noy2tics nomy2tics nocbtics nomcbtics +set grid layerdefault lt 0 linewidth 0.500 dashtype solid, lt 0 linewidth 0.500 dashtype solid +set raxis +set style parallel front lt black linewidth 2.000 dashtype solid +set key title "" center +set key inside right top vertical Right noreverse enhanced autotitle nobox +set key noinvert samplen 4 spacing 1 width 0 height 0 +set key maxcolumns 0 maxrows 0 +set key noopaque +unset label +unset arrow +set style increment default +unset style line +unset style arrow +set style histogram clustered gap 2 title textcolor lt -1 +unset object +set style textbox transparent margins 1.0, 1.0 border +unset logscale +set logscale x 10 +set offsets 0, 0, 0, 0 +set pointsize 1 +set pointintervalbox 1 +set encoding default +unset polar +unset parametric +unset decimalsign +set view 60, 30, 1, 1 +set samples 100, 100 +set isosamples 10, 10 +set surface +unset contour +set cntrlabel format '%8.3g' font '' start 5 interval 20 +set mapping cartesian +set datafile separator whitespace +unset hidden3d +set cntrparam order 4 +set cntrparam linear +set cntrparam levels auto 5 +set cntrparam points 5 +set size ratio 0 1,1 +set origin 0,0 +set style data points +set style function lines +unset xzeroaxis +unset yzeroaxis +unset zzeroaxis +unset x2zeroaxis +unset y2zeroaxis +set xyplane relative 0.5 +set tics scale 1, 0.5, 1, 1, 1 +set mxtics default +set mytics default +set mztics default +set mx2tics default +set my2tics default +set mcbtics default +set mrtics default +set xtics border in scale 1,0.5 mirror norotate autojustify +set xtics norangelimit autofreq +set ytics border in scale 1,0.5 mirror norotate autojustify +set ytics norangelimit autofreq +set ztics border in scale 1,0.5 nomirror norotate autojustify +set ztics norangelimit autofreq +unset x2tics +unset y2tics +set cbtics border in scale 1,0.5 mirror norotate autojustify +set cbtics norangelimit autofreq +set rtics axis in scale 1,0.5 nomirror norotate autojustify +set rtics norangelimit autofreq +unset paxis 1 tics +unset paxis 2 tics +unset paxis 3 tics +unset paxis 4 tics +unset paxis 5 tics +unset paxis 6 tics +unset paxis 7 tics +set title "" +set title font "" norotate +set timestamp bottom +set timestamp "" +set timestamp font "" norotate +set rrange [ * : * ] noreverse nowriteback +set trange [ * : * ] noreverse nowriteback +set urange [ * : * ] noreverse nowriteback +set vrange [ * : * ] noreverse nowriteback +set xlabel "Number of determinants" +set xlabel font "" textcolor lt -1 norotate +set x2label "" +set x2label font "" textcolor lt -1 norotate +set xrange [ * : * ] noreverse nowriteback +set x2range [ * : * ] noreverse nowriteback +set ylabel "Energy (au)" +set ylabel font "" textcolor lt -1 rotate by -270 +set y2label "" +set y2label font "" textcolor lt -1 rotate by -270 +set yrange [ * : * ] noreverse nowriteback +set y2range [ * : * ] noreverse nowriteback +set zlabel "" +set zlabel font "" textcolor lt -1 norotate +set zrange [ * : * ] noreverse nowriteback +set cblabel "" +set cblabel font "" textcolor lt -1 rotate by -270 +set cbrange [ * : * ] noreverse nowriteback +set paxis 1 range [ * : * ] noreverse nowriteback +set paxis 2 range [ * : * ] noreverse nowriteback +set paxis 3 range [ * : * ] noreverse nowriteback +set paxis 4 range [ * : * ] noreverse nowriteback +set paxis 5 range [ * : * ] noreverse nowriteback +set paxis 6 range [ * : * ] noreverse nowriteback +set paxis 7 range [ * : * ] noreverse nowriteback +set zero 1e-08 +set lmargin -1 +set bmargin -1 +set rmargin -1 +set tmargin -1 +set pm3d explicit at s +set pm3d scansautomatic +set pm3d interpolate 1,1 flush begin noftriangles noborder corners2color mean +set palette positive nops_allcF maxcolors 0 gamma 1.5 color model RGB +set palette rgbformulae 7, 5, 15 +set colorbox default +set colorbox vertical origin screen 0.9, 0.2, 0 size screen 0.05, 0.6, 0 front bdefault +set style boxplot candles range 1.50 outliers pt 7 separation 1 labels auto unsorted +set loadpath +set fontpath +set psdir +set fit brief errorvariables nocovariancevariables errorscaling prescale nowrap v5 +GNUTERM = "qt" +x = 0.0 +plot 'energy.dat' index 0 u 1:3 w lp title 'Excited state, E_{var}', 'energy.dat' index 0 u 1:2 w lp title 'Ground state, E_{var}', 'energy.dat' index 1 u 1:3 w errorlines title 'Excited state, E_{var}+PT2', 'energy.dat' index 1 u 1:2 w errorlines title 'Ground state, E_{var}+PT2' +# EOF diff --git a/docs/source/perf/pt2.dat b/docs/source/perf/pt2.dat new file mode 100644 index 00000000..3995121a --- /dev/null +++ b/docs/source/perf/pt2.dat @@ -0,0 +1,20 @@ +# Olympe (9356952-det wf) +1 7883.7443 9829.1905 +5 1629.0612 2022.3561 +10 832.8855 1029.9124 +20 440.7613 537.3718 +30 303.3079 378.6931 +40 246.1153 296.3148 +50 201.8428 241.5513 + + +# Irene (9356952-det and 42959496-det wf) +1 4935.8080 6152.2916 24586.619 37440.5884 +10 525.9498 652.2256 2458.6619 3086.1905 +25 237.4663 286.0648 1041.6905 1295.4332 +50 144.3955 174.1193 588.3493 724.2491 +75 109.1328 129.1688 446.7379 537.5888 +100 100.7537 103.4340 367.2123 450.3191 +150 82.0399 91.7735 298.6261 358.2490 +200 75.6173 85.2520 268.9631 312.2254 + diff --git a/docs/source/perf/pt2.plt b/docs/source/perf/pt2.plt new file mode 100644 index 00000000..58e75d8c --- /dev/null +++ b/docs/source/perf/pt2.plt @@ -0,0 +1,176 @@ +#!/usr/bin/gnuplot -persist +# +# +# G N U P L O T +# Version 5.0 patchlevel 3 last modified 2016-02-21 +# +# Copyright (C) 1986-1993, 1998, 2004, 2007-2016 +# Thomas Williams, Colin Kelley and many others +# +# gnuplot home: http://www.gnuplot.info +# faq, bugs, etc: type "help FAQ" +# immediate help: type "help" (plot window: hit 'h') +set terminal png +set output '../_static/speedup_pt2.png' +unset clip points +set clip one +unset clip two +set bar 1.000000 front +set border 31 front lt black linewidth 1.000 dashtype solid +set zdata +set ydata +set xdata +set y2data +set x2data +set boxwidth +set style fill empty border +set style rectangle back fc bgnd fillstyle solid 1.00 border lt -1 +set style circle radius graph 0.02, first 0.00000, 0.00000 +set style ellipse size graph 0.05, 0.03, first 0.00000 angle 0 units xy +set dummy x, y +set format x "% h" +set format y "% h" +set format x2 "% h" +set format y2 "% h" +set format z "% h" +set format cb "% h" +set format r "% h" +set timefmt "%d/%m/%y,%H:%M" +set angles radians +set tics back +unset grid +set raxis +set style parallel front lt black linewidth 2.000 dashtype solid +set key title "" center +set key inside right top vertical Right noreverse enhanced autotitle nobox +set key noinvert samplen 4 spacing 1 width 0 height 0 +set key maxcolumns 0 maxrows 0 +set key noopaque +unset label +unset arrow +set style increment default +unset style line +unset style arrow +set style histogram clustered gap 2 title textcolor lt -1 +unset object +set style textbox transparent margins 1.0, 1.0 border +unset logscale +set offsets 0, 0, 0, 0 +set pointsize 1 +set pointintervalbox 1 +set encoding default +unset polar +unset parametric +unset decimalsign +set view 60, 30, 1, 1 +set samples 100, 100 +set isosamples 10, 10 +set surface +unset contour +set cntrlabel format '%8.3g' font '' start 5 interval 20 +set mapping cartesian +set datafile separator whitespace +unset hidden3d +set cntrparam order 4 +set cntrparam linear +set cntrparam levels auto 5 +set cntrparam points 5 +set size ratio 0 1,1 +set origin 0,0 +set style data points +set style function lines +unset xzeroaxis +unset yzeroaxis +unset zzeroaxis +unset x2zeroaxis +unset y2zeroaxis +set xyplane relative 0.5 +set tics scale 1, 0.5, 1, 1, 1 +set mxtics default +set mytics default +set mztics default +set mx2tics default +set my2tics default +set mcbtics default +set mrtics default +set xtics border in scale 1,0.5 mirror norotate autojustify +set xtics norangelimit autofreq +set ytics border in scale 1,0.5 mirror norotate autojustify +set ytics norangelimit autofreq +set ztics border in scale 1,0.5 nomirror norotate autojustify +set ztics norangelimit autofreq +unset x2tics +unset y2tics +set cbtics border in scale 1,0.5 mirror norotate autojustify +set cbtics norangelimit autofreq +set rtics axis in scale 1,0.5 nomirror norotate autojustify +set rtics norangelimit autofreq +unset paxis 1 tics +unset paxis 2 tics +unset paxis 3 tics +unset paxis 4 tics +unset paxis 5 tics +unset paxis 6 tics +unset paxis 7 tics +set title "" +set title font "" norotate +set timestamp bottom +set timestamp "" +set timestamp font "" norotate +set rrange [ * : * ] noreverse nowriteback +set trange [ * : * ] noreverse nowriteback +set urange [ * : * ] noreverse nowriteback +set vrange [ * : * ] noreverse nowriteback +set xlabel "Number of 36-core nodes" +set xlabel font "" textcolor lt -1 norotate +set x2label "" +set x2label font "" textcolor lt -1 norotate +set xrange [ * : * ] noreverse nowriteback +set x2range [ * : * ] noreverse nowriteback +set ylabel "Speedup" +set ylabel font "" textcolor lt -1 rotate by -270 +set y2label "" +set y2label font "" textcolor lt -1 rotate by -270 +set yrange [ * : * ] noreverse nowriteback +set y2range [ * : * ] noreverse nowriteback +set zlabel "" +set zlabel font "" textcolor lt -1 norotate +set zrange [ * : * ] noreverse nowriteback +set cblabel "" +set cblabel font "" textcolor lt -1 rotate by -270 +set cbrange [ * : * ] noreverse nowriteback +set paxis 1 range [ * : * ] noreverse nowriteback +set paxis 2 range [ * : * ] noreverse nowriteback +set paxis 3 range [ * : * ] noreverse nowriteback +set paxis 4 range [ * : * ] noreverse nowriteback +set paxis 5 range [ * : * ] noreverse nowriteback +set paxis 6 range [ * : * ] noreverse nowriteback +set paxis 7 range [ * : * ] noreverse nowriteback +set zero 1e-08 +set lmargin -1 +set bmargin -1 +set rmargin -1 +set tmargin -1 +set pm3d explicit at s +set pm3d scansautomatic +set pm3d interpolate 1,1 flush begin noftriangles noborder corners2color mean +set palette positive nops_allcF maxcolors 0 gamma 1.5 color model RGB +set palette rgbformulae 7, 5, 15 +set colorbox default +set colorbox vertical origin screen 0.9, 0.2, 0 size screen 0.05, 0.6, 0 front bdefault +set style boxplot candles range 1.50 outliers pt 7 separation 1 labels auto unsorted +set loadpath +set fontpath +set psdir +set fit brief errorvariables nocovariancevariables errorscaling prescale nowrap v5 +set key bottom +set grid +GNUTERM = "qt" +x = 0.0 + + + +plot 'pt2.dat' index 0 u 1:(7883.7443/$2) w lp title '9 356 952 dets, Olympe', \ + 'pt2.dat' index 1 u 1:(4935.8080/$2) w lp title '9 356 952 dets, Irene', \ + 'pt2.dat' index 1 u 1:(24586.6194/$4) w lp title '42 959 496 dets, Irene', \ +# EOF