mirror of
https://github.com/fjosw/pyerrors.git
synced 2025-05-14 19:43:41 +02:00
Documentation updated
This commit is contained in:
fjosw
parent
8095130228
commit
a3a1df9027
1 changed files with 202 additions and 192 deletions
|
|
@ -1157,105 +1157,110 @@
|
|||
</span><span id="L-1061"><a href="#L-1061"><span class="linenos">1061</span></a> <span class="n">found</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="L-1062"><a href="#L-1062"><span class="linenos">1062</span></a> <span class="n">files</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="L-1063"><a href="#L-1063"><span class="linenos">1063</span></a> <span class="n">names</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="L-1064"><a href="#L-1064"><span class="linenos">1064</span></a> <span class="k">for</span> <span class="p">(</span><span class="n">dirpath</span><span class="p">,</span> <span class="n">dirnames</span><span class="p">,</span> <span class="n">filenames</span><span class="p">)</span> <span class="ow">in</span> <span class="n">os</span><span class="o">.</span><span class="n">walk</span><span class="p">(</span><span class="n">path</span> <span class="o">+</span> <span class="s2">"/"</span><span class="p">):</span>
|
||||
</span><span id="L-1065"><a href="#L-1065"><span class="linenos">1065</span></a> <span class="n">found</span><span class="o">.</span><span class="n">extend</span><span class="p">(</span><span class="n">filenames</span><span class="p">)</span>
|
||||
</span><span id="L-1066"><a href="#L-1066"><span class="linenos">1066</span></a> <span class="k">break</span>
|
||||
</span><span id="L-1064"><a href="#L-1064"><span class="linenos">1064</span></a>
|
||||
</span><span id="L-1065"><a href="#L-1065"><span class="linenos">1065</span></a> <span class="k">if</span> <span class="s2">"names"</span> <span class="ow">in</span> <span class="n">kwargs</span><span class="p">:</span>
|
||||
</span><span id="L-1066"><a href="#L-1066"><span class="linenos">1066</span></a> <span class="n">names</span> <span class="o">=</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="s2">"names"</span><span class="p">)</span>
|
||||
</span><span id="L-1067"><a href="#L-1067"><span class="linenos">1067</span></a>
|
||||
</span><span id="L-1068"><a href="#L-1068"><span class="linenos">1068</span></a> <span class="k">for</span> <span class="n">f</span> <span class="ow">in</span> <span class="n">found</span><span class="p">:</span>
|
||||
</span><span id="L-1069"><a href="#L-1069"><span class="linenos">1069</span></a> <span class="k">if</span> <span class="n">fnmatch</span><span class="o">.</span><span class="n">fnmatch</span><span class="p">(</span><span class="n">f</span><span class="p">,</span> <span class="n">prefix</span> <span class="o">+</span> <span class="s2">"*.ms5_xsf_"</span> <span class="o">+</span> <span class="n">qc</span> <span class="o">+</span> <span class="s2">".dat"</span><span class="p">):</span>
|
||||
</span><span id="L-1070"><a href="#L-1070"><span class="linenos">1070</span></a> <span class="n">files</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">f</span><span class="p">)</span>
|
||||
</span><span id="L-1071"><a href="#L-1071"><span class="linenos">1071</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="n">sep</span> <span class="o">==</span> <span class="s2">""</span><span class="p">:</span>
|
||||
</span><span id="L-1072"><a href="#L-1072"><span class="linenos">1072</span></a> <span class="n">names</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">prefix</span> <span class="o">+</span> <span class="s2">"|r"</span> <span class="o">+</span> <span class="n">f</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="s2">"."</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="n">sep</span><span class="p">)[</span><span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="L-1073"><a href="#L-1073"><span class="linenos">1073</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-1074"><a href="#L-1074"><span class="linenos">1074</span></a> <span class="n">names</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">prefix</span><span class="p">)</span>
|
||||
</span><span id="L-1075"><a href="#L-1075"><span class="linenos">1075</span></a> <span class="n">files</span> <span class="o">=</span> <span class="nb">sorted</span><span class="p">(</span><span class="n">files</span><span class="p">)</span>
|
||||
</span><span id="L-1076"><a href="#L-1076"><span class="linenos">1076</span></a>
|
||||
</span><span id="L-1077"><a href="#L-1077"><span class="linenos">1077</span></a> <span class="k">if</span> <span class="s2">"names"</span> <span class="ow">in</span> <span class="n">kwargs</span><span class="p">:</span>
|
||||
</span><span id="L-1078"><a href="#L-1078"><span class="linenos">1078</span></a> <span class="n">names</span> <span class="o">=</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="s2">"names"</span><span class="p">)</span>
|
||||
</span><span id="L-1079"><a href="#L-1079"><span class="linenos">1079</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-1080"><a href="#L-1080"><span class="linenos">1080</span></a> <span class="n">names</span> <span class="o">=</span> <span class="nb">sorted</span><span class="p">(</span><span class="n">names</span><span class="p">)</span>
|
||||
</span><span id="L-1081"><a href="#L-1081"><span class="linenos">1081</span></a>
|
||||
</span><span id="L-1082"><a href="#L-1082"><span class="linenos">1082</span></a> <span class="n">cnfgs</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="L-1083"><a href="#L-1083"><span class="linenos">1083</span></a> <span class="n">realsamples</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="L-1084"><a href="#L-1084"><span class="linenos">1084</span></a> <span class="n">imagsamples</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="L-1085"><a href="#L-1085"><span class="linenos">1085</span></a> <span class="n">repnum</span> <span class="o">=</span> <span class="mi">0</span>
|
||||
</span><span id="L-1086"><a href="#L-1086"><span class="linenos">1086</span></a> <span class="k">for</span> <span class="n">file</span> <span class="ow">in</span> <span class="n">files</span><span class="p">:</span>
|
||||
</span><span id="L-1087"><a href="#L-1087"><span class="linenos">1087</span></a> <span class="k">with</span> <span class="nb">open</span><span class="p">(</span><span class="n">path</span> <span class="o">+</span> <span class="s2">"/"</span> <span class="o">+</span> <span class="n">file</span><span class="p">,</span> <span class="s2">"rb"</span><span class="p">)</span> <span class="k">as</span> <span class="n">fp</span><span class="p">:</span>
|
||||
</span><span id="L-1088"><a href="#L-1088"><span class="linenos">1088</span></a>
|
||||
</span><span id="L-1089"><a href="#L-1089"><span class="linenos">1089</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="L-1090"><a href="#L-1090"><span class="linenos">1090</span></a> <span class="n">kappa</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1091"><a href="#L-1091"><span class="linenos">1091</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="L-1092"><a href="#L-1092"><span class="linenos">1092</span></a> <span class="n">csw</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1093"><a href="#L-1093"><span class="linenos">1093</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="L-1094"><a href="#L-1094"><span class="linenos">1094</span></a> <span class="n">dF</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1095"><a href="#L-1095"><span class="linenos">1095</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="L-1096"><a href="#L-1096"><span class="linenos">1096</span></a> <span class="n">zF</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1097"><a href="#L-1097"><span class="linenos">1097</span></a>
|
||||
</span><span id="L-1098"><a href="#L-1098"><span class="linenos">1098</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">4</span><span class="p">)</span>
|
||||
</span><span id="L-1099"><a href="#L-1099"><span class="linenos">1099</span></a> <span class="n">tmax</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'i'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1100"><a href="#L-1100"><span class="linenos">1100</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">4</span><span class="p">)</span>
|
||||
</span><span id="L-1101"><a href="#L-1101"><span class="linenos">1101</span></a> <span class="n">bnd</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'i'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1068"><a href="#L-1068"><span class="linenos">1068</span></a> <span class="k">for</span> <span class="p">(</span><span class="n">dirpath</span><span class="p">,</span> <span class="n">dirnames</span><span class="p">,</span> <span class="n">filenames</span><span class="p">)</span> <span class="ow">in</span> <span class="n">os</span><span class="o">.</span><span class="n">walk</span><span class="p">(</span><span class="n">path</span> <span class="o">+</span> <span class="s2">"/"</span><span class="p">):</span>
|
||||
</span><span id="L-1069"><a href="#L-1069"><span class="linenos">1069</span></a> <span class="n">found</span><span class="o">.</span><span class="n">extend</span><span class="p">(</span><span class="n">filenames</span><span class="p">)</span>
|
||||
</span><span id="L-1070"><a href="#L-1070"><span class="linenos">1070</span></a> <span class="k">break</span>
|
||||
</span><span id="L-1071"><a href="#L-1071"><span class="linenos">1071</span></a>
|
||||
</span><span id="L-1072"><a href="#L-1072"><span class="linenos">1072</span></a> <span class="k">for</span> <span class="n">f</span> <span class="ow">in</span> <span class="n">found</span><span class="p">:</span>
|
||||
</span><span id="L-1073"><a href="#L-1073"><span class="linenos">1073</span></a> <span class="k">if</span> <span class="n">fnmatch</span><span class="o">.</span><span class="n">fnmatch</span><span class="p">(</span><span class="n">f</span><span class="p">,</span> <span class="n">prefix</span> <span class="o">+</span> <span class="s2">"*.ms5_xsf_"</span> <span class="o">+</span> <span class="n">qc</span> <span class="o">+</span> <span class="s2">".dat"</span><span class="p">):</span>
|
||||
</span><span id="L-1074"><a href="#L-1074"><span class="linenos">1074</span></a> <span class="n">files</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">f</span><span class="p">)</span>
|
||||
</span><span id="L-1075"><a href="#L-1075"><span class="linenos">1075</span></a> <span class="k">if</span> <span class="s2">"names"</span> <span class="ow">not</span> <span class="ow">in</span> <span class="n">kwargs</span><span class="p">:</span>
|
||||
</span><span id="L-1076"><a href="#L-1076"><span class="linenos">1076</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="n">sep</span> <span class="o">==</span> <span class="s2">""</span><span class="p">:</span>
|
||||
</span><span id="L-1077"><a href="#L-1077"><span class="linenos">1077</span></a> <span class="n">se</span> <span class="o">=</span> <span class="n">f</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="s2">"."</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1078"><a href="#L-1078"><span class="linenos">1078</span></a> <span class="k">for</span> <span class="n">s</span> <span class="ow">in</span> <span class="n">f</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="s2">"."</span><span class="p">)[</span><span class="mi">1</span><span class="p">:</span><span class="o">-</span><span class="mi">1</span><span class="p">]:</span>
|
||||
</span><span id="L-1079"><a href="#L-1079"><span class="linenos">1079</span></a> <span class="n">se</span> <span class="o">+=</span> <span class="s2">"."</span> <span class="o">+</span> <span class="n">s</span>
|
||||
</span><span id="L-1080"><a href="#L-1080"><span class="linenos">1080</span></a> <span class="n">names</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">se</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="n">sep</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span> <span class="o">+</span> <span class="s2">"|r"</span> <span class="o">+</span> <span class="n">se</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="n">sep</span><span class="p">)[</span><span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="L-1081"><a href="#L-1081"><span class="linenos">1081</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-1082"><a href="#L-1082"><span class="linenos">1082</span></a> <span class="n">names</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">prefix</span><span class="p">)</span>
|
||||
</span><span id="L-1083"><a href="#L-1083"><span class="linenos">1083</span></a>
|
||||
</span><span id="L-1084"><a href="#L-1084"><span class="linenos">1084</span></a> <span class="n">names</span> <span class="o">=</span> <span class="nb">sorted</span><span class="p">(</span><span class="n">names</span><span class="p">)</span>
|
||||
</span><span id="L-1085"><a href="#L-1085"><span class="linenos">1085</span></a> <span class="n">files</span> <span class="o">=</span> <span class="nb">sorted</span><span class="p">(</span><span class="n">files</span><span class="p">)</span>
|
||||
</span><span id="L-1086"><a href="#L-1086"><span class="linenos">1086</span></a>
|
||||
</span><span id="L-1087"><a href="#L-1087"><span class="linenos">1087</span></a> <span class="n">cnfgs</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="L-1088"><a href="#L-1088"><span class="linenos">1088</span></a> <span class="n">realsamples</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="L-1089"><a href="#L-1089"><span class="linenos">1089</span></a> <span class="n">imagsamples</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="L-1090"><a href="#L-1090"><span class="linenos">1090</span></a> <span class="n">repnum</span> <span class="o">=</span> <span class="mi">0</span>
|
||||
</span><span id="L-1091"><a href="#L-1091"><span class="linenos">1091</span></a> <span class="k">for</span> <span class="n">file</span> <span class="ow">in</span> <span class="n">files</span><span class="p">:</span>
|
||||
</span><span id="L-1092"><a href="#L-1092"><span class="linenos">1092</span></a> <span class="k">with</span> <span class="nb">open</span><span class="p">(</span><span class="n">path</span> <span class="o">+</span> <span class="s2">"/"</span> <span class="o">+</span> <span class="n">file</span><span class="p">,</span> <span class="s2">"rb"</span><span class="p">)</span> <span class="k">as</span> <span class="n">fp</span><span class="p">:</span>
|
||||
</span><span id="L-1093"><a href="#L-1093"><span class="linenos">1093</span></a>
|
||||
</span><span id="L-1094"><a href="#L-1094"><span class="linenos">1094</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="L-1095"><a href="#L-1095"><span class="linenos">1095</span></a> <span class="n">kappa</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1096"><a href="#L-1096"><span class="linenos">1096</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="L-1097"><a href="#L-1097"><span class="linenos">1097</span></a> <span class="n">csw</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1098"><a href="#L-1098"><span class="linenos">1098</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="L-1099"><a href="#L-1099"><span class="linenos">1099</span></a> <span class="n">dF</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1100"><a href="#L-1100"><span class="linenos">1100</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="L-1101"><a href="#L-1101"><span class="linenos">1101</span></a> <span class="n">zF</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1102"><a href="#L-1102"><span class="linenos">1102</span></a>
|
||||
</span><span id="L-1103"><a href="#L-1103"><span class="linenos">1103</span></a> <span class="n">placesBI</span> <span class="o">=</span> <span class="p">[</span><span class="s2">"gS"</span><span class="p">,</span> <span class="s2">"gP"</span><span class="p">,</span>
|
||||
</span><span id="L-1104"><a href="#L-1104"><span class="linenos">1104</span></a> <span class="s2">"gA"</span><span class="p">,</span> <span class="s2">"gV"</span><span class="p">,</span>
|
||||
</span><span id="L-1105"><a href="#L-1105"><span class="linenos">1105</span></a> <span class="s2">"gVt"</span><span class="p">,</span> <span class="s2">"lA"</span><span class="p">,</span>
|
||||
</span><span id="L-1106"><a href="#L-1106"><span class="linenos">1106</span></a> <span class="s2">"lV"</span><span class="p">,</span> <span class="s2">"lVt"</span><span class="p">,</span>
|
||||
</span><span id="L-1107"><a href="#L-1107"><span class="linenos">1107</span></a> <span class="s2">"lT"</span><span class="p">,</span> <span class="s2">"lTt"</span><span class="p">]</span>
|
||||
</span><span id="L-1108"><a href="#L-1108"><span class="linenos">1108</span></a> <span class="n">placesBB</span> <span class="o">=</span> <span class="p">[</span><span class="s2">"g1"</span><span class="p">,</span> <span class="s2">"l1"</span><span class="p">]</span>
|
||||
</span><span id="L-1109"><a href="#L-1109"><span class="linenos">1109</span></a>
|
||||
</span><span id="L-1110"><a href="#L-1110"><span class="linenos">1110</span></a> <span class="c1"># the chunks have the following structure:</span>
|
||||
</span><span id="L-1111"><a href="#L-1111"><span class="linenos">1111</span></a> <span class="c1"># confignumber, 10x timedependent complex correlators as doubles, 2x timeindependent complex correlators as doubles</span>
|
||||
</span><span id="L-1112"><a href="#L-1112"><span class="linenos">1112</span></a>
|
||||
</span><span id="L-1113"><a href="#L-1113"><span class="linenos">1113</span></a> <span class="n">chunksize</span> <span class="o">=</span> <span class="mi">4</span> <span class="o">+</span> <span class="p">(</span><span class="mi">8</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="mi">10</span><span class="p">)</span> <span class="o">+</span> <span class="p">(</span><span class="mi">8</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="mi">2</span><span class="p">)</span>
|
||||
</span><span id="L-1114"><a href="#L-1114"><span class="linenos">1114</span></a> <span class="n">packstr</span> <span class="o">=</span> <span class="s1">'=i'</span> <span class="o">+</span> <span class="p">(</span><span class="s1">'d'</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="mi">10</span><span class="p">)</span> <span class="o">+</span> <span class="p">(</span><span class="s1">'d'</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="mi">2</span><span class="p">)</span>
|
||||
</span><span id="L-1115"><a href="#L-1115"><span class="linenos">1115</span></a> <span class="n">cnfgs</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="L-1116"><a href="#L-1116"><span class="linenos">1116</span></a> <span class="n">realsamples</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="L-1117"><a href="#L-1117"><span class="linenos">1117</span></a> <span class="n">imagsamples</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="L-1118"><a href="#L-1118"><span class="linenos">1118</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">tmax</span><span class="p">):</span>
|
||||
</span><span id="L-1119"><a href="#L-1119"><span class="linenos">1119</span></a> <span class="n">realsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="L-1120"><a href="#L-1120"><span class="linenos">1120</span></a> <span class="n">imagsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="L-1121"><a href="#L-1121"><span class="linenos">1121</span></a>
|
||||
</span><span id="L-1122"><a href="#L-1122"><span class="linenos">1122</span></a> <span class="k">while</span> <span class="kc">True</span><span class="p">:</span>
|
||||
</span><span id="L-1123"><a href="#L-1123"><span class="linenos">1123</span></a> <span class="n">cnfgt</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="n">chunksize</span><span class="p">)</span>
|
||||
</span><span id="L-1124"><a href="#L-1124"><span class="linenos">1124</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="n">cnfgt</span><span class="p">:</span>
|
||||
</span><span id="L-1125"><a href="#L-1125"><span class="linenos">1125</span></a> <span class="k">break</span>
|
||||
</span><span id="L-1126"><a href="#L-1126"><span class="linenos">1126</span></a> <span class="n">asascii</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="n">packstr</span><span class="p">,</span> <span class="n">cnfgt</span><span class="p">)</span>
|
||||
</span><span id="L-1127"><a href="#L-1127"><span class="linenos">1127</span></a> <span class="n">cnfg</span> <span class="o">=</span> <span class="n">asascii</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1128"><a href="#L-1128"><span class="linenos">1128</span></a> <span class="n">cnfgs</span><span class="p">[</span><span class="n">repnum</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">cnfg</span><span class="p">)</span>
|
||||
</span><span id="L-1129"><a href="#L-1129"><span class="linenos">1129</span></a>
|
||||
</span><span id="L-1130"><a href="#L-1130"><span class="linenos">1130</span></a> <span class="k">if</span> <span class="n">corr</span> <span class="ow">not</span> <span class="ow">in</span> <span class="n">placesBB</span><span class="p">:</span>
|
||||
</span><span id="L-1131"><a href="#L-1131"><span class="linenos">1131</span></a> <span class="n">tmpcorr</span> <span class="o">=</span> <span class="n">asascii</span><span class="p">[</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="n">placesBI</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">):</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="n">placesBI</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span><span class="p">]</span>
|
||||
</span><span id="L-1132"><a href="#L-1132"><span class="linenos">1132</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-1133"><a href="#L-1133"><span class="linenos">1133</span></a> <span class="n">tmpcorr</span> <span class="o">=</span> <span class="n">asascii</span><span class="p">[</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="nb">len</span><span class="p">(</span><span class="n">placesBI</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">placesBB</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">):</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="nb">len</span><span class="p">(</span><span class="n">placesBI</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">placesBB</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span><span class="p">]</span>
|
||||
</span><span id="L-1103"><a href="#L-1103"><span class="linenos">1103</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">4</span><span class="p">)</span>
|
||||
</span><span id="L-1104"><a href="#L-1104"><span class="linenos">1104</span></a> <span class="n">tmax</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'i'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1105"><a href="#L-1105"><span class="linenos">1105</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">4</span><span class="p">)</span>
|
||||
</span><span id="L-1106"><a href="#L-1106"><span class="linenos">1106</span></a> <span class="n">bnd</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'i'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1107"><a href="#L-1107"><span class="linenos">1107</span></a>
|
||||
</span><span id="L-1108"><a href="#L-1108"><span class="linenos">1108</span></a> <span class="n">placesBI</span> <span class="o">=</span> <span class="p">[</span><span class="s2">"gS"</span><span class="p">,</span> <span class="s2">"gP"</span><span class="p">,</span>
|
||||
</span><span id="L-1109"><a href="#L-1109"><span class="linenos">1109</span></a> <span class="s2">"gA"</span><span class="p">,</span> <span class="s2">"gV"</span><span class="p">,</span>
|
||||
</span><span id="L-1110"><a href="#L-1110"><span class="linenos">1110</span></a> <span class="s2">"gVt"</span><span class="p">,</span> <span class="s2">"lA"</span><span class="p">,</span>
|
||||
</span><span id="L-1111"><a href="#L-1111"><span class="linenos">1111</span></a> <span class="s2">"lV"</span><span class="p">,</span> <span class="s2">"lVt"</span><span class="p">,</span>
|
||||
</span><span id="L-1112"><a href="#L-1112"><span class="linenos">1112</span></a> <span class="s2">"lT"</span><span class="p">,</span> <span class="s2">"lTt"</span><span class="p">]</span>
|
||||
</span><span id="L-1113"><a href="#L-1113"><span class="linenos">1113</span></a> <span class="n">placesBB</span> <span class="o">=</span> <span class="p">[</span><span class="s2">"g1"</span><span class="p">,</span> <span class="s2">"l1"</span><span class="p">]</span>
|
||||
</span><span id="L-1114"><a href="#L-1114"><span class="linenos">1114</span></a>
|
||||
</span><span id="L-1115"><a href="#L-1115"><span class="linenos">1115</span></a> <span class="c1"># the chunks have the following structure:</span>
|
||||
</span><span id="L-1116"><a href="#L-1116"><span class="linenos">1116</span></a> <span class="c1"># confignumber, 10x timedependent complex correlators as doubles, 2x timeindependent complex correlators as doubles</span>
|
||||
</span><span id="L-1117"><a href="#L-1117"><span class="linenos">1117</span></a>
|
||||
</span><span id="L-1118"><a href="#L-1118"><span class="linenos">1118</span></a> <span class="n">chunksize</span> <span class="o">=</span> <span class="mi">4</span> <span class="o">+</span> <span class="p">(</span><span class="mi">8</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="mi">10</span><span class="p">)</span> <span class="o">+</span> <span class="p">(</span><span class="mi">8</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="mi">2</span><span class="p">)</span>
|
||||
</span><span id="L-1119"><a href="#L-1119"><span class="linenos">1119</span></a> <span class="n">packstr</span> <span class="o">=</span> <span class="s1">'=i'</span> <span class="o">+</span> <span class="p">(</span><span class="s1">'d'</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="mi">10</span><span class="p">)</span> <span class="o">+</span> <span class="p">(</span><span class="s1">'d'</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="mi">2</span><span class="p">)</span>
|
||||
</span><span id="L-1120"><a href="#L-1120"><span class="linenos">1120</span></a> <span class="n">cnfgs</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="L-1121"><a href="#L-1121"><span class="linenos">1121</span></a> <span class="n">realsamples</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="L-1122"><a href="#L-1122"><span class="linenos">1122</span></a> <span class="n">imagsamples</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="L-1123"><a href="#L-1123"><span class="linenos">1123</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">tmax</span><span class="p">):</span>
|
||||
</span><span id="L-1124"><a href="#L-1124"><span class="linenos">1124</span></a> <span class="n">realsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="L-1125"><a href="#L-1125"><span class="linenos">1125</span></a> <span class="n">imagsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="L-1126"><a href="#L-1126"><span class="linenos">1126</span></a>
|
||||
</span><span id="L-1127"><a href="#L-1127"><span class="linenos">1127</span></a> <span class="k">while</span> <span class="kc">True</span><span class="p">:</span>
|
||||
</span><span id="L-1128"><a href="#L-1128"><span class="linenos">1128</span></a> <span class="n">cnfgt</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="n">chunksize</span><span class="p">)</span>
|
||||
</span><span id="L-1129"><a href="#L-1129"><span class="linenos">1129</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="n">cnfgt</span><span class="p">:</span>
|
||||
</span><span id="L-1130"><a href="#L-1130"><span class="linenos">1130</span></a> <span class="k">break</span>
|
||||
</span><span id="L-1131"><a href="#L-1131"><span class="linenos">1131</span></a> <span class="n">asascii</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="n">packstr</span><span class="p">,</span> <span class="n">cnfgt</span><span class="p">)</span>
|
||||
</span><span id="L-1132"><a href="#L-1132"><span class="linenos">1132</span></a> <span class="n">cnfg</span> <span class="o">=</span> <span class="n">asascii</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1133"><a href="#L-1133"><span class="linenos">1133</span></a> <span class="n">cnfgs</span><span class="p">[</span><span class="n">repnum</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">cnfg</span><span class="p">)</span>
|
||||
</span><span id="L-1134"><a href="#L-1134"><span class="linenos">1134</span></a>
|
||||
</span><span id="L-1135"><a href="#L-1135"><span class="linenos">1135</span></a> <span class="n">corrres</span> <span class="o">=</span> <span class="p">[[],</span> <span class="p">[]]</span>
|
||||
</span><span id="L-1136"><a href="#L-1136"><span class="linenos">1136</span></a> <span class="k">for</span> <span class="n">i</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)):</span>
|
||||
</span><span id="L-1137"><a href="#L-1137"><span class="linenos">1137</span></a> <span class="n">corrres</span><span class="p">[</span><span class="n">i</span> <span class="o">%</span> <span class="mi">2</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">[</span><span class="n">i</span><span class="p">])</span>
|
||||
</span><span id="L-1138"><a href="#L-1138"><span class="linenos">1138</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">int</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)</span> <span class="o">/</span> <span class="mi">2</span><span class="p">)):</span>
|
||||
</span><span id="L-1139"><a href="#L-1139"><span class="linenos">1139</span></a> <span class="n">realsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">][</span><span class="n">t</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">corrres</span><span class="p">[</span><span class="mi">0</span><span class="p">][</span><span class="n">t</span><span class="p">])</span>
|
||||
</span><span id="L-1140"><a href="#L-1140"><span class="linenos">1140</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">int</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)</span> <span class="o">/</span> <span class="mi">2</span><span class="p">)):</span>
|
||||
</span><span id="L-1141"><a href="#L-1141"><span class="linenos">1141</span></a> <span class="n">imagsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">][</span><span class="n">t</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">corrres</span><span class="p">[</span><span class="mi">1</span><span class="p">][</span><span class="n">t</span><span class="p">])</span>
|
||||
</span><span id="L-1142"><a href="#L-1142"><span class="linenos">1142</span></a> <span class="n">repnum</span> <span class="o">+=</span> <span class="mi">1</span>
|
||||
</span><span id="L-1143"><a href="#L-1143"><span class="linenos">1143</span></a>
|
||||
</span><span id="L-1144"><a href="#L-1144"><span class="linenos">1144</span></a> <span class="n">s</span> <span class="o">=</span> <span class="s2">"Read correlator "</span> <span class="o">+</span> <span class="n">corr</span> <span class="o">+</span> <span class="s2">" from "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">repnum</span><span class="p">)</span> <span class="o">+</span> <span class="s2">" replika with "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">realsamples</span><span class="p">[</span><span class="mi">0</span><span class="p">][</span><span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="L-1145"><a href="#L-1145"><span class="linenos">1145</span></a> <span class="k">for</span> <span class="n">rep</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span> <span class="n">repnum</span><span class="p">):</span>
|
||||
</span><span id="L-1146"><a href="#L-1146"><span class="linenos">1146</span></a> <span class="n">s</span> <span class="o">+=</span> <span class="s2">", "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">realsamples</span><span class="p">[</span><span class="n">rep</span><span class="p">][</span><span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="L-1147"><a href="#L-1147"><span class="linenos">1147</span></a> <span class="n">s</span> <span class="o">+=</span> <span class="s2">" samples"</span>
|
||||
</span><span id="L-1148"><a href="#L-1148"><span class="linenos">1148</span></a> <span class="nb">print</span><span class="p">(</span><span class="n">s</span><span class="p">)</span>
|
||||
</span><span id="L-1149"><a href="#L-1149"><span class="linenos">1149</span></a> <span class="nb">print</span><span class="p">(</span><span class="s2">"Asserted run parameters:</span><span class="se">\n</span><span class="s2"> T:"</span><span class="p">,</span> <span class="n">tmax</span><span class="p">,</span> <span class="s2">"kappa:"</span><span class="p">,</span> <span class="n">kappa</span><span class="p">,</span> <span class="s2">"csw:"</span><span class="p">,</span> <span class="n">csw</span><span class="p">,</span> <span class="s2">"dF:"</span><span class="p">,</span> <span class="n">dF</span><span class="p">,</span> <span class="s2">"zF:"</span><span class="p">,</span> <span class="n">zF</span><span class="p">,</span> <span class="s2">"bnd:"</span><span class="p">,</span> <span class="n">bnd</span><span class="p">)</span>
|
||||
</span><span id="L-1150"><a href="#L-1150"><span class="linenos">1150</span></a>
|
||||
</span><span id="L-1151"><a href="#L-1151"><span class="linenos">1151</span></a> <span class="c1"># we have the data now... but we need to re format the whole thing and put it into Corr objects.</span>
|
||||
</span><span id="L-1152"><a href="#L-1152"><span class="linenos">1152</span></a>
|
||||
</span><span id="L-1153"><a href="#L-1153"><span class="linenos">1153</span></a> <span class="n">compObs</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="L-1154"><a href="#L-1154"><span class="linenos">1154</span></a>
|
||||
</span><span id="L-1155"><a href="#L-1155"><span class="linenos">1155</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">int</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)</span> <span class="o">/</span> <span class="mi">2</span><span class="p">)):</span>
|
||||
</span><span id="L-1156"><a href="#L-1156"><span class="linenos">1156</span></a> <span class="n">compObs</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">CObs</span><span class="p">(</span><span class="n">Obs</span><span class="p">([</span><span class="n">realsamples</span><span class="p">[</span><span class="n">rep</span><span class="p">][</span><span class="n">t</span><span class="p">]</span> <span class="k">for</span> <span class="n">rep</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">repnum</span><span class="p">)],</span> <span class="n">names</span><span class="o">=</span><span class="n">names</span><span class="p">,</span> <span class="n">idl</span><span class="o">=</span><span class="n">cnfgs</span><span class="p">),</span>
|
||||
</span><span id="L-1157"><a href="#L-1157"><span class="linenos">1157</span></a> <span class="n">Obs</span><span class="p">([</span><span class="n">imagsamples</span><span class="p">[</span><span class="n">rep</span><span class="p">][</span><span class="n">t</span><span class="p">]</span> <span class="k">for</span> <span class="n">rep</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">repnum</span><span class="p">)],</span> <span class="n">names</span><span class="o">=</span><span class="n">names</span><span class="p">,</span> <span class="n">idl</span><span class="o">=</span><span class="n">cnfgs</span><span class="p">)))</span>
|
||||
</span><span id="L-1158"><a href="#L-1158"><span class="linenos">1158</span></a>
|
||||
</span><span id="L-1159"><a href="#L-1159"><span class="linenos">1159</span></a> <span class="k">if</span> <span class="nb">len</span><span class="p">(</span><span class="n">compObs</span><span class="p">)</span> <span class="o">==</span> <span class="mi">1</span><span class="p">:</span>
|
||||
</span><span id="L-1160"><a href="#L-1160"><span class="linenos">1160</span></a> <span class="k">return</span> <span class="n">compObs</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1161"><a href="#L-1161"><span class="linenos">1161</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-1162"><a href="#L-1162"><span class="linenos">1162</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">compObs</span><span class="p">)</span>
|
||||
</span><span id="L-1135"><a href="#L-1135"><span class="linenos">1135</span></a> <span class="k">if</span> <span class="n">corr</span> <span class="ow">not</span> <span class="ow">in</span> <span class="n">placesBB</span><span class="p">:</span>
|
||||
</span><span id="L-1136"><a href="#L-1136"><span class="linenos">1136</span></a> <span class="n">tmpcorr</span> <span class="o">=</span> <span class="n">asascii</span><span class="p">[</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="n">placesBI</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">):</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="n">placesBI</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span><span class="p">]</span>
|
||||
</span><span id="L-1137"><a href="#L-1137"><span class="linenos">1137</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-1138"><a href="#L-1138"><span class="linenos">1138</span></a> <span class="n">tmpcorr</span> <span class="o">=</span> <span class="n">asascii</span><span class="p">[</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="nb">len</span><span class="p">(</span><span class="n">placesBI</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">placesBB</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">):</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="nb">len</span><span class="p">(</span><span class="n">placesBI</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">placesBB</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span><span class="p">]</span>
|
||||
</span><span id="L-1139"><a href="#L-1139"><span class="linenos">1139</span></a>
|
||||
</span><span id="L-1140"><a href="#L-1140"><span class="linenos">1140</span></a> <span class="n">corrres</span> <span class="o">=</span> <span class="p">[[],</span> <span class="p">[]]</span>
|
||||
</span><span id="L-1141"><a href="#L-1141"><span class="linenos">1141</span></a> <span class="k">for</span> <span class="n">i</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)):</span>
|
||||
</span><span id="L-1142"><a href="#L-1142"><span class="linenos">1142</span></a> <span class="n">corrres</span><span class="p">[</span><span class="n">i</span> <span class="o">%</span> <span class="mi">2</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">[</span><span class="n">i</span><span class="p">])</span>
|
||||
</span><span id="L-1143"><a href="#L-1143"><span class="linenos">1143</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">int</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)</span> <span class="o">/</span> <span class="mi">2</span><span class="p">)):</span>
|
||||
</span><span id="L-1144"><a href="#L-1144"><span class="linenos">1144</span></a> <span class="n">realsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">][</span><span class="n">t</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">corrres</span><span class="p">[</span><span class="mi">0</span><span class="p">][</span><span class="n">t</span><span class="p">])</span>
|
||||
</span><span id="L-1145"><a href="#L-1145"><span class="linenos">1145</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">int</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)</span> <span class="o">/</span> <span class="mi">2</span><span class="p">)):</span>
|
||||
</span><span id="L-1146"><a href="#L-1146"><span class="linenos">1146</span></a> <span class="n">imagsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">][</span><span class="n">t</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">corrres</span><span class="p">[</span><span class="mi">1</span><span class="p">][</span><span class="n">t</span><span class="p">])</span>
|
||||
</span><span id="L-1147"><a href="#L-1147"><span class="linenos">1147</span></a> <span class="n">repnum</span> <span class="o">+=</span> <span class="mi">1</span>
|
||||
</span><span id="L-1148"><a href="#L-1148"><span class="linenos">1148</span></a>
|
||||
</span><span id="L-1149"><a href="#L-1149"><span class="linenos">1149</span></a> <span class="n">s</span> <span class="o">=</span> <span class="s2">"Read correlator "</span> <span class="o">+</span> <span class="n">corr</span> <span class="o">+</span> <span class="s2">" from "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">repnum</span><span class="p">)</span> <span class="o">+</span> <span class="s2">" replika with "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">realsamples</span><span class="p">[</span><span class="mi">0</span><span class="p">][</span><span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="L-1150"><a href="#L-1150"><span class="linenos">1150</span></a> <span class="k">for</span> <span class="n">rep</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span> <span class="n">repnum</span><span class="p">):</span>
|
||||
</span><span id="L-1151"><a href="#L-1151"><span class="linenos">1151</span></a> <span class="n">s</span> <span class="o">+=</span> <span class="s2">", "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">realsamples</span><span class="p">[</span><span class="n">rep</span><span class="p">][</span><span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="L-1152"><a href="#L-1152"><span class="linenos">1152</span></a> <span class="n">s</span> <span class="o">+=</span> <span class="s2">" samples"</span>
|
||||
</span><span id="L-1153"><a href="#L-1153"><span class="linenos">1153</span></a> <span class="nb">print</span><span class="p">(</span><span class="n">s</span><span class="p">)</span>
|
||||
</span><span id="L-1154"><a href="#L-1154"><span class="linenos">1154</span></a> <span class="nb">print</span><span class="p">(</span><span class="s2">"Asserted run parameters:</span><span class="se">\n</span><span class="s2"> T:"</span><span class="p">,</span> <span class="n">tmax</span><span class="p">,</span> <span class="s2">"kappa:"</span><span class="p">,</span> <span class="n">kappa</span><span class="p">,</span> <span class="s2">"csw:"</span><span class="p">,</span> <span class="n">csw</span><span class="p">,</span> <span class="s2">"dF:"</span><span class="p">,</span> <span class="n">dF</span><span class="p">,</span> <span class="s2">"zF:"</span><span class="p">,</span> <span class="n">zF</span><span class="p">,</span> <span class="s2">"bnd:"</span><span class="p">,</span> <span class="n">bnd</span><span class="p">)</span>
|
||||
</span><span id="L-1155"><a href="#L-1155"><span class="linenos">1155</span></a>
|
||||
</span><span id="L-1156"><a href="#L-1156"><span class="linenos">1156</span></a> <span class="c1"># we have the data now... but we need to re format the whole thing and put it into Corr objects.</span>
|
||||
</span><span id="L-1157"><a href="#L-1157"><span class="linenos">1157</span></a>
|
||||
</span><span id="L-1158"><a href="#L-1158"><span class="linenos">1158</span></a> <span class="n">compObs</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="L-1159"><a href="#L-1159"><span class="linenos">1159</span></a>
|
||||
</span><span id="L-1160"><a href="#L-1160"><span class="linenos">1160</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">int</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)</span> <span class="o">/</span> <span class="mi">2</span><span class="p">)):</span>
|
||||
</span><span id="L-1161"><a href="#L-1161"><span class="linenos">1161</span></a> <span class="n">compObs</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">CObs</span><span class="p">(</span><span class="n">Obs</span><span class="p">([</span><span class="n">realsamples</span><span class="p">[</span><span class="n">rep</span><span class="p">][</span><span class="n">t</span><span class="p">]</span> <span class="k">for</span> <span class="n">rep</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">repnum</span><span class="p">)],</span> <span class="n">names</span><span class="o">=</span><span class="n">names</span><span class="p">,</span> <span class="n">idl</span><span class="o">=</span><span class="n">cnfgs</span><span class="p">),</span>
|
||||
</span><span id="L-1162"><a href="#L-1162"><span class="linenos">1162</span></a> <span class="n">Obs</span><span class="p">([</span><span class="n">imagsamples</span><span class="p">[</span><span class="n">rep</span><span class="p">][</span><span class="n">t</span><span class="p">]</span> <span class="k">for</span> <span class="n">rep</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">repnum</span><span class="p">)],</span> <span class="n">names</span><span class="o">=</span><span class="n">names</span><span class="p">,</span> <span class="n">idl</span><span class="o">=</span><span class="n">cnfgs</span><span class="p">)))</span>
|
||||
</span><span id="L-1163"><a href="#L-1163"><span class="linenos">1163</span></a>
|
||||
</span><span id="L-1164"><a href="#L-1164"><span class="linenos">1164</span></a> <span class="k">if</span> <span class="nb">len</span><span class="p">(</span><span class="n">compObs</span><span class="p">)</span> <span class="o">==</span> <span class="mi">1</span><span class="p">:</span>
|
||||
</span><span id="L-1165"><a href="#L-1165"><span class="linenos">1165</span></a> <span class="k">return</span> <span class="n">compObs</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="L-1166"><a href="#L-1166"><span class="linenos">1166</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="L-1167"><a href="#L-1167"><span class="linenos">1167</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">compObs</span><span class="p">)</span>
|
||||
</span></pre></div>
|
||||
|
||||
|
||||
|
|
@ -2349,105 +2354,110 @@ projection to the topological charge sector defined by target</li>
|
|||
</span><span id="read_ms5_xsf-1062"><a href="#read_ms5_xsf-1062"><span class="linenos">1062</span></a> <span class="n">found</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="read_ms5_xsf-1063"><a href="#read_ms5_xsf-1063"><span class="linenos">1063</span></a> <span class="n">files</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="read_ms5_xsf-1064"><a href="#read_ms5_xsf-1064"><span class="linenos">1064</span></a> <span class="n">names</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="read_ms5_xsf-1065"><a href="#read_ms5_xsf-1065"><span class="linenos">1065</span></a> <span class="k">for</span> <span class="p">(</span><span class="n">dirpath</span><span class="p">,</span> <span class="n">dirnames</span><span class="p">,</span> <span class="n">filenames</span><span class="p">)</span> <span class="ow">in</span> <span class="n">os</span><span class="o">.</span><span class="n">walk</span><span class="p">(</span><span class="n">path</span> <span class="o">+</span> <span class="s2">"/"</span><span class="p">):</span>
|
||||
</span><span id="read_ms5_xsf-1066"><a href="#read_ms5_xsf-1066"><span class="linenos">1066</span></a> <span class="n">found</span><span class="o">.</span><span class="n">extend</span><span class="p">(</span><span class="n">filenames</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1067"><a href="#read_ms5_xsf-1067"><span class="linenos">1067</span></a> <span class="k">break</span>
|
||||
</span><span id="read_ms5_xsf-1065"><a href="#read_ms5_xsf-1065"><span class="linenos">1065</span></a>
|
||||
</span><span id="read_ms5_xsf-1066"><a href="#read_ms5_xsf-1066"><span class="linenos">1066</span></a> <span class="k">if</span> <span class="s2">"names"</span> <span class="ow">in</span> <span class="n">kwargs</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1067"><a href="#read_ms5_xsf-1067"><span class="linenos">1067</span></a> <span class="n">names</span> <span class="o">=</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="s2">"names"</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1068"><a href="#read_ms5_xsf-1068"><span class="linenos">1068</span></a>
|
||||
</span><span id="read_ms5_xsf-1069"><a href="#read_ms5_xsf-1069"><span class="linenos">1069</span></a> <span class="k">for</span> <span class="n">f</span> <span class="ow">in</span> <span class="n">found</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1070"><a href="#read_ms5_xsf-1070"><span class="linenos">1070</span></a> <span class="k">if</span> <span class="n">fnmatch</span><span class="o">.</span><span class="n">fnmatch</span><span class="p">(</span><span class="n">f</span><span class="p">,</span> <span class="n">prefix</span> <span class="o">+</span> <span class="s2">"*.ms5_xsf_"</span> <span class="o">+</span> <span class="n">qc</span> <span class="o">+</span> <span class="s2">".dat"</span><span class="p">):</span>
|
||||
</span><span id="read_ms5_xsf-1071"><a href="#read_ms5_xsf-1071"><span class="linenos">1071</span></a> <span class="n">files</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">f</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1072"><a href="#read_ms5_xsf-1072"><span class="linenos">1072</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="n">sep</span> <span class="o">==</span> <span class="s2">""</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1073"><a href="#read_ms5_xsf-1073"><span class="linenos">1073</span></a> <span class="n">names</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">prefix</span> <span class="o">+</span> <span class="s2">"|r"</span> <span class="o">+</span> <span class="n">f</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="s2">"."</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="n">sep</span><span class="p">)[</span><span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="read_ms5_xsf-1074"><a href="#read_ms5_xsf-1074"><span class="linenos">1074</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1075"><a href="#read_ms5_xsf-1075"><span class="linenos">1075</span></a> <span class="n">names</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">prefix</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1076"><a href="#read_ms5_xsf-1076"><span class="linenos">1076</span></a> <span class="n">files</span> <span class="o">=</span> <span class="nb">sorted</span><span class="p">(</span><span class="n">files</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1077"><a href="#read_ms5_xsf-1077"><span class="linenos">1077</span></a>
|
||||
</span><span id="read_ms5_xsf-1078"><a href="#read_ms5_xsf-1078"><span class="linenos">1078</span></a> <span class="k">if</span> <span class="s2">"names"</span> <span class="ow">in</span> <span class="n">kwargs</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1079"><a href="#read_ms5_xsf-1079"><span class="linenos">1079</span></a> <span class="n">names</span> <span class="o">=</span> <span class="n">kwargs</span><span class="o">.</span><span class="n">get</span><span class="p">(</span><span class="s2">"names"</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1080"><a href="#read_ms5_xsf-1080"><span class="linenos">1080</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1081"><a href="#read_ms5_xsf-1081"><span class="linenos">1081</span></a> <span class="n">names</span> <span class="o">=</span> <span class="nb">sorted</span><span class="p">(</span><span class="n">names</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1082"><a href="#read_ms5_xsf-1082"><span class="linenos">1082</span></a>
|
||||
</span><span id="read_ms5_xsf-1083"><a href="#read_ms5_xsf-1083"><span class="linenos">1083</span></a> <span class="n">cnfgs</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="read_ms5_xsf-1084"><a href="#read_ms5_xsf-1084"><span class="linenos">1084</span></a> <span class="n">realsamples</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="read_ms5_xsf-1085"><a href="#read_ms5_xsf-1085"><span class="linenos">1085</span></a> <span class="n">imagsamples</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="read_ms5_xsf-1086"><a href="#read_ms5_xsf-1086"><span class="linenos">1086</span></a> <span class="n">repnum</span> <span class="o">=</span> <span class="mi">0</span>
|
||||
</span><span id="read_ms5_xsf-1087"><a href="#read_ms5_xsf-1087"><span class="linenos">1087</span></a> <span class="k">for</span> <span class="n">file</span> <span class="ow">in</span> <span class="n">files</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1088"><a href="#read_ms5_xsf-1088"><span class="linenos">1088</span></a> <span class="k">with</span> <span class="nb">open</span><span class="p">(</span><span class="n">path</span> <span class="o">+</span> <span class="s2">"/"</span> <span class="o">+</span> <span class="n">file</span><span class="p">,</span> <span class="s2">"rb"</span><span class="p">)</span> <span class="k">as</span> <span class="n">fp</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1089"><a href="#read_ms5_xsf-1089"><span class="linenos">1089</span></a>
|
||||
</span><span id="read_ms5_xsf-1090"><a href="#read_ms5_xsf-1090"><span class="linenos">1090</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1091"><a href="#read_ms5_xsf-1091"><span class="linenos">1091</span></a> <span class="n">kappa</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1092"><a href="#read_ms5_xsf-1092"><span class="linenos">1092</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1093"><a href="#read_ms5_xsf-1093"><span class="linenos">1093</span></a> <span class="n">csw</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1094"><a href="#read_ms5_xsf-1094"><span class="linenos">1094</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1095"><a href="#read_ms5_xsf-1095"><span class="linenos">1095</span></a> <span class="n">dF</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1096"><a href="#read_ms5_xsf-1096"><span class="linenos">1096</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1097"><a href="#read_ms5_xsf-1097"><span class="linenos">1097</span></a> <span class="n">zF</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1098"><a href="#read_ms5_xsf-1098"><span class="linenos">1098</span></a>
|
||||
</span><span id="read_ms5_xsf-1099"><a href="#read_ms5_xsf-1099"><span class="linenos">1099</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">4</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1100"><a href="#read_ms5_xsf-1100"><span class="linenos">1100</span></a> <span class="n">tmax</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'i'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1101"><a href="#read_ms5_xsf-1101"><span class="linenos">1101</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">4</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1102"><a href="#read_ms5_xsf-1102"><span class="linenos">1102</span></a> <span class="n">bnd</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'i'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1069"><a href="#read_ms5_xsf-1069"><span class="linenos">1069</span></a> <span class="k">for</span> <span class="p">(</span><span class="n">dirpath</span><span class="p">,</span> <span class="n">dirnames</span><span class="p">,</span> <span class="n">filenames</span><span class="p">)</span> <span class="ow">in</span> <span class="n">os</span><span class="o">.</span><span class="n">walk</span><span class="p">(</span><span class="n">path</span> <span class="o">+</span> <span class="s2">"/"</span><span class="p">):</span>
|
||||
</span><span id="read_ms5_xsf-1070"><a href="#read_ms5_xsf-1070"><span class="linenos">1070</span></a> <span class="n">found</span><span class="o">.</span><span class="n">extend</span><span class="p">(</span><span class="n">filenames</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1071"><a href="#read_ms5_xsf-1071"><span class="linenos">1071</span></a> <span class="k">break</span>
|
||||
</span><span id="read_ms5_xsf-1072"><a href="#read_ms5_xsf-1072"><span class="linenos">1072</span></a>
|
||||
</span><span id="read_ms5_xsf-1073"><a href="#read_ms5_xsf-1073"><span class="linenos">1073</span></a> <span class="k">for</span> <span class="n">f</span> <span class="ow">in</span> <span class="n">found</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1074"><a href="#read_ms5_xsf-1074"><span class="linenos">1074</span></a> <span class="k">if</span> <span class="n">fnmatch</span><span class="o">.</span><span class="n">fnmatch</span><span class="p">(</span><span class="n">f</span><span class="p">,</span> <span class="n">prefix</span> <span class="o">+</span> <span class="s2">"*.ms5_xsf_"</span> <span class="o">+</span> <span class="n">qc</span> <span class="o">+</span> <span class="s2">".dat"</span><span class="p">):</span>
|
||||
</span><span id="read_ms5_xsf-1075"><a href="#read_ms5_xsf-1075"><span class="linenos">1075</span></a> <span class="n">files</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">f</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1076"><a href="#read_ms5_xsf-1076"><span class="linenos">1076</span></a> <span class="k">if</span> <span class="s2">"names"</span> <span class="ow">not</span> <span class="ow">in</span> <span class="n">kwargs</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1077"><a href="#read_ms5_xsf-1077"><span class="linenos">1077</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="n">sep</span> <span class="o">==</span> <span class="s2">""</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1078"><a href="#read_ms5_xsf-1078"><span class="linenos">1078</span></a> <span class="n">se</span> <span class="o">=</span> <span class="n">f</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="s2">"."</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1079"><a href="#read_ms5_xsf-1079"><span class="linenos">1079</span></a> <span class="k">for</span> <span class="n">s</span> <span class="ow">in</span> <span class="n">f</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="s2">"."</span><span class="p">)[</span><span class="mi">1</span><span class="p">:</span><span class="o">-</span><span class="mi">1</span><span class="p">]:</span>
|
||||
</span><span id="read_ms5_xsf-1080"><a href="#read_ms5_xsf-1080"><span class="linenos">1080</span></a> <span class="n">se</span> <span class="o">+=</span> <span class="s2">"."</span> <span class="o">+</span> <span class="n">s</span>
|
||||
</span><span id="read_ms5_xsf-1081"><a href="#read_ms5_xsf-1081"><span class="linenos">1081</span></a> <span class="n">names</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">se</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="n">sep</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span> <span class="o">+</span> <span class="s2">"|r"</span> <span class="o">+</span> <span class="n">se</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="n">sep</span><span class="p">)[</span><span class="mi">1</span><span class="p">])</span>
|
||||
</span><span id="read_ms5_xsf-1082"><a href="#read_ms5_xsf-1082"><span class="linenos">1082</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1083"><a href="#read_ms5_xsf-1083"><span class="linenos">1083</span></a> <span class="n">names</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">prefix</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1084"><a href="#read_ms5_xsf-1084"><span class="linenos">1084</span></a>
|
||||
</span><span id="read_ms5_xsf-1085"><a href="#read_ms5_xsf-1085"><span class="linenos">1085</span></a> <span class="n">names</span> <span class="o">=</span> <span class="nb">sorted</span><span class="p">(</span><span class="n">names</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1086"><a href="#read_ms5_xsf-1086"><span class="linenos">1086</span></a> <span class="n">files</span> <span class="o">=</span> <span class="nb">sorted</span><span class="p">(</span><span class="n">files</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1087"><a href="#read_ms5_xsf-1087"><span class="linenos">1087</span></a>
|
||||
</span><span id="read_ms5_xsf-1088"><a href="#read_ms5_xsf-1088"><span class="linenos">1088</span></a> <span class="n">cnfgs</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="read_ms5_xsf-1089"><a href="#read_ms5_xsf-1089"><span class="linenos">1089</span></a> <span class="n">realsamples</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="read_ms5_xsf-1090"><a href="#read_ms5_xsf-1090"><span class="linenos">1090</span></a> <span class="n">imagsamples</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="read_ms5_xsf-1091"><a href="#read_ms5_xsf-1091"><span class="linenos">1091</span></a> <span class="n">repnum</span> <span class="o">=</span> <span class="mi">0</span>
|
||||
</span><span id="read_ms5_xsf-1092"><a href="#read_ms5_xsf-1092"><span class="linenos">1092</span></a> <span class="k">for</span> <span class="n">file</span> <span class="ow">in</span> <span class="n">files</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1093"><a href="#read_ms5_xsf-1093"><span class="linenos">1093</span></a> <span class="k">with</span> <span class="nb">open</span><span class="p">(</span><span class="n">path</span> <span class="o">+</span> <span class="s2">"/"</span> <span class="o">+</span> <span class="n">file</span><span class="p">,</span> <span class="s2">"rb"</span><span class="p">)</span> <span class="k">as</span> <span class="n">fp</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1094"><a href="#read_ms5_xsf-1094"><span class="linenos">1094</span></a>
|
||||
</span><span id="read_ms5_xsf-1095"><a href="#read_ms5_xsf-1095"><span class="linenos">1095</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1096"><a href="#read_ms5_xsf-1096"><span class="linenos">1096</span></a> <span class="n">kappa</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1097"><a href="#read_ms5_xsf-1097"><span class="linenos">1097</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1098"><a href="#read_ms5_xsf-1098"><span class="linenos">1098</span></a> <span class="n">csw</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1099"><a href="#read_ms5_xsf-1099"><span class="linenos">1099</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1100"><a href="#read_ms5_xsf-1100"><span class="linenos">1100</span></a> <span class="n">dF</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1101"><a href="#read_ms5_xsf-1101"><span class="linenos">1101</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">8</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1102"><a href="#read_ms5_xsf-1102"><span class="linenos">1102</span></a> <span class="n">zF</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'d'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1103"><a href="#read_ms5_xsf-1103"><span class="linenos">1103</span></a>
|
||||
</span><span id="read_ms5_xsf-1104"><a href="#read_ms5_xsf-1104"><span class="linenos">1104</span></a> <span class="n">placesBI</span> <span class="o">=</span> <span class="p">[</span><span class="s2">"gS"</span><span class="p">,</span> <span class="s2">"gP"</span><span class="p">,</span>
|
||||
</span><span id="read_ms5_xsf-1105"><a href="#read_ms5_xsf-1105"><span class="linenos">1105</span></a> <span class="s2">"gA"</span><span class="p">,</span> <span class="s2">"gV"</span><span class="p">,</span>
|
||||
</span><span id="read_ms5_xsf-1106"><a href="#read_ms5_xsf-1106"><span class="linenos">1106</span></a> <span class="s2">"gVt"</span><span class="p">,</span> <span class="s2">"lA"</span><span class="p">,</span>
|
||||
</span><span id="read_ms5_xsf-1107"><a href="#read_ms5_xsf-1107"><span class="linenos">1107</span></a> <span class="s2">"lV"</span><span class="p">,</span> <span class="s2">"lVt"</span><span class="p">,</span>
|
||||
</span><span id="read_ms5_xsf-1108"><a href="#read_ms5_xsf-1108"><span class="linenos">1108</span></a> <span class="s2">"lT"</span><span class="p">,</span> <span class="s2">"lTt"</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1109"><a href="#read_ms5_xsf-1109"><span class="linenos">1109</span></a> <span class="n">placesBB</span> <span class="o">=</span> <span class="p">[</span><span class="s2">"g1"</span><span class="p">,</span> <span class="s2">"l1"</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1110"><a href="#read_ms5_xsf-1110"><span class="linenos">1110</span></a>
|
||||
</span><span id="read_ms5_xsf-1111"><a href="#read_ms5_xsf-1111"><span class="linenos">1111</span></a> <span class="c1"># the chunks have the following structure:</span>
|
||||
</span><span id="read_ms5_xsf-1112"><a href="#read_ms5_xsf-1112"><span class="linenos">1112</span></a> <span class="c1"># confignumber, 10x timedependent complex correlators as doubles, 2x timeindependent complex correlators as doubles</span>
|
||||
</span><span id="read_ms5_xsf-1113"><a href="#read_ms5_xsf-1113"><span class="linenos">1113</span></a>
|
||||
</span><span id="read_ms5_xsf-1114"><a href="#read_ms5_xsf-1114"><span class="linenos">1114</span></a> <span class="n">chunksize</span> <span class="o">=</span> <span class="mi">4</span> <span class="o">+</span> <span class="p">(</span><span class="mi">8</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="mi">10</span><span class="p">)</span> <span class="o">+</span> <span class="p">(</span><span class="mi">8</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="mi">2</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1115"><a href="#read_ms5_xsf-1115"><span class="linenos">1115</span></a> <span class="n">packstr</span> <span class="o">=</span> <span class="s1">'=i'</span> <span class="o">+</span> <span class="p">(</span><span class="s1">'d'</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="mi">10</span><span class="p">)</span> <span class="o">+</span> <span class="p">(</span><span class="s1">'d'</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="mi">2</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1116"><a href="#read_ms5_xsf-1116"><span class="linenos">1116</span></a> <span class="n">cnfgs</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="read_ms5_xsf-1117"><a href="#read_ms5_xsf-1117"><span class="linenos">1117</span></a> <span class="n">realsamples</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="read_ms5_xsf-1118"><a href="#read_ms5_xsf-1118"><span class="linenos">1118</span></a> <span class="n">imagsamples</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="read_ms5_xsf-1119"><a href="#read_ms5_xsf-1119"><span class="linenos">1119</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">tmax</span><span class="p">):</span>
|
||||
</span><span id="read_ms5_xsf-1120"><a href="#read_ms5_xsf-1120"><span class="linenos">1120</span></a> <span class="n">realsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="read_ms5_xsf-1121"><a href="#read_ms5_xsf-1121"><span class="linenos">1121</span></a> <span class="n">imagsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="read_ms5_xsf-1122"><a href="#read_ms5_xsf-1122"><span class="linenos">1122</span></a>
|
||||
</span><span id="read_ms5_xsf-1123"><a href="#read_ms5_xsf-1123"><span class="linenos">1123</span></a> <span class="k">while</span> <span class="kc">True</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1124"><a href="#read_ms5_xsf-1124"><span class="linenos">1124</span></a> <span class="n">cnfgt</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="n">chunksize</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1125"><a href="#read_ms5_xsf-1125"><span class="linenos">1125</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="n">cnfgt</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1126"><a href="#read_ms5_xsf-1126"><span class="linenos">1126</span></a> <span class="k">break</span>
|
||||
</span><span id="read_ms5_xsf-1127"><a href="#read_ms5_xsf-1127"><span class="linenos">1127</span></a> <span class="n">asascii</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="n">packstr</span><span class="p">,</span> <span class="n">cnfgt</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1128"><a href="#read_ms5_xsf-1128"><span class="linenos">1128</span></a> <span class="n">cnfg</span> <span class="o">=</span> <span class="n">asascii</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1129"><a href="#read_ms5_xsf-1129"><span class="linenos">1129</span></a> <span class="n">cnfgs</span><span class="p">[</span><span class="n">repnum</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">cnfg</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1130"><a href="#read_ms5_xsf-1130"><span class="linenos">1130</span></a>
|
||||
</span><span id="read_ms5_xsf-1131"><a href="#read_ms5_xsf-1131"><span class="linenos">1131</span></a> <span class="k">if</span> <span class="n">corr</span> <span class="ow">not</span> <span class="ow">in</span> <span class="n">placesBB</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1132"><a href="#read_ms5_xsf-1132"><span class="linenos">1132</span></a> <span class="n">tmpcorr</span> <span class="o">=</span> <span class="n">asascii</span><span class="p">[</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="n">placesBI</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">):</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="n">placesBI</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1133"><a href="#read_ms5_xsf-1133"><span class="linenos">1133</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1134"><a href="#read_ms5_xsf-1134"><span class="linenos">1134</span></a> <span class="n">tmpcorr</span> <span class="o">=</span> <span class="n">asascii</span><span class="p">[</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="nb">len</span><span class="p">(</span><span class="n">placesBI</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">placesBB</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">):</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="nb">len</span><span class="p">(</span><span class="n">placesBI</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">placesBB</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1104"><a href="#read_ms5_xsf-1104"><span class="linenos">1104</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">4</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1105"><a href="#read_ms5_xsf-1105"><span class="linenos">1105</span></a> <span class="n">tmax</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'i'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1106"><a href="#read_ms5_xsf-1106"><span class="linenos">1106</span></a> <span class="n">t</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="mi">4</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1107"><a href="#read_ms5_xsf-1107"><span class="linenos">1107</span></a> <span class="n">bnd</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="s1">'i'</span><span class="p">,</span> <span class="n">t</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1108"><a href="#read_ms5_xsf-1108"><span class="linenos">1108</span></a>
|
||||
</span><span id="read_ms5_xsf-1109"><a href="#read_ms5_xsf-1109"><span class="linenos">1109</span></a> <span class="n">placesBI</span> <span class="o">=</span> <span class="p">[</span><span class="s2">"gS"</span><span class="p">,</span> <span class="s2">"gP"</span><span class="p">,</span>
|
||||
</span><span id="read_ms5_xsf-1110"><a href="#read_ms5_xsf-1110"><span class="linenos">1110</span></a> <span class="s2">"gA"</span><span class="p">,</span> <span class="s2">"gV"</span><span class="p">,</span>
|
||||
</span><span id="read_ms5_xsf-1111"><a href="#read_ms5_xsf-1111"><span class="linenos">1111</span></a> <span class="s2">"gVt"</span><span class="p">,</span> <span class="s2">"lA"</span><span class="p">,</span>
|
||||
</span><span id="read_ms5_xsf-1112"><a href="#read_ms5_xsf-1112"><span class="linenos">1112</span></a> <span class="s2">"lV"</span><span class="p">,</span> <span class="s2">"lVt"</span><span class="p">,</span>
|
||||
</span><span id="read_ms5_xsf-1113"><a href="#read_ms5_xsf-1113"><span class="linenos">1113</span></a> <span class="s2">"lT"</span><span class="p">,</span> <span class="s2">"lTt"</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1114"><a href="#read_ms5_xsf-1114"><span class="linenos">1114</span></a> <span class="n">placesBB</span> <span class="o">=</span> <span class="p">[</span><span class="s2">"g1"</span><span class="p">,</span> <span class="s2">"l1"</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1115"><a href="#read_ms5_xsf-1115"><span class="linenos">1115</span></a>
|
||||
</span><span id="read_ms5_xsf-1116"><a href="#read_ms5_xsf-1116"><span class="linenos">1116</span></a> <span class="c1"># the chunks have the following structure:</span>
|
||||
</span><span id="read_ms5_xsf-1117"><a href="#read_ms5_xsf-1117"><span class="linenos">1117</span></a> <span class="c1"># confignumber, 10x timedependent complex correlators as doubles, 2x timeindependent complex correlators as doubles</span>
|
||||
</span><span id="read_ms5_xsf-1118"><a href="#read_ms5_xsf-1118"><span class="linenos">1118</span></a>
|
||||
</span><span id="read_ms5_xsf-1119"><a href="#read_ms5_xsf-1119"><span class="linenos">1119</span></a> <span class="n">chunksize</span> <span class="o">=</span> <span class="mi">4</span> <span class="o">+</span> <span class="p">(</span><span class="mi">8</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="mi">10</span><span class="p">)</span> <span class="o">+</span> <span class="p">(</span><span class="mi">8</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="mi">2</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1120"><a href="#read_ms5_xsf-1120"><span class="linenos">1120</span></a> <span class="n">packstr</span> <span class="o">=</span> <span class="s1">'=i'</span> <span class="o">+</span> <span class="p">(</span><span class="s1">'d'</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="mi">10</span><span class="p">)</span> <span class="o">+</span> <span class="p">(</span><span class="s1">'d'</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="mi">2</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1121"><a href="#read_ms5_xsf-1121"><span class="linenos">1121</span></a> <span class="n">cnfgs</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="read_ms5_xsf-1122"><a href="#read_ms5_xsf-1122"><span class="linenos">1122</span></a> <span class="n">realsamples</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="read_ms5_xsf-1123"><a href="#read_ms5_xsf-1123"><span class="linenos">1123</span></a> <span class="n">imagsamples</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="read_ms5_xsf-1124"><a href="#read_ms5_xsf-1124"><span class="linenos">1124</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">tmax</span><span class="p">):</span>
|
||||
</span><span id="read_ms5_xsf-1125"><a href="#read_ms5_xsf-1125"><span class="linenos">1125</span></a> <span class="n">realsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="read_ms5_xsf-1126"><a href="#read_ms5_xsf-1126"><span class="linenos">1126</span></a> <span class="n">imagsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">([])</span>
|
||||
</span><span id="read_ms5_xsf-1127"><a href="#read_ms5_xsf-1127"><span class="linenos">1127</span></a>
|
||||
</span><span id="read_ms5_xsf-1128"><a href="#read_ms5_xsf-1128"><span class="linenos">1128</span></a> <span class="k">while</span> <span class="kc">True</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1129"><a href="#read_ms5_xsf-1129"><span class="linenos">1129</span></a> <span class="n">cnfgt</span> <span class="o">=</span> <span class="n">fp</span><span class="o">.</span><span class="n">read</span><span class="p">(</span><span class="n">chunksize</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1130"><a href="#read_ms5_xsf-1130"><span class="linenos">1130</span></a> <span class="k">if</span> <span class="ow">not</span> <span class="n">cnfgt</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1131"><a href="#read_ms5_xsf-1131"><span class="linenos">1131</span></a> <span class="k">break</span>
|
||||
</span><span id="read_ms5_xsf-1132"><a href="#read_ms5_xsf-1132"><span class="linenos">1132</span></a> <span class="n">asascii</span> <span class="o">=</span> <span class="n">struct</span><span class="o">.</span><span class="n">unpack</span><span class="p">(</span><span class="n">packstr</span><span class="p">,</span> <span class="n">cnfgt</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1133"><a href="#read_ms5_xsf-1133"><span class="linenos">1133</span></a> <span class="n">cnfg</span> <span class="o">=</span> <span class="n">asascii</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1134"><a href="#read_ms5_xsf-1134"><span class="linenos">1134</span></a> <span class="n">cnfgs</span><span class="p">[</span><span class="n">repnum</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">cnfg</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1135"><a href="#read_ms5_xsf-1135"><span class="linenos">1135</span></a>
|
||||
</span><span id="read_ms5_xsf-1136"><a href="#read_ms5_xsf-1136"><span class="linenos">1136</span></a> <span class="n">corrres</span> <span class="o">=</span> <span class="p">[[],</span> <span class="p">[]]</span>
|
||||
</span><span id="read_ms5_xsf-1137"><a href="#read_ms5_xsf-1137"><span class="linenos">1137</span></a> <span class="k">for</span> <span class="n">i</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)):</span>
|
||||
</span><span id="read_ms5_xsf-1138"><a href="#read_ms5_xsf-1138"><span class="linenos">1138</span></a> <span class="n">corrres</span><span class="p">[</span><span class="n">i</span> <span class="o">%</span> <span class="mi">2</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">[</span><span class="n">i</span><span class="p">])</span>
|
||||
</span><span id="read_ms5_xsf-1139"><a href="#read_ms5_xsf-1139"><span class="linenos">1139</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">int</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)</span> <span class="o">/</span> <span class="mi">2</span><span class="p">)):</span>
|
||||
</span><span id="read_ms5_xsf-1140"><a href="#read_ms5_xsf-1140"><span class="linenos">1140</span></a> <span class="n">realsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">][</span><span class="n">t</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">corrres</span><span class="p">[</span><span class="mi">0</span><span class="p">][</span><span class="n">t</span><span class="p">])</span>
|
||||
</span><span id="read_ms5_xsf-1141"><a href="#read_ms5_xsf-1141"><span class="linenos">1141</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">int</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)</span> <span class="o">/</span> <span class="mi">2</span><span class="p">)):</span>
|
||||
</span><span id="read_ms5_xsf-1142"><a href="#read_ms5_xsf-1142"><span class="linenos">1142</span></a> <span class="n">imagsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">][</span><span class="n">t</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">corrres</span><span class="p">[</span><span class="mi">1</span><span class="p">][</span><span class="n">t</span><span class="p">])</span>
|
||||
</span><span id="read_ms5_xsf-1143"><a href="#read_ms5_xsf-1143"><span class="linenos">1143</span></a> <span class="n">repnum</span> <span class="o">+=</span> <span class="mi">1</span>
|
||||
</span><span id="read_ms5_xsf-1144"><a href="#read_ms5_xsf-1144"><span class="linenos">1144</span></a>
|
||||
</span><span id="read_ms5_xsf-1145"><a href="#read_ms5_xsf-1145"><span class="linenos">1145</span></a> <span class="n">s</span> <span class="o">=</span> <span class="s2">"Read correlator "</span> <span class="o">+</span> <span class="n">corr</span> <span class="o">+</span> <span class="s2">" from "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">repnum</span><span class="p">)</span> <span class="o">+</span> <span class="s2">" replika with "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">realsamples</span><span class="p">[</span><span class="mi">0</span><span class="p">][</span><span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="read_ms5_xsf-1146"><a href="#read_ms5_xsf-1146"><span class="linenos">1146</span></a> <span class="k">for</span> <span class="n">rep</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span> <span class="n">repnum</span><span class="p">):</span>
|
||||
</span><span id="read_ms5_xsf-1147"><a href="#read_ms5_xsf-1147"><span class="linenos">1147</span></a> <span class="n">s</span> <span class="o">+=</span> <span class="s2">", "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">realsamples</span><span class="p">[</span><span class="n">rep</span><span class="p">][</span><span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="read_ms5_xsf-1148"><a href="#read_ms5_xsf-1148"><span class="linenos">1148</span></a> <span class="n">s</span> <span class="o">+=</span> <span class="s2">" samples"</span>
|
||||
</span><span id="read_ms5_xsf-1149"><a href="#read_ms5_xsf-1149"><span class="linenos">1149</span></a> <span class="nb">print</span><span class="p">(</span><span class="n">s</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1150"><a href="#read_ms5_xsf-1150"><span class="linenos">1150</span></a> <span class="nb">print</span><span class="p">(</span><span class="s2">"Asserted run parameters:</span><span class="se">\n</span><span class="s2"> T:"</span><span class="p">,</span> <span class="n">tmax</span><span class="p">,</span> <span class="s2">"kappa:"</span><span class="p">,</span> <span class="n">kappa</span><span class="p">,</span> <span class="s2">"csw:"</span><span class="p">,</span> <span class="n">csw</span><span class="p">,</span> <span class="s2">"dF:"</span><span class="p">,</span> <span class="n">dF</span><span class="p">,</span> <span class="s2">"zF:"</span><span class="p">,</span> <span class="n">zF</span><span class="p">,</span> <span class="s2">"bnd:"</span><span class="p">,</span> <span class="n">bnd</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1151"><a href="#read_ms5_xsf-1151"><span class="linenos">1151</span></a>
|
||||
</span><span id="read_ms5_xsf-1152"><a href="#read_ms5_xsf-1152"><span class="linenos">1152</span></a> <span class="c1"># we have the data now... but we need to re format the whole thing and put it into Corr objects.</span>
|
||||
</span><span id="read_ms5_xsf-1153"><a href="#read_ms5_xsf-1153"><span class="linenos">1153</span></a>
|
||||
</span><span id="read_ms5_xsf-1154"><a href="#read_ms5_xsf-1154"><span class="linenos">1154</span></a> <span class="n">compObs</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="read_ms5_xsf-1155"><a href="#read_ms5_xsf-1155"><span class="linenos">1155</span></a>
|
||||
</span><span id="read_ms5_xsf-1156"><a href="#read_ms5_xsf-1156"><span class="linenos">1156</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">int</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)</span> <span class="o">/</span> <span class="mi">2</span><span class="p">)):</span>
|
||||
</span><span id="read_ms5_xsf-1157"><a href="#read_ms5_xsf-1157"><span class="linenos">1157</span></a> <span class="n">compObs</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">CObs</span><span class="p">(</span><span class="n">Obs</span><span class="p">([</span><span class="n">realsamples</span><span class="p">[</span><span class="n">rep</span><span class="p">][</span><span class="n">t</span><span class="p">]</span> <span class="k">for</span> <span class="n">rep</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">repnum</span><span class="p">)],</span> <span class="n">names</span><span class="o">=</span><span class="n">names</span><span class="p">,</span> <span class="n">idl</span><span class="o">=</span><span class="n">cnfgs</span><span class="p">),</span>
|
||||
</span><span id="read_ms5_xsf-1158"><a href="#read_ms5_xsf-1158"><span class="linenos">1158</span></a> <span class="n">Obs</span><span class="p">([</span><span class="n">imagsamples</span><span class="p">[</span><span class="n">rep</span><span class="p">][</span><span class="n">t</span><span class="p">]</span> <span class="k">for</span> <span class="n">rep</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">repnum</span><span class="p">)],</span> <span class="n">names</span><span class="o">=</span><span class="n">names</span><span class="p">,</span> <span class="n">idl</span><span class="o">=</span><span class="n">cnfgs</span><span class="p">)))</span>
|
||||
</span><span id="read_ms5_xsf-1159"><a href="#read_ms5_xsf-1159"><span class="linenos">1159</span></a>
|
||||
</span><span id="read_ms5_xsf-1160"><a href="#read_ms5_xsf-1160"><span class="linenos">1160</span></a> <span class="k">if</span> <span class="nb">len</span><span class="p">(</span><span class="n">compObs</span><span class="p">)</span> <span class="o">==</span> <span class="mi">1</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1161"><a href="#read_ms5_xsf-1161"><span class="linenos">1161</span></a> <span class="k">return</span> <span class="n">compObs</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1162"><a href="#read_ms5_xsf-1162"><span class="linenos">1162</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1163"><a href="#read_ms5_xsf-1163"><span class="linenos">1163</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">compObs</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1136"><a href="#read_ms5_xsf-1136"><span class="linenos">1136</span></a> <span class="k">if</span> <span class="n">corr</span> <span class="ow">not</span> <span class="ow">in</span> <span class="n">placesBB</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1137"><a href="#read_ms5_xsf-1137"><span class="linenos">1137</span></a> <span class="n">tmpcorr</span> <span class="o">=</span> <span class="n">asascii</span><span class="p">[</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="n">placesBI</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">):</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="n">placesBI</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1138"><a href="#read_ms5_xsf-1138"><span class="linenos">1138</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1139"><a href="#read_ms5_xsf-1139"><span class="linenos">1139</span></a> <span class="n">tmpcorr</span> <span class="o">=</span> <span class="n">asascii</span><span class="p">[</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="nb">len</span><span class="p">(</span><span class="n">placesBI</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">placesBB</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">):</span><span class="mi">1</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">tmax</span> <span class="o">*</span> <span class="nb">len</span><span class="p">(</span><span class="n">placesBI</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">placesBB</span><span class="o">.</span><span class="n">index</span><span class="p">(</span><span class="n">corr</span><span class="p">)</span> <span class="o">+</span> <span class="mi">2</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1140"><a href="#read_ms5_xsf-1140"><span class="linenos">1140</span></a>
|
||||
</span><span id="read_ms5_xsf-1141"><a href="#read_ms5_xsf-1141"><span class="linenos">1141</span></a> <span class="n">corrres</span> <span class="o">=</span> <span class="p">[[],</span> <span class="p">[]]</span>
|
||||
</span><span id="read_ms5_xsf-1142"><a href="#read_ms5_xsf-1142"><span class="linenos">1142</span></a> <span class="k">for</span> <span class="n">i</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)):</span>
|
||||
</span><span id="read_ms5_xsf-1143"><a href="#read_ms5_xsf-1143"><span class="linenos">1143</span></a> <span class="n">corrres</span><span class="p">[</span><span class="n">i</span> <span class="o">%</span> <span class="mi">2</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">[</span><span class="n">i</span><span class="p">])</span>
|
||||
</span><span id="read_ms5_xsf-1144"><a href="#read_ms5_xsf-1144"><span class="linenos">1144</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">int</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)</span> <span class="o">/</span> <span class="mi">2</span><span class="p">)):</span>
|
||||
</span><span id="read_ms5_xsf-1145"><a href="#read_ms5_xsf-1145"><span class="linenos">1145</span></a> <span class="n">realsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">][</span><span class="n">t</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">corrres</span><span class="p">[</span><span class="mi">0</span><span class="p">][</span><span class="n">t</span><span class="p">])</span>
|
||||
</span><span id="read_ms5_xsf-1146"><a href="#read_ms5_xsf-1146"><span class="linenos">1146</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">int</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)</span> <span class="o">/</span> <span class="mi">2</span><span class="p">)):</span>
|
||||
</span><span id="read_ms5_xsf-1147"><a href="#read_ms5_xsf-1147"><span class="linenos">1147</span></a> <span class="n">imagsamples</span><span class="p">[</span><span class="n">repnum</span><span class="p">][</span><span class="n">t</span><span class="p">]</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">corrres</span><span class="p">[</span><span class="mi">1</span><span class="p">][</span><span class="n">t</span><span class="p">])</span>
|
||||
</span><span id="read_ms5_xsf-1148"><a href="#read_ms5_xsf-1148"><span class="linenos">1148</span></a> <span class="n">repnum</span> <span class="o">+=</span> <span class="mi">1</span>
|
||||
</span><span id="read_ms5_xsf-1149"><a href="#read_ms5_xsf-1149"><span class="linenos">1149</span></a>
|
||||
</span><span id="read_ms5_xsf-1150"><a href="#read_ms5_xsf-1150"><span class="linenos">1150</span></a> <span class="n">s</span> <span class="o">=</span> <span class="s2">"Read correlator "</span> <span class="o">+</span> <span class="n">corr</span> <span class="o">+</span> <span class="s2">" from "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="n">repnum</span><span class="p">)</span> <span class="o">+</span> <span class="s2">" replika with "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">realsamples</span><span class="p">[</span><span class="mi">0</span><span class="p">][</span><span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="read_ms5_xsf-1151"><a href="#read_ms5_xsf-1151"><span class="linenos">1151</span></a> <span class="k">for</span> <span class="n">rep</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span> <span class="n">repnum</span><span class="p">):</span>
|
||||
</span><span id="read_ms5_xsf-1152"><a href="#read_ms5_xsf-1152"><span class="linenos">1152</span></a> <span class="n">s</span> <span class="o">+=</span> <span class="s2">", "</span> <span class="o">+</span> <span class="nb">str</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">realsamples</span><span class="p">[</span><span class="n">rep</span><span class="p">][</span><span class="n">t</span><span class="p">]))</span>
|
||||
</span><span id="read_ms5_xsf-1153"><a href="#read_ms5_xsf-1153"><span class="linenos">1153</span></a> <span class="n">s</span> <span class="o">+=</span> <span class="s2">" samples"</span>
|
||||
</span><span id="read_ms5_xsf-1154"><a href="#read_ms5_xsf-1154"><span class="linenos">1154</span></a> <span class="nb">print</span><span class="p">(</span><span class="n">s</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1155"><a href="#read_ms5_xsf-1155"><span class="linenos">1155</span></a> <span class="nb">print</span><span class="p">(</span><span class="s2">"Asserted run parameters:</span><span class="se">\n</span><span class="s2"> T:"</span><span class="p">,</span> <span class="n">tmax</span><span class="p">,</span> <span class="s2">"kappa:"</span><span class="p">,</span> <span class="n">kappa</span><span class="p">,</span> <span class="s2">"csw:"</span><span class="p">,</span> <span class="n">csw</span><span class="p">,</span> <span class="s2">"dF:"</span><span class="p">,</span> <span class="n">dF</span><span class="p">,</span> <span class="s2">"zF:"</span><span class="p">,</span> <span class="n">zF</span><span class="p">,</span> <span class="s2">"bnd:"</span><span class="p">,</span> <span class="n">bnd</span><span class="p">)</span>
|
||||
</span><span id="read_ms5_xsf-1156"><a href="#read_ms5_xsf-1156"><span class="linenos">1156</span></a>
|
||||
</span><span id="read_ms5_xsf-1157"><a href="#read_ms5_xsf-1157"><span class="linenos">1157</span></a> <span class="c1"># we have the data now... but we need to re format the whole thing and put it into Corr objects.</span>
|
||||
</span><span id="read_ms5_xsf-1158"><a href="#read_ms5_xsf-1158"><span class="linenos">1158</span></a>
|
||||
</span><span id="read_ms5_xsf-1159"><a href="#read_ms5_xsf-1159"><span class="linenos">1159</span></a> <span class="n">compObs</span> <span class="o">=</span> <span class="p">[]</span>
|
||||
</span><span id="read_ms5_xsf-1160"><a href="#read_ms5_xsf-1160"><span class="linenos">1160</span></a>
|
||||
</span><span id="read_ms5_xsf-1161"><a href="#read_ms5_xsf-1161"><span class="linenos">1161</span></a> <span class="k">for</span> <span class="n">t</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">int</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">tmpcorr</span><span class="p">)</span> <span class="o">/</span> <span class="mi">2</span><span class="p">)):</span>
|
||||
</span><span id="read_ms5_xsf-1162"><a href="#read_ms5_xsf-1162"><span class="linenos">1162</span></a> <span class="n">compObs</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">CObs</span><span class="p">(</span><span class="n">Obs</span><span class="p">([</span><span class="n">realsamples</span><span class="p">[</span><span class="n">rep</span><span class="p">][</span><span class="n">t</span><span class="p">]</span> <span class="k">for</span> <span class="n">rep</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">repnum</span><span class="p">)],</span> <span class="n">names</span><span class="o">=</span><span class="n">names</span><span class="p">,</span> <span class="n">idl</span><span class="o">=</span><span class="n">cnfgs</span><span class="p">),</span>
|
||||
</span><span id="read_ms5_xsf-1163"><a href="#read_ms5_xsf-1163"><span class="linenos">1163</span></a> <span class="n">Obs</span><span class="p">([</span><span class="n">imagsamples</span><span class="p">[</span><span class="n">rep</span><span class="p">][</span><span class="n">t</span><span class="p">]</span> <span class="k">for</span> <span class="n">rep</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">repnum</span><span class="p">)],</span> <span class="n">names</span><span class="o">=</span><span class="n">names</span><span class="p">,</span> <span class="n">idl</span><span class="o">=</span><span class="n">cnfgs</span><span class="p">)))</span>
|
||||
</span><span id="read_ms5_xsf-1164"><a href="#read_ms5_xsf-1164"><span class="linenos">1164</span></a>
|
||||
</span><span id="read_ms5_xsf-1165"><a href="#read_ms5_xsf-1165"><span class="linenos">1165</span></a> <span class="k">if</span> <span class="nb">len</span><span class="p">(</span><span class="n">compObs</span><span class="p">)</span> <span class="o">==</span> <span class="mi">1</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1166"><a href="#read_ms5_xsf-1166"><span class="linenos">1166</span></a> <span class="k">return</span> <span class="n">compObs</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span>
|
||||
</span><span id="read_ms5_xsf-1167"><a href="#read_ms5_xsf-1167"><span class="linenos">1167</span></a> <span class="k">else</span><span class="p">:</span>
|
||||
</span><span id="read_ms5_xsf-1168"><a href="#read_ms5_xsf-1168"><span class="linenos">1168</span></a> <span class="k">return</span> <span class="n">Corr</span><span class="p">(</span><span class="n">compObs</span><span class="p">)</span>
|
||||
</span></pre></div>
|
||||
|
||||
|
||||
|
|
|
|||
Loading…
Add table
Add a link
Reference in a new issue