Ignore:
Timestamp:
Jun 29, 2015 9:06:28 AM (7 years ago)
Author:
vondreele
Message:

more twin stuff - nonmerohedral mods
refl[3] now twin id not mul.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrMath.py

    r1910 r1912  
    853853        refl = refDict['RefList'][iBeg:iFin]    #array(blkSize,nItems)
    854854        H = refl.T[:3]                          #array(blkSize,3)
    855         H = np.squeeze(np.inner(H.T,TwinLaw))   #maybe array(blkSize,3,nTwins) or (blkSize,3)
     855        H = np.squeeze(np.inner(H.T,TwinLaw))   #maybe array(blkSize,nTwins,3) or (blkSize,3)
     856        TwMask = np.any(H,axis=2)
     857        if TwinLaw.shape[0] > 1 and TwinDict:
     858            for ir in range(blkSize):
     859                iref = ir+iBeg
     860                if iref in TwinDict:
     861                    for i in TwinDict[iref]:
     862                        H[ir][i] = TwinDict[iref][i]
     863            TwMask = np.any(H,axis=2)
    856864        SQ = 1./(2.*refl.T[4])**2               #array(blkSize)
    857865        SQfactor = 4.0*SQ*twopisq               #ditto prev.
Note: See TracChangeset for help on using the changeset viewer.