Python CoordinateSystem.zetaの例

コード例 #1

ファイル: angle_test.py プロジェクト: dials/dials_scratch

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

コード例 #2

ファイル: angle_test.py プロジェクト: dials/dials_scratch

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