Changeset 3577


Ignore:
Timestamp:
Aug 30, 2018 4:49:24 PM (5 years ago)
Author:
vondreele
Message:

delete MAXMAGN - not using it any more
fix (maybe) mag constraints; works for orthogonal axes
fix bug in ReadPowerInstprm?
spelling of Levenberg!
fixes to mag structure from Bilbao stuff
set dmin to be min 1A for UnitCell? show reflections
fix Latt2text to show '-' & Trans2text
fix ext test in GetHistogramPhaseData? for mag structures

Location:
trunk
Files:
1 deleted
7 edited

Legend:

Unmodified
Added
Removed

  • trunk/GSASIIconstrGUI.py

    r3567 r3577  
    12891289    detTrans = np.abs(nl.det(Trans))
    12901290    invTrans = nl.inv(Trans)
    1291 #    print 'invTrans',invTrans
    12921291    nAcof = G2lat.cell2A(newPhase['General']['Cell'][1:7])
    12931292   
    … …  
    13231322        if Nop < 0:         #inversion
    13241323            Opr *= -1
    1325         XOpr = np.inner(Opr,invTrans)
    1326 #        for ix,name in enumerate(xnames):
     1324        XOpr = np.inner(Opr,Trans)
    13271325        for ix in list(set(CSX[0])):
    13281326            if not ix:
    … …  
    13351333                    DepCons.append([opval,G2obj.G2VarObj('%d::%s:%s'%(opId,xnames[iop],iat))])
    13361334            if len(DepCons) == 1:
    1337                 constraints['Phase'].append([IndpCon,DepCons[0],None,None,'e'])
     1335                constraints['Phase'].append([DepCons[0],IndpCon,None,None,'e'])
    13381336            elif len(DepCons) > 1:
    13391337                for Dep in DepCons:
    … …  
    13431341            IndpCon = [1.0,G2obj.G2VarObj('%d::%s:%d'%(npId,name,ia))]
    13441342            DepCons = [1.0,G2obj.G2VarObj('%d::%s:%s'%(opId,name,iat))]
    1345             constraints['Phase'].append([IndpCon,DepCons,None,None,'e'])
     1343            constraints['Phase'].append([DepCons,IndpCon,None,None,'e'])
    13461344#        for iu,Uid in enumerate(Uids):
    13471345        for iu in list(set(CSU[0])):
    … …  
    13591357                        DepCons.append([Trans[ibu,iau]/detTrans,G2obj.G2VarObj(parm)])
    13601358            if len(DepCons) == 1:
    1361                 constraints['Phase'].append([IndpCon,DepCons[0],None,None,'e'])
     1359                constraints['Phase'].append([DepCons[0],IndpCon,None,None,'e'])
    13621360            elif len(DepCons) > 1:       
    13631361                for Dep in DepCons:
    … …  
    13791377        for iat in range(3):
    13801378            if nSGData['SGLaue'] in ['-1','2/m']:       #set holds
    1381                 if (abs(nAcof[iA]) < 1.e-8) and (abs(Trans[Aid[0],Aid[1]]) < 1.e-8):
     1379                if (abs(nAcof[iA]) < 1.e-8) and (abs(invTrans[Aid[0],Aid[1]]) < 1.e-8):
    13821380                    if Axes[iat] != nSGData['SGUniq'] and nSGData['SGLaue'] != oSGData['SGLaue']:
    13831381                        HoldObj = G2obj.G2VarObj('%d::%s'%(npId,Aid[2]))
    … …  
    13871385                            continue
    13881386            for ibt in range(3):
    1389                 if abs(Trans[Aid[0],iat]) > 1.e-4 and abs(Trans[Aid[1],ibt]) > 1.e-4 and abs(oAcof[iA]) > 1.e-8:
     1387                if abs(Trans[Aid[0],iat]) > 1.e-4 and abs(invTrans[Aid[1],ibt]) > 1.e-4 and abs(oAcof[iA]) > 1.e-8:
    13901388                    parm = SetUniqAj(opId,Anames[ibt][iat],nSGData['SGLaue'])
    13911389                    parmDict[parm] = oAcof[As.index(Anames[ibt][iat])]
    13921390                    if not parm in varyList:
    13931391                        varyList.append(parm)
    1394                     DepCons.append([Trans[ibt,iat],G2obj.G2VarObj(parm)])
     1392                    DepCons.append([invTrans[ibt,iat]*invTrans[iat,ibt],G2obj.G2VarObj(parm)])
    13951393        if len(DepCons) == 1:
    13961394            constraints['Phase'].append([IndpCon,DepCons[0],None,None,'e'])
    … …  
    14101408        IndpCon = [1.0,G2obj.G2VarObj(ohapkey+'Scale')]
    14111409        DepCons = [detTrans,G2obj.G2VarObj(nhapkey+'Scale')]
    1412         constraints['HAP'].append([IndpCon,DepCons,None,None,'e'])
     1410        constraints['HAP'].append([DepCons,IndpCon,None,None,'e'])
    14131411        for name in ['Size;i','Mustrain;i']:
    14141412            IndpCon = [1.0,G2obj.G2VarObj(ohapkey+name)]

  • trunk/GSASIIdataGUI.py

    r3571 r3577  
    12451245            if len(banklist) > 1: # yes, the user must make a selection
    12461246                choices = [str(i) for i in banklist]
    1247                 bank = int(G2G.ItemSelector(choices,self.G2frame,multiple=False))
     1247                bank = int(G2G.ItemSelector(choices,self,multiple=False))
    12481248            else:
    12491249                bank = 1

  • trunk/GSASIImath.py

    r3568 r3577  
    158158    nfev = 0
    159159    if Print:
    160         print (' Hessian Levenburg-Marquardt SVD refinement on %d variables:'%(n))
     160        print (' Hessian Levenberg-Marquardt SVD refinement on %d variables:'%(n))
    161161    Lam = np.zeros((n,n))
    162162    while icycle < maxcyc:

  • trunk/GSASIIphsGUI.py

    r3575 r3577  
    308308        self.Mtrans = False
    309309        self.kvec = [0.,0.,0.]
    310         self.Bilbao = ''
    311310        self.Draw()
    312311
    … …  
    527526        self.Phase['General']['SGData'] = self.SGData
    528527        if self.ifMag:
    529             self.Phase['General']['Name'] += ' mag: '+self.Bilbao
     528            self.Phase['General']['Name'] += ' mag: '
    530529        else:
    531530            self.Phase['General']['Name'] += ' %s'%(self.Common)
    … …  
    23692368                if dlg.ShowModal() == wx.ID_OK:
    23702369                    newPhase,Trans,Uvec,Vvec,ifMag,ifConstr,Common = dlg.GetSelection()
     2370                    newPhase['ranId'] = ran.randint(0,sys.maxsize),
    23712371                    if ifMag:
    23722372                        BNSlatt = newPhase['General']['SGData']['BNSlattsym'][0]
    … …  
    24282428        G2frame.GPXtree.SetItemPyData(sub,newPhase)
    24292429        newPhase['Drawing'] = []
    2430        
     2430        print(newPhase)
    24312431        if ifConstr:
    24322432            G2cnstG.TransConstraints(G2frame,data,newPhase,Trans,Vvec,atCodes)     #data is old phase
    … …  
    24542454            sel = dlg.GetSelection()
    24552455            magchoice = magKeep[sel]
    2456             phaseName = '%s mag (%d)'%(data['General']['Name'],sel)
     2456            phaseName = '%s mag_%d'%(data['General']['Name'],sel)
    24572457            newPhase = copy.deepcopy(data)
     2458            newPhase['ranId'] = ran.randint(0,sys.maxsize),
    24582459            del newPhase['magPhases']
    24592460            generalData = newPhase['General']
    24602461            generalData['SGData'] = copy.deepcopy(magchoice['SGData'])           
    2461             generalData['Cell'][1:] = magchoice['Cell']
     2462            generalData['Cell'][1:] = magchoice['Cell'][:]
    24622463            SGData = generalData['SGData']
    24632464            vvec = np.array([0.,0.,0.])

  • trunk/GSASIIpwdGUI.py

    r3575 r3577  
    27652765    else:
    27662766        difC = Inst['difC'][1]
    2767         dmin = G2lat.Pos2dsp(Inst,Limits[0])
     2767        dmin = max(1.0,G2lat.Pos2dsp(Inst,Limits[0]))
    27682768   
    27692769    def SetLattice(controls):
    … …  
    30693069        Cell = Phase['General']['Cell']
    30703070        SGData = Phase['General']['SGData']
     3071        if 'mono' in SGData['SGSys']:
     3072            SpGrp = G2spc.fixMono(SGData['SpGrp'])
     3073            if SpGrp == None:
     3074                wx.MessageBox('Unusable space group',caption='Monoclinic '+SGData['SpGrp']+' not usable here',style=wx.ICON_EXCLAMATION)
     3075                return
     3076        controls[13] = SpGrp
    30713077        controls[4] = 1
    30723078        controls[5] = (SGData['SGLatt']+SGData['SGLaue']).replace('-','')
     3079        if controls[5][1:] == 'm3': controls[5] += 'm'
    30733080        if 'R' in controls[5]: controls[5] = 'R3-H'
    30743081        controls[6:13] = Cell[1:8]
    3075         controls[13] = SGData['SpGrp']
     3082       
    30763083        if 'N' in Inst['Type'][0]:
    30773084            G2frame.dataWindow.RunSubGroupsMag.Enable(True)
    … …  
    30813088       
    30823089    def ImportUnitCell(event):
    3083         controls,bravais,cells,dminx,ssopt = G2frame.GPXtree.GetItemPyData(UnitCellsId)
     3090        controls,bravais,cells,dminx,ssopt = G2frame.GPXtree.GetItemPyData(UnitCellsId)[:5]
    30843091        reqrdr = G2frame.dataWindow.ReImportMenuId.get(event.GetId())
    30853092        rdlist = G2frame.OnImportGeneric(reqrdr,
    … …  
    30893096        Cell = rd.Phase['General']['Cell']
    30903097        SGData = rd.Phase['General']['SGData']
     3098        if '1 1' in SGData['SpGrp']:
     3099            wx.MessageBox('Unusable space group',caption='Monoclinic '+SGData['SpGrp']+' not usable here',style=wx.ICON_EXCLAMATION)
     3100            return
    30913101        controls[4] = 1
    30923102        controls[5] = (SGData['SGLatt']+SGData['SGLaue']).replace('-','')
     3103        if controls[5][1:] == 'm3': controls[5] += 'm'
    30933104        if 'R' in controls[5]: controls[5] = 'R3-H'
    30943105        controls[6:13] = Cell[1:8]

  • trunk/GSASIIspc.py

    r3575 r3577  
    594594                txt += '0,'
    595595                continue
     596            if icen < 0:
     597                txt += '-'
     598                icen *= -1
    596599            for mul,prod,frac in zip(mulList,prodList,fracList):
    597600                if abs(icen*mul-prod) < 1.e-5:
    … …  
    10861089    return ' '.join(magSym)
    10871090
     1091def fixMono(SpGrp):
     1092    'fixes b-unique monoclinics in e.g. P 1 2/1c 1 --> P 21/c '
     1093    Flds = SpGrp.split()
     1094    if len(Flds) == 4:
     1095        if Flds[2] != '1':
     1096            return '%s %s'%(Flds[0],Flds[2])
     1097        else:
     1098            return None
     1099    else:
     1100        return SpGrp
     1101
    10881102def Trans2Text(Trans):
    10891103    "from transformation matrix to text"
    … …  
    10911105    Text = ''
    10921106    for row in Trans:
     1107        Fld = ''
    10931108        for i in [0,1,2]:
    10941109            if row[i]:
    1095                 Text += str(int(row[i]))+cells[i]
     1110                if Fld and row[i] > 0.:
     1111                    Fld += '+'
     1112                Fld += '%3.1f'%(row[i])+cells[i]
     1113        Text += Fld
    10961114        Text += ','
    1097         Text = Text.replace('1','')
     1115        Text = Text.replace('1.0','').replace('.0','').replace('0.5','1/2')
    10981116    return Text[:-1]
    10991117
    … …  
    36033621'''
    36043622altSettingsOrtho = {
    3605         'P 2 2 21' :{'abc':'P 2 2 21','cab':'P 21 2 2','bca':'P 2 21 2','a-cb':'P 2 21 2','ba-c':'P 2 2 21','-cba':'P 21 2 2'},
    3606         'P 21 21 2':{'abc':'P 2121 2','cab':'P 2 21 21','bca':'P 21 2 21','a-cb':'P 21 2 21','ba-c':'P 21 21 2','-cba':'P 2 21 21'},
    3607         'C 2 2 21':{'abc':'C 2 2 21','cab':'A 21 2 2','bca':'B 2 21 2','a-cb':'B 2 21 2','ba-c':'C 2 2 21','-cba':'A 21 2 2'},
    3608         'C 2 2 2':{'abc':'C 2 2 2','cab':'A 2 2 2','bca':'B 2 2 2','a-cb':'B 2 2 2','ba-c':'C 2 2 2','-cba':'A 2 2 2'},
    3609         'P m m 2':{'abc':'P m m 2','cab':'P 2 m m','bca':'P m 2 m','a-cb':'P m 2 m','ba-c':'P m m 2','-cba':'P 2 m m'},
    3610         'P m c 21':{'abc':'P m c 21','cab':'P 21 m a','bca':'P b 21 m','a-cb':'P m 21 b','ba-c':'P c m 21','-cba':'P 21 a m'},
    3611         'P c c 2':{'abc':'P c c 2','cab':'P 2 a a','bca':'P b 2 b','a-cb':'P b 2 b','ba-c':'P c c 2','-cba':'P 2 a a'},
    3612         'P m a 2':{'abc':'P m a 2','cab':'P 2 m b','bca':'P c 2 m','a-cb':'P m 2 a','ba-c':'P b m 2','-cba':'P 2 c m'},
    3613         'P c a 21':{'abc':'P c a 21','cab':'P 21 a b','bca':'P c 21 b','a-cb':'P b 21 a','ba-c':'P b c 21','-cba':'P 21 c a'},
    3614         'P n c 2':{'abc':'P n c 2','cab':'P 2 n a','bca':'P b 2 n','a-cb':'P n 2 b','ba-c':'P c n 2','-cba':'P 2 a n'},
    3615         'P m n 21':{'abc':'P m n 21','cab':'P 21 m n','bca':'P n 21 m','a-cb':'P m 21 n','ba-c':'P n m 21','-cba':'P 21 n m'},
    3616         'P b a 2':{'abc':'P b a 2','cab':'P 2 c b','bca':'P c 2 a','a-cb':'P c 2 a','ba-c':'P b a 2','-cba':'P 2 c b'},
    3617         'P n a 21':{'abc':'P n a 21','cab':'P 21 n b','bca':'P c 21 n','a-cb':'P n 21 a','ba-c':'P b n 21','-cba':'P 21 c n'},
    3618         'P n n 2':{'abc':'P n n 2','cab':'P 2 n n','bca':'P n 2 n','a-cb':'P n 2 n','ba-c':'P n n 2','-cba':'P 2 n n'},
    3619         'C m m 2':{'abc':'C m m 2','cab':'A 2 m m','bca':'B m 2 m','a-cb':'B m 2 m','ba-c':'C m m 2','-cba':'A 2 m m'},
    3620         'C m c 21':{'abc':'C m c 21','cab':'A 21 m a','bca':'B b 21 m','a-cb':'B m 21 b','ba-c':'C c m 21','-cba':'A 21 a m'},
    3621         'C c c 2':{'abc':'C c c 2','cab':'A 2 a a','bca':'B b 2 b','a-cb':'B b 2 b','ba-c':'C c c 2','-cba':'A 2 a a'},
    3622         'A m m 2':{'abc':'A m m 2','cab':'B 2 m m','bca':'C m 2 m','a-cb':'A m 2 m','ba-c':'B m m 2','-cba':'C 2 m m'},
    3623         'A b m 2':{'abc':'A b m 2','cab':'B 2 c m','bca':'C m 2 a','a-cb':'A c 2 m','ba-c':'B m a 2','-cba':'C 2 m b'},
    3624         'A m a 2':{'abc':'A m a 2','cab':'B 2 m b','bca':'C c 2 m','a-cb':'A m 2 a','ba-c':'B b m 2','-cba':'C 2 c m'},
    3625         'A b a 2':{'abc':'A b a 2','cab':'B 2 c b','bca':'C c 2 a','a-cb':'A c 2 a','ba-c':'B b a 2','-cba':'C 2 c b'},
    3626         'F m m 2':{'abc':'F m m 2','cab':'F 2 m m','bca':'F m 2 m','a-cb':'F m 2 m','ba-c':'F m m 2','-cba':'F 2 m m'},
    3627         'F d d 2':{'abc':'F d d 2','cab':'F 2 d d','bca':'F d 2 d','a-cb':'F d 2 d','ba-c':'F d d 2','-cba':'F 2 d d'},
    3628         'I m m 2':{'abc':'I m m 2','cab':'I 2 m m','bca':'I m 2 m','a-cb':'I m 2 m','ba-c':'I m m 2','-cba':'I 2 m m'},
    3629         'I b a 2':{'abc':'I b a 2','cab':'I 2 c b','bca':'I c 2 a','a-cb':'I c 2 a','ba-c':'I b a 2','-cba':'I 2 c b'},
    3630         'I m a 2':{'abc':'I m a 2','cab':'I 2 m b','bca':'I c 2 m','a-cb':'I m 2 a','ba-c':'I b m 2','-cba':'I 2 c m'},
    3631         'P c c m':{'abc':'P c c m','cab':'P m a a','bca':'P b m b','a-cb':'P b m b','ba-c':'P c c m','-cba':'P m a a'},
    3632         'P b a n':{'abc':'P b a n','cab':'P n c b','bca':'P c n a','a-cb':'P c n a','ba-c':'P b a n','-cba':'P n c b'},
    3633         'P m m a':{'abc':'P m m a','cab':'P b m m','bca':'P m c m','a-cb':'P m a m','ba-c':'P m m b','-cba':'P c m m'},
    3634         'P n n a':{'abc':'P n n a','cab':'P b n n','bca':'P n c n','a-cb':'P n a n','ba-c':'P n n b','-cba':'P c n n'},
    3635         'P m n a':{'abc':'P m n a','cab':'P b m n','bca':'P n c m','a-cb':'P m a n','ba-c':'P n m b','-cba':'P c n m'},
    3636         'P c c a':{'abc':'P c c a','cab':'P b a a','bca':'P b c b','a-cb':'P b a b','ba-c':'P c c b','-cba':'P c a a'},
    3637         'P b a m':{'abc':'P b a m','cab':'P m c b','bca':'P c m a','a-cb':'P c m a','ba-c':'P b a m','-cba':'P m c b'},
    3638         'P c c n':{'abc':'P c c n','cab':'P n a a','bca':'P b n b','a-cb':'P b n b','ba-c':'P c c n','-cba':'P n a a'},
    3639         'P b c m':{'abc':'P b c m','cab':'P m c a','bca':'P b m a','a-cb':'P c m b','ba-c':'P c a m','-cba':'P m a b'},
    3640         'P n n m':{'abc':'P n n m','cab':'P m n n','bca':'P n m n','a-cb':'P n m n','ba-c':'P n n m','-cba':'P m n n'},
    3641         'P m m n':{'abc':'P m m n','cab':'P n m m','bca':'P m n m','a-cb':'P m n m','ba-c':'P m m n','-cba':'P n m m'},
    3642         'P b c n':{'abc':'P b c n','cab':'P n c a','bca':'P b n a','a-cb':'P c n b','ba-c':'P c a n','-cba':'P n a b'},
    3643         'P b c a':{'abc':'P b c a','cab':'P b c a','bca':'P b c a','a-cb':'P c a b','ba-c':'P c a b','-cba':'P c a b'},
    3644         'P n m a':{'abc':'P n m a','cab':'P b n m','bca':'P m c n','a-cb':'P n a m','ba-c':'P m n b','-cba':'P c m n'},
    3645         'C m c m':{'abc':'C m c m','cab':'A m m a','bca':'B b m m','a-cb':'B m m b','ba-c':'C c m m','-cba':'A m a m'},
    3646         'C m c a':{'abc':'C m c a','cab':'A b m a','bca':'B b c m','a-cb':'B m a b','ba-c':'C c m b','-cba':'A c a m'},
    3647         'C m m m':{'abc':'C m m m','cab':'A m m m','bca':'B m m m','a-cb':'B m m m','ba-c':'C m m m','-cba':'A m m m'},
    3648         'C c c m':{'abc':'C c c m','cab':'A m a a','bca':'B b m b','a-cb':'B b m b','ba-c':'C c c m','-cba':'A m a a'},
    3649         'C m m a':{'abc':'C m m a','cab':'A b m m','bca':'B m c m','a-cb':'B m a m','ba-c':'C m m b','-cba':'A c m m'},
    3650         'C c c a':{'abc':'C c a a','cab':'A b a a','bca':'B b c b','a-cb':'B b a b','ba-c':'C c c b','-cba':'A c a a'},
    3651         'I b a m':{'abc':'I b a m','cab':'I m c b','bca':'I c m a','a-cb':'I c m a','ba-c':'I b a m','-cba':'I m c b'},
    3652         'I b c a':{'abc':'I b c a','cab':'I b c a','bca':'I b c a','a-cb':'I c a b','ba-c':'I c a b','-cba':'I c a b'},
    3653         'I m m a':{'abc':'I m m a','cab':'I b m m','bca':'I m c m','a-cb':'I m a m','ba-c':'I m m  b','-cba':'I c m m'},
     3623        'P 2 2 21' :{'abc':'P 2 2 21','cab':'P 21 2 2','bca':'P 2 21 2','acb':'P 2 21 2','bac':'P 2 2 21','cba':'P 21 2 2'},
     3624        'P 21 21 2':{'abc':'P 2121 2','cab':'P 2 21 21','bca':'P 21 2 21','acb':'P 21 2 21','bac':'P 21 21 2','cba':'P 2 21 21'},
     3625        'C 2 2 21':{'abc':'C 2 2 21','cab':'A 21 2 2','bca':'B 2 21 2','acb':'B 2 21 2','bac':'C 2 2 21','cba':'A 21 2 2'},
     3626        'C 2 2 2':{'abc':'C 2 2 2','cab':'A 2 2 2','bca':'B 2 2 2','acb':'B 2 2 2','bac':'C 2 2 2','cba':'A 2 2 2'},
     3627        'P m m 2':{'abc':'P m m 2','cab':'P 2 m m','bca':'P m 2 m','acb':'P m 2 m','bac':'P m m 2','cba':'P 2 m m'},
     3628        'P m c 21':{'abc':'P m c 21','cab':'P 21 m a','bca':'P b 21 m','acb':'P m 21 b','bac':'P c m 21','cba':'P 21 a m'},
     3629        'P c c 2':{'abc':'P c c 2','cab':'P 2 a a','bca':'P b 2 b','acb':'P b 2 b','bac':'P c c 2','cba':'P 2 a a'},
     3630        'P m a 2':{'abc':'P m a 2','cab':'P 2 m b','bca':'P c 2 m','acb':'P m 2 a','bac':'P b m 2','cba':'P 2 c m'},
     3631        'P c a 21':{'abc':'P c a 21','cab':'P 21 a b','bca':'P c 21 b','acb':'P b 21 a','bac':'P b c 21','cba':'P 21 c a'},
     3632        'P n c 2':{'abc':'P n c 2','cab':'P 2 n a','bca':'P b 2 n','acb':'P n 2 b','bac':'P c n 2','cba':'P 2 a n'},
     3633        'P m n 21':{'abc':'P m n 21','cab':'P 21 m n','bca':'P n 21 m','acb':'P m 21 n','bac':'P n m 21','cba':'P 21 n m'},
     3634        'P b a 2':{'abc':'P b a 2','cab':'P 2 c b','bca':'P c 2 a','acb':'P c 2 a','bac':'P b a 2','cba':'P 2 c b'},
     3635        'P n a 21':{'abc':'P n a 21','cab':'P 21 n b','bca':'P c 21 n','acb':'P n 21 a','bac':'P b n 21','cba':'P 21 c n'},
     3636        'P n n 2':{'abc':'P n n 2','cab':'P 2 n n','bca':'P n 2 n','acb':'P n 2 n','bac':'P n n 2','cba':'P 2 n n'},
     3637        'C m m 2':{'abc':'C m m 2','cab':'A 2 m m','bca':'B m 2 m','acb':'B m 2 m','bac':'C m m 2','cba':'A 2 m m'},
     3638        'C m c 21':{'abc':'C m c 21','cab':'A 21 m a','bca':'B b 21 m','acb':'B m 21 b','bac':'C c m 21','cba':'A 21 a m'},
     3639        'C c c 2':{'abc':'C c c 2','cab':'A 2 a a','bca':'B b 2 b','acb':'B b 2 b','bac':'C c c 2','cba':'A 2 a a'},
     3640        'A m m 2':{'abc':'A m m 2','cab':'B 2 m m','bca':'C m 2 m','acb':'A m 2 m','bac':'B m m 2','cba':'C 2 m m'},
     3641        'A b m 2':{'abc':'A b m 2','cab':'B 2 c m','bca':'C m 2 a','acb':'A c 2 m','bac':'B m a 2','cba':'C 2 m b'},
     3642        'A m a 2':{'abc':'A m a 2','cab':'B 2 m b','bca':'C c 2 m','acb':'A m 2 a','bac':'B b m 2','cba':'C 2 c m'},
     3643        'A b a 2':{'abc':'A b a 2','cab':'B 2 c b','bca':'C c 2 a','acb':'A c 2 a','bac':'B b a 2','cba':'C 2 c b'},
     3644        'F m m 2':{'abc':'F m m 2','cab':'F 2 m m','bca':'F m 2 m','acb':'F m 2 m','bac':'F m m 2','cba':'F 2 m m'},
     3645        'F d d 2':{'abc':'F d d 2','cab':'F 2 d d','bca':'F d 2 d','acb':'F d 2 d','bac':'F d d 2','cba':'F 2 d d'},
     3646        'I m m 2':{'abc':'I m m 2','cab':'I 2 m m','bca':'I m 2 m','acb':'I m 2 m','bac':'I m m 2','cba':'I 2 m m'},
     3647        'I b a 2':{'abc':'I b a 2','cab':'I 2 c b','bca':'I c 2 a','acb':'I c 2 a','bac':'I b a 2','cba':'I 2 c b'},
     3648        'I m a 2':{'abc':'I m a 2','cab':'I 2 m b','bca':'I c 2 m','acb':'I m 2 a','bac':'I b m 2','cba':'I 2 c m'},
     3649        'P c c m':{'abc':'P c c m','cab':'P m a a','bca':'P b m b','acb':'P b m b','bac':'P c c m','cba':'P m a a'},
     3650        'P b a n':{'abc':'P b a n','cab':'P n c b','bca':'P c n a','acb':'P c n a','bac':'P b a n','cba':'P n c b'},
     3651        'P m m a':{'abc':'P m m a','cab':'P b m m','bca':'P m c m','acb':'P m a m','bac':'P m m b','cba':'P c m m'},
     3652        'P n n a':{'abc':'P n n a','cab':'P b n n','bca':'P n c n','acb':'P n a n','bac':'P n n b','cba':'P c n n'},
     3653        'P m n a':{'abc':'P m n a','cab':'P b m n','bca':'P n c m','acb':'P m a n','bac':'P n m b','cba':'P c n m'},
     3654        'P c c a':{'abc':'P c c a','cab':'P b a a','bca':'P b c b','acb':'P b a b','bac':'P c c b','cba':'P c a a'},
     3655        'P b a m':{'abc':'P b a m','cab':'P m c b','bca':'P c m a','acb':'P c m a','bac':'P b a m','cba':'P m c b'},
     3656        'P c c n':{'abc':'P c c n','cab':'P n a a','bca':'P b n b','acb':'P b n b','bac':'P c c n','cba':'P n a a'},
     3657        'P b c m':{'abc':'P b c m','cab':'P m c a','bca':'P b m a','acb':'P c m b','bac':'P c a m','cba':'P m a b'},
     3658        'P n n m':{'abc':'P n n m','cab':'P m n n','bca':'P n m n','acb':'P n m n','bac':'P n n m','cba':'P m n n'},
     3659        'P m m n':{'abc':'P m m n','cab':'P n m m','bca':'P m n m','acb':'P m n m','bac':'P m m n','cba':'P n m m'},
     3660        'P b c n':{'abc':'P b c n','cab':'P n c a','bca':'P b n a','acb':'P c n b','bac':'P c a n','cba':'P n a b'},
     3661        'P b c a':{'abc':'P b c a','cab':'P b c a','bca':'P b c a','acb':'P c a b','bac':'P c a b','cba':'P c a b'},
     3662        'P n m a':{'abc':'P n m a','cab':'P b n m','bca':'P m c n','acb':'P n a m','bac':'P m n b','cba':'P c m n'},
     3663        'C m c m':{'abc':'C m c m','cab':'A m m a','bca':'B b m m','acb':'B m m b','bac':'C c m m','cba':'A m a m'},
     3664        'C m c a':{'abc':'C m c a','cab':'A b m a','bca':'B b c m','acb':'B m a b','bac':'C c m b','cba':'A c a m'},
     3665        'C m m m':{'abc':'C m m m','cab':'A m m m','bca':'B m m m','acb':'B m m m','bac':'C m m m','cba':'A m m m'},
     3666        'C c c m':{'abc':'C c c m','cab':'A m a a','bca':'B b m b','acb':'B b m b','bac':'C c c m','cba':'A m a a'},
     3667        'C m m a':{'abc':'C m m a','cab':'A b m m','bca':'B m c m','acb':'B m a m','bac':'C m m b','cba':'A c m m'},
     3668        'C c c a':{'abc':'C c a a','cab':'A b a a','bca':'B b c b','acb':'B b a b','bac':'C c c b','cba':'A c a a'},
     3669        'I b a m':{'abc':'I b a m','cab':'I m c b','bca':'I c m a','acb':'I c m a','bac':'I b a m','cba':'I m c b'},
     3670        'I b c a':{'abc':'I b c a','cab':'I b c a','bca':'I b c a','acb':'I c a b','bac':'I c a b','cba':'I c a b'},
     3671        'I m m a':{'abc':'I m m a','cab':'I b m m','bca':'I m c m','acb':'I m a m','bac':'I m m  b','cba':'I c m m'},
    36543672        }
    36553673spg2origins = {}

  • trunk/GSASIIstrIO.py

    r3573 r3577  
    23872387                        for h,k,l,d in HKLd:
    23882388                            ext,mul,uniq,phi = G2spc.GenHKLf([h,k,l],SGData)
    2389                             if 'N' in inst['Type'][0] and  'MagSpGrp' in SGData:
     2389                            if ext and 'N' in inst['Type'][0] and  'MagSpGrp' in SGData:
    23902390                                ext = G2spc.checkMagextc([h,k,l],SGData)
    23912391                            mul *= 2      # for powder overlap of Friedel pairs
Note: See TracChangeset for help on using the changeset viewer.