Ignore:
Timestamp:
Aug 5, 2021 8:52:37 AM (6 months ago)
Author:
vondreele
Message:

add 'det2thet' to image importers data dict.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrMath.py

    r5001 r5009  
    16181618            MF = refDict['FF']['MF'][iBeg:iFin].T[Tindx].T   #Nref,Natm
    16191619            TMcorr = 0.539*(np.reshape(Tiso,Tuij.shape)*Tuij)[:,0,:]*Mdata*Fdata*MF/(2*Nops)     #Nref,Natm
    1620             HM = np.inner(Bmat,HP.T)                #put into cartesian space X||H,Z||H*L;
     1620            HM = np.inner(Bmat.T,HP.T)                #put into cartesian space X||H,Z||H*L;
    16211621            eM = (HM*refl.T[5]).T                   # normalize HP by d*    Nref,hkl=Unit vectors || Q
    16221622
     
    16271627            fams = TMcorr[:,nxs,nxs,:,nxs]*SGData['MagMom'][nxs,nxs,:,nxs,nxs]*np.array([np.where(H[3,i]!=0,(
    16281628                (MmodAR+H[3,i]*MmodBR)*cosm[i,nxs,:,:,nxs]+   
    1629                 GamI[nxs,:,nxs,nxs]*(-MmodAI+H[3,i]*MmodBI)*sinm[i,nxs,:,:,nxs]),
     1629                GamI[nxs,:,nxs,nxs]*(MmodAI-H[3,i]*MmodBI)*sinm[i,nxs,:,:,nxs]),
    16301630                0.) for i in range(mRef)])/2.          #Nref,Ntau,Nops,Natm,Mxyz
    16311631                       
    16321632            fbms = TMcorr[:,nxs,nxs,:,nxs]*SGData['MagMom'][nxs,nxs,:,nxs,nxs]*np.array([np.where(H[3,i]!=0,(
    16331633                (MmodAR+H[3,i]*MmodBR)*sinm[i,nxs,:,:,nxs]+   
    1634                 GamI[nxs,:,nxs,nxs]*(MmodAI-H[3,i]*MmodBI)*cosm[i,nxs,:,:,nxs]),
     1634                GamI[nxs,:,nxs,nxs]*(-MmodAI+H[3,i]*MmodBI)*cosm[i,nxs,:,:,nxs]),
    16351635                0.) for i in range(mRef)])/2.          #Nref,Ntau,Nops,Natm,Mxyz
    16361636           
     
    16401640                               
    16411641#sum ops & atms                               
    1642             fasm = np.sum(np.sum(fams,axis=-2),axis=-2)    #Nref,Mxyz; sum ops & atoms
     1642            fasm = np.sum(np.sum(fams,axis=-2),axis=-2)    #Nref,Ntau,Mxyz; sum ops & atoms
    16431643            fbsm = np.sum(np.sum(fbms,axis=-2),axis=-2)
    16441644# #put into cartesian space
    1645             facm = np.inner(fasm,uBmat)       #uBmat best fit for DyMnGe
    1646             fbcm = np.inner(fbsm,uBmat)
     1645            facm = np.inner(fasm,uAmat)       #uBmat best fit for DyMnGe
     1646            fbcm = np.inner(fbsm,uAmat)
    16471647#form e.F dot product
    16481648            eDotFa = np.sum(eM[:,nxs,:]*facm,axis=-1)    #Nref,Ntau       
Note: See TracChangeset for help on using the changeset viewer.