Changeset 3596


Ignore:
Timestamp:
Sep 13, 2018 12:24:59 PM (3 years ago)
Author:
vondreele
Message:

fix OnTransform? for magnetic phases to give correct atom sites/symmetry

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIlattice.py

    r3594 r3596  
    356356    abc = np.eye(3)
    357357    cba = np.rot90(np.eye(3))
    358     cba[0,2] *= -1      #makes -cba
     358    cba[1,1] *= -1      #makes -cba
    359359    Mats = {'abc':abc,'cab':np.roll(abc,1,1),'bca':np.roll(abc,2,1),
    360360            'acb':np.roll(cba,1,1),'bac':np.roll(cba,2,1),'cba':cba}        #ok
  • trunk/GSASIIphsGUI.py

    r3594 r3596  
    23812381                if dlg.ShowModal() == wx.ID_OK:
    23822382                    newPhase,Trans,Uvec,Vvec,ifMag,ifConstr,Common = dlg.GetSelection()
    2383                     newPhase['ranId'] = ran.randint(0,sys.maxsize),
     2383                    newPhase['ranId'] = ran.randint(0,sys.maxsize)
     2384                    SGData = newPhase['General']['SGData']
    23842385                    if ifMag:
    2385                         BNSlatt = newPhase['General']['SGData']['BNSlattsym'][0]
     2386                        BNSlatt = SGData['BNSlattsym'][0]
     2387                       
     2388                    if not '_' in BNSlatt:
     2389                        SGData['SGSpin'] = G2spc.GetSGSpin(SGData,SGData['MagSpGrp'])
     2390                    SGData['GenSym'],SGData['GenFlg'],BNSsym = G2spc.GetGenSym(SGData)
     2391                    if '_' in BNSlatt:
     2392                        SGData['BNSlattsym'] = [BNSlatt,BNSsym[BNSlatt]]
     2393                        G2spc.ApplyBNSlatt(SGData,SGData['BNSlattsym'])
     2394                    SGData['SpnFlp'] = G2spc.GenMagOps(SGData)[1]
     2395                    SGData['MagSpGrp'] = G2spc.MagSGSym(SGData)
     2396                       
     2397                       
     2398                       
    23862399                else:
    23872400                    return
     
    24022415            else:
    24032416                phaseName = newPhase['General']['Name']
    2404                 newPhase,atCodes = G2lat.TransformPhase(data,newPhase,Trans.T,Uvec,Vvec,ifMag)
     2417               
     2418               
     2419                newPhase,atCodes = G2lat.TransformPhase(data,newPhase,Trans,Uvec,Vvec,ifMag)
    24052420                detTrans = np.abs(nl.det(Trans))
    24062421                generalData = newPhase['General']
     
    24102425                if ifMag:
    24112426                    atMxyz = []                   
    2412                     G2spc.ApplyBNSlatt(SGData,SGData['BNSlattsym'])
    2413                     SGData['GenSym'],SGData['GenFlg'],BNSsym = G2spc.GetGenSym(SGData)
    2414                     SGData['OprNames'],SGData['SpnFlp'] = G2spc.GenMagOps(SGData)
    2415                     SGData['MagSpGrp'] = G2spc.MagSGSym(SGData)
    24162427                    for atom in Atoms:
    24172428                        SytSym,Mul,Nop,dupDir = G2spc.SytSym(atom[3:6],SGData)
Note: See TracChangeset for help on using the changeset viewer.