Changeset 2718


Ignore:
Timestamp:
Feb 21, 2017 3:51:44 PM (5 years ago)
Author:
vondreele
Message:

implement LaBail? - not finished.

Location:
trunk
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIddataGUI.py

    r2574 r2718  
    816816            Obj = event.GetEventObject()
    817817            UseList[G2frame.hist]['Use'] = Obj.GetValue()
    818        
     818
     819        def OnLeBail(event):
     820            Obj = event.GetEventObject()
     821            UseList[G2frame.hist]['LeBail'] = Obj.GetValue()
     822
    819823        def OnResetSize(event):
    820824            Obj = event.GetEventObject()
     
    853857        if 'Use' not in UseList[G2frame.hist]:      #patch
    854858            UseList[G2frame.hist]['Use'] = True
     859        if 'LeBail' not in UseList[G2frame.hist]:
     860            UseList[G2frame.hist]['LeBail'] = False
    855861        if 'Babinet' not in UseList[G2frame.hist]:
    856862            UseList[G2frame.hist]['Babinet'] = {'BabA':[0.0,False],'BabU':[0.0,False]}
    857863        bottomSizer = wx.BoxSizer(wx.VERTICAL)
     864        useBox = wx.BoxSizer(wx.HORIZONTAL)
    858865        useData = wx.CheckBox(DData,wx.ID_ANY,label='Use Histogram: '+G2frame.hist+' ?')
    859866        useData.Bind(wx.EVT_CHECKBOX, OnUseData)
    860867        useData.SetValue(UseList[G2frame.hist]['Use'])
    861         bottomSizer.Add(useData,0,WACV|wx.TOP|wx.BOTTOM|wx.LEFT,5)
     868        useBox.Add(useData,0,WACV)
     869        lebail = wx.CheckBox(DData,wx.ID_ANY,label='Do LeBail extraction?')
     870        lebail.Bind(wx.EVT_CHECKBOX, OnLeBail)
     871        lebail.SetValue(UseList[G2frame.hist]['LeBail'])
     872        useBox.Add(lebail,0,WACV)
     873        bottomSizer.Add(useBox,0,WACV|wx.TOP|wx.BOTTOM|wx.LEFT,5)
    862874       
    863875        bottomSizer.Add(ScaleSizer(),0,WACV|wx.BOTTOM,5)
  • trunk/GSASIIphsGUI.py

    r2693 r2718  
    206206            generalData['Flip']['testHKL'] = [[0,0,2],[2,0,0],[1,1,1],[0,2,0],[1,2,3]]
    207207        if 'doPawley' not in generalData:
    208             generalData['doPawley'] = False
     208            generalData['doPawley'] = False     #ToDo: change to ''
    209209        if 'Pawley dmin' not in generalData:
    210210            generalData['Pawley dmin'] = 1.0
     
    828828            pawleySizer.Add(wx.StaticText(General,label=' Pawley controls: '),0,WACV)
    829829            pawlRef = wx.CheckBox(General,-1,label=' Do Pawley refinement?')
     830            #ToDo: change parameter to ComboBox of blank, Pawley, LeBail
    830831            pawlRef.SetValue(generalData['doPawley'])
    831832            pawlRef.Bind(wx.EVT_CHECKBOX,OnPawleyRef)
     
    55505551            refDict,reflData = G2frame.PatternTree.GetItemPyData(Id)
    55515552            UseList[histoName] = {'Histogram':histoName,'Show':False,'Scale':[1.0,True],
    5552                 'Babinet':{'BabA':[0.0,False],'BabU':[0.0,False]},
     5553                'Babinet':{'BabA':[0.0,False],'BabU':[0.0,False]},'LeBail':False,
    55535554                'Extinction':['Lorentzian','None',
    55545555                {'Tbar':0.1,'Cos2TM':0.955,'Eg':[1.e-7,False],'Es':[1.e-7,False],'Ep':[1.e-7,False]},],
     
    56075608            copyNames = ['Scale','Extinction','Babinet','Flack','Twins']
    56085609        else:  #PWDR 
    5609             copyNames = ['Scale','Pref.Ori.','Size','Mustrain','HStrain','Extinction','Babinet']
     5610            copyNames = ['Scale','Pref.Ori.','Size','Mustrain','HStrain','Extinction','Babinet','LeBail']
    56105611        copyDict = {}
    56115612        for name in copyNames:
     
    57025703            copyNames = ['Scale','Extinction','Babinet','Flack','Twins']
    57035704        else:  #PWDR 
    5704             copyNames = ['Scale','Pref.Ori.','Size','Mustrain','HStrain','Extinction','Babinet']
     5705            copyNames = ['Scale','Pref.Ori.','Size','Mustrain','HStrain','Extinction','Babinet','LeBail']
    57055706        dlg = G2G.G2MultiChoiceDialog(G2frame.dataFrame,'Select which parameters to copy',
    57065707            'Select phase data parameters', copyNames)
     
    57515752                    for histoName in newList:
    57525753                        Id = G2gd.GetPatternTreeItemId(G2frame,G2frame.root,histoName)
    5753                         UseList[histoName] = {'Histogram':histoName,'Show':False,
     5754                        UseList[histoName] = {'Histogram':histoName,'Show':False,'LeBail':False,
    57545755                            'Scale':[1.0,False],'Pref.Ori.':['MD',1.0,False,[0,0,1],0,{},['',],0.1],
    57555756                            'Size':['isotropic',[1.,1.,1.],[False,False,False],[0,0,1],
     
    73777378            im = 1
    73787379        HistoNames = filter(lambda a:Histograms[a]['Use']==True,Histograms.keys())
    7379         PatternId = G2gd.GetPatternTreeItemId(G2frame,G2frame.root,HistoNames[0])
     7380        PatternId = G2gd.GetPatternTreeItemId(G2frame,G2frame.root,HistoNames[0])       #only use 1st histogram
    73807381        xdata = G2frame.PatternTree.GetItemPyData(PatternId)[1]
    73817382        Inst = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame,PatternId,'Instrument Parameters'))[0]
     
    74187419    def OnPawleyUpdate(event):
    74197420        '''This is the place for any reflection modification trick
    7420         Patterson squared, leBail extraction, etc.
     7421        Patterson squared, etc.
    74217422        '''
    74227423        try:
  • trunk/GSASIIstrIO.py

    r2715 r2718  
    22332233                elif 'T' in inst['Type'][0]:
    22342234                    dmin = limits[0]/inst['difC'][1]
     2235                pfx = str(pId)+':'+str(hId)+':'
    22352236                if Phases[phase]['General']['doPawley']:
     2237                    hapDict[pfx+'LeBail'] = False           #Pawley supercedes LeBail
    22362238                    Tmin = G2lat.Dsp2pos(inst,dmin)
    22372239                    if 'C' in inst['Type'][1]:
     
    22392241                    else:
    22402242                        limits[0] = max(limits[0],Tmin)
     2243                else:
     2244                    hapDict[pfx+'LeBail'] = hapData.get('LeBail',False)
    22412245                if Phases[phase]['General']['Type'] == 'magnetic':
    22422246                    dmin = max(dmin,Phases[phase]['General']['MagDmin'])
    2243                 pfx = str(pId)+':'+str(hId)+':'
    22442247                for item in ['Scale','Extinction']:
    22452248                    hapDict[pfx+item] = hapData[item][0]
    2246                     if hapData[item][1]:
     2249                    if hapData[item][1] and not hapDict[pfx+'LeBail']:
    22472250                        hapVary.append(pfx+item)
    22482251                names = G2spc.HStrainNames(SGData)
     
    22532256                    if hapData['HStrain'][1][i]:
    22542257                        hapVary.append(pfx+name)
    2255                 controlDict[pfx+'poType'] = hapData['Pref.Ori.'][0]
    2256                 if hapData['Pref.Ori.'][0] == 'MD':
    2257                     hapDict[pfx+'MD'] = hapData['Pref.Ori.'][1]
    2258                     controlDict[pfx+'MDAxis'] = hapData['Pref.Ori.'][3]
    2259                     if hapData['Pref.Ori.'][2]:
    2260                         hapVary.append(pfx+'MD')
    2261                 else:                           #'SH' spherical harmonics
    2262                     controlDict[pfx+'SHord'] = hapData['Pref.Ori.'][4]
    2263                     controlDict[pfx+'SHncof'] = len(hapData['Pref.Ori.'][5])
    2264                     controlDict[pfx+'SHnames'] = G2lat.GenSHCoeff(SGData['SGLaue'],'0',controlDict[pfx+'SHord'],False)
    2265                     controlDict[pfx+'SHhkl'] = []
    2266                     try: #patch for old Pref.Ori. items
    2267                         controlDict[pfx+'SHtoler'] = 0.1
    2268                         if hapData['Pref.Ori.'][6][0] != '':
    2269                             controlDict[pfx+'SHhkl'] = [eval(a.replace(' ',',')) for a in hapData['Pref.Ori.'][6]]
    2270                         controlDict[pfx+'SHtoler'] = hapData['Pref.Ori.'][7]
    2271                     except IndexError:
    2272                         pass
    2273                     for item in hapData['Pref.Ori.'][5]:
    2274                         hapDict[pfx+item] = hapData['Pref.Ori.'][5][item]
    2275                         if hapData['Pref.Ori.'][2]:
    2276                             hapVary.append(pfx+item)
     2258                    controlDict[pfx+'poType'] = hapData['Pref.Ori.'][0]
     2259                    if hapData['Pref.Ori.'][0] == 'MD':
     2260                        hapDict[pfx+'MD'] = hapData['Pref.Ori.'][1]
     2261                        controlDict[pfx+'MDAxis'] = hapData['Pref.Ori.'][3]
     2262                        if hapData['Pref.Ori.'][2] and not hapDict[pfx+'LeBail']:
     2263                            hapVary.append(pfx+'MD')
     2264                    else:                           #'SH' spherical harmonics
     2265                        controlDict[pfx+'SHord'] = hapData['Pref.Ori.'][4]
     2266                        controlDict[pfx+'SHncof'] = len(hapData['Pref.Ori.'][5])
     2267                        controlDict[pfx+'SHnames'] = G2lat.GenSHCoeff(SGData['SGLaue'],'0',controlDict[pfx+'SHord'],False)
     2268                        controlDict[pfx+'SHhkl'] = []
     2269                        try: #patch for old Pref.Ori. items
     2270                            controlDict[pfx+'SHtoler'] = 0.1
     2271                            if hapData['Pref.Ori.'][6][0] != '':
     2272                                controlDict[pfx+'SHhkl'] = [eval(a.replace(' ',',')) for a in hapData['Pref.Ori.'][6]]
     2273                            controlDict[pfx+'SHtoler'] = hapData['Pref.Ori.'][7]
     2274                        except IndexError:
     2275                            pass
     2276                        for item in hapData['Pref.Ori.'][5]:
     2277                            hapDict[pfx+item] = hapData['Pref.Ori.'][5][item]
     2278                            if hapData['Pref.Ori.'][2] and not hapDict[pfx+'LeBail']:
     2279                                hapVary.append(pfx+item)
    22772280                for item in ['Mustrain','Size']:
    22782281                    controlDict[pfx+item+'Type'] = hapData[item][0]
     
    23062309                    for bab in ['BabA','BabU']:
    23072310                        hapDict[pfx+bab] = hapData['Babinet'][bab][0]
    2308                         if hapData['Babinet'][bab][1]:
     2311                        if hapData['Babinet'][bab][1] and not hapDict[pfx+'LeBail']:
    23092312                            hapVary.append(pfx+bab)
    23102313                               
     
    23122315                    print >>pFile,'\n Phase: ',phase,' in histogram: ',histogram
    23132316                    print >>pFile,135*'='
    2314                     print >>pFile,' Phase fraction  : %10.4f'%(hapData['Scale'][0]),' Refine?',hapData['Scale'][1]
    2315                     print >>pFile,' Extinction coeff: %10.4f'%(hapData['Extinction'][0]),' Refine?',hapData['Extinction'][1]
    2316                     if hapData['Pref.Ori.'][0] == 'MD':
    2317                         Ax = hapData['Pref.Ori.'][3]
    2318                         print >>pFile,' March-Dollase PO: %10.4f'%(hapData['Pref.Ori.'][1]),' Refine?',hapData['Pref.Ori.'][2], \
    2319                             ' Axis: %d %d %d'%(Ax[0],Ax[1],Ax[2])
    2320                     else: #'SH' for spherical harmonics
    2321                         PrintSHPO(hapData['Pref.Ori.'])
    2322                         print >>pFile,' Penalty hkl list: '+str(controlDict[pfx+'SHhkl'])+' tolerance: %.2f'%(controlDict[pfx+'SHtoler'])
     2317                    if hapDict[pfx+'LeBail']:
     2318                        print >>pFile,' Perform Lebail extraction'                       
     2319                    else:
     2320                        print >>pFile,' Phase fraction  : %10.4f'%(hapData['Scale'][0]),' Refine?',hapData['Scale'][1]
     2321                        print >>pFile,' Extinction coeff: %10.4f'%(hapData['Extinction'][0]),' Refine?',hapData['Extinction'][1]
     2322                        if hapData['Pref.Ori.'][0] == 'MD':
     2323                            Ax = hapData['Pref.Ori.'][3]
     2324                            print >>pFile,' March-Dollase PO: %10.4f'%(hapData['Pref.Ori.'][1]),' Refine?',hapData['Pref.Ori.'][2], \
     2325                                ' Axis: %d %d %d'%(Ax[0],Ax[1],Ax[2])
     2326                        else: #'SH' for spherical harmonics
     2327                            PrintSHPO(hapData['Pref.Ori.'])
     2328                            print >>pFile,' Penalty hkl list: '+str(controlDict[pfx+'SHhkl'])+' tolerance: %.2f'%(controlDict[pfx+'SHtoler'])
    23232329                    PrintSize(hapData['Size'])
    23242330                    PrintMuStrain(hapData['Mustrain'],SGData)
    23252331                    PrintHStrain(hapData['HStrain'],SGData)
    2326                     if Phases[phase]['General']['Type'] != 'magnetic':
     2332                    if Phases[phase]['General']['Type'] != 'magnetic' and not hapDict[pfx+'LeBail']:
    23272333                        if hapData['Babinet']['BabA'][0]:
    23282334                            PrintBabinet(hapData['Babinet'])                       
    2329                 if resetRefList:
     2335                if resetRefList and not hapDict[pfx+'LeBail']:
    23302336                    refList = []
    23312337                    Uniq = []
     
    23432349                                    pos = G2lat.Dsp2pos(inst,d)
    23442350                                    if limits[0] < pos < limits[1]:
    2345                                         refList.append([h,k,l,m,mul,d, pos,0.0,0.0,0.0,0.0, 0.0,0.0,1.0,1.0,1.0])
     2351                                        refList.append([h,k,l,m,mul,d, pos,0.0,0.0,0.0,100., 0.0,0.0,1.0,1.0,1.0])
    23462352                                        #... sig,gam,fotsq,fctsq, phase,icorr,prfo,abs,ext
    23472353                                        Uniq.append(uniq)
     
    23512357                                    if limits[0] < pos < limits[1]:
    23522358                                        wave = inst['difC'][1]*d/(252.816*inst['fltPath'][0])
    2353                                         refList.append([h,k,l,m,mul,d, pos,0.0,0.0,0.0,0.0, 0.0,0.0,0.0,0.0,wave, 1.0,1.0,1.0])
     2359                                        refList.append([h,k,l,m,mul,d, pos,0.0,0.0,0.0,100., 0.0,0.0,0.0,0.0,wave, 1.0,1.0,1.0])
    23542360                                        # ... sig,gam,fotsq,fctsq, phase,icorr,alp,bet,wave, prfo,abs,ext
    23552361                                        Uniq.append(uniq)
     
    23672373                                pos = G2lat.Dsp2pos(inst,d)
    23682374                                if limits[0] < pos < limits[1]:
    2369                                     refList.append([h,k,l,mul,d, pos,0.0,0.0,0.0,0.0, 0.0,0.0,1.0,1.0,1.0])
     2375                                    refList.append([h,k,l,mul,d, pos,0.0,0.0,0.0,100., 0.0,0.0,1.0,1.0,1.0])
    23702376                                    #... sig,gam,fotsq,fctsq, phase,icorr,prfo,abs,ext
    23712377                                    Uniq.append(uniq)
     
    23752381                                if limits[0] < pos < limits[1]:
    23762382                                    wave = inst['difC'][1]*d/(252.816*inst['fltPath'][0])
    2377                                     refList.append([h,k,l,mul,d, pos,0.0,0.0,0.0,0.0, 0.0,0.0,0.0,0.0,wave, 1.0,1.0,1.0])
     2383                                    refList.append([h,k,l,mul,d, pos,0.0,0.0,0.0,100., 0.0,0.0,0.0,0.0,wave, 1.0,1.0,1.0])
    23782384                                    # ... sig,gam,fotsq,fctsq, phase,icorr,alp,bet,wave, prfo,abs,ext
    23792385                                    Uniq.append(uniq)
     
    26682674                wtFrSum[hId] = 0.
    26692675            if 'PWDR' in histogram:
    2670                 for item in ['Scale','Extinction']:
    2671                     hapData[item][0] = parmDict[pfx+item]
    2672                     if pfx+item in sigDict:
    2673                         PhFrExtPOSig.update({pfx+item:sigDict[pfx+item],})
    2674                 wtFrSum[hId] += hapData['Scale'][0]*General['Mass']
    2675                 if hapData['Pref.Ori.'][0] == 'MD':
    2676                     hapData['Pref.Ori.'][1] = parmDict[pfx+'MD']
    2677                     if pfx+'MD' in sigDict:
    2678                         PhFrExtPOSig.update({pfx+'MD':sigDict[pfx+'MD'],})
    2679                 else:                           #'SH' spherical harmonics
    2680                     for item in hapData['Pref.Ori.'][5]:
    2681                         hapData['Pref.Ori.'][5][item] = parmDict[pfx+item]
     2676                if not parmDict[pfx+'LeBail']:
     2677                    for item in ['Scale','Extinction']:
     2678                        hapData[item][0] = parmDict[pfx+item]
    26822679                        if pfx+item in sigDict:
    26832680                            PhFrExtPOSig.update({pfx+item:sigDict[pfx+item],})
     2681                    wtFrSum[hId] += hapData['Scale'][0]*General['Mass']
     2682                    if hapData['Pref.Ori.'][0] == 'MD':
     2683                        hapData['Pref.Ori.'][1] = parmDict[pfx+'MD']
     2684                        if pfx+'MD' in sigDict:
     2685                            PhFrExtPOSig.update({pfx+'MD':sigDict[pfx+'MD'],})
     2686                    else:                           #'SH' spherical harmonics
     2687                        for item in hapData['Pref.Ori.'][5]:
     2688                            hapData['Pref.Ori.'][5][item] = parmDict[pfx+item]
     2689                            if pfx+item in sigDict:
     2690                                PhFrExtPOSig.update({pfx+item:sigDict[pfx+item],})
    26842691                SizeMuStrSig.update({pfx+'Mustrain':[[0,0,0],[0 for i in range(len(hapData['Mustrain'][4]))]],
    26852692                    pfx+'Size':[[0,0,0],[0 for i in range(len(hapData['Size'][4]))]],
     
    27142721                    if pfx+name in sigDict:
    27152722                        SizeMuStrSig[pfx+'HStrain'][name] = sigDict[pfx+name]
    2716                 if Phases[phase]['General']['Type'] != 'magnetic':
     2723                if Phases[phase]['General']['Type'] != 'magnetic' and not parmDict[pfx+'LeBail']:
    27172724                    for name in ['BabA','BabU']:
    27182725                        hapData['Babinet'][name][0] = parmDict[pfx+name]
     
    27762783                    print >>pFile,' Durbin-Watson statistic = %.3f'%(Histogram['Residuals']['Durbin-Watson'])
    27772784                    print >>pFile,' Bragg intensity sum = %.3g'%(Histogram['Residuals'][pfx+'sumInt'])
    2778                
    2779                     if pfx+'Scale' in PhFrExtPOSig:
    2780                         wtFr = hapData['Scale'][0]*General['Mass']/wtFrSum[hId]
    2781                         sigwtFr = PhFrExtPOSig[pfx+'Scale']*wtFr/hapData['Scale'][0]
    2782                         print >>pFile,' Phase fraction  : %10.5f, sig %10.5f Weight fraction  : %8.5f, sig %10.5f' \
    2783                             %(hapData['Scale'][0],PhFrExtPOSig[pfx+'Scale'],wtFr,sigwtFr)
    2784                     if pfx+'Extinction' in PhFrExtPOSig:
    2785                         print >>pFile,' Extinction coeff: %10.4f, sig %10.4f'%(hapData['Extinction'][0],PhFrExtPOSig[pfx+'Extinction'])
    2786                     if hapData['Pref.Ori.'][0] == 'MD':
    2787                         if pfx+'MD' in PhFrExtPOSig:
    2788                             print >>pFile,' March-Dollase PO: %10.4f, sig %10.4f'%(hapData['Pref.Ori.'][1],PhFrExtPOSig[pfx+'MD'])
     2785                   
     2786                    if parmDict[pfx+'LeBail']:
     2787                        print >>pFile,' Performed LeBail extraction for phase '+phase+' in histogram '+histogram
    27892788                    else:
    2790                         PrintSHPOAndSig(pfx,hapData['Pref.Ori.'],PhFrExtPOSig)
     2789                        if pfx+'Scale' in PhFrExtPOSig:
     2790                            wtFr = hapData['Scale'][0]*General['Mass']/wtFrSum[hId]
     2791                            sigwtFr = PhFrExtPOSig[pfx+'Scale']*wtFr/hapData['Scale'][0]
     2792                            print >>pFile,' Phase fraction  : %10.5f, sig %10.5f Weight fraction  : %8.5f, sig %10.5f' \
     2793                                %(hapData['Scale'][0],PhFrExtPOSig[pfx+'Scale'],wtFr,sigwtFr)
     2794                        if pfx+'Extinction' in PhFrExtPOSig:
     2795                            print >>pFile,' Extinction coeff: %10.4f, sig %10.4f'%(hapData['Extinction'][0],PhFrExtPOSig[pfx+'Extinction'])
     2796                        if hapData['Pref.Ori.'][0] == 'MD':
     2797                            if pfx+'MD' in PhFrExtPOSig:
     2798                                print >>pFile,' March-Dollase PO: %10.4f, sig %10.4f'%(hapData['Pref.Ori.'][1],PhFrExtPOSig[pfx+'MD'])
     2799                        else:
     2800                            PrintSHPOAndSig(pfx,hapData['Pref.Ori.'],PhFrExtPOSig)
    27912801                    PrintSizeAndSig(hapData['Size'],SizeMuStrSig[pfx+'Size'])
    27922802                    PrintMuStrainAndSig(hapData['Mustrain'],SizeMuStrSig[pfx+'Mustrain'],SGData)
    27932803                    PrintHStrainAndSig(hapData['HStrain'],SizeMuStrSig[pfx+'HStrain'],SGData)
    2794                     if Phases[phase]['General']['Type'] != 'magnetic':
     2804                    if Phases[phase]['General']['Type'] != 'magnetic' and not parmDict[pfx+'LeBail']:
    27952805                        if len(BabSig):
    27962806                            PrintBabinetAndSig(pfx,hapData['Babinet'],BabSig)
  • trunk/GSASIIstrMath.py

    r2715 r2718  
    29662966                                    sumInt += refl[11+im]*refl[9+im]*kRatio
    29672967                            refl[8+im] = np.sum(np.where(ratio[iBeg:iFin2]>0.,yp[iBeg:iFin2]*ratio[iBeg:iFin2]/(refl[11+im]*(1.+kRatio)),0.0))
     2968                            if parmDict[phfx+'LeBail']:
     2969                                refl[9+im] = refl[8+im]
    29682970                               
    29692971                    elif 'T' in calcControls[hfx+'histType']:
     
    29832985                            yp[iBeg:iFin] = refl[11+im]*refl[9+im]*G2pwd.getEpsVoigt(refl[5+im],refl[12+im],refl[13+im],refl[6+im],refl[7+im],ma.getdata(x[iBeg:iFin]))  #>90% of time spent here
    29842986                            refl[8+im] = np.sum(np.where(ratio[iBeg:iFin]>0.,yp[iBeg:iFin]*ratio[iBeg:iFin]/refl[11+im],0.0))
     2987                            if parmDict[phfx+'LeBail']:
     2988                                refl[9+im] = refl[8+im]
    29852989                            sumInt += refl[11+im]*refl[9+im]
    29862990                    Fo = np.sqrt(np.abs(refl[8+im]))
     
    30753079        GA,GB = G2lat.Gmat2AB(G)    #Orthogonalization matricies
    30763080        Vst = np.sqrt(nl.det(G))    #V*
    3077         if not Phase['General'].get('doPawley'):
     3081        if not Phase['General'].get('doPawley') and not parmDict[phfx+'LeBail']:
    30783082            if im:
    30793083                SStructureFactor(refDict,G,hfx,pfx,SGData,SSGData,calcControls,parmDict)
     
    32503254        G,g = G2lat.A2Gmat(A)       #recip & real metric tensors
    32513255        GA,GB = G2lat.Gmat2AB(G)    #Orthogonalization matricies
    3252         if not Phase['General'].get('doPawley'):
     3256        if not Phase['General'].get('doPawley')  and not parmDict[phfx+'LeBail']:
    32533257            if im:
    32543258                dFdvDict = SStructureFactorDerv(refDict,im,G,hfx,pfx,SGData,SSGData,calcControls,parmDict)
     
    34453449                        if Ka2 and iFin2-iBeg2:
    34463450                            depDerivDict[phfx+name][iBeg2:iFin2] += parmDict[phfx+'Scale']*dFdvDict[phfx+name][iref]*dervDict2['int']/refl[9+im]                 
    3447             if not Phase['General'].get('doPawley'):
     3451            if not Phase['General'].get('doPawley')  and not parmDict[phfx+'LeBail']:
    34483452                #do atom derivatives -  for RB,F,X & U so far - how do I scale mixed phase constraints?
    34493453                corr = 0.
Note: See TracChangeset for help on using the changeset viewer.