Changeset 1615 for trunk/GSASIIspc.py


Ignore:
Timestamp:
Dec 21, 2014 8:12:54 AM (8 years ago)
Author:
vondreele
Message:

more work on incommensurate wave input to LS
fix FillAtomLookup? problem in G2restrGUI

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIspc.py

    r1613 r1615  
    14231423   
    14241424def GetSSfxuinel(XYZ,UIJ,SGData,SSGData):
    1425     CSI = {'Sfrac':[[0,1],[1.,1.]],'Spos':[[0,1,2, 0,1,2],[1.,1.,1., 1.,1.,1.]],    #sin & cos
    1426         'Sadp':[[0,1,2,3,4,5, 0,1,2,3,4,5],[1.,1.,1.,1.,1.,1., 1.,1.,1.,1.,1.,1.]],
    1427         'Smag':[[0,1,2, 0,1,2],[1.,1.,1., 1.,1.,1.]]}
     1425    CSI = {'Sfrac':[[1,2],[1.,1.]],'Spos':[[1,2,3, 4,5,6],[1.,1.,1., 1.,1.,1.]],    #sin & cos
     1426        'Sadp':[[1,2,3,4,5,6, 7,8,9,10,11,12],[1.,1.,1.,1.,1.,1., 1.,1.,1.,1.,1.,1.]],
     1427        'Smag':[[1,2,3, 4,5,6],[1.,1.,1., 1.,1.,1.]]}
    14281428    deltx = np.ones((3,4))*.01
    14291429    deltx[:3,:3] = np.eye((3))*.001
     
    14531453        return CSI
    14541454    elif siteSym == '-1':   #"-1" site symmetry
    1455         CSI['Spos'][0] = [0,1,2, -1,-1,-1]
    1456         CSI['Sadp'][0] = [-1,-1,-1,-1,-1,-1, 0,1,2,3,4,5,]
     1455        CSI['Spos'][0] = [1,2,3, 0,0,0]
     1456        CSI['Sadp'][0] = [0,0,0,0,0,0, 1,2,3,4,5,6]
    14571457        return CSI       
    14581458    print siteSym,OpText,SSOptext   
     
    14621462            if 's' in SSGData['SSpGrp'].split('(')[1]:
    14631463                if siteSym == 'm':
    1464                     CSI['Spos'][0] = [0,-1,1, -1,1,-1]
    1465                     CSI['Sadp'][0] = [0,1,2,-1,3,-1, -1,-1,-1,4,-1,5]
     1464                    CSI['Spos'][0] = [1,0,2, 0,3,0]
     1465                    CSI['Sadp'][0] = [1,2,3,0,4,0, 0,0,0,5,0,6]
    14661466                elif siteSym == '2/m':
    1467                     CSI['Spos'][0] = [0,-1,1, -1,-1,-1]
    1468                     CSI['Sadp'][0] = [-1,-1,-1,-1,-1,-1, -1,-1,-1,0,-1,1]
     1467                    CSI['Spos'][0] = [1,0,2, 0,0,0]
     1468                    CSI['Sadp'][0] = [0,0,0,0,0,0, 0,0,0,1,0,2]
    14691469                elif siteSym == '2':
    1470                     CSI['Spos'][0] = [0,-1,1, 2,-1,3]               
    1471                     CSI['Sadp'][0] = [-1,-1,-1,0,-1,1, -1,-1,-1,2,-1,3]
     1470                    CSI['Spos'][0] = [1,0,2, 3,0,4]               
     1471                    CSI['Sadp'][0] = [0,0,0,1,0,2, 0,0,0,3,0,4]
    14721472            elif '1/2' in SSGData['modSymb']:    #e.g. (0b1/2)
    14731473                if siteSym == 'm':
    1474                     CSI['Spos'][0] = [-1,-1,-1, -1,-1,-1]
    1475                     CSI['Sadp'][0] = [-1,-1,-1,-1,-1,-1, -1,-1,-1,-1,-1,-1]
     1474                    CSI['Spos'][0] = [0,0,0, 0,0,0]
     1475                    CSI['Sadp'][0] = [0,0,0,0,0,0, 0,0,0,0,0,0]
    14761476                elif siteSym == '2/m':
    1477                     CSI['Spos'][0] = [0,-1,1, -1,-1,-1]
    1478                     CSI['Sadp'][0] = [-1,-1,-1,-1,-1,-1, 0,1,2,-1,3,-1]
     1477                    CSI['Spos'][0] = [1,0,2, 0,0,0]
     1478                    CSI['Sadp'][0] = [0,0,0,0,0,0, 1,2,3,0,4,0]
    14791479                elif siteSym == '2':
    1480                     CSI['Spos'][0] = [0,-1,1, 2,-1,3]
    1481                     CSI['Sadp'][0] = [-1,-1,-1,0,-1,1, -1,-1,-1,2,-1,3]
     1480                    CSI['Spos'][0] = [1,0,2, 3,0,4]
     1481                    CSI['Sadp'][0] = [0,0,0,1,0,2, 0,0,0,3,0,4]
    14821482            else:
    14831483                if siteSym == 'm':
    1484                     CSI['Spos'][0] = [-1,0,-1, 1,-1,2]
    1485                     CSI['Sadp'][0] = [-1,-1,-1,0,-1,1, 2,3,4,-1,5,-1]
     1484                    CSI['Spos'][0] = [0,1,0, 2,0,3]
     1485                    CSI['Sadp'][0] = [0,0,0,1,0,2, 3,4,5,0,6,0]
    14861486                elif siteSym == '2/m':
    1487                     CSI['Spos'][0] = [-1,0,-1, -1,-1,-1]
    1488                     CSI['Sadp'][0] = [-1,-1,-1,-1,-1,-1, 0,1,2,-1,3,-1]
     1487                    CSI['Spos'][0] = [0,1,0, 0,0,0]
     1488                    CSI['Sadp'][0] = [0,0,0,0,0,0, 1,2,3,0,4,0]
    14891489                elif siteSym == '2':
    1490                     CSI['Spos'][0] = [-1,0,-1, -1,1,-1]
    1491                     CSI['Sadp'][0] = [0,1,2,-1,3,-1, 4,5,6,-1,7,-1]
     1490                    CSI['Spos'][0] = [0,1,0, 0,2,0]
     1491                    CSI['Sadp'][0] = [1,2,3,0,4,0, 5,6,7,0,8,0]
    14921492           
    14931493        else:   #e.g. (a0g)
    14941494            if 's' in SSGData['SSpGrp'].split('(')[1]:
    14951495                if siteSym == 'm':
    1496                     CSI['Spos'][0] = [-1,0,-1, -1,1,-1]
    1497                     CSI['Sadp'][0] = [-1,-1,-1,0,-1,1, -1,-1,-1,2,-1,3]
     1496                    CSI['Spos'][0] = [0,1,0, 0,2,0]
     1497                    CSI['Sadp'][0] = [0,0,0,1,0,2, 0,0,0,3,0,4]
    14981498                elif siteSym == '2/m':
    1499                     CSI['Spos'][0] = [-1,0,-1, -1,-1,-1]
    1500                     CSI['Sadp'][0] = [-1,-1,-1,-1,-1,-1, -1,-1,-1,0,-1,1]
     1499                    CSI['Spos'][0] = [0,1,0, 0,0,0]
     1500                    CSI['Sadp'][0] = [0,0,0,0,0,0, 0,0,0,1,0,2]
    15011501                elif siteSym == '2':
    1502                     CSI['Spos'][0] = [-1,0,-1, 1,-1,2]
    1503                     CSI['Sadp'][0] = [0,1,2,-1,3,-1, -1,-1,-1,4,-1,5]
     1502                    CSI['Spos'][0] = [0,1,0, 2,0,3]
     1503                    CSI['Sadp'][0] = [1,2,3,0,4,0, 0,0,0,5,0,6]
    15041504            else:
    15051505                if siteSym == 'm':
    1506                     CSI['Spos'][0] = [0,-1,1, 2,-1,3]
    1507                     CSI['Sadp'][0] = [0,1,2,-1,3,-1, 4,5,6,-1,7,-1]
     1506                    CSI['Spos'][0] = [1,0,2, 3,0,4]
     1507                    CSI['Sadp'][0] = [1,2,3,0,4,0, 5,6,7,0,8,0]
    15081508                elif siteSym == '2/m':
    1509                     CSI['Spos'][0] = [0,-1,1, -1,-1,-1]
    1510                     CSI['Sadp'][0] = [-1,-1,-1,-1,-1,-1, 0,1,2,-1,3,-1]
     1509                    CSI['Spos'][0] = [1,0,2, 0,0,0]
     1510                    CSI['Sadp'][0] = [0,0,0,0,0,0, 1,2,3,0,4,0]
    15111511                elif siteSym == '2':
    1512                     CSI['Spos'][0] = [0,-1,1, -1,2,-1]
    1513                     CSI['Sadp'][0] = [-1,-1,-1,0,-1,1, 2,3,4,-1,5,-1]
     1512                    CSI['Spos'][0] = [1,0,2, 0,3,0]
     1513                    CSI['Sadp'][0] = [0,0,0,1,0,2, 3,4,5,0,6,0]
    15141514        Sx = CSI['Spos'][0]
    15151515        Su = CSI['Sadp'][0]
Note: See TracChangeset for help on using the changeset viewer.