Ignore:
Timestamp:
May 7, 2018 8:08:45 AM (4 years ago)
Author:
vondreele
Message:

further mag derivatives - getting closer

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrMath.py

    r3372 r3374  
    11641164        dqdk = np.sum(HM[:,:,nxs,nxs]*deDotK-dG.T[:,nxs,:,:],axis=3)     #Nref,Nops,Mxyyz
    11651165        NQ = np.where(np.abs(Q)>0.,1./np.abs(Q),0.)     #this sort of works esp for 1 axis moments
    1166 #        dqdk = (HM*HM)-1.
    11671166        dqdm = dqdk[:,:,:,nxs]*dkdm[:,nxs,:,:]
    1168         dmx = Q*dMdm[:,nxs,nxs,:]
    1169         dmx += dqdm*Mag[nxs,nxs,:,:]
     1167        dmx = NQ*Q*dMdm[:,nxs,nxs,:]
     1168        dmx += NQ*dqdm*Mag[nxs,nxs,:,:]
    11701169       
    11711170        fam = Q*TMcorr[nxs,:,nxs,:]*cosm[nxs,:,:,:]*Mag[nxs,nxs,:,:]    #Mxyz,Nref,Nop,Natm
Note: See TracChangeset for help on using the changeset viewer.