Ignore:
Timestamp:
Sep 23, 2015 3:28:12 PM (7 years ago)
Author:
vondreele
Message:

More modulation structure factor stuff, begin SSF derivatives

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrMath.py

    r1976 r1977  
    983983    waveTypes,FSSdata,XSSdata,USSdata,MSSdata = GetAtomSSFXU(pfx,calcControls,parmDict)
    984984    SStauM = list(GetSSTauM(SGData['SGOps'],SSGData['SSGOps'],pfx,calcControls,Xdata))
    985     SST = SSGT[:,3]
    986985    if SGData['SGInv']:
    987986        SStauM[0] = np.hstack((SStauM[0],SStauM[0]))
    988987        SStauM[1] = np.hstack((SStauM[1],SStauM[1]))
    989         SST = np.hstack((SST,-SST))
    990988    modQ = np.array([parmDict[pfx+'mV0'],parmDict[pfx+'mV1'],parmDict[pfx+'mV2']])
    991989    FF = np.zeros(len(Tdata))
     
    10361034            SSPhi = np.hstack((SSPhi,-SSPhi))
    10371035#        GSASIIpath.IPyBreak()
    1038         GfpuA = G2mth.Modulation(waveTypes,SSUniq,SST,FSSdata,XSSdata,USSdata,Mast)
    1039         phase = twopi*(np.inner(Uniq,(dXdata.T+Xdata.T))+Phi[:,np.newaxis])
     1036        GfpuA = G2mth.Modulation(waveTypes,SSUniq,FSSdata,XSSdata,USSdata,Mast)
     1037        phase = twopi*(np.inner(Uniq,(dXdata.T+Xdata.T))+SSPhi[:,np.newaxis])
    10401038        sinp = np.sin(phase)
    10411039        cosp = np.cos(phase)
     
    11031101        Phi = np.inner(H[:3],SGT)
    11041102        SSPhi = np.inner(H,SSGT)
    1105         GfpuA,GfpuB = G2mth.Modulation(waveTypes,SSUniq,SSPhi,FSSdata,XSSdata,USSdata,SStauM,Mast)
    1106         dGAdk,dGBdk = G2mth.ModulationDerv(waveTypes,SSUniq,SSPhi,FSSdata,XSSdata,USSdata,SStauM,Mast)
     1103        GfpuA = G2mth.Modulation(waveTypes,SSUniq,FSSdata,XSSdata,USSdata,Mast)
     1104        dGAdk = G2mth.ModulationDerv(waveTypes,SSUniq,FSSdata,XSSdata,USSdata,Mast)
    11071105        #need ModulationDerv here dGAdXsin, etc 
    1108         phase = twopi*(np.inner((dXdata.T+Xdata.T),Uniq)+Phi[np.newaxis,:])
     1106        phase = twopi*(np.inner((dXdata.T+Xdata.T),Uniq)+SSPhi[np.newaxis,:])
    11091107        sinp = np.sin(phase)
    11101108        cosp = np.cos(phase)
     
    11191117        fot = (FF+FP-Bab)*occ*Tcorr
    11201118        fotp = FPP*occ*Tcorr
    1121         fa = np.array([fot[:,np.newaxis]*cosp,fotp[:,np.newaxis]*cosp])       #non positions
    1122         fb = np.array([fot[:,np.newaxis]*sinp,-fotp[:,np.newaxis]*sinp])
    1123         fa = fa*GfpuA[:,:,np.newaxis]-fb*GfpuB[:,:,np.newaxis]
    1124         fb = fb*GfpuA[:,:,np.newaxis]+fa*GfpuB[:,:,np.newaxis]
     1119        fa *= GfpuA
     1120        fb *= GfpuA       
    11251121       
    11261122        fas = np.sum(np.sum(fa,axis=1),axis=1)
Note: See TracChangeset for help on using the changeset viewer.