Changeset 1999


Ignore:
Timestamp:
Oct 12, 2015 9:12:45 AM (6 years ago)
Author:
vondreele
Message:
 
Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIlattice.py

    r1947 r1999  
    273273    """
    274274    U = np.array([
    275         [U6[0],  U6[3]/2.,  U6[4]/2.],
    276         [U6[3]/2.,  U6[1],  U6[5]/2.],
    277         [U6[4]/2.,  U6[5]/2.,  U6[2]]])
     275        [U6[0],  U6[3],  U6[4]],
     276        [U6[3],  U6[1],  U6[5]],
     277        [U6[4],  U6[5],  U6[2]]])
    278278    return U
    279279
     
    282282    NB: there is a non numpy version in GSASIIspc: Uij2U
    283283    """
    284     U6 = np.array([U[0][0],U[1][1],U[2][2],U[0][1]*2.,U[0][2]*2.,U[1][2]*2.])
     284    U6 = np.array([U[0][0],U[1][1],U[2][2],U[0][1],U[0][2],U[1][2]])
    285285    return U6
    286286
  • trunk/GSASIIstrMath.py

    r1996 r1999  
    754754        else:
    755755            if len(TwinLaw) > 1:
    756                 refl.T[9] = np.sum(fas[:,:,0],axis=0)**2+np.sum(fbs[:,:,0],axis=0)**2   #FcT from primary twin element
     756                refl.T[9] = np.sum(fas[:,:,0]**2,axis=0)+np.sum(fbs[:,:,0]**2,axis=0)   #FcT from primary twin element
    757757                refl.T[7] = np.sum(TwinFr*np.sum(TwMask[np.newaxis,:,:]*fas,axis=0)**2,axis=-1)+   \
    758758                    np.sum(TwinFr*np.sum(TwMask[np.newaxis,:,:]*fbs,axis=0)**2,axis=-1)                        #Fc sum over twins
    759759                refl.T[10] = atan2d(fbs[0].T[0],fas[0].T[0])  #ignore f' & f"
    760760            else:
    761                 refl.T[9] = np.sum(fas,axis=0)**2+np.sum(fbs,axis=0)**2
     761                refl.T[9] = np.sum(fas**2,axis=0)+np.sum(fbs**2,axis=0)
    762762                refl.T[7] = np.copy(refl.T[9])               
    763763                refl.T[10] = atan2d(fbs[0],fas[0])  #ignore f' & f"
     
    930930            dFdvDict[pfx+'AU22:'+str(i)] = np.sum(dFdua.T[1][i]*TwinFr[:,np.newaxis],axis=0)
    931931            dFdvDict[pfx+'AU33:'+str(i)] = np.sum(dFdua.T[2][i]*TwinFr[:,np.newaxis],axis=0)
    932             dFdvDict[pfx+'AU12:'+str(i)] = np.sum(0.5*dFdua.T[3][i]*TwinFr[:,np.newaxis],axis=0)
    933             dFdvDict[pfx+'AU13:'+str(i)] = np.sum(0.5*dFdua.T[4][i]*TwinFr[:,np.newaxis],axis=0)
    934             dFdvDict[pfx+'AU23:'+str(i)] = np.sum(0.5*dFdua.T[5][i]*TwinFr[:,np.newaxis],axis=0)
     932            dFdvDict[pfx+'AU12:'+str(i)] = np.sum(dFdua.T[3][i]*TwinFr[:,np.newaxis],axis=0)
     933            dFdvDict[pfx+'AU13:'+str(i)] = np.sum(dFdua.T[4][i]*TwinFr[:,np.newaxis],axis=0)
     934            dFdvDict[pfx+'AU23:'+str(i)] = np.sum(dFdua.T[5][i]*TwinFr[:,np.newaxis],axis=0)
    935935    else:
    936936        for i in range(len(Mdata)):
     
    943943            dFdvDict[pfx+'AU22:'+str(i)] = dFdua.T[1][i]
    944944            dFdvDict[pfx+'AU33:'+str(i)] = dFdua.T[2][i]
    945             dFdvDict[pfx+'AU12:'+str(i)] = 0.5*dFdua.T[3][i]
    946             dFdvDict[pfx+'AU13:'+str(i)] = 0.5*dFdua.T[4][i]
    947             dFdvDict[pfx+'AU23:'+str(i)] = 0.5*dFdua.T[5][i]
     945            dFdvDict[pfx+'AU12:'+str(i)] = dFdua.T[3][i]
     946            dFdvDict[pfx+'AU13:'+str(i)] = dFdua.T[4][i]
     947            dFdvDict[pfx+'AU23:'+str(i)] = dFdua.T[5][i]
    948948        dFdvDict[phfx+'Flack'] = 4.*dFdfl.T
    949949    dFdvDict[phfx+'BabA'] = dFdbab.T[0]
     
    10751075        else:
    10761076            if len(TwinLaw) > 1:
    1077                 refl.T[10] = np.sum(fas[:,:,0],axis=0)**2+np.sum(fbs[:,:,0],axis=0)**2   #FcT from primary twin element
     1077                refl.T[10] = np.sum(fas[:,:,0]**2,axis=0)+np.sum(fbs[:,:,0]**2,axis=0)   #FcT from primary twin element
    10781078                refl.T[8] = np.sum(TwinFr*np.sum(TwMask[np.newaxis,:,:]*fas,axis=0)**2,axis=-1)+   \
    10791079                    np.sum(TwinFr*np.sum(TwMask[np.newaxis,:,:]*fbs,axis=0)**2,axis=-1)                        #Fc sum over twins
    10801080                refl.T[11] = atan2d(fbs[0].T[0],fas[0].T[0])  #ignore f' & f"
    10811081            else:
    1082                 refl.T[10] = np.sum(fas,axis=0)**2+np.sum(fbs,axis=0)**2
     1082                refl.T[10] = np.sum(fas**2,axis=0)+np.sum(fbs**2,axis=0)
    10831083                refl.T[8] = np.copy(refl.T[10])               
    10841084                refl.T[11] = atan2d(fbs[0],fas[0])  #ignore f' & f"
     
    12801280        dFdvDict[pfx+'AU22:'+str(i)] = dFdua.T[1][i]
    12811281        dFdvDict[pfx+'AU33:'+str(i)] = dFdua.T[2][i]
    1282         dFdvDict[pfx+'AU12:'+str(i)] = .5*dFdua.T[3][i]
    1283         dFdvDict[pfx+'AU13:'+str(i)] = .5*dFdua.T[4][i]
    1284         dFdvDict[pfx+'AU23:'+str(i)] = .5*dFdua.T[5][i]
     1282        dFdvDict[pfx+'AU12:'+str(i)] = dFdua.T[3][i]
     1283        dFdvDict[pfx+'AU13:'+str(i)] = dFdua.T[4][i]
     1284        dFdvDict[pfx+'AU23:'+str(i)] = dFdua.T[5][i]
    12851285        for j in range(XSSdata.shape[1]):
    12861286            dFdvDict[pfx+'Xsin:'+str(i)+':'+str(j)] = dFdGx.T[0][j][i]
Note: See TracChangeset for help on using the changeset viewer.