Changeset 4051


Ignore:
Timestamp:
Jul 2, 2019 3:22:24 PM (3 years ago)
Author:
vondreele
Message:

Correct error in PWDR extinction derivatives

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrMath.py

    r4050 r4051  
    25712571    if 0 < xfac <= 1.:
    25722572        xn = np.array([i*flv2*xfac**i for i in [1,2,3,4,5,6]])
    2573         dlde = np.sum(xn*coef)
     2573        dlde = np.sum(xn*coef)/xfac
    25742574    elif xfac > 1.:
    25752575        xfac2 = 1./np.sqrt(xfac)
    2576         dlde = flv2*pi2*xfac2*(-1./xfac+0.375/xfac**2)
     2576        dlde = 0.5*flv2*pi2*xfac2*(-1./xfac+0.375/xfac**2)
    25772577       
    25782578    return dbde*sth2+dlde*(1.-sth2)
Note: See TracChangeset for help on using the changeset viewer.