Changeset 1995
- Timestamp:
- Oct 9, 2015 10:07:14 AM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/GSASIImath.py
r1993 r1995 1092 1092 1093 1093 def posFourier(tau,psin,pcos,smul): 1094 A = np.array([ps[:,np.newaxis]*np.sin(2*np.pi*(i+1)*tau) for i,ps in enumerate(psin)]) *smul1094 A = np.array([ps[:,np.newaxis]*np.sin(2*np.pi*(i+1)*tau) for i,ps in enumerate(psin)]) #*smul 1095 1095 B = np.array([pc[:,np.newaxis]*np.cos(2*np.pi*(i+1)*tau) for i,pc in enumerate(pcos)]) 1096 1096 return np.sum(A,axis=0)+np.sum(B,axis=0) … … 2322 2322 h,k,l,m = -hkl+Hmax 2323 2323 Fhkl[h,k,l,m] = dF*phasem 2324 SSrho = fft.fftn(fft.fftshift(Fhkl))/ (10.*cell[6])#4D map2324 SSrho = fft.fftn(fft.fftshift(Fhkl))/cell[6] #4D map 2325 2325 rho = fft.fftn(fft.fftshift(Fhkl[:,:,:,maxM+1]))/cell[6] #3D map 2326 2326 map4DData['rho'] = np.real(SSrho)
Note: See TracChangeset
for help on using the changeset viewer.