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))
