Ignore:
Timestamp:
Mar 9, 2017 8:54:25 AM (5 years ago)
Author:
vondreele
Message:

fix (again) a LeBail? bug on loading a new phase

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrMath.py

    r2750 r2752  
    849849    if not SGData['SGInv'] and 'S' in calcControls[hfx+'histType'] and phfx+'Flack' in parmDict:
    850850        Flack = 1.-2.*parmDict[phfx+'Flack']
    851     time0 = time.time()
     851#    time0 = time.time()
    852852#reflection processing begins here - big arrays!
    853853    iBeg = 0
     
    941941#        GSASIIpath.IPyBreak()
    942942        iBeg += blkSize
    943     print ' %d derivative time %.4f\r'%(nRef,time.time()-time0)
     943#    print ' %d derivative time %.4f\r'%(nRef,time.time()-time0)
    944944        #loop over atoms - each dict entry is list of derivatives for all the reflections
    945945    for i in range(len(Mdata)):
     
    10781078        #accumulate derivatives   
    10791079        dFdfr[iBeg:iFin] = 2.*np.sum((fams[:,:,nxs]*dfadfr+fbms[:,:,nxs]*dfbdfr)*Mdata/Nops,axis=0) #ok
    1080         dFdx[iBeg:iFin] =  2.*np.sum(fams[:,:,nxs,nxs]*dfadx+fbms[:,:,nxs,nxs]*dfbdx,axis=0)        #ok
     1080        dFdx[iBeg:iFin] = 2.*np.sum(fams[:,:,nxs,nxs]*dfadx+fbms[:,:,nxs,nxs]*dfbdx,axis=0)         #ok
    10811081        dFdMx[:,iBeg:iFin,:] = 2.*(fams[:,:,nxs]*dfadmx+fbms[:,:,nxs]*dfbdmx)                       #problems
    10821082        dFdui[iBeg:iFin] = 2.*np.sum(fams[:,:,nxs]*dfadui+fbms[:,:,nxs]*dfbdui,axis=0)              #ok
Note: See TracChangeset for help on using the changeset viewer.