Changeset 712


Ignore:
Timestamp:
Aug 16, 2012 9:26:36 AM (10 years ago)
Author:
vondreele
Message:

fis Mac crash after PWDR delete

Location:
trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASII.py

    r709 r712  
    16731673                    for item in DelItemList:
    16741674                        self.PatternTree.Delete(item)
    1675                     G2plt.PlotPatterns(self,True)                        #so plot gets updated
     1675                    self.PickId = 0
     1676                    wx.CallAfter(G2plt.PlotPatterns,self,True)                        #so plot gets updated
    16761677            finally:
    16771678                dlg.Destroy()
  • trunk/GSASIIplot.py

    r711 r712  
    627627    Ymax = 1.0
    628628    lenX = 0
    629     if G2frame.PatternTree.GetItemText(PickId) in ['Reflection Lists']:
    630         Phases = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame,PatternId,'Reflection Lists'))
    631         HKL = []
    632         if Phases:
    633             for peak in Phases[G2frame.RefList]:
    634                 HKL.append(peak[:6])
    635             HKL = np.array(HKL)
    636     else:
    637         HKL = np.array(G2frame.HKL)
     629    if PickId:
     630        if G2frame.PatternTree.GetItemText(PickId) in ['Reflection Lists']:
     631            Phases = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame,PatternId,'Reflection Lists'))
     632            HKL = []
     633            if Phases:
     634                for peak in Phases[G2frame.RefList]:
     635                    HKL.append(peak[:6])
     636                HKL = np.array(HKL)
     637        else:
     638            HKL = np.array(G2frame.HKL)
    638639    for Pattern in PlotList:
    639640        xye = Pattern[1]
  • trunk/fsource/powsubs/psvoigt.for

    r210 r712  
    2525      REAL*4        ACOFG(7),ACOFL(7)   
    2626      REAL*4        GNORM               !Gaussian Normalization constant
    27       REAL*4        COFT(6),COFN(3)    
     27      REAL*4        COFT(6),COFN(3)
    2828
    2929!SUBROUTINES CALLED:
     
    3737      DATA COFT/1.0,2.69269,2.42843,4.47163,0.07842,1.0/
    3838      DATA COFN/1.36603,-0.47719,0.11116/
    39       DATA ITYPE/1/
    4039
    4140!CODE:
     
    8483      RETURN
    8584      END
     85
     86      SUBROUTINE PSVOIGT2(DX,SIG,GAM,FUNC,DFDX,DFDS,DFDG)
     87
     88!PURPOSE: Compute function & derivatives pseudovoigt
     89!pseudo Voigt W.I.F. David, J. Appl. Cryst. 19, 63-64 (1986)
     90
     91      INCLUDE       '../INCLDS/COPYRIGT.FOR'
     92
     93!PSEUDOCODE:
     94
     95!CALLING ARGUMENTS:
     96
     97      REAL*4        DX                  !Delta-x from center
     98      REAL*4        SIG                 !Gaussian variance
     99      REAL*4        GAM                 !Lorentzian FWHM
     100      REAL*4        FUNC                !Value of pseudo-Voigt at DX
     101      REAL*4        DFDX                !dF/dx
     102      REAL*4        DFDS                !dF/ds
     103      REAL*4        DFDG                !dF/dg
     104
     105!INCLUDE STATEMENTS:
     106
     107!LOCAL VARIABLES:
     108
     109      REAL*4        GNORM               !Gaussian Normalization constant
     110      REAL*4        COFEG(7),COFEL(7),COFGG(6),COFGL(6)   
     111
     112!SUBROUTINES CALLED:
     113
     114!FUNCTION DEFINITIONS:
     115
     116!DATA STATEMENTS:
     117
     118      DATA STOFW/2.35482005/            !2*SQRT(2LN2)
     119      DATA SQ2PI/2.506628275/            !SQRT(2PI)
     120      DATA ITYPE/1/
     121      DATA COFEG/0.00268,0.75458,2.88898,-3.85144,-.55765,3.03824,
     122     1  -1.27539/
     123      DATA COFEL/1.35248,0.41168,-2.18731,6.42452,-10.29036,6.88093,
     124     1  -1.59194/
     125      DATA COFGG/-.50734,-.22744,1.63804,-2.28532,1.31943,0.0/
     126      DATA COFGL/-.99725,1.14594,2.56150,-6.52088,5.82647,-1.91086/
     127
     128!CODE:
     129
     130      SQSG = MAX(SQRT(SIG),0.001)
     131      FWHG = STOFW*SQSG
     132      FW = FWHG+GAM
     133      R1 = FWHG/FW
     134      R17 = R1
     135      R2 = GAM/FWHG
     136      R27 = R2
     137      DSDL = 0.0
     138      DSDG = 0.0
     139      ETAG = 0.0
     140      ETAL = 0.0
     141      DO ITRM=1,7
     142        ETAG = ETAG+R17*COFEG(ITRM)
     143        ETAL = ETAL+R27*COFEL(ITRM)
     144        R17 = R17*R1
     145        R27 = R27*R2
     146      END DO
     147      WG = 1.0
     148      WL = 1.0
     149      R16 = R1
     150      R26 = R2
     151      DO ITRM=1,6
     152        WG = WG+R26*COFGG(ITRM)
     153        WL = WL+R16*COFGL(ITRM)
     154        R16 = R16*R1
     155        R26 = R26*R2
     156      END DO
     157      CALL LORENTZ(DX,WL,TL,DTLDT,DTLDFW)
     158      SIGP = (WG/STOFW)**2
     159      EX = MAX(-20.0,-0.5*DX**2/SIGP)
     160      TG = STOFW*EXP(EX)/(SQ2PI*FWHM)
     161      FUNC = ETAL*TL+ETAG*TG
     162
     163      RETURN
     164      END
Note: See TracChangeset for help on using the changeset viewer.