Changeset 2230 for trunk/GSASIIphsGUI.py
- Timestamp:
- Apr 29, 2016 8:30:58 AM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/GSASIIphsGUI.py
r2229 r2230 2109 2109 for ind in indx: 2110 2110 XYZ = np.array(atomData[ind][cx:cx+3]) 2111 for unit in Unit: 2112 XYZ += unit 2113 XYZ -= T 2114 XYZ = np.inner(A,XYZ) #to Cartesian 2115 XYZ = np.inner(M,XYZ) #rotate 2116 XYZ = np.inner(B,XYZ)+T #back to crystal & translate 2117 if np.all(XYZ>=0.) and np.all(XYZ<1.0): 2118 atom = atomData[ind] 2119 atom[cx:cx+3] = XYZ 2120 atom[css:css+2] = G2spc.SytSym(XYZ,SGData) 2121 break 2111 XYZS = XYZ+Unit 2112 XYZS -= T 2113 XYZS = np.inner(A,XYZS).T #to Cartesian 2114 XYZS = np.inner(M,XYZS).T #rotate 2115 XYZS = np.inner(B,XYZS).T+T #back to crystal & translate 2116 GSASIIpath.IPyBreak() 2117 atomData[ind][cx:cx+3] = XYZ 2118 # for unit in Unit: 2119 # XYZ = np.copy(np.array(atomData[ind][cx:cx+3])) 2120 # XYZ += unit 2121 # XYZ -= T 2122 # XYZ = np.inner(A,XYZ) #to Cartesian 2123 # XYZ = np.inner(M,XYZ) #rotate 2124 # XYZ = np.inner(B,XYZ)+T #back to crystal & translate 2125 # if np.all(XYZ>=0.) and np.all(XYZ<1.0): 2126 # atomData[ind][cx:cx+3] = XYZ 2127 ## atom[css:css+2] = G2spc.SytSym(XYZ,SGData) 2128 # break 2122 2129 finally: 2123 2130 dlg.Destroy()
Note: See TracChangeset
for help on using the changeset viewer.