Ignore:
Timestamp:
Sep 11, 2018 4:13:58 PM (3 years ago)
Author:
vondreele
Message:

work on FindNonstandard? for mag space groups
enhance UseMagAtomDialog? to show mag site symmetry

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIlattice.py

    r3590 r3591  
    357357    cba = np.rot90(np.eye(3))
    358358    cba[0,2] *= -1      #makes -cba
    359     Mats = {'abc':abc,'cab':np.roll(abc,2,1),'bca':np.roll(abc,1,1),
    360             'acb':np.roll(cba,1,1),'bac':np.roll(cba,2,1),'cba':cba}
    361     BNS = {'A':{'abc':'A','cab':'C','bca':'B','acb':'B','bac':'C','cba':'A'},
     359    Mats = {'abc':abc,'cab':np.roll(abc,1,1),'bca':np.roll(abc,2,1),
     360            'acb':np.roll(cba,1,1),'bac':np.roll(cba,2,1),'cba':cba}        #ok
     361    BNS = {'A':{'abc':'A','cab':'C','bca':'B','acb':'A','bac':'B','cba':'C'},   
    362362           'B':{'abc':'B','cab':'A','bca':'C','acb':'C','bac':'A','cba':'B'},
    363            'C':{'abc':'C','cab':'B','bca':'A','acb':'A','bac':'B','cba':'C'},
    364            'a':{'abc':'a','cab':'c','bca':'b','acb':'b','bac':'c','cba':'a'},
     363           'C':{'abc':'C','cab':'B','bca':'A','acb':'B','bac':'C','cba':'A'},
     364           'a':{'abc':'a','cab':'c','bca':'b','acb':'a','bac':'b','cba':'c'},   #Ok
    365365           'b':{'abc':'b','cab':'a','bca':'c','acb':'c','bac':'a','cba':'b'},
    366            'c':{'abc':'c','cab':'b','bca':'a','acb':'a','bac':'b','cba':'c'},
     366           'c':{'abc':'c','cab':'b','bca':'a','acb':'b','bac':'c','cba':'a'},
    367367           'S':{'abc':'S','cab':'S','bca':'S','acb':'S','bac':'S','cba':'S'},
    368368           }
    369     Fives = {'ababc':'abc','bcbca':'cba','acacb':'acb'}
    370369    Trans = Phase['Trans']
    371370    Uvec = Phase['Uvec']
     
    377376    spn = SGData['SGSpin']
    378377    if 'ortho' in SGData['SGSys']:
    379 #        transText = G2spc.Trans2Text(nl.inv(Trans.T))
    380         transText = G2spc.Trans2Text(Trans.T)
    381         lattSym = ''
    382         for fld in transText.split(','):
    383             if 'a' in fld: lattSym += 'a'
    384             if 'b' in fld: lattSym += 'b'
    385             if 'c' in fld: lattSym += 'c'
    386         if len(lattSym) == 5:
    387             print(transText,lattSym)
    388             lattSym = Fives[lattSym]
    389 #            return None
     378        lattSym = G2spc.getlattSym(Trans)
    390379        SpGrp = SGData['SpGrp']
    391         NUvec = np.inner(np.abs(Mats[lattSym]),Uvec)
    392         NTrans = np.inner(Mats[lattSym],Trans.T)
    393         spn[1:4] = np.inner(Mats[lattSym],spn[1:4])
    394         if lattSym != 'abc' and SpGrp in G2spc.altSettingOrtho:
    395             NSG = G2spc.altSettingOrtho[SpGrp].get(lattSym,SpGrp).replace("'",'').split(' ')
     380        NUvec = np.inner(Uvec,Mats[lattSym])    #ok
     381        NTrans = np.inner(Mats[lattSym],Trans.T)        #ok
     382        spn[1:4] = np.inner(np.abs(Mats[lattSym]),spn[1:4])         #ok
     383        SGsym = G2spc.getlattSym(Mats[lattSym])
     384       
     385        if lattSym != 'abc':
     386            NSG = G2spc.altSettingOrtho[SpGrp].get(SGsym,SpGrp).replace("'",'').split(' ')
    396387            Bns = ''
    397388            if bns:
Note: See TracChangeset for help on using the changeset viewer.