Ignore:
Timestamp:
Mar 23, 2021 2:02:01 PM (7 months ago)
Author:
vondreele
Message:

add A-centered monoclinic to list of indexable lattices.
Fix a number of monoclinic indexing issues.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIlattice.py

    r4826 r4862  
    868868    ''' convert d-spacing to powder pattern position (2-theta or TOF, musec)
    869869    '''
    870     if 'T' in Inst['Type'][0] or 'PKS' in Inst['Type'][0]:
     870    if 'T' in Inst['Type'][0]:
    871871        pos = Inst['difC'][1]*dsp+Inst['Zero'][1]+Inst['difA'][1]*dsp**2+Inst.get('difB',[0,0,False])[1]/dsp
    872872    else:   #'C' or 'B'
     
    11851185            * 12 P orthorhombic
    11861186            * 13 I monoclinic
    1187             * 14 C monoclinic
    1188             * 15 P monoclinic
    1189             * 16 P triclinic
     1187            = 14 A monoclinic
     1188            * 15 C monoclinic
     1189            * 16 P monoclinic
     1190            * 17 P triclinic
    11901191           
    11911192    :param A: reciprocal metric tensor elements as [G11,G22,G33,2*G12,2*G13,2*G23]
     
    12031204                    cctbx_args['sg_type'], cctbx_args['uctbx_unit_cell'], cctbx_args['miller_index_generator'])
    12041205   
    1205     if Bravais in [9,]:
     1206    if Bravais in [9,14]:
    12061207        Cent = 'A'
    12071208    elif Bravais in [10,]:
    12081209        Cent = 'B'
    1209     elif Bravais in [11,14]:
     1210    elif Bravais in [11,15]:
    12101211        Cent = 'C'
    12111212    elif Bravais in [1,5,8,13]:
     
    12201221    dminsq = 1./(dmin**2)
    12211222    HKL = []
    1222     if Bravais == 16:                       #triclinic
     1223    if Bravais == 17:                       #triclinic
    12231224        for l in range(-Hmax[2],Hmax[2]+1):
    12241225            for k in range(-Hmax[1],Hmax[1]+1):
     
    12311232                    if 0 < rdsq <= dminsq:
    12321233                        HKL.append([h,k,l,rdsq2d(rdsq,6),-1])
    1233     elif Bravais in [13,14,15]:                #monoclinic - b unique
     1234    elif Bravais in [13,14,15,16]:                #monoclinic - b unique
    12341235        Hmax = SwapIndx(2,Hmax)
    12351236        for h in range(Hmax[0]+1):
Note: See TracChangeset for help on using the changeset viewer.