r25518 - /branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/ -- September 01, 2014

Author: tlinnet
Date: Mon Sep  1 23:25:34 2014
New Revision: 25518

URL: http://svn.gna.org/viewcvs/relax?rev=25518&view=rev
Log:
Fix for data-stripping script, looping over the wrong indexes.

task #7824(https://gna.org/task/index.php?7824): Model parameter ERROR 
estimation from Jacobian and Co-variance matrix of dispersion models.

Modified:
    
branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/final_results_mc_strip.bz2
    
branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/final_state_mc_strip_data.py

Modified: 
branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/final_results_mc_strip.bz2
URL: 
http://svn.gna.org/viewcvs/relax/branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/final_results_mc_strip.bz2?rev=25518&r1=25517&r2=25518&view=diff
==============================================================================
--- 
branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/final_results_mc_strip.bz2
        (original)
+++ 
branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/final_results_mc_strip.bz2
        Mon Sep  1 23:25:34 2014
@@ -1,18 +1,17 @@
-BZh91AY&SY³>}¿ÀÿTtDPïÿ÷½Ëü
-¿ïÿú`]ô qíáÄ7· Ó ¡ í=·AªÙ¶Û[8I&ÄÐSô) P=2hm#!ê
6 FÒ&M0Ð11L¦mTÔÉ  
$Ò1Dôõ4òê4Ñê= Ä zã&M&Lbi¦ T "!M

F§´0
4ÓLÛéIÎòY Rª0ðfbgôø$Oµ$©">÷ÊZ^I5hX $F$"$D ÒöÙ    
ADdD"!@X
-°ðXI"¬R@çõúCæzyôØeUT$TUYDrÝå}Çn
-¯mäßEýû²¯¾º7Vuö´Éè¬
,nÁc|\!#÷k73L\´þ?Íã¦ý%á?¬öÍîZ1ÌäÀZÎPcco±O
½qÜ2Ûöb"vÝ¤7omÕJÈL°XÈÙïu5¯wÂµÃ¡Ýðõg¹EOñÂ&-J*J`Dhc¡JÑ´¨bÑ¢UK
-a(â´kkV[iE¶ª ÚÑÃF4Exú8ú¡ý}l$($àâàèæÉîà
a,ùùûÆïß¡ÖÈ^T!JIhvr%$ÈI$ÈY©&B%)A    $
-%&B^´ Y&AiÕ©M
-µnÄÒ,ÖèêSI»mÑµ4«tuVèì)¤±ÑÆ&4Ý¼5a!Â`×dvËmâPM[7bôÓe×]ÃÊÂÒKï¿+7j/²(¢RIp«p
 âI©Ë.8ÓBË©²í®ÍeÝÖíËkìÛ*XÍF÷½/
SR®\©0AX±BG"¤<gVÀHÍvf®·¥&÷â.($¸XrÌ,Vï[>7µf1
pêe³+4Ü8ßMa7Q(m³jc.Rz>ïÖäGªC¥Ü»"Ê]vÖhE¢E0ð 
ÙÛ@ú¿7»ì|_;B|t#Tèù,·OiôI-¢H0$Iâ'@1D28ù0@ÖãÀ)¨B¬Ô$à3Hwd;¼SGöD¹%wlTú/n08a`2Â%°ll!e0»YK0ÂÓ¤èTK"GôílåZc¶Þ4ÔÃS^UXe¶¹(¦³,ç
$Ãá»glD5±Tâ¤sc-$ß°Ì§f^Ëµ-eÝw`Á¥êzÂ$ªJ¶ö6PÉ½£ÕN×^ï4ÙPÎË3uY-Ï§  
ÂêIK4§7rµ4ºaáÂçJ¬µp:
-J¥;ûWeE)U)Nä/}²¤,Ü¤
-D$°xVK´º
-R8uÇø{=>/g=½ú.ÇE¼¬§±=^¯k|øÝksì¤õêJ3J¤;Q©
-RÊ0¥]];Â©©»:±eRîfSÑLE3RU#
-9rK½Ë-/3cËL^bY§¥#!I2*¤*
-JÕîºlöÈ!K®ÇjeTì=ÄFçõ)+M5iTÑif4**'¥BIY¯Õ 
ÀL².BDËL¡4Ð¦$)KaX_VaÇ8Ãå$ßy&Ë£¿ì¦´¦ë¼,®Þ-_^´¶ë0ËÃn\ºUo¿ë[¥¹Óq$ÓgìÐl[6Èeec`Ö²Ýkñ±gºÖ;¥'-©å¦£XRåfzð½Õ<3§³w
±©YÝÂë-³Í¼a%Úë,©w±)
ùIà¿ªÌzó¿k3-´ÉÒud0µ"l¤¬å[ÐÞêgsÑ/^e²HØó$¦yd2Âüf2¤UdéÙºl¾VLt£¨ROZÓtNÏ4WòÊÝ;íjÕjõ¿é1mévöY$au!u:»2àbrX6EÕÖj+8Í)f]dL0m51,`¦A°"R1·ÞøÅHËO"6e8TªyÃ¦ëåÃp·7é&ðªJIEEIBz{%gf¯v½z÷å¶Ì±¶7bÌiáÄ¸1½xh-Ãwa¦^Ýu
-H¹ÔHPiu[cB¯KÆyA:SÄÎºÙúÅZrÐÌÝá³Õ<$Ã
.í¼t$ÂÌ<ït6õñeÝnÏoEL³}Õ:jíâxÃi$Mü+iwk-ðß¶7iÌ¥Ú[®!§1      , 
¤4à8oR@B·Ë¯f"'¥ï bd×'FØ°¼+'¡ñ|ð"´í)&Ìv¶í¶£"'*þöÃË$¼¶Ãu
-s¸©b¤c¯â_`!ÉÆ2gÙû"ûò3£QÝÅ6zÑÔDRÒâDd@ÒÌ'C2ôöÔ°n~¤/¯¢ú@C¤
Ýíq#AÒÄF3ÉØÎ¸k}JCà»úþZ!kâçýÿý!SçæÇÖD@íákI7)çh
-,Pd*H,"*(   ;8:íµ¥µµ¶Ûmmh(¢(¢ ª 
4ùnM³oßæú2ÓM±bÔs!!ôÀHÀHÀB0>¦0¨RT(
'ÙDá$%áIPúà¨Y?þ?Fñ¹é·k^èñûóòü×«_åüzÉÉøþýºÿ;ØÓ²»ÑJ6v31®×(XÊ(þ3_K²ú¼\¯[]Àsöo];½ìiV
Íi6º[0UicØ{»×\À     ï¹ 
i*ÚTÀûætryA5ª¬òÞ¤Tµ©ºT¢:¶8÷¹æ)ámKÆáµÎ,Hw$Ýlº% 
ÎP²;ÚJÂ"^ÏØÑ©QrSkI-Hª|¸2¾ö^hòå@0P¶ó50PKEdF5·¶/Ò¦p2àê#Ó<ðcby|8P¢Z,¼£·s4²F.RÙJÍ
± "/ï®ptÏ¦Î33pÌg]ìO*RÎJXqÌ,MjAGü×³ObâX& 
d|XcqBÖ±Û'mkGëCÀØÑõkp\¹:ïFÀWUoÈRfµKÆIRçe=*<iKi$8}÷»·ºÔ½íMA:£NôCJrm
Ûø81¤5Ì8åµ
-Á 
%Å×@ßjFÁM,Ç´oV;ØOx!Ö\¸SSÂîyõamÀ3Ó¯^O{KB-k7E%£[-çÊI¼kü®öÄ$9üÇ;üýÏUWMVx)4ñ4ûLTèJñ´ê³D¤chájæõ\oî'C¨|®ö_kþ+¾õþãã«ZÊÌíl^^^^^zîÏ5ÔÕVdª´ºb¼â6yV¢à6åPõ9n|Q½ÁÒ7\bU+Ìc
 Ä¤ x@Î.;    ñ9bJÃÜÛáêaËz¯öë¦¤}_ô@ÖÜLÇ1« .fÖ<U¤±hA`v2¤¬Õ'Ö
-ÃÂnDÖ³*t9Ày^>>*ð9ÿrE8P³>}¿
+BZh91AY&SY4   ÀÿTtDPïÿ÷½Ëü
+¿ïÿú`_< pàîîô=ÓÛP;
+ -f»k°$Bdh'¤ÔÑ¦Óõ4ÚPiµ i
+Ð
ESýª42hdC ¦Ma&¢§y&H¡êhõ
PÆ¦Ð=@q&CL       ¦Á14Ó@*HB
+z§¤õ dzÇ?*þ0¤ú¥-$,µUSù        nÊó0²AõüäO¹$©">ïá-   /$´
+¬ÈA
+#D@"B/K!(0$<ID"!@X
+°ðXI"¬R@çõúCåzyCôØ±UU UV@µyßQß§ëßzöäQõù×6êÎ¿<Õ±Í
-Ö,o$
aæiÇó´tÛ¤¢¼'õÌÞ+F9AÊ
cìmÖ)ãîÙÙ£-!'5ÙûnéäÞZP6!2Á`\Rkw¿ñ5¯6.5+
kCÅÕ¤M-?ÇX´b%(¨a)¢%B+FÒ¢UDJV-,)£ÑYm¥ÚRªTVkF¶·
ÑãçãêöH1õ° WwG.H·wlàÁÛÛñÎ^þçÙ¸É&B´!TÔìåZII¶RL
i&JÒ¢ª©2I2ZL·Òd%$èËDnKOÃ>Æe¥ÑÁÄ!¡fZVe¦Èh¹iµï/Æ¯<¨
E.ß4®qÏ,SnbôÓ×]ÃÊÂÒKÞúY¶¢ü!Ê%$
+·*d²çm0d,º).ÐZèpÖX
°²ÚÐ]²Úå¶m,f£wÎ3L±&s$#Ef
+X÷JÍ3¬h,fË3\æ"óÎ«" Á¢Ì É£&  ËÎB¼ânµ¥¨%À/©
¬ÓaÎôÖj%8qLeÒOsÏÏÿ{ÙÍ"FÕ!Úî"G²]Å`¨L<avxÄHw§O/¬x¼Ì&Â]þw¸È;ÒN@Q$
$ñ ¢H|$Ú
 lP0Ú¶\è TÁúÑ#ùq$¯ÅI9°(¯£µíÎ,$la)°@ØØKÊav%[²a
§i%Ð©*D¯ÂÙÊ´ÇtÔÃS^_<Ì4ÞW¹(§wDéfém"F#Ïã¶*©ØËI7à3)ÙrúTµwâÁKÔøB$ªJ¶ìp¡vR^ïTáPÏ3uY-×·#
 ÊîIK4§Wtµ4ºaåÊçj¬µp;
+J¥<}Ë²¢Å*§'ß±d.ÙHI`ò¬*et,¥H»çÿßµðß~c¢ÞVSÝ=Ïyøíkuð¤÷ÔfHxQòÖT
)eIR®¯8U!u9VtbË,¥Ý,Í§ºf¢fF4rôkZf4±eÃ&/3eÌD
+¤©
+BµföºpùH!K®Ç+jeTð
ÄFÏAï
+RmuÝ³v;YtEªË(9ãã|Û¾1`-&XJLHpà
(Ä¥°¬/«0ß8Ãå$Þä.ü;²BÒ]áexUöÊùö¥¸Y^råÓVº
çä¿kqÃ9iÃüôÃ2YXà5¬¶µùà³ãk'¨R§ôËSQÕÏL)OY
Ü¯ucO,ééï·,±©YÛÖ[«yÃRKµÞYRîWbSõ$È=ÿÆÌ{ë~gi´îÈajE.áI/YÊ¸¡»©³ÙO}¹ÃN¢DÝRJk¦C,/Ð8c*EVNÞ
]áÂùY1áGp$¤õ¦Ñ<©/é»xðµjÍ5zß§¤íQmÒìì²I2Âë,BêwveÀÄè°p«¼Õ[<ÖqRÌ0p 
RÁ¸p¥2A$
2!½ñ{ß©bIèG²SÏ*OXbÚùrØunÜbI¸U%$¢¢¤¡E=¾µ·»^ûéÇ±Æ¶ÅÓËq8cuå 
·-½
2öpÆ×mO.2îÓð»ß;Ýôçn±ÃqáVÓj^Ye)æIgm»áûÅZtÐÌÛËºyI*2]á¸ì9I#zÝÐ÷óÊîöÏ²Ö¦Y¾Õ;jíÄ$ó.ÒHò®%ÞÃ~`ÛN¥.ÒÝö©
:Ie!§!ÌA!   
ßäÛ·»"'¥î±0Ä2É/t"Â÷6'õÂ½0"²À¤S2Å¶ìµçlb"o_T=0Ä2É.ûaªI©ÔT1DR
1ÙòMßp!ÉÊ2giÕæDO£ØXÎbÜWuUQòàElPÖÌ'C2îõT°n>·ð¾¾»é7gÁÄÊ
f"1~fvC]îPBR@$$À>_ïøD&ºíp«#Á! Ð!{øBñP)&¥;Z T!Y
+%HO7eµ¶´µµµ¶¶¶ÛmEQEPDTè=Âé³6þ_«-4Ùc£ Ø             
FØÀ¡FFF $é!!êRbFÄ
+öÇ÷àS ¹þu{¢GÕÑ/G¦
ðôtÐûÇEøfe`HÃú]P±0Qüæ¾eéñHAr*»Úíã¿rê<ÙócB°fl[A°m²Ö    
 
+«H££ãÚ¹æOuÉH×UnÂ¦Hvß3a  Ug¤VêEKYZJöGvÇ?×XO;i\5¶ýpbC¹&çeË)R
ÜÒV÷ï¶4jAT\âÒKR*¢0\2¾ö^hòä ÇÊ-M,f¬Æ¶ú¢üÁg>ôThe£ÃÇé¾
Íeät(Í·¡Ê1pØºTn(mmxo|D$]sC|ÖpøûdÌ3:íbyéK9°à#:È,MJACþkÆY'Ô±i@,P³L18!ëXì¶¥zuàkhÁ{58.NGSsÈÖ"ªê£ïúÊLÔ©xÉíô\ðSÚ£Î©I¤¶HõÞ|¼>Vª|Í
 
ÅhÓÁÆ!«o4¹+Ôûx÷ÃN/wÁÛ¾îû=Hù¦¤øÒ½îa=â^h9rßM/¹JÓ§máZ³]:;¹`<îÐ(>
     $Zf ØBA#mÀP¼sÛçâ÷wUUtÖ³ÃI§Ýi   B¿+O,%#³GWªáx¯q (v
à\Ä÷¤ëOLê¼C,EäÄ¯~ýûùë«<×Y¬2ÉUitÅyÈ$kñ
E¼mªê8rm|Q¹ÁÒ7\bU+Êc Ä¤ 
t6$\qGûN°ù8þæ7Uø;í©!F_géöÜaÛó&4çæmHA½2ÅZKS*JÍ2}@¬=ÓpÏRÈª9@ýDàs[ÛíöþîcÿrE8P4
        

Modified: 
branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/final_state_mc_strip_data.py
URL: 
http://svn.gna.org/viewcvs/relax/branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/final_state_mc_strip_data.py?rev=25518&r1=25517&r2=25518&view=diff
==============================================================================
--- 
branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/final_state_mc_strip_data.py
      (original)
+++ 
branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/final_state_mc_strip_data.py
      Mon Sep  1 23:25:34 2014
@@ -56,7 +56,9 @@
     # Generate spin string.
     spin_string = generate_spin_string(spin=cur_spin, mol_name=mol_name, 
res_num=resi, res_name=resn)
 
-    for exp_type, frq, offset, point, ei, mi, oi, di in 
loop_exp_frq_offset_point(return_indices=True):
+    #for exp_type, frq, offset, point, ei, mi, oi, di in 
loop_exp_frq_offset_point(return_indices=True):
+    for exp_type, frq, ei, mi in loop_exp_frq(return_indices=True):
+
         # Generate the param_key.
         #param_key = return_param_key_from_data(exp_type=exp_type, frq=frq, 
offset=offset, point=point)
         param_key = generate_r20_key(exp_type=exp_type, frq=frq)
@@ -102,6 +104,10 @@
         my_dic[spin_id][param_key]['dw_err_sim'] = dw_sim_err
         my_dic[spin_id][param_key]['k_AB_err_sim'] = k_AB_sim_err
 
+        # Check for correct size.
+        if len(r2a_sim_l) != 10000:
+            print asd
+
 
 ## Loop over spin attributes, to delete them.
 for cur_spin, mol_name, resi, resn, spin_id in spin_loop(full_info=True, 
return_id=True, skip_desel=True):
r25518 - /branches/est_par_error/test_suite/shared_data/dispersion/estimate_par_err/tsmfk01/

Header

Content

Related Messages
