Page History
...
Code Block |
---|
xrot2.shape: (13, 9) row: -0.0504 -0.0106 0.0292 0.0690 0.1088 0.1485 0.1883 0.2281 0.2679 row: -0.0685 -0.0287 0.0111 0.0508 0.0906 0.1304 0.1702 0.2100 0.2498 row: -0.0866 -0.0469 -0.0071 0.0327 0.0725 0.1123 0.1521 0.1919 0.2316 row: -0.1048 -0.0650 -0.0252 0.0146 0.0544 0.0942 0.1339 0.1737 0.2135 row: -0.1229 -0.0831 -0.0433 -0.0035 0.0363 0.0760 0.1158 0.1556 0.1954 row: -0.1410 -0.1012 -0.0614 -0.0217 0.0181 0.0579 0.0977 0.1375 0.1773 row: -0.1591 -0.1194 -0.0796 -0.0398 0.0000 0.0398 0.0796 0.1194 0.1591 row: -0.1773 -0.1375 -0.0977 -0.0579 -0.0181 0.0217 0.0614 0.1012 0.1410 row: -0.1954 -0.1556 -0.1158 -0.0760 -0.0363 0.0035 0.0433 0.0831 0.1229 row: -0.2135 -0.1737 -0.1339 -0.0942 -0.0544 -0.0146 0.0252 0.0650 0.1048 row: -0.2316 -0.1919 -0.1521 -0.1123 -0.0725 -0.0327 0.0071 0.0469 0.0866 row: -0.2498 -0.2100 -0.1702 -0.1304 -0.0906 -0.0508 -0.0111 0.0287 0.0685 row: -0.2679 -0.2281 -0.1883 -0.1485 -0.1088 -0.0690 -0.0292 0.0106 0.0504 yrot1.shape: (13, 9) row: -0.3409 -0.3044 -0.2680 -0.2315 -0.1950 -0.1585 -0.1221 -0.0856 -0.0491 row: -0.3084 -0.2719 -0.2355 -0.1990 -0.1625 -0.1260 -0.0896 -0.0531 -0.0166 row: -0.2759 -0.2394 -0.2030 -0.1665 -0.1300 -0.0935 -0.0571 -0.0206 0.0159 row: -0.2434 -0.2069 -0.1705 -0.1340 -0.0975 -0.0610 -0.0245 0.0119 0.0484 row: -0.2109 -0.1744 -0.1380 -0.1015 -0.0650 -0.0285 0.0080 0.0444 0.0809 row: -0.1784 -0.1419 -0.1055 -0.0690 -0.0325 0.0040 0.0405 0.0769 0.1134 row: -0.1459 -0.1094 -0.0730 -0.0365 0.0000 0.0365 0.0730 0.1094 0.1459 row: -0.1134 -0.0769 -0.0405 -0.0040 0.0325 0.0690 0.1055 0.1419 0.1784 row: -0.0809 -0.0444 -0.0080 0.0285 0.0650 0.1015 0.1380 0.1744 0.2109 row: -0.0484 -0.0119 0.0245 0.0610 0.0975 0.1340 0.1705 0.2069 0.2434 row: -0.0159 0.0206 0.0571 0.0935 0.1300 0.1665 0.2030 0.2394 0.2759 row: 0.0166 0.0531 0.0896 0.1260 0.1625 0.1990 0.2355 0.2719 0.3084 row: 0.0491 0.0856 0.1221 0.1585 0.1950 0.2315 0.2680 0.3044 0.3409 zrot2.shape: (13, 9) row: -0.0135 -0.0028 0.0078 0.0185 0.0291 0.0398 0.0505 0.0611 0.0718 row: -0.0184 -0.0077 0.0030 0.0136 0.0243 0.0349 0.0456 0.0563 0.0669 row: -0.0232 -0.0126 -0.0019 0.0088 0.0194 0.0301 0.0407 0.0514 0.0621 row: -0.0281 -0.0174 -0.0068 0.0039 0.0146 0.0252 0.0359 0.0466 0.0572 row: -0.0329 -0.0223 -0.0116 -0.0009 0.0097 0.0204 0.0310 0.0417 0.0524 row: -0.0378 -0.0271 -0.0165 -0.0058 0.0049 0.0155 0.0262 0.0368 0.0475 row: -0.0426 -0.0320 -0.0213 -0.0107 0.0000 0.0107 0.0213 0.0320 0.0426 row: -0.0475 -0.0368 -0.0262 -0.0155 -0.0049 0.0058 0.0165 0.0271 0.0378 row: -0.0524 -0.0417 -0.0310 -0.0204 -0.0097 0.0009 0.0116 0.0223 0.0329 row: -0.0572 -0.0466 -0.0359 -0.0252 -0.0146 -0.0039 0.0068 0.0174 0.0281 row: -0.0621 -0.0514 -0.0407 -0.0301 -0.0194 -0.0088 0.0019 0.0126 0.0232 row: -0.0669 -0.0563 -0.0456 -0.0349 -0.0243 -0.0136 -0.0030 0.0077 0.0184 row: -0.0718 -0.0611 -0.0505 -0.0398 -0.0291 -0.0185 -0.0078 0.0028 0.0135 |
Look-up table
Look-up table is generated in method make_lookup_table
called as
lut = make_lookup_table(b1, b2, b3, hmax, kmax, lmax, np.float32, evald_rad, sigma_q, fout, bpq, bpomega, bpbeta)
for each crystal orientation in omega_deg and beta_deg method
dr, qv, qh = radial_distance(xrot2, yrot1, zrot2, evald_rad)
evaluates
dr [1/A] - distance between lattice node and Evald's sphere
qv, qh [1/A] - vertical and horizontal components of scattering vector q
in a simple code:
Code Block |
---|
def radial_distance(X, Y, Z, evald_rad=0.5) :
DX = X + evald_rad
L = np.sqrt(DX*DX + Y*Y + Z*Z)
dr = L - evald_rad
qv = evald_rad * Z/L
ql = evald_rad * (DX/L-1)
qy = evald_rad * Y/L
qh = np.sqrt(ql*ql + qy*qy) * np.select([Y<0], [-1], default=1)
return dr, qv, qh |
Finally, in method make_lookup_table
txt = str_omega_drhkl(ind, beta_deg, omega_deg, dr, r, qv, qh, h, k, l, sigma_q)
generates text record for file with look-up table, and
lut[iomega,:] += fill_row(dr, qv, qh, h, k, l, sigma_q, bpq)
generates 2-d array, which presents look-up table in plot.