Changeset 1797


Ignore:
Timestamp:
Apr 18, 2015 9:18:31 AM (7 years ago)
Author:
vondreele
Message:

flag HKLF space group extinct reflections (mul=0) & not count them

Location:
trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASII.py

    r1792 r1797  
    699699                        ref[4+Super] = np.sqrt(1./G2lat.calc_rDsq2(H,G))
    700700                        iabsnt,ref[3+Super],Uniq,phi = G2spc.GenHKLf(H,SGData)
     701                        if iabsnt:  #flag space gp. absences
     702                            ref[3+Super] = 0
    701703                    UseList[histoName] = SetDefaultDData(reflData['Type'],histoName)
    702704                elif histoName in PWDRlist:
     
    841843                    ref[4+Super] = np.sqrt(1./G2lat.calc_rDsq2(H,G))
    842844                    iabsnt,ref[3+Super],Uniq,phi = G2spc.GenHKLf(H,SGData)
     845                    if iabsnt:  #flag space gp. absences
     846                        ref[3+Super] = 0
    843847        wx.EndBusyCursor()
    844848       
  • trunk/GSASIIpwdGUI.py

    r1787 r1797  
    30113011            for r in range(G2frame.refTable[phaseName].GetNumberRows()):
    30123012                if HKLF:
    3013                     if float(G2frame.refTable[phaseName].GetCellValue(r,3+im)) < 0.:
     3013                    if float(G2frame.refTable[phaseName].GetCellValue(r,3+im)) <= 0.:
    30143014                        G2frame.refTable[phaseName].SetCellBackgroundColour(r,3+im,wx.RED)
    30153015                    Fosq = float(G2frame.refTable[phaseName].GetCellValue(r,5+im))
     
    30303030        G2frame.dataFrame.SetLabel('Reflection List for '+phaseName)
    30313031        if HKLF:
    3032             Status.SetStatusText('abs(Fo-Fc)/sig > 10 in red; > 3 in yellow; mul < 0 (user rejected) in red')
     3032            Status.SetStatusText('abs(DF)/sig > 10 red; > 3 yellow; mul < 0 (user rejected) red; mul=0 (sp. gp. absent) red')
    30333033        else:
    30343034            Status.SetStatusText('Prfo < 0. in red')
  • trunk/GSASIIstrMath.py

    r1793 r1797  
    27262726                        ref[7+im] = parmDict[phfx+'Scale']*ref[9+im]*ref[11+im]  #correct Fc^2 for extinction
    27272727                        ref[8+im] = ref[5+im]/(parmDict[phfx+'Scale']*ref[11+im])
    2728                         if UserRejectHKL(ref,im,calcControls['UsrReject']):
     2728                        if UserRejectHKL(ref,im,calcControls['UsrReject']) and ref[3+im]:    #skip sp.gp. absences (mul=0)
    27292729                            ref[3+im] = abs(ref[3+im])      #mark as allowed
    27302730                            Fo = np.sqrt(ref[5+im])
     
    27482748                        Fc = np.sqrt(ref[7+im])
    27492749                        w = 2.0*(Fo/ref[6+im])**2    # 1/sig(F)?
    2750                         if UserRejectHKL(ref,im,calcControls['UsrReject']):
     2750                        if UserRejectHKL(ref,im,calcControls['UsrReject']) and ref[3+im]:    #skip sp.gp. absences (mul=0)
    27512751                            ref[3+im] = abs(ref[3+im])      #mark as allowed
    27522752                            sumFo += Fo
Note: See TracChangeset for help on using the changeset viewer.