def rt_random2(T1, T2): theta = (random.random() - 0.5) * math.pi R = AtomMath.rmatrixu(random_vec(), theta) T1R = numpy.dot(numpy.dot(R, T1), numpy.transpose(R)) T2R = numpy.dot(numpy.dot(R, T2), numpy.transpose(R)) return T1, T2
def rt_random(T): theta = (random.random() - 0.5) * math.pi R = AtomMath.rmatrixu(random_vec(), theta) return numpy.dot(numpy.dot(R, T), numpy.transpose(R))