Changeset 5015 for trunk/GSASIImath.py


Ignore:
Timestamp:
Aug 23, 2021 7:51:22 AM (3 months ago)
Author:
vondreele
Message:

changes to comments

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIImath.py

    r5013 r5015  
    16651665    ngl = 36                    #selected for integer steps for 1/6,1/4,1/3...
    16661666    glTau,glWt = pwd.pygauleg(0.,1.,ngl)         #get Gauss-Legendre intervals & weights
    1667     mglTau = np.arange(0.,1.,1./ngl)
    16681667    Ax = np.array(XSSdata[:3]).T   #atoms x waves x sin pos mods
    16691668    Bx = np.array(XSSdata[3:]).T   #...cos pos mods
     
    17141713        Umod = 1.0
    17151714    if nWaves[3]:
    1716         tauM = np.arange(1.,nWaves[3]+1-nx)[:,nxs]*mglTau  #Mwaves x ngl
     1715        tauM = np.arange(1.,nWaves[3]+1-nx)[:,nxs]*glTau  #Mwaves x ngl
    17171716        MmodA = Am[:,:,:,nxs]*np.sin(twopi*tauM[nxs,:,nxs,:]) #atoms X waves X 3 X tau
    17181717        MmodB = Bm[:,:,:,nxs]*np.cos(twopi*tauM[nxs,:,nxs,:]) #ditto
Note: See TracChangeset for help on using the changeset viewer.