Ignore:
Timestamp:
Aug 12, 2013 11:22:30 AM (9 years ago)
Author:
vondreele
Message:

fix excluded region problems

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrMath.py

    r1017 r1025  
    13161316                elif not iBeg-iFin:     #peak above high limit - done
    13171317                    break
    1318                 yc[iBeg:iFin] += refl[13]*refl[9]*G2pwd.getFCJVoigt3(refl[5],refl[6],refl[7],shl,x[iBeg:iFin])    #>90% of time spent here
     1318                   
     1319                print refl[5],ma.count(x[iBeg:iFin]),iFin-iBeg,refl[13],refl[9]
     1320                yc[iBeg:iFin] += refl[13]*refl[9]*G2pwd.getFCJVoigt3(refl[5],refl[6],refl[7],shl,ma.getdata(x[iBeg:iFin]))    #>90% of time spent here
    13191321                if Ka2:
    13201322                    pos2 = refl[5]+lamRatio*tand(refl[5]/2.0)       # + 360/pi * Dlam/lam * tan(th)
     
    13261328                    elif not iBeg-iFin:     #peak above high limit - done
    13271329                        return yc,yb
    1328                     yc[iBeg:iFin] += refl[13]*refl[9]*kRatio*G2pwd.getFCJVoigt3(pos2,refl[6],refl[7],shl,x[iBeg:iFin])        #and here
     1330                    yc[iBeg:iFin] += refl[13]*refl[9]*kRatio*G2pwd.getFCJVoigt3(pos2,refl[6],refl[7],shl,ma.getdata(x[iBeg:iFin]))        #and here
    13291331            elif 'T' in calcControls[hfx+'histType']:
    13301332                print 'TOF Undefined at present'
     
    14341436                lenBF = iFin-iBeg
    14351437                dMdpk = np.zeros(shape=(6,lenBF))
    1436                 dMdipk = G2pwd.getdFCJVoigt3(refl[5],refl[6],refl[7],shl,x[iBeg:iFin])
     1438                dMdipk = G2pwd.getdFCJVoigt3(refl[5],refl[6],refl[7],shl,ma.getdata(x[iBeg:iFin]))
    14371439                for i in range(5):
    14381440                    dMdpk[i] += 100.*cw[iBeg:iFin]*refl[13]*refl[9]*dMdipk[i]
     
    14451447                        lenBF2 = iFin2-iBeg2
    14461448                        dMdpk2 = np.zeros(shape=(6,lenBF2))
    1447                         dMdipk2 = G2pwd.getdFCJVoigt3(pos2,refl[6],refl[7],shl,x[iBeg2:iFin2])
     1449                        dMdipk2 = G2pwd.getdFCJVoigt3(pos2,refl[6],refl[7],shl,ma.getdata(x[iBeg2:iFin2]))
    14481450                        for i in range(5):
    14491451                            dMdpk2[i] = 100.*cw[iBeg2:iFin2]*refl[13]*refl[9]*kRatio*dMdipk2[i]
Note: See TracChangeset for help on using the changeset viewer.