From 61747745f44aa28e361ebebd42ab08623db0fd09 Mon Sep 17 00:00:00 2001 From: Manuel Zingl Date: Fri, 21 Aug 2015 12:24:54 +0200 Subject: [PATCH] New test added for build_sigma_from_txt.py --- python/build_sigma_from_txt.py | 5 +++++ test/CMakeLists.txt | 1 + test/sigma_from_file.output.h5 | Bin 0 -> 51312 bytes test/sigma_from_file.py | 30 ++++++++++++++++++++++++++++++ 4 files changed, 36 insertions(+) create mode 100644 test/sigma_from_file.output.h5 create mode 100644 test/sigma_from_file.py diff --git a/python/build_sigma_from_txt.py b/python/build_sigma_from_txt.py index 387d810b..5ce91a55 100644 --- a/python/build_sigma_from_txt.py +++ b/python/build_sigma_from_txt.py @@ -1,3 +1,7 @@ +import numpy +import string +from pytriqs.gf.local import * + def read_fortran_file (filename): """ Returns a generator that yields all numbers in the Fortran file as float, one by one""" import os.path @@ -17,6 +21,7 @@ def constr_Sigma_real_axis(self, filename, hdf=True, hdf_dataset='SigmaReFreq',n if (len(ol)==1): # if blocks are of size one Fname = filename+'_'+bl+'.dat' else: + print 'TEST' Fname = filename+'_'+bl+'/'+str(ol[0])+'_'+str(ol[0])+'.dat' R = read_fortran_file(Fname) diff --git a/test/CMakeLists.txt b/test/CMakeLists.txt index 62fa3f94..0c2c2d11 100644 --- a/test/CMakeLists.txt +++ b/test/CMakeLists.txt @@ -7,3 +7,4 @@ triqs_add_test_hdf(hk_convert " -p 1.e-6" ) triqs_add_test_hdf(sumkdft_basic " -d 1.e-6" ) triqs_add_test_hdf(srvo3_Gloc " -d 1.e-6" ) triqs_add_test_hdf(srvo3_transp " -d 1.e-6" ) +triqs_add_test_hdf(sigma_from_file " -d 1.e-6" ) diff --git a/test/sigma_from_file.output.h5 b/test/sigma_from_file.output.h5 new file mode 100644 index 0000000000000000000000000000000000000000..47dee8e441915cc84633f83f40da6b72ae78e241 GIT binary patch literal 51312 zcmeI530xM{`o{-FB{#s`ye3`K{FK?Ne@hYTye=r-yb5YArM@7bARsS`ikV6-NqWnj zB(I5<>EGS1CaBltMo76_%QW|0bGa%8Dk8{z-*@KW^gf(3!|{%NK(|HESF*K){-IPb>1-u=$qg!u@~O&zpKm9B`vwnGwi?D$mNL#D8KdJu zV=W;;@j*%jj7u{#5G;A)*RP`7QJ69=wZ`w2XXN(A_!Zm%Pe^m#_)SvA4{T=N@DWi# zmhdR6Wvrv>e;vPld-d*WAHy*2yp(YZ{o-*9aTuj@23(|yPUkTXKM#yL!* zt#S*P7fkZ1=3hJyn5sGVGf-JiKyhex{FF7|%>(Azaw1AQ4{WF-w|Bnle#(;h6mKq` z2jDnPRn8BP8n~s63m-LNbYxI$c>Fkr$u=K^*%!h`g@gy&<_4G#U_JpJm-0MNvG{3E zU<7XAybjl^?cL>i@Dr%4C!+J2HPX4=LS;S8F`GH}U7pwAc^1w?b*a5$vs|S5wtFu( zsb8NS0evmDhi$Jpo=#Zeg2O_iLM5rTEs^%zLO#0>?-SZBHZ%ru;v!X%ZAt?d$xrUr zRj#mg;rJ-Gb+)?{`2&wdwCojaPv-1Bd3q7$9?>kdttPcCQaY++tF_xJRH)!^=%%!< z=#p1g_H~eY&*A>ge27 zwszFo|FCC7dEB}xUL+}1nQx0Rv;fCJYJb@Pegc)mB#w~ly?{gQvWCir7D_!tmDNkG zhyJpab&mk(sBD+h6Sw*I{Y)}9_$9+F`eBXwPU_jt(-bVGbQtT z?e=ZXuOS`}n>Taffas_l^UYTE3m7&*GJlnuIytA=9^w3`^St~(T#4zJ|7Akw`x8b< z=JkK7Q*%b=OyPLgJnRqfi%V+9d_usT%Y7`8*-!HG+!Xwsa6EeR`Tj6;0=8rc>9%V8~=8VqMzw$l*FX4FDJRKi+ zxZS3DlQdB4VWvZLH3^RbhShn(uUTR1;#o^HQZelQ=g ztk#5Hg+-)Y`i9)CNxOvOVe>HG@#Y_D$NceG5z{-IC?e4lJja}!vr{-8Hec*dId9-B zQ9I^8ZQUdwHLZwvAG@~s#|d;({WQk@lAV+%?;U<6-kAPF^5S)Q)-g z0sB@Qizy;w+D_}8(QTV>JnDR+MqG*Bn140Qz2Vxw7Lmw7`&*B{nl2oV-uzGd>($9` zQA8djRL);NKTS9uy?L*G8Kzwoi%9J&G3FtDTZH4$o1ed6=*5Sp9+Fmzg1by|{YE$* zHc!VZj9}-I+A%+J-rJKmZ+S?&_v!}Wc=YDqoY~>$mf3|QWyz#%$KF^c9FID$J$^7BvZX`%;><$QBeUn) z(=FEs$D=ntWr%P8iGMF7)yAcqsQ=0;;ds<}?fhVVMA51$`(p~p@{6^XEPZpOa6IZf z?|KYo$Ha8Zhs6!Ak<_h_96t4Njm1G<3CE+(Yv%{^meEUQu6?$U6s_IS%#{44a6EeR zdwab+)8}RZ8J%#hVA{!L!tt!zB8MxdA{=l5}8pUbaK=T;dnT_#(0*B?U)aUn=-ev`vbD4fzKuP z4?h%+hr=uWti^tv^8@i%(=k6HtlH)EC-O%T6D-Pl9J% zSvseCxNtmr^VjRe1|LekPjXj$a-;p{!NT#d`Li1Pi`p?iyMEJ83y0h%H{9mFwAaHT z9FN}of-%Y42EA~f%xoMmc5~BrgyUiJ&fm?J96{8M`Pg9#x?a4RM`mqkwf?zJ1`5Z+ z<_{|4nKwUEJLcc1c4l$Px;!%YPGFD7k>dH$n~yK*`pTHtJhJn{_pZ^0PmdGmK~yxwBmJ#s9-=iS36#t6s5=IQx> z6HoE@!hC*+_u9Gb?-BpEtGrpU$pqne^yV#X;ttomcbA+!|K72;I*MP9>CKOAvCAiA z(_NBqsH#`cIPvRMb)GjL!d|qdV}5VvUq5ISbC=9=ZM6Kko->8_m)?Bv>2<+-Tizx4 zM}4Z*9Wq-u9=-XZE#ccEf4xJ}*6bR#HDsP}JZ#>>8}%UY>WzTb=_~Bux5TI39JLcfNqV zX!XXt^@rVyFYU=CbGDy5WF9U4ypGL#X^bCg$9(nIKQPBn$R!cq%saE~#23Q(QRlVi zAI$${OK@^t>s)f^xm$0%GFtrkGnC)+wpXQK#eX^!odx(FZqBlR-bgrRKj~ud~akKv4Cao9F51Tjfp0g+i7PVvk-RCcU zJ@fu=Ua^&wsD{jSQ}G;L^>Z;@{`$%^zDk&ZqK#-$;(nfrPNZ zn}zq6-u%{<4==QF{f(T;h)6s#Ui|k7Y~G~IhrHvR+A%-9U(F8&Y`9HY``4McYKQpm zNz{4m^9AOcO)nTa_TAg0jmM4`hc_4h{g2-K&WF!zNv?F8TzNeubI~U8-$$wQ+WEo! zf4oxmY+aX4mXmjU$8Hw?J)1hO9S`PHF7&*4wSP9*w)Oq6kD7`9{*cX^HReNV$NZqV zK}iD&ZjnCqeOr!LE&h8^b)I*86<`0PV?Og^w}e_fZjtHtYPBioF8=#oz4_anKAhHY z{!Ox}@8xqH3o?c0LpD$M2V}yzq;|~zG}raJsi$s`nm0^A&A%7_y}I7~uv-!Tj&62? z9B9;i+>Z4>2DUisU7ntmS6D9kN6J>X!66z+8g%?$HV67{&3z;utx2eUvPKG zpSy3oP8wfowA$yu0pWP`=6lZX&~Vb7Yh=!eUKuY&FMPe|9}^>gQ+6^@6^Lmc22&ZAILJLV@v*E{;z z%je0pfF1!q&N(L>kKQ~_rwhNAsin)Yu#dd$@V+x-!YH4P4B+?iy_C&1q4#>Dqrx3; zvPT8QIo_Oy_w(U>V2vLplt`wFR2jYJEssFx&;F*qec8iS?dqVt9Dn)mXJ;vIEX%DO zjctGM``MB-8%$-Y1bu=AN>thX)Wd`Kvd{KaZwJS9655~YRQa`*99JzI8{;Cs&R?t- z0J;GsOVM$+PZ@VFI3=y*xU)Lh|ELu$Jq1|^RVtAe&{NQ|Kf0bmwvxBhO25wq&YRwP z3cT|s);-|Oi+)Pf5VJBZN>M=xj;gru_Z&(KFpr7Pi}nah>c?XJ2D@IY{n*y}58^aR zBLhYTj0_kVD2o}$R_-UW3FsG@TDP(NYgfIqE8>aG_PpN8;+7jRJz55|?(6e{(GG__+$qpQ(y2B{+@@9UgCAjtw6XW?z8wmf2qU2%nP2sqERsMee!eb3j=K)GoJa zZ^1uaS1QLV^dJXUKF)^~igiWvQ2ud>;du2Qi7ij}`r{i{tIOrM(#i)sSBA2pZ3Nnn zFS{IYTIGgt94eSmJS7>x=S`Ec-A9qW zID9Xwmu@=fePccICtvT|>bU-pB>K*{J>~ekuJ`v-SnqqbGEeY^c?-pmwMkx(uVy4E zZMVFB^jo>~k_*YCY!$RKeJvaO6&~i7D=~6CwDX9S*Uf8{j>>q!`<@SPht|BCwNBnp zcpgy3H`@qL=M3=15k78sZn5q=uN|+R&6}0qi6(M>*Guh~$GY!4zr1=jZ|1}Sy-+*m zvF|bx)c%ENgJ@$`v-+3Ibp3TGY%%i9s z^Lp#X^CIBYWB*w9oyYO&**x8THb?E4$GY!qo?px6VZP(dKh%zStozRQ%(k%kVt>kc z1GChQdA)Vxv1hGy>>umCYw_65=1rWu08Q;W+Uir1yKh}NclWdDRuRVS+ueWYIJ8FI{_K$Vn`6Sz-&THof^H}$t z?Ui4v&hxIvpg%)*9x}lzv2H!qedlq!dN%K+@jQXrF^_fMd475IY##hsIoD&Dd$ub$0&XpA3f$GqOU@w~Ws_1M4Oy74@}yn5_k zZ{2vFUtT@-k9FU99Iu|u)A7OPs2%fI_npo2YdO5ec$SLon8&*Be9vqPhgbYri~Sop=%~k9l>%MDtfEuXtyz>QgV(898CY(#EBW&o-J7)&3VcmMH`;Oadt+RR1 zndZ^dj(NRx<9QMA>al;U`_AKd^=#hxHv%PDYR5d*eJ{zX*RXl2E6=8>9rJqY#_o+A*)UZamvRzZUz~ zTQ{EXnQg)TvFS0XXwsDCfKDl-O!z9X9jAqZoS^R@!+SlCH9YX-%HyM zv}5xo-g6cOqjt>ets4(!=;`Soa z9sHEGWb%Nz+7urgl4=f%hl=%IUp*#QkY_OKGZoS^R@hqQeb=W`FeOF^y zHa73CF>a_G^Lp#Xvtm}O!~XTwjaT!{vSI&N_npP4)vmi&3j%^N^Q7B~Q2qYR5d*eaEEMI-7?$z%TBF+A*)UZans^wT}Jc_t3qR4g%48 z!DU}B{2N72;>pkp_u8nwC*i$rj}7YW;63Lwv_I9U@@p-rS0cYr|M>L6w<_b#&QarvZKbQ9AKh z_XGSJrI7(614agn3>X TRIQS -> TXT -> TRIQS successful!'