Changeset 738 for trunk/GSASIImath.py


Ignore:
Timestamp:
Aug 26, 2012 12:54:17 PM (11 years ago)
Author:
vondreele
Message:

put wx.BusyCursor? in various places

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIImath.py

    r736 r738  
    938938        rho[rMM[0]:rMP[0],rMM[1]:rMP[1],rMM[2]:rMP[2]] = peakFunc(x1,rX,rY,rZ,rhoPeak,res,SGData['SGLaue'])
    939939        rho = np.roll(np.roll(np.roll(rho,-rMI[2],axis=2),-rMI[1],axis=1),-rMI[0],axis=0)
    940     if SGData['SGInv']:                 #check origin location & fix it if needed - centrosymmetric only
    941         Ocheck = np.zeros_like(rho)
    942         for ipeak in peaks:
    943             for opeak in peaks:           
    944                 dx = ((opeak-ipeak)*incre)%incre
    945                 if np.any(dx):      #avoid self vector
    946                     Ocheck[dx[0],dx[1],dx[2]] += 1
    947         dxMax = np.unravel_index(np.argmax(Ocheck),rho.shape)
    948         print 'Inversion at:',dxMax,' shift by ;',dxMax-mapHalf,' map size:',rho.shape
    949         rho = np.roll(np.roll(np.roll(rho,dxMax[2],axis=2),dxMax[1],axis=1),dxMax[0],axis=0)
    950         for peak in peaks:
    951             peak = (peak-(dxMax+mapHalf)/incre)%1.0
     940#    if SGData['SGInv']:                 #check origin location & fix it if needed - centrosymmetric only
     941#        Ocheck = np.zeros_like(rho)
     942#        for ipeak in peaks:
     943#            for opeak in peaks:           
     944#                dx = ((opeak-ipeak)*incre)%incre
     945#                if np.any(dx):      #avoid self vector
     946#                    Ocheck[dx[0],dx[1],dx[2]] += 1
     947#        dxMax = np.unravel_index(np.argmax(Ocheck),rho.shape)
     948#        print 'Inversion at:',dxMax,' shift by ;',dxMax-mapHalf,' map size:',rho.shape
     949#        rho = np.roll(np.roll(np.roll(rho,dxMax[2],axis=2),dxMax[1],axis=1),dxMax[0],axis=0)
     950#        for peak in peaks:
     951#            peak = (peak-(dxMax+mapHalf)/incre)%1.0
    952952       
    953953    return np.array(peaks),np.array([mags,]).T,np.array([dzeros,]).T
Note: See TracChangeset for help on using the changeset viewer.