Ignore:
Timestamp:
May 26, 2021 7:43:09 AM (5 months ago)
Author:
vondreele
Message:

add some doc strings for profile functions in G2pwd
rename getPowderProfileDervMP as getPowderProfileDerv since it is more generally used than just multiprocessing

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrMath.py

    r4914 r4915  
    30613061                            profArgs[iref%G2mp.ncores].append((refl,iref))
    30623062                        else:
    3063                             icod= G2mp.ComputeFobsSqCW(refl,iref)
     3063                            icod = G2mp.ComputeFobsSqCW(refl,iref)
    30643064                            if type(icod) is tuple:
    30653065                                refl[8+im] = icod[0]
     
    30873087                            profArgs[iref%G2mp.ncores].append((refl,iref))
    30883088                        else:
    3089                             icod= G2mp.ComputeFobsSqTOF(refl,iref)
     3089                            icod = G2mp.ComputeFobsSqTOF(refl,iref)
    30903090                            if type(icod) is tuple:
    30913091                                refl[8+im] = icod[0]
     
    31133113                            profArgs[iref%G2mp.ncores].append((refl,iref))
    31143114                        else:
    3115                             icod= G2mp.ComputeFobsSqPink(refl,iref)
     3115                            icod = G2mp.ComputeFobsSqPink(refl,iref)
    31163116                            if type(icod) is tuple:
    31173117                                refl[8+im] = icod[0]
     
    31413141                for iref,refl in enumerate(refDict['RefList']):
    31423142                    Fo = np.sqrt(np.abs(refl[8+im]))
    3143                     Fc = np.sqrt(np.abs(refl[9]+im))
     3143                    Fc = np.sqrt(np.abs(refl[9+im]))
    31443144                    sumFo += Fo
    31453145                    sumFosq += refl[8+im]**2
     
    34223422    return yc,yb
    34233423   
    3424 def getPowderProfileDervMP(args):
     3424def getPowderProfileDerv(args):
    34253425    '''Computes the derivatives of the computed powder pattern with respect to all
    34263426    refined parameters.
    3427     Multiprocessing version.
     3427    Used for single processor & Multiprocessor versions
    34283428    '''
    34293429    import pytexture as ptx
     
    39273927            xB = np.searchsorted(x,Limits[0])
    39283928            xF = np.searchsorted(x,Limits[1])+1
    3929             dMdv,depDerivDict = getPowderProfileDervMP([parmDict,x[xB:xF],
     3929            dMdv,depDerivDict = getPowderProfileDerv([parmDict,x[xB:xF],
    39303930                varylist,Histogram,Phases,rigidbodyDict,calcControls,pawleyLookup,dependentVars])
    39313931            G2mv.Dict2Deriv(varylist,depDerivDict,dMdv)
     
    39993999#                    (parmDict,x[xB:xF],varylist,Histogram,Phases,rigidbodyDict,calcControls,pawleyLookup,dependentVars,
    40004000#                     i,ncores,histogram) for i in range(ncores)]
    4001 #                for dmdv,depDerivs in MPpool.imap_unordered(getPowderProfileDervMP,profArgs):
     4001#                for dmdv,depDerivs in MPpool.imap_unordered(getPowderProfileDerv,profArgs):
    40024002                # better, use a generator so arg is created as used
    40034003                profGenArgs = (
    40044004                    (parmDict,x[xB:xF],varylist,Histogram,Phases,rigidbodyDict,calcControls,pawleyLookup,dependentVars,
    40054005                     i,ncores,histogram) for i in range(ncores))
    4006                 for dmdv,depDerivs in MPpool.imap_unordered(getPowderProfileDervMP,profGenArgs):
     4006                for dmdv,depDerivs in MPpool.imap_unordered(getPowderProfileDerv,profGenArgs):
    40074007                    if dMdvh is None:
    40084008                       dMdvh = dmdv
     
    40134013                MPpool.terminate()
    40144014            else:
    4015                 dMdvh,depDerivDict = getPowderProfileDervMP([parmDict,x[xB:xF],
     4015                dMdvh,depDerivDict = getPowderProfileDerv([parmDict,x[xB:xF],
    40164016                    varylist,Histogram,Phases,rigidbodyDict,calcControls,pawleyLookup,dependentVars,0,1,histogram])
    40174017                #dMdvh = getPowderProfileDerv(parmDict,x[xB:xF],
Note: See TracChangeset for help on using the changeset viewer.