Changeset 4138 for trunk/GSASIIstrMath.py
- Timestamp:
- Sep 5, 2019 4:08:51 PM (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/GSASIIstrMath.py
r4137 r4138 1510 1510 1511 1511 mXYZ = np.array([[xyz[0] for xyz in list(G2spc.GenAtom(xyz,SGData,All=True,Move=True))] for xyz in (Xdata+dXdata).T])%1. #Natn,Nop,xyz 1512 # TmagA,TmagB = G2mth.MagMod2(mXYZ,modQ,MSSdata,SGData,SSGData) #Nops,Natm,Mxyz-TmagA+TmagB matches drawing moments @ tau=0.1513 1512 MmodA,MmodB = G2mth.MagMod(glTau,mXYZ,modQ,MSSdata,SGData,SSGData) #Ntau,Nops,Natm,Mxyz sum matches drawing 1514 1513 1515 1514 if not SGData['SGGray']: #for fixed Mx,My,Mz 1516 # Mmod += Gdata.T[:,nxs,:]1517 1515 GSdata = np.inner(Gdata.T,np.swapaxes(SGMT,1,2)) #apply sym. ops.--> Natm,Nops,Nxyz 1518 1516 if SGData['SGInv'] and not SGData['SGFixed']: #inversion if any … … 1613 1611 1614 1612 if not SGData['SGGray']: 1613 fams *= 0.5 1614 fbms *= 0.5 1615 1615 fams += fam0[:,nxs,:,:,:] 1616 1616 fbms += fbm0[:,nxs,:,:,:]
Note: See TracChangeset
for help on using the changeset viewer.