Changeset 4130


Ignore:
Timestamp:
Sep 2, 2019 9:09:57 AM (2 years ago)
Author:
vondreele
Message:

a useful fix to incomm. mag. str.. fctr.. calc.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrMath.py

    r4129 r4130  
    15941594            sinm = np.sin(phasem)
    15951595            MF = refDict['FF']['MF'][iBeg:iFin].T[Tindx].T   #Nref,Natm
    1596             TMcorr = 0.539*(np.reshape(Tiso,Tuij.shape)*Tuij)[:,0,:]*Fdata*Mdata*MF/(2.*Nops)     #Nref,Natm
     1596            TMcorr = 0.539*(np.reshape(Tiso,Tuij.shape)*Tuij)[:,0,:]*Mdata*MF/(2.*Nops)     #Nref,Natm
    15971597#                     
    15981598            HM = np.inner(Bmat,HP.T)                            #put into cartesian space
     
    16191619            fbmqs = np.sum(np.sum(fbms,axis=-2),axis=-2)
    16201620           
    1621             fass = np.sum(famqs,axis=-1)**2-np.sum(eM.T[:,nxs,:]*famqs,axis=-1)**2      #mag intensity calc F^2-(e.F)^2
    1622             fbss = np.sum(fbmqs,axis=-1)**2-np.sum(eM.T[:,nxs,:]*fbmqs,axis=-1)**2
     1621            famcs = np.swapaxes(np.inner(Amat,famqs).T,0,1)     #convert to cartesian
     1622            fbmcs = np.swapaxes(np.inner(Amat,fbmqs).T,0,1)     # as Nref,Ntau,Mxyz
     1623           
     1624            famcs /= np.sqrt(np.sum(famcs**2,axis=-1))[:,:,nxs]     #normalize
     1625            fbmcs /= np.sqrt(np.sum(fbmcs**2,axis=-1))[:,:,nxs]
     1626           
     1627            famcs = np.nan_to_num(famcs)    #nan --> 0.0
     1628            fbmcs = np.nan_to_num(fbmcs)
     1629           
     1630            fass = np.sum(famqs**2,axis=-1)*(1.-np.sum(eM.T[:,nxs,:]*famcs,axis=-1)**2)      #mag intensity calc F^2-(e.F)^2
     1631            fbss = np.sum(fbmqs**2,axis=-1)*(1.-np.sum(eM.T[:,nxs,:]*fbmcs,axis=-1)**2)
    16231632           
    16241633            fas = np.sum(glWt*fass,axis=1)
Note: See TracChangeset for help on using the changeset viewer.