Changeset 1862 for trunk/GSASIIstrMath.py
- Timestamp:
- May 21, 2015 8:52:52 AM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/GSASIIstrMath.py
r1860 r1862 2473 2473 if ref[5+im] > 0.: 2474 2474 dervDict,dervCor = SCExtinction(ref,im,phfx,hfx,pfx,calcControls,parmDict,varylist+dependentVars)[1:] 2475 print ref[:3],dervCor,ref[11+im]2476 2475 Fo = np.sqrt(ref[5+im]) 2477 2476 Fc = np.sqrt(ref[7+im]) 2478 2477 w = 1.0/ref[6+im] 2479 2478 if ref[3+im] > 0: 2480 wdf[iref] = 2.0*F c*w*(Fo-Fc)2479 wdf[iref] = 2.0*Fo*w*(Fo-Fc) 2481 2480 for j,var in enumerate(varylist): 2482 2481 if var in dFdvDict: … … 2488 2487 dMdvh[varylist.index(phfx+'Scale')][iref] = w*ref[9+im]*ref[11+im] 2489 2488 elif phfx+'Scale' in dependentVars: 2490 depDerivDict[phfx+'Scale'][iref] = w*ref[9+im]*ref[11+im] 2489 depDerivDict[phfx+'Scale'][iref] = w*ref[9+im]*ref[11+im] 2491 2490 for item in ['Ep','Es','Eg']: 2492 2491 if phfx+item in varylist and phfx+item in dervDict: 2493 dMdvh[varylist.index(phfx+item)][iref] = w*dervDict[phfx+item] *dervCor/ref[11+im]2492 dMdvh[varylist.index(phfx+item)][iref] = w*dervDict[phfx+item]/ref[11+im] 2494 2493 elif phfx+item in dependentVars and phfx+item in dervDict: 2495 depDerivDict[phfx+item][iref] = w*dervDict[phfx+item] *dervCor/ref[11+im]2494 depDerivDict[phfx+item][iref] = w*dervDict[phfx+item]/ref[11+im] 2496 2495 for item in ['BabA','BabU']: 2497 2496 if phfx+item in varylist:
Note: See TracChangeset
for help on using the changeset viewer.