[cctbxbb] cctbx.miller.complete_array() bug?
Luc Bourhis
luc_j_bourhis at mac.com
Thu Jul 12 08:37:46 PDT 2012
Hi Jan,
> As you can see the index (0, 1, 2) is added with zero amplitude where instead
> the index (1, 0, 2) should be renamed to (0, 1, 2) as those two are
> symmetrically equivalent for any tetragonal spacegroup (like P4). So to me it
> looks like the complete_array method does not honor symmetrically equivalent
> indices at all, where it should imho. What do you think on this issue?
Considering the current implementation of method complete_array, it seems to me that the problem lies in method complete_set of class miller.set, that replaces (1, 0, 2) with (0, 1, 2) in your example for no particularly good reason I can see. The implementation of complete_array clearly assumes that complete_set should not do that. But then, the way complete_set is implemented, such surprises are bound to happen. A more involved algorithm is necessary there.
I put that on the todo list.
Best wishes,
Luc Bourhis
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://phenix-online.org/pipermail/cctbxbb/attachments/20120712/8d645d0b/attachment.htm>
More information about the cctbxbb
mailing list