from dials.algorithms.reflection_basis import FromRotationAngleFast # s0 = matrix.col((0.00801639379479156, -3.8514801707506e-14, -1.0208975723528189) ) # m2 = matrix.col((1.0, -1.5919306617286774e-16, -6.904199434387693e-16)) # s1 = matrix.col((0.502248600762, -0.0181543707198, -0.888657908118)) # s0 = matrix.col((0.0, 0.0, -1.0)) # m2 = matrix.col((1.0, 0.0, 0.0)) # s1 = matrix.col((1.0, 0.05, 0.0)).normalize() phi = 0 # c3 = -0.013439 n = 5 c3 = -n * 0.157 * pi / 180.0 cs = CoordinateSystem(m2, s0, s1, phi) print(cs.zeta()) # e1 = matrix.col(cs.e1_axis()) # e3 = matrix.col(cs.e3_axis()) # ps = matrix.col(cs.p_star()) # from math import sqrt # m2e1 = m2.dot(e1) # m2e3 = m2.dot(e3) # m2ps = m2.dot(ps.normalize()) # print tuple(e1) # print tuple(e3) # print tuple(ps) # print "m2.e1", m2e1 # print "m2.e3", m2e3 # print "m2.ps", m2ps
from dials.algorithms.reflection_basis import FromRotationAngleFast #s0 = matrix.col((0.00801639379479156, -3.8514801707506e-14, -1.0208975723528189) ) #m2 = matrix.col((1.0, -1.5919306617286774e-16, -6.904199434387693e-16)) #s1 = matrix.col((0.502248600762, -0.0181543707198, -0.888657908118)) #s0 = matrix.col((0.0, 0.0, -1.0)) #m2 = matrix.col((1.0, 0.0, 0.0)) #s1 = matrix.col((1.0, 0.05, 0.0)).normalize() phi = 0 #c3 = -0.013439 n = 5 c3 = -n * 0.157 * pi / 180.0 cs = CoordinateSystem(m2, s0, s1, phi) print cs.zeta() #e1 = matrix.col(cs.e1_axis()) #e3 = matrix.col(cs.e3_axis()) #ps = matrix.col(cs.p_star()) #from math import sqrt #m2e1 = m2.dot(e1) #m2e3 = m2.dot(e3) #m2ps = m2.dot(ps.normalize()) #print tuple(e1) #print tuple(e3) #print tuple(ps) #print "m2.e1", m2e1 #print "m2.e3", m2e3 #print "m2.ps", m2ps