Changeset 1601 for trunk/GSASII.py


Ignore:
Timestamp:
Dec 9, 2014 3:46:07 PM (7 years ago)
Author:
vondreele
Message:

remove some more (3+2) & (3+3) supersymmetry stuff
begin implementation of 4D Fourier maps & their display
implement super lattice for cif reflection import
fix space groups imported from cif

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASII.py

    r1594 r1601  
    655655            Super = generalData.get('Super',0)
    656656            SuperVec = []
    657             for i in range(Super):
    658                 SuperVec.append(generalData['SuperVec'][i][0])
    659                 SuperVec = np.array(SuperVec)
     657            if Super:
     658                SuperVec = np.array(generalData['SuperVec'][0])
    660659            UseList = data['Histograms']
    661660            NShkl = len(G2spc.MustrainNames(SGData))
     
    671670                    for iref,ref in enumerate(reflData['RefList']):
    672671                        hkl = ref[:3]
    673                         H = list(hkl+np.sum(SuperVec*ref[3:3+Super],axis=0))
     672                        if Super:
     673                            H = list(hkl+SuperVec*ref[3])
     674                        else:
     675                            H = hkl
    674676                        ref[4+Super] = np.sqrt(1./G2lat.calc_rDsq2(H,G))
    675677                        iabsnt,ref[3+Super],Uniq,phi = G2spc.GenHKLf(H,SGData)
     
    799801            Super = generalData.get('Super',0)
    800802            SuperVec = []
    801             for i in range(Super):
    802                 SuperVec.append(generalData['SuperVec'][i][0])
    803                 SuperVec = np.array(SuperVec)
     803            if Super:
     804                SuperVec = np.array(generalData['SuperVec'][0])
    804805            UseList = data['Histograms']
    805806            for histoName in newHistList:
     
    811812                for iref,ref in enumerate(reflData['RefList']):
    812813                    hkl = ref[:3]
    813                     H = list(hkl+np.sum(SuperVec*ref[3:3+Super],axis=0))
     814                    if Super:
     815                        H = list(hkl+SuperVec*ref[3])
     816                    else:
     817                        H = hkl
    814818                    ref[4+Super] = np.sqrt(1./G2lat.calc_rDsq2(H,G))
    815819                    iabsnt,ref[3+Super],Uniq,phi = G2spc.GenHKLf(H,SGData)
Note: See TracChangeset for help on using the changeset viewer.