Closedown

2018-11-20 10:43:04 +01:00
parent 53ba051cb6
commit ba528066ca
1 changed files with 24 additions and 8 deletions
@@ -3,6 +3,12 @@
 ###################################################################################################
 set_geometry("fourcv", True)

+###################################################################################################\
+#Constraints
+###################################################################################################
+hkl.con('a_eq_b')
+
+
 ###################################################################################################\
 #Orientation
 ###################################################################################################
@@ -42,15 +48,13 @@ ub.newub('test')
 ub.setub([[ 0.81374,    0.02138,    0.00228], [-0.02199,    0.81543,    0.00123], [ 0.00784,   -0.00100,    0.81347]])
 ub.checkub()

-###################################################################################################\
-#Constraints
-###################################################################################################
-hkl.con('a_eq_b')
 print hkl_to_angles(0, 1, 2)[0]
 print hkl_to_angles(2, 2, 2)[0]
 print hkl_to_angles(2, 0, 2)[0]

-for r in get_ub_matrix(): print r
+orig_ub = get_ub_matrix()
+
+



@@ -67,12 +71,24 @@ en = 9.5
 ub.addref([0, 1, 2], [ 9.7212  ,  9.7488 ,  19.5860, 93.1295], en)
 ub.addref([2, 2, 2], [ 9.8097 ,  27.7194 ,  20.9065 , 57.5840], en)
 ub.addref([2, 0, 2], [9.8217  , 19.5132 ,  20.2488 ,  8.4510], en)
-#ub.ub()
-#ub.checkub()
+ub.ub()
+ub.checkub()


-for r in get_ub_matrix(): print r
+calc_ub =  get_ub_matrix()
 print hkl_to_angles(0, 1, 2)[0]
 print hkl_to_angles(2, 2, 2)[0]
 print hkl_to_angles(2, 0, 2)[0]

+
+
+for r in orig_ub: print r
+for r in calc_ub: print r
+print "---"
+
+print e
+for i in range(3):
+    for j in range(3):
+        e[i][j] = abs(orig_ub[i][j] - calc_ub[i][j])
+
+for r in e: print r