Changeset 3821


Ignore:
Timestamp:
Feb 13, 2019 4:11:20 PM (3 years ago)
Author:
vondreele
Message:

more changes for incommensurate mag str fctr
remove GUI update after changes to PDF file items - wasn't needed. Does recompute PDF though.

Location:
trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIImath.py

    r3818 r3821  
    12981298    ngl = 36                    #selected for integer steps for 1/6,1/4,1/3...
    12991299    glTau,glWt = pwd.pygauleg(0.,1.,ngl)         #get Gauss-Legendre intervals & weights
     1300    mglTau = np.arange(0.,1.,1./ngl)
    13001301    Ax = np.array(XSSdata[:3]).T   #atoms x waves x sin pos mods
    13011302    Bx = np.array(XSSdata[3:]).T   #...cos pos mods
     
    13461347        Umod = 1.0
    13471348    if nWaves[3]:
    1348         tauM = np.arange(1.,nWaves[3]+1-nx)[:,nxs]*glTau  #Mwaves x ngl
     1349        tauM = np.arange(1.,nWaves[3]+1-nx)[:,nxs]*mglTau  #Mwaves x ngl
    13491350        MmodA = Am[:,:,:,nxs]*np.sin(twopi*tauM[nxs,:,nxs,:]) #atoms X waves X 3 X ngl
    13501351        MmodB = Bm[:,:,:,nxs]*np.cos(twopi*tauM[nxs,:,nxs,:]) #ditto
  • trunk/GSASIIpwdGUI.py

    r3810 r3821  
    64096409                data[fileKey]['Mult'] = GetExposure(value)
    64106410                ResetFlatBkg()
    6411                 wx.CallLater(100,UpdatePDFGrid,G2frame,data)
    64126411                wx.CallAfter(OnComputePDF,None)
    64136412               
     
    64166415                mult.SetValue(data[key]['Mult'])
    64176416                multSpin.SetValue(0)
    6418 #                ResetFlatBkg()
    6419                 wx.CallLater(100,UpdatePDFGrid,G2frame,data)
    64206417                wx.CallAfter(OnComputePDF,None)
    64216418                           
    64226419            def OnMult(invalid,value,tc):
    64236420                if invalid: return
    6424 #                ResetFlatBkg()
    6425                 wx.CallLater(100,UpdatePDFGrid,G2frame,data)
    64266421                wx.CallAfter(OnComputePDF,None)
    64276422               
  • trunk/GSASIIstrMath.py

    r3818 r3821  
    15221522        GSdata = GSdata[:,:,nxs,:]+np.swapaxes(MSmod,0,1)         #Natm,Nops,Ntau,Mxyz
    15231523        GSdata = SGData['MagMom'][nxs,:,nxs,nxs]*GSdata   #flip vectors according to spin flip * det(opM)
    1524        
    1525         Kdata = np.inner(GSdata,uAmat).T     #Cartesian unit vectors
     1524        mXYZ = np.array([[xyz[0] for xyz in list(G2spc.GenAtom(xyz,SGData,All=True,Move=True))] for xyz in (Xdata+dXdata).T])
     1525        mXYZ = np.array(np.inner(mXYZ,modQ)*ngl,dtype=int)
     1526        RGSdata = np.array([[np.roll(GSdata[i,j],mXYZ[i,j],0) for j in range(mXYZ.shape[1])] for i in range(mXYZ.shape[0])])
     1527        Kdata = np.inner(RGSdata,uAmat).T     #Cartesian unit vectors
    15261528        SMag = np.sqrt(np.sum(Kdata**2,axis=0))
    15271529        Kdata /= SMag      #mxyz,ntau,nops,natm
     
    16071609            fbm = (Q*TMcorr[nxs,:,nxs,nxs,:]*sinm[nxs,:,nxs,:,:]*SMag[nxs,nxs,:,:,:])
    16081610           
    1609             fas = np.sum(np.sum(fam,axis=-1),axis=-1)      #xyz,Nref,ntau; sum ops & atoms
    1610             fbs = np.sum(np.sum(fbm,axis=-1),axis=-1)      #ditto
     1611            fas = np.sum(np.sum(fam,axis=-1)**2,axis=-1)      #xyz,Nref,ntau; sum ops & atoms
     1612            fbs = np.sum(np.sum(fbm,axis=-1)**2,axis=-1)      #ditto
    16111613           
    1612             refl.T[10] = np.sum(np.sum(fas**2,axis=0)*glWt[nxs,nxs,:],axis=-1)+     \
    1613                 np.sum(np.sum(fbs**2,axis=0)*glWt[nxs,nxs,:],axis=-1)    #square of sums
    1614             refl.T[11] = atan2d(fbs[:,0],fas[:,0])  #ignore f' & f"
     1614            refl.T[10] = np.sum(np.sum(fas,axis=0)/ngl,axis=-1)+np.sum(np.sum(fbs,axis=0)/ngl,axis=-1)    #square of sums
     1615#            refl.T[11] = mphase[:,0,0]  #ignore f' & f"
    16151616           
    16161617        else:
Note: See TracChangeset for help on using the changeset viewer.