Changeset 730


Ignore:
Timestamp:
Aug 22, 2012 6:52:30 PM (10 years ago)
Author:
vondreele
Message:

a try to fix the neg pawley peak problem

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstruct.py

    r709 r730  
    26232623                        pInd =pfx+'PWLref:%d'%(pawleyLookup[pfx+'%d,%d,%d'%(h,k,l)])
    26242624                        parmDict[pInd] = max(parmDict[pInd]/2.,parmDict[pInd])       
    2625                         refl[9] = parmDict[pInd]
     2625                        refl[9] = abs(parmDict[pInd])
    26262626                    except KeyError:
    26272627#                        print ' ***Error %d,%d,%d missing from Pawley reflection list ***'%(h,k,l)
     
    27732773                        dMdpw[iBeg:iFin] = dervDict['int']/refl[9]
    27742774                        if parmDict[pIdx] < 0.:
    2775                             dMdpw[iBeg:iFin] = 2.*dervDict['int']/refl[9]
     2775                            dMdpw[iBeg:iFin] = -dervDict['int']/refl[9]
    27762776                        if Ka2:
    27772777                            dMdpw[iBeg2:iFin2] += dervDict2['int']/refl[9]
    27782778                            if parmDict[pIdx] < 0.:
    2779                                 dMdpw[iBeg2:iFin2] += 2.*dervDict['int']/refl[9]
     2779                                dMdpw[iBeg2:iFin2] -= dervDict['int']/refl[9]
    27802780                        dMdv[idx] = dMdpw
    27812781                    except: # ValueError:
Note: See TracChangeset for help on using the changeset viewer.