Changeset 2078


Ignore:
Timestamp:
Dec 3, 2015 9:32:46 AM (6 years ago)
Author:
vondreele
Message:

fix initial refinement of HKLF Scale problem; now estimates Scale if large discrepancy.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrMath.py

    r2076 r2078  
    29872987            sumFo = 0
    29882988            sumFo2 = 0
     2989            sumFc2 = 0
    29892990            sumdF = 0
    29902991            sumdF2 = 0
     
    30133014                            sumFo += Fo
    30143015                            sumFo2 += ref[5+im]
     3016                            sumFc2 += ref[7+im]
    30153017                            sumdF += abs(Fo-np.sqrt(ref[7+im]))
    30163018                            sumdF2 += abs(ref[5+im]-ref[7+im])
     
    30473049                            sumFo += Fo
    30483050                            sumFo2 += ref[5+im]
     3051                            sumFc2 += ref[7+im]
    30493052                            sumdF += abs(Fo-Fc)
    30503053                            sumdF2 += abs(ref[5+im]-ref[7+im])
     
    30683071                            else:   #sp.gp.extinct
    30693072                                next += 1
     3073            Scale = sumFo2/sumFc2
     3074            if (Scale < 0.8 or Scale > 1.2) and phfx+'Scale' in varylist:
     3075                print 'New scale:',Scale*parmDict[phfx+'Scale']
     3076                indx = varylist.index(phfx+'Scale')
     3077                values[indx] = Scale*parmDict[phfx+'Scale']             
    30703078            Histogram['Residuals']['Nobs'] = nobs
    30713079            Histogram['Residuals']['sumwYo'] = sumwYo
Note: See TracChangeset for help on using the changeset viewer.