<div dir="ltr"><div>Hi Bert,</div><div><br></div><div>actually, my "first" thought solves this problem. If you turn all your points into an xray_structure, and get the asu_mappings, this gives you the coordinates of the equivalent site in the asu (returned by the "mapped_site" method). However, it won't make a list of equivalent sites unique, i.e. if you have site1 = ( 0.1, 0.1, 0.1 ) and site2 = ( 1.1, 0.1, 0.1 ) (in fractional coordinates), you get a mapping for the same site twice, so perhaps the most efficient method is just to transform your scalar field into the unit cell, and then discard points that are outside the asu. If the grid is not congruent with the unit cell, you are probably better off using asu_mappings, since you will be able to interpolate on a finer grid.</div><div><br></div><div>BW, Gabor<br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 21, 2016 at 2:02 PM, Tijskens Engelbert <span dir="ltr"><<a href="mailto:Engelbert.Tijskens@uantwerpen.be" target="_blank">Engelbert.Tijskens@uantwerpen.be</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">> Message: 3<br>
> Date: Mon, 21 Nov 2016 13:38:21 +0000<br>
> From: Gabor Bunkoczi <<a href="mailto:gabor.bunkoczi@googlemail.com">gabor.bunkoczi@googlemail.com</a><wbr>><br>
> To: cctbx mailing list <<a href="mailto:cctbxbb@phenix-online.org">cctbxbb@phenix-online.org</a>><br>
> Subject: Re: [cctbxbb] convert an arbitrary point to one inside the<br>
> asymmetric unit<br>
> Message-ID:<br>
> <CABYxhoUxqFsc0f_<wbr>TjWwn4oC4hKFk8AEbzwx5=9=_+<a href="mailto:gwW3xmkpg@mail.gmail.com">gwW3<wbr>xmkpg@mail.gmail.com</a>><br>
> Content-Type: text/plain; charset="utf-8"<br>
><br>
> On second thought, the simplest solution is to use the direct_space_asu<br>
> object. This is available for all space groups:<br>
><br>
>>>> from cctbx import sgtbx<br>
>>>> sgi = sgtbx.space_group_info( "P2" )<br>
>>>> asu = sgi.direct_space_asu()<br>
>>>> asu.is_inside( ( 0, 0, 0 ) )<br>
> True<br>
>>>> asu.is_inside( (0, 0, 0.51 ) )<br>
> False<br>
><br>
> BW, Gabor<br>
><br>
Thanks for your quick reply, Gabor.<br>
Unfortunately, it does not solve my question. For the example you give above, the second point is not inside the asu, and my question is how i can the corresponding point inside the asu, For a simple space group as P2 the answer is easily figured out by hand. For a space group such as I m -3 m, it becomes complicated. I would like to have a (computationally efficient) method for finding the point inside the asu given an arbitrary point which works for every space group.<br>
Any “third” thoughts?<br>
<br>
kindest regards, bert<br>
<br>
______________________________<wbr>_________________<br>
cctbxbb mailing list<br>
<a href="mailto:cctbxbb@phenix-online.org">cctbxbb@phenix-online.org</a><br>
<a href="http://phenix-online.org/mailman/listinfo/cctbxbb" target="_blank" rel="noreferrer">http://phenix-online.org/<wbr>mailman/listinfo/cctbxbb</a><br>
</blockquote></div><br></div>