Changeset 2501 for trunk/GSASIIphsGUI.py


Ignore:
Timestamp:
Oct 22, 2016 1:21:40 PM (6 years ago)
Author:
vondreele
Message:

fix (again!) mag moment drawings
clean up testDeriv

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIphsGUI.py

    r2497 r2501  
    23172317            atomData = data['Atoms']
    23182318            SGData = generalData['SGData']
    2319             SpnFlp = SGData.get('SpnFlp',[])
    23202319            dlg = G2gd.SymOpDialog(G2frame,SGData,True,True)
    23212320            New = False
     
    43394338                                opNum = G2spc.GetOpNum(OprNum,SGData)
    43404339                                mom = np.inner(np.array(atom[cmx:cmx+3]),Bmat)
    4341                                 print OprNum,newOp,opNum,SpnFlp
    4342                                 if Inv:
    4343                                     atom[cmx:cmx+3] = np.inner(np.inner(mom,-M),Amat)*nl.det(M)*SpnFlp[opNum-1]
    4344                                 else:
    4345                                     atom[cmx:cmx+3] = np.inner(np.inner(mom,M),Amat)*nl.det(M)*SpnFlp[opNum-1]
     4340#                                print OprNum,newOp,opNum,SpnFlp
     4341                                atom[cmx:cmx+3] = np.inner(np.inner(mom,M),Amat)*nl.det(M)*SpnFlp[opNum-1]
    43464342                            if atom[cui] == 'A':
    43474343                                Uij = atom[cuij:cuij+6]
     
    43934389                                opNum = G2spc.GetOpNum(item[2],SGData)
    43944390                                mom = np.inner(np.array(atom[cmx:cmx+3]),Bmat)
    4395                                 if item[2] < 0:
    4396                                     atom[cmx:cmx+3] = np.inner(np.inner(mom,-M),Amat)*nl.det(M)*SpnFlp[opNum-1]
    4397                                 else:
    4398                                     atom[cmx:cmx+3] = np.inner(np.inner(mom,M),Amat)*nl.det(M)*SpnFlp[opNum-1]
     4391                                atom[cmx:cmx+3] = np.inner(np.inner(mom,M),Amat)*nl.det(M)*SpnFlp[opNum-1]
    43994392                            atom[cs-1] = str(item[2])+'+'
    44004393                            atom[cuij:cuij+6] = item[1]
     
    44524445                            opNum = G2spc.GetOpNum(OprNum,SGData)
    44534446                            mom = np.inner(np.array(atom[cmx:cmx+3]),Bmat)
    4454                             if Inv:
    4455                                 atom[cmx:cmx+3] = np.inner(np.inner(mom,-M),Amat)*nl.det(M)*SpnFlp[opNum-1]
    4456                             else:
    4457                                 atom[cmx:cmx+3] = np.inner(np.inner(mom,M),Amat)*nl.det(M)*SpnFlp[opNum-1]
     4447                            atom[cmx:cmx+3] = np.inner(np.inner(mom,M),Amat)*nl.det(M)*SpnFlp[opNum-1]
    44584448                        atomOp = atom[cs-1]
    44594449                        newOp = str(((Opr+1)+100*Cent)*(1-2*Inv))+'+'+ \
     
    45374527            SGData = generalData['SGData']
    45384528            SpnFlp = SGData.get('SpnFlp',[])
     4529            MagMom = SGData.get('MagMom',[])
    45394530            wx.BeginBusyCursor()
    45404531            try:
     
    45524543                            opNum = G2spc.GetOpNum(item[2],SGData)
    45534544                            mom = np.inner(np.array(atom[cmx:cmx+3]),Bmat)
    4554                             if item[2] < 0:
    4555                                 atom[cmx:cmx+3] = np.inner(np.inner(mom,-M),Amat)*nl.det(M)*SpnFlp[opNum-1]
    4556                             else:
    4557                                 atom[cmx:cmx+3] = np.inner(np.inner(mom,M),Amat)*nl.det(M)*SpnFlp[opNum-1]
     4545                            atom[cmx:cmx+3] = np.inner(np.inner(mom,M),Amat)*nl.det(M)*SpnFlp[opNum-1]
    45584546                        atom[cs-1] = str(item[2])+'+' \
    45594547                            +str(item[3][0])+','+str(item[3][1])+','+str(item[3][2])
Note: See TracChangeset for help on using the changeset viewer.