Changeset 4077 for trunk/GSASIIsasd.py


Ignore:
Timestamp:
Aug 4, 2019 3:28:46 PM (2 years ago)
Author:
vondreele
Message:

add SHAPES capability & all the interface stuff

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIsasd.py

    r4076 r4077  
    10841084
    10851085def PairDistFxn(Profile,ProfDict,Limits,Sample,data):
    1086     print('create P(R) fit to data - TBD')
    10871086   
    10881087    def CalcMoore():
    1089        
    1090 #        def MooreRg(cw,q):      #lines 1400-1409
    1091 #            n = 0
    1092 #            nmax = len(cw)-5
    1093 #            POR = 0.
    1094 #            dmax = cw[nmax+2]
    1095 #            while n < nmax:
    1096 #                POR += cw[n]*((-1**n)*((np.pi*(n+1))**2-6)/((n+1)**3))
    1097 #                n += 1
    1098 #            POR *= 4*(dmax**4)/((np.pi**2)*cw[nmax+3])
    1099 #            return POR
    11001088       
    11011089        def MoorePOR(cw,r,dmax):    #lines 1417-1428
     
    11181106                n += 1
    11191107            return POR
    1120        
    1121 #        def MooreINaught(cw,q,dmax):         #lines 1457-1466
    1122 #            n = 0
    1123 #            nmax = len(cw)
    1124 #            POR = 0.
    1125 #            while n < nmax:
    1126 #                POR += cw[n]*8*(dmax**2)*((-1.)**n)/(n+1)
    1127 #                n += 1
    1128 #            return POR
    11291108       
    11301109        def calcSASD(values,Q,Io,wt,Ifb,dmax,ifBack):
     
    11501129            N += 1
    11511130        MPV = np.zeros(N)
    1152         MPS = np.zeros(N)
    11531131        MPV[0] = Q[Ibeg]
    11541132        dmax = pairData['MaxRadius']
     
    11641142        chisq = np.sum(result[2]['fvec']**2)
    11651143        Ic[Ibeg:Ifin] = MooreIOREFF(MPVR,Q[Ibeg:Ifin],dmax)+Ifb+Back
    1166         covM = result[1]
     1144        ncalc = result[2]['nfev']
    11671145        GOF = chisq/(Ifin-Ibeg-N)
    1168         MPS = np.sqrt(np.diag(covM)*GOF)
     1146        Rwp = np.sqrt(chisq/np.sum(wt[Ibeg:Ifin]*Io[Ibeg:Ifin]**2))*100.      #to %
     1147        print (' Results of small angle data modelling fit of P(R):')
     1148        print ('Number of function calls: %d Number of observations: %d Number of parameters: %d'%(ncalc,Ifin-Ibeg,N))
     1149        print ('Rwp = %7.2f%%, chi**2 = %12.6g, reduced chi**2 = %6.2f'%(Rwp,chisq,GOF))
    11691150        BinMag = MoorePOR(MPVR,Bins,dmax)/2.
    11701151        return Bins,Dbins,BinMag
Note: See TracChangeset for help on using the changeset viewer.