<div dir="ltr"><div class="gmail_extra">On Sun, Oct 27, 2013 at 6:36 AM, Keitaro Yamashita <span dir="ltr"><<a href="mailto:k.yamashita@spring8.or.jp" target="_blank">k.yamashita@spring8.or.jp</a>></span> wrote:<br><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">When I changed the cctbx code to sort reflections after mapping to<br>
a.s.u., CC1/2 values by phenix.merging_statistics became closer to<br>
XSCALE.<br>
<br>
In merging_stats.__init__() in iotbx/merging_statistics.py, I inserted a line<br>
> � �array = array.sort("packed_indices")<br>
after<br>
> � �array = array.customized_copy(anomalous_flag=anomalous).map_to_asu()<br>
<br>
but I think there would be more correct way to fix this unexpected behavior.<br></blockquote><div><br></div><div style>If you can think of a better method, feel free to implement it, but I think the above is quite appropriate. �Assuming this passes the regression test (iotbx/regression/tst_merging_statistics.py), it's fine to check in.��The only drawback is that it's probably inefficient (especially for arrays that are already sorted), but the entire program isn't very efficient to begin with, and a more comprehensive solution would require a ground-up approach mostly coded in C++. �This will definitely become necessary at some point but none of us have time to spare right now.</div>
<div style><br></div><div style>-Nat</div></div></div></div>