def test_correction_on_random_coordinates(model):
convert = ParallaxCorrectedPxMmStrategy(model["mu"], model["t0"])
for i in range(10000):
xy = matrix.col((random.uniform(-1000, 1000), random.uniform(-1000, 1000)))
# Do the forward and reverse corrections
xy_corr_gold = matrix.col(correct_gold(model, xy))
xy_corr = matrix.col(convert.to_pixel(model["detector"], xy))
xy_corr_inv = matrix.col(convert.to_millimeter(model["detector"], xy_corr))
# Check the values
assert abs(xy_corr_gold - xy_corr) < 1e-7
assert abs(xy_corr_inv - xy) < 1e-3
def correct_inv(self, xy):
from dxtbx.model import ParallaxCorrectedPxMmStrategy
convert = ParallaxCorrectedPxMmStrategy(self.mu, self.t0)
return convert.to_millimeter(self.detector[0], xy)