def tst_centroid_points3d(self): from dials.algorithms.image.centroid import centroid_points from scitbx import matrix centroid = centroid_points(self.pixels3d.as_1d(), self.points3d.as_1d()) assert(abs(self.gold3d - matrix.col(centroid.mean())) < self.EPS) assert(abs(self.gold3dvar - matrix.col(centroid.variance())) < self.EPS) assert(abs(self.gold3dubvar - matrix.col(centroid.unbiased_variance())) < self.EPS) print 'OK'
def test_centroid_points3d(self): from dials.algorithms.image.centroid import centroid_points from scitbx import matrix centroid = centroid_points(self.pixels3d.as_1d(), self.points3d.as_1d()) assert abs(self.gold3d - matrix.col(centroid.mean())) < self.EPS assert abs(self.gold3dvar - matrix.col(centroid.variance())) < self.EPS assert ( abs(self.gold3dubvar - matrix.col(centroid.unbiased_variance())) < self.EPS )