Changeset 4080 for trunk/GSASIIpwd.py
- Timestamp:
- Aug 7, 2019 11:35:30 AM (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/GSASIIpwd.py
r4048 r4080 1112 1112 1113 1113 def ellipseSize(H,Sij,GB): 1114 ' needs a doc string'1114 'Implements r=1/sqrt(sum((1/S)*(q.v)^2) per note from Alexander Brady' 1115 1115 HX = np.inner(H.T,GB) 1116 1116 lenHX = np.sqrt(np.sum(HX**2)) 1117 1117 Esize,Rsize = nl.eigh(G2lat.U6toUij(Sij)) 1118 R = np.inner(HX/lenHX,Rsize)* Esize #want column length for hkl in crystal1119 lenR = np.sqrt(np.sum(R**2))1118 R = np.inner(HX/lenHX,Rsize)**2*Esize #want column length for hkl in crystal 1119 lenR = 1./np.sqrt(np.sum(R)) 1120 1120 return lenR 1121 1121
Note: See TracChangeset
for help on using the changeset viewer.