Changeset 656


Ignore:
Timestamp:
Jun 26, 2012 9:55:52 AM (9 years ago)
Author:
vondreele
Message:

fix Debye background model

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIpwd.py

    r655 r656  
    752752            cqr = np.cos(q*dbR)
    753753            temp = np.exp(-dbU*q**2)
    754             dyddb[3*iD] = 100.*dx*ff*sqr*temp
    755             dyddb[3*iD+1] = 100.*dx*ff*dbA*temp*(cqr-sqr)/dbR
    756             dyddb[3*iD+2] = -100.*dx*ff*dbA*sqr*temp*q**2
    757             iD += 1       
     754            dyddb[3*iD] = ff*sqr*temp/(np.pi*dx)
     755            dyddb[3*iD+1] = ff*dbA*temp*(cqr-sqr)/(np.pi*dbR*dx)
     756            dyddb[3*iD+2] = -ff*dbA*sqr*temp*q**2/(np.pi*dx)
     757            iD += 1       #ff*dbA*np.sin(q*dbR)*np.exp(-dbU*q**2)/(q*dbR)
    758758        except KeyError:
    759759            break
  • trunk/GSASIIstruct.py

    r655 r656  
    15831583                line += '%8s refine?'%(parm)
    15841584            print line
    1585             for j,term in enumerate(Debye['debyeTerms']):
     1585            for j,term in enumerate(DebyePeaks['debyeTerms']):
    15861586                line = ' term'+'%2d'%(j)+':'
    15871587                for i in range(3):
    1588                     line += '%10.4g %5s'%(term[2*i],bool(term[2*i+1]))                   
     1588                    line += '%10.3f %5s'%(term[2*i],bool(term[2*i+1]))                   
    15891589                print line
    15901590        if DebyePeaks['nPeaks']:
     
    17701770        if DebyePeaks['nDebye']:
    17711771            ifAny = False
    1772             ptfmt = "%12.5f"
     1772            ptfmt = "%12.3f"
    17731773            names =  ' names :'
    17741774            ptstr =  ' values:'
Note: See TracChangeset for help on using the changeset viewer.