Changeset 4791 for trunk/GSASIIindex.py


Ignore:
Timestamp:
Feb 2, 2021 9:37:33 AM (10 months ago)
Author:
toby
Message:

docs chapter numbering

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIindex.py

    r4466 r4791  
    787787    return len(HKL),M20,X20,Aref,Vref,Z
    788788   
    789 def refinePeaks(peaks,ibrav,A,ifX20=True):
     789def refinePeaks(peaks,ibrav,A,ifX20=True,sg_type=None):
    790790    'needs a doc string'
    791791    dmin = getDmin(peaks)
     
    795795    OK = False
    796796    tries = 0
    797     HKL = G2lat.GenHBravais(dmin,ibrav,A)
     797    sg_type = G2lat.make_sgtype(ibrav)
     798    HKL = G2lat.GenHBravais(dmin,ibrav,A,sg_type)
    798799    while len(HKL) > 2 and IndexPeaks(peaks,HKL)[0]:
    799800        Pwr = pwr - (tries % 2)
     
    810811            continue
    811812        try:
    812             HKL = G2lat.GenHBravais(dmin,ibrav,A)
     813            HKL = G2lat.GenHBravais(dmin,ibrav,A,sg_type)
    813814        except FloatingPointError:
    814815            A = oldA
     
    919920    return A
    920921
    921 def DoIndexPeaks(peaks,controls,bravais,dlg,ifX20=True):
     922def DoIndexPeaks(peaks,controls,bravais,dlg,ifX20=True,
     923            timeout=None,M20_min=2.0,X20_max=None,return_Nc=False):
    922924    'needs a doc string'
    923925   
     
    963965                        N2 = 0
    964966                        while N2 < N2s[ibrav]:                                  #Table 2 step (iii)               
     967                            if timeout and time.time() - begin > timeout:
     968                                GoOn = False
     969                                break
     970
    965971                            if ibrav > 2:
    966972                                if not N2:
     
    993999                                    a,b,c,alp,bet,gam = G2lat.A2cell(A)
    9941000                                    V = G2lat.calc_V(A)
    995                                     if M20 >= 2.0:
     1001                                    if (
     1002                                        (M20 >= M20_min) and
     1003                                        (X20_max is None or X20 <= X20_max)
     1004                                    ):
    9961005                                        cell = [M20,X20,ibrav,a,b,c,alp,bet,gam,V,False,False]
     1006                                        if return_Nc: cell.append(Nc)
    9971007                                        newcell = np.array(cell[3:10])
    9981008                                        if not np.allclose(newcell,lastcell):
    999                                             print ("%10.3f %3d %3d %10.5f %10.5f %10.5f %10.3f %10.3f %10.3f %10.2f %10.2f"  \
    1000                                                 %(M20,X20,Nc,a,b,c,alp,bet,gam,V,V1))
     1009                                            print ("%10.3f %3d %3d %10.5f %10.5f %10.5f %10.3f %10.3f %10.3f %10.2f %10.2f %s"
     1010                                                %(M20,X20,Nc,a,b,c,alp,bet,gam,V,V1,bravaisNames[ibrav]))
    10011011                                            cells.append(cell)
    10021012                                        lastcell = np.array(cell[3:10])
     
    10181028                                    V1 = cells[0][9]
    10191029                                else:
    1020                                     V1 = 25
     1030                                    V1 = controls[3]
    10211031                                ncMax += Nobs
    10221032                                cycle += 1
Note: See TracChangeset for help on using the changeset viewer.