Changeset 2195 for trunk/GSASIIpwd.py
- Timestamp:
- Apr 3, 2016 2:30:21 PM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/GSASIIpwd.py
r2194 r2195 40 40 import pydiffax as pyx 41 41 except ImportError: 42 print ' the pydiffax library could not be loaded'42 print 'pydiffax is not available for this platform - under develpment' 43 43 44 44 … … 1916 1916 df.close() 1917 1917 time0 = time.time() 1918 subp.call(DIFFaX) 1918 try: 1919 subp.call(DIFFaX) 1920 except OSError: 1921 print 'DIFFax.exe is not available for this platform - under development' 1919 1922 print 'DIFFaX time = %.2fs'%(time.time()-time0) 1920 1923 if os.path.exists('GSASII-DIFFaX.spc'): … … 1979 1982 planeId = ['h0l','0kl','hhl','h-hl'].index(Layers['Sadp']['Plane'])+1 1980 1983 lmax = int(Layers['Sadp']['Lmax']) 1981 mult = int(Layers['Sadp']['Mult'])1984 mult = 1 1982 1985 # Sequences 1983 1986 StkType = ['recursive','explicit'].index(Layers['Stacking'][0]) … … 2046 2049 hkLim,Incr,Nblk = pyx.pygetsadp(controls,Nspec,spec) 2047 2050 Sapd = np.zeros((256,256)) 2048 maxInt = np.max(spec[1:])2049 Scale = mult*32767./maxInt2050 2051 iB = 0 2051 2052 for i in range(hkLim): … … 2064 2065 Sapd[:,p2] = spec[iB:iF] 2065 2066 iB += Nblk 2066 Sapd *= Scale2067 Sapd = np.where(Sapd<32767.,Sapd,32767.)2068 2067 Layers['Sadp']['Img'] = Sapd 2069 2068 print 'GETSAD time = %.2fs'%(time.time()-time0)
Note: See TracChangeset
for help on using the changeset viewer.