Changeset 378


Ignore:
Timestamp:
Sep 20, 2011 3:33:06 PM (10 years ago)
Author:
vondreele
Message:

a Pawley intensity fix, format fixes & Sph. Harm PO correction

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIlattice.py

    r375 r378  
    718718    for iord in [2*i+2 for i in range(L/2)]:
    719719        for m in [i-iord for i in range(2*iord+1)]:
    720             if OdfChk(SamSym,iord,m):
     720            if SamSym and OdfChk(SamSym,iord,m):
    721721                for n in [i-iord for i in range(2*iord+1)]:
    722722                    if OdfChk(SGLaue,iord,n):
    723723                        coeffNames.append('C(%d,%d,%d)'%(iord,m,n))
    724 #            else:                  #what's this for?
    725 #                for n in [i-iord for i in range(2*iord+1)]:
    726 #                    if OdfChk(SGLaue,iord,n):
    727 #                        coeffNames.append('C(%d,%d)'%(iord,n))
     724            else:                  #use for powder sample PO when SamSym = None
     725                for n in [i-iord for i in range(2*iord+1)]:
     726                    if OdfChk(SGLaue,iord,n):
     727                        coeffNames.append('C(%d,%d)'%(iord,n))
    728728    return coeffNames
    729729
  • trunk/GSASIIstruct.py

    r377 r378  
    503503        line = '\n Size model    : '+hapData[0]
    504504        if hapData[0] in ['isotropic','uniaxial']:
    505             line += ' equatorial:'+'%12.2f'%(hapData[1][0])+' Refine? '+str(hapData[2][0])
     505            line += ' equatorial:'+'%12.3f'%(hapData[1][0])+' Refine? '+str(hapData[2][0])
    506506            if hapData[0] == 'uniaxial':
    507                 line += ' axial:'+'%12.2f'%(hapData[1][1])+' Refine? '+str(hapData[2][1])
     507                line += ' axial:'+'%12.3f'%(hapData[1][1])+' Refine? '+str(hapData[2][1])
    508508            print line
    509509        else:
     
    524524        line = '\n Mustrain model: '+hapData[0]
    525525        if hapData[0] in ['isotropic','uniaxial']:
    526             line += ' equatorial:'+'%12.4f'%(hapData[1][0])+' Refine? '+str(hapData[2][0])
     526            line += ' equatorial:'+'%12.1f'%(hapData[1][0])+' Refine? '+str(hapData[2][0])
    527527            if hapData[0] == 'uniaxial':
    528                 line += ' axial:'+'%12.4f'%(hapData[1][1])+' Refine? '+str(hapData[2][1])
     528                line += ' axial:'+'%12.1f'%(hapData[1][1])+' Refine? '+str(hapData[2][1])
    529529            print line
    530530        else:
     
    661661   
    662662    def PrintSizeAndSig(hapData,sizeSig):
    663         line = '\n Size model: '+hapData[0]
     663        line = '\n Size model:     '+hapData[0]
    664664        if hapData[0] in ['isotropic','uniaxial']:
    665             line += ' equatorial:%12.2f'%(hapData[1][0])
     665            line += ' equatorial:%12.3f'%(hapData[1][0])
    666666            if sizeSig[0][0]:
    667                 line += ', sig: %8.2f'%(sizeSig[0][0])
     667                line += ', sig: %8.3f'%(sizeSig[0][0])
    668668            if hapData[0] == 'uniaxial':
    669                 line += ' axial:%12.2f'%(hapData[1][1])
     669                line += ' axial:%12.3f'%(hapData[1][1])
    670670                if sizeSig[0][1]:
    671                     line += ', sig: %8.2f'%(sizeSig[0][1])
     671                    line += ', sig: %8.3f'%(sizeSig[0][1])
    672672            print line
    673673        else:
     
    691691        line = '\n Mustrain model: '+hapData[0]
    692692        if hapData[0] in ['isotropic','uniaxial']:
    693             line += ' equatorial:%12.4f'%(hapData[1][0])
     693            line += ' equatorial:%12.1f'%(hapData[1][0])
    694694            if mustrainSig[0][0]:
    695                 line += ',sig: %12.4f'%(mustrainSig[0][0])
     695                line += ', sig: %8.1f'%(mustrainSig[0][0])
    696696            if hapData[0] == 'uniaxial':
    697                 line += ' axial:%12.4f'%(hapData[1][1])
     697                line += ' axial:%12.1f'%(hapData[1][1])
    698698                if mustrainSig[0][1]:
    699                      line += ',sig: %12.4f'%(mustrainSig[0][1])
     699                     line += ', sig: %8.1f'%(mustrainSig[0][1])
    700700            print line
    701701        else:
     
    748748                hapData[item][0] = parmDict[pfx+item]
    749749                if pfx+item in sigDict:
    750                     PhFrExtPOSig[item] = sigDict[pfx+item]           
     750                    PhFrExtPOSig[item] = sigDict[pfx+item]
    751751            if hapData['Pref.Ori.'][0] == 'MD':
    752752                hapData['Pref.Ori.'][1] = parmDict[pfx+'MD']
    753                 if pfx+item in sigDict:
    754                     PhFrExtPOSig[item] = sigDict[pfx+item]
     753                if pfx+'MD' in sigDict:
     754                    PhFrExtPOSig['MD'] = sigDict[pfx+'MD']
    755755            else:                           #'SH' spherical harmonics
    756756                for item in hapData['Pref.Ori.'][5]:
    757                     hapData['Pref.Ori.'][5][item] = parmDict[pfx+item]
     757                    hapData['Pref.Ori.'][5][item] = parmDict[pfx+'MD']
    758758                    if pfx+item in sigDict:
    759759                        PhFrExtPOSig[item] = sigDict[pfx+item]
    760 #            print '\n Phase fraction  : %10.4f, sig %10.4f'%(hapData['Scale'][0],PhFrExtPOSig['Scale'])
    761 #            print ' Extinction coeff: %10.4f, sig %10.4f'%(hapData['Extinction'][0],PhFrExtPOSig['Extinction'])
    762 #            if hapData['Pref.Ori.'][0] == 'MD':
    763 #                Ax = hapData['Pref.Ori.'][3]
    764 #                print ' March-Dollase PO: %10.4f'%(hapData['Pref.Ori.'][1]),' Refine?',hapData['Pref.Ori.'][2], \
    765 #                    ' Axis: %d %d %d'%(Ax[0],Ax[1],Ax[2])
    766                
     760            print '\n'
     761            if 'Scale' in PhFrExtPOSig:
     762                print ' Phase fraction  : %10.4f, sig %10.4f'%(hapData['Scale'][0],PhFrExtPOSig['Scale'])
     763            if 'Extinction' in PhFrExtPOSig:
     764                print ' Extinction coeff: %10.4f, sig %10.4f'%(hapData['Extinction'][0],PhFrExtPOSig['Extinction'])
     765            if 'MD' in PhFrExtPOSig:
     766                print ' March-Dollase PO: %10.4f, sig %10.4f'%(hapData['Pref.Ori.'][1],PhFrExtPOSig['MD'])
     767             
    767768            SizeMuStrSig = {'Mustrain':[[0,0],[0 for i in range(len(hapData['Mustrain'][4]))]],
    768769                'Size':[[0,0],[0 for i in range(len(hapData['Size'][4]))]],
     
    15261527                try:
    15271528                    idx = varylist.index(pfx+'PWLref:'+str(iref))
    1528                     dMdv[idx] = dervDict['int']
     1529                    dMdv[idx] = dervDict['int']/refl[8]
    15291530                except ValueError:
    15301531                    pass
Note: See TracChangeset for help on using the changeset viewer.