U12 = U21 = atom_obj[u-1]["u[0][1]"]
U13 = U31 = atom_obj[u-1]["u[0][2]"]
U23 = U32 = atom_obj[u-1]["u[1][2]"]
U = tlsvld.array([[U11,U12,U13],[U21,U22,U23],[U31,U32,U33]])
V11 = atom_obj[u]["u[0][0]"]
V22 = atom_obj[u]["u[1][1]"]
V33 = atom_obj[u]["u[2][2]"]
V12 = V21 = atom_obj[u]["u[0][1]"]
V13 = V31 = atom_obj[u]["u[0][2]"]
V23 = V32 = atom_obj[u]["u[1][2]"]
V = tlsvld.array([[V11,V12,V13],[V21,V22,V23],[V31,V32,V33]])
## Calculate all of the residuals
Suij = rosenfeld = diff_trace_UV = sum_square_diff = 0.0
cc_UV = tlsvld.calc_ccuij(U, V)
#cc_UV = AtomMath.calc_CCuij([[U11,U12,U13],
# [U21,U22,U23],
# [U31,U32,U33]],
# [[V11,V12,V13],
# [V21,V22,V23],
# [V31,V32,V33]])
rosenfeld = tlsvld.rosenfeld(Upos, Vpos, U, V)
sum_square_diff = tlsvld.sum_square_diff(U, V)
diff_trace_UV = tlsvld.diff_trace_uv(U, V)
## The actual output of this entire script.
print "%s\t%s:%s-%s:%s\t%s\t%.10f\t%.10f\t%.10f\t%.10f\t%.10f" % (
bond_pos, cid1, res_name1, res_num1, bond, tls_group,
cc_UV, Suij, rosenfeld/10000,
diff_trace_UV/10000, sum_square_diff/10000**2)
#print "="*80 #U = scipy.array([[0.37471362,-0.07411912,-0.0379023],[-0.07411912,0.39010298,-0.03240841],[-0.0379023,-0.03240841,0.33057791]]) #V = scipy.array([[0.38484156,-0.06933809,-0.04994903],[-0.06933809,0.41376374,-0.03275701],[-0.04994903,-0.03275701,0.33967632]]) #V = scipy.array([[100.38484156,-0.06933809,-0.04994903],[-0.06933809,0.41376374,-0.03275701],[-0.04994903,-0.03275701,0.33967632]]) #print U #print V #print "="*80 #U = tlsvld.array([[7,7,7],[7,7,7],[7,7,7]]) #V = tlsvld.array([[4,4,4],[4,4,4],[4,4,4]]) W = tlsvld.array([[1.0,1.0,1.0],[1.0,1.0,1.0],[1.0,1.0,1.0]]) X = tlsvld.array([[4.0,4.0,4.0],[4.0,4.0,4.0],[4.0,4.0,4.0]]) print "SQD = ",tlsvld.sum_square_diff(W,X) #print "trace(U) = ", tlsvld.trace3x3(X) print "dtUV = ",tlsvld.diff_trace_uv(W,X) print tlsvld.mul3x3(U,V) print tlsvld.add3x3(U,V) print "SUM(W) = ",tlsvld.sum3x3(X) print tlsvld.subtract3x3(U,V) #print tlsvld.trn3x3(U) #print tlsvld.inv3x3(U) #W = scipy.array([[0,0,0],[0,0,0],[0,0,0]]) #det, invU = tlsvld.inv3x3(U) #print "PYTHON: det(U) = %s" % det #print "PYTHON: inv(U) = %s" % invU #print "="*80 #X = scipy.array([[1,0,0],[0,2,0],[0,0,3]]) #Y = scipy.array([[1,0,0],[0,2,0],[0,0,3]])
continue
#U = Ueq[u-1]
#V = Ueq[u]
U11 = atom_obj[u-1]["u[0][0]"]
U22 = atom_obj[u-1]["u[1][1]"]
U33 = atom_obj[u-1]["u[2][2]"]
U12 = U21 = atom_obj[u-1]["u[0][1]"]
U13 = U31 = atom_obj[u-1]["u[0][2]"]
U23 = U32 = atom_obj[u-1]["u[1][2]"]
U = tlsvld.array([[U11,U12,U13],[U21,U22,U23],[U31,U32,U33]])
V11 = atom_obj[u]["u[0][0]"]
V22 = atom_obj[u]["u[1][1]"]
V33 = atom_obj[u]["u[2][2]"]
V12 = V21 = atom_obj[u]["u[0][1]"]
V13 = V31 = atom_obj[u]["u[0][2]"]
V23 = V32 = atom_obj[u]["u[1][2]"]
V = tlsvld.array([[V11,V12,V13],[V21,V22,V23],[V31,V32,V33]])
Suij = rosenfeld = diff_trace_UV = sum_square_diff = 0.0
cc_UV = tlsvld.calc_ccuij(U, V)
#cc_UV = AtomMath.calc_CCuij([[U11,U12,U13],[U21,U22,U23],[U31,U32,U33]], [[V11,V12,V13],[V21,V22,V23],[V31,V32,V33]])
rosenfeld = tlsvld.rosenfeld(Upos, Vpos, U, V)
sum_square_diff = tlsvld.sum_square_diff(U, V)
diff_trace_UV = tlsvld.diff_trace_uv(U, V)
print "%s\t%s:%s-%s:%s\t%s\t%.10f\t%.10f\t%.10f\t%.10f\t%.10f" % (
bond_pos, cid1, res_name1, res_num1, bond, tls_group,
cc_UV, Suij, rosenfeld/10000,
diff_trace_UV/10000, sum_square_diff/10000**2)
