Changeset 910


Ignore:
Timestamp:
May 16, 2013 1:29:11 PM (9 years ago)
Author:
vondreele
Message:

closer to RB
fix plot problem

Location:
trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIImath.py

    r905 r910  
    204204            UIJ = atom[cia+2:cia+8]
    205205               
    206 def TLS2Uij(xyz,g,Amat,rbObj):
     206def TLS2Uij(xyz,g,Amat,rbObj):    #not used anywhere, but could be?
    207207    TLStype,TLS = rbObj['ThermalMotion'][:2]
    208208    Tmat = np.zeros((3,3))
     
    296296    if 'S' in TLStype:
    297297        S = np.array(TLS[12:])*(np.pi/180.)
    298     g = np.inner(Bmat,Bmat.T)
     298    g = np.inner(Bmat,Bmat)
    299299    gvec = 1./np.sqrt(np.array([g[0][0]**2,g[1][1]**2,g[2][2]**2,
    300300        g[0][0]*g[1][1],g[0][0]*g[2][2],g[1][1]*g[2][2]]))
     
    317317                S[0]*X[1]-S[1]*X[2]+S[7]*X[0]
    318318            Umat = G2lat.U6toUij(U)
    319             beta = np.inner(np.inner(Bmat,Umat),Bmat.T)
     319            beta = np.inner(np.inner(Bmat.T,Umat),Bmat)
    320320            Uout.append(['A',0.0,]+list(G2lat.UijtoU6(beta)*gvec))
    321321        else:
  • trunk/GSASIIplot.py

    r885 r910  
    26292629        flipData = generalData['Flip']                       
    26302630        flipData['mapRoll'] = [0,0,0]
    2631     cx,ct,cs,ci = drawingData['atomPtrs']
    26322631    Wt = np.array([255,255,255])
    26332632    Rd = np.array([255,0,0])
     
    26942693            key = str(event.key).upper()
    26952694        indx = drawingData['selectedAtoms']
     2695        cx,ct = drawingData['atomPtrs'][:2]
    26962696        if key in ['C']:
    26972697            drawingData['viewPoint'] = [[.5,.5,.5],[0,0]]
     
    27722772                        SetSelectedAtoms(i,Add)
    27732773        else:
     2774            cx = drawingData['atomPtrs'][0]
    27742775            for i,atom in enumerate(drawAtoms):
    27752776                x,y,z = atom[cx:cx+3]
  • trunk/GSASIIstruct.py

    r906 r910  
    745745    rpd = np.pi/180.
    746746    rpd2 = rpd**2
    747     g = np.inner(Bmat,Bmat.T)
     747    g = np.inner(Bmat,Bmat)
    748748    gvec = np.sqrt(np.array([g[0][0]**2,g[1][1]**2,g[2][2]**2,
    749749        g[0][0]*g[1][1],g[0][0]*g[2][2],g[1][1]*g[2][2]]))
     
    760760        dXdv = []
    761761        for iv in range(len(VModel['VectMag'])):
    762             dXdv.append(np.inner(Bmat,VModel['rbVect'][iv]).T)
     762            dCdv = []
     763            for vec in VModel['rbVect'][iv]:
     764                dCdv.append(G2mth.prodQVQ(Q,vec))
     765            dXdv.append(np.inner(Bmat,np.array(dCdv)).T)
    763766        XYZ,Cart = G2mth.UpdateRBXYZ(Bmat,RBObj,RBData,'Vector')
    764767        for ia,atId in enumerate(RBObj['Ids']):
     
    783786            dFdu = G2lat.U6toUij(dFdu.T)
    784787            dFdu = np.tensordot(Amat,np.tensordot(Amat,dFdu,([1,0])),([0,1]))
    785 #            dFdu = np.tensordot(QM,np.tensordot(QM,dFdu,([1,0])),([0,1]))
     788           
    786789            dFdu = G2lat.UijtoU6(dFdu)
    787790            atNum = AtLookup[atId]
Note: See TracChangeset for help on using the changeset viewer.