Changeset 1597 for trunk/GSASIIstrIO.py


Ignore:
Timestamp:
Dec 6, 2014 3:12:32 PM (7 years ago)
Author:
vondreele
Message:

Incommensurate Pawley refinement works, including refinement of modulation vector.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrIO.py

    r1596 r1597  
    562562        print >>pFile,' Process histograms in reverse order: ',Controls['Reverse Seq']
    563563   
    564 def GetPawleyConstr(SGLaue,PawleyRef,pawleyVary):
     564def GetPawleyConstr(SGLaue,PawleyRef,im,pawleyVary):
    565565    'needs a doc string'
    566566#    if SGLaue in ['-1','2/m','mmm']:
     
    571571        refI = int(varyI.split(':')[-1])
    572572        ih,ik,il = PawleyRef[refI][:3]
    573         dspI = PawleyRef[refI][4]
     573        dspI = PawleyRef[refI][4+im]
    574574        for varyJ in pawleyVary[i+1:]:
    575575            refJ = int(varyJ.split(':')[-1])
    576576            jh,jk,jl = PawleyRef[refJ][:3]
    577             dspJ = PawleyRef[refJ][4]
     577            dspJ = PawleyRef[refJ][4+im]
    578578            if SGLaue in ['4/m','4/mmm']:
    579579                isum = ih**2+ik**2
     
    631631def modVary(pfx,SSGData):
    632632    vary = []
    633     for i,item in SSGData['modSymb']:
     633    for i,item in enumerate(SSGData['modSymb']):
    634634        if item in ['a','b','g']:
    635635            vary.append(pfx+'mV%d'%(i))
     
    11481148                if refl[5+im]:
    11491149                    pawleyVary.append(pfx+'PWLref:'+str(i))
    1150             GetPawleyConstr(SGData['SGLaue'],PawleyRef,pawleyVary)      #does G2mv.StoreEquivalence
     1150            GetPawleyConstr(SGData['SGLaue'],PawleyRef,im,pawleyVary)      #does G2mv.StoreEquivalence
    11511151            phaseVary += pawleyVary
    11521152               
     
    16071607                for var in ['mV0','mV1','mV2']:
    16081608                    namstr += '%12s'%(pfx+var)
    1609                     ptstr += '%12.4f'%(parmDict[pfx+var])
     1609                    ptstr += '%12.6f'%(parmDict[pfx+var])
    16101610                    if pfx+var in sigDict:
    1611                         sigstr += '%12.4f'%(sigDict[pfx+var])
     1611                        sigstr += '%12.6f'%(sigDict[pfx+var])
    16121612                    else:
    16131613                        sigstr += 12*' '
Note: See TracChangeset for help on using the changeset viewer.