Changeset 4147


Ignore:
Timestamp:
Sep 13, 2019 3:19:10 PM (2 years ago)
Author:
vondreele
Message:

mag. stuff update

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrMath.py

    r4145 r4147  
    16241624                fbms += fbm0[:,nxs,:,:,:]
    16251625# do sum on ops, atms 1st                       
    1626             cos2 = np.sqrt(1.-eDotK**2)
     1626            sinsq = np.sqrt(1.-eDotK**2)        #projection - Nref,Ntau,Nops,Natm
    16271627           
    1628             famqs = np.sum(np.sum(fams*cos2[:,:,:,:,nxs],axis=-2),axis=-2)      #Nref,Ntau,Mxyz; sum ops & atoms
    1629             fbmqs = np.sum(np.sum(fbms*cos2[:,:,:,:,nxs],axis=-2),axis=-2)
     1628            famqs = np.sum(np.sum(fams*sinsq[:,:,:,:,nxs],axis=-2),axis=-2)      #Nref,Ntau,Mxyz; sum ops & atoms
     1629            fbmqs = np.sum(np.sum(fbms*sinsq[:,:,:,:,nxs],axis=-2),axis=-2)
    16301630           
    16311631            fass = np.sum(famqs**2,axis=-1)      #mag intensity calc F^2-(e.F)^2
Note: See TracChangeset for help on using the changeset viewer.