- Timestamp:
- Jul 19, 2020 10:04:30 AM (3 years ago)
- Location:
- trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/GSASIImath.py
r4519 r4526 1503 1503 MmodA = TA[nxs,:,:,nxs]*(Am[nxs,nxs,:,:]*pcos[:,:,:,nxs]-Bm[nxs,nxs,:,:]*psin[:,:,:,nxs])/2. #Nref,Nops,Natm,Mxyz 1504 1504 MmodB = TA[nxs,:,:,nxs]*(Am[nxs,nxs,:,:]*psin[:,:,:,nxs]+Bm[nxs,nxs,:,:]*pcos[:,:,:,nxs])/2. #Nref,Nops,Natm,Mxyz 1505 MmodA = np.sum(SGMT[nxs,:,nxs,:,:]*MmodA[:,:,:,nxs,:],axis=-1)*SGData[' MagMom'][nxs,:,nxs,nxs]1506 MmodB = np.sum(SGMT[nxs,:,nxs,:,:]*MmodB[:,:,:,nxs,:],axis=-1)*SGData[' MagMom'][nxs,:,nxs,nxs]1505 MmodA = np.sum(SGMT[nxs,:,nxs,:,:]*MmodA[:,:,:,nxs,:],axis=-1)*SGData['SpnFlp'][nxs,:,nxs,nxs] 1506 MmodB = np.sum(SGMT[nxs,:,nxs,:,:]*MmodB[:,:,:,nxs,:],axis=-1)*SGData['SpnFlp'][nxs,:,nxs,nxs] 1507 1507 return MmodA,MmodB #Nref,Nops,Natm,Mxyz; cos & sin parts 1508 1508 -
trunk/GSASIIspc.py
r4510 r4526 3569 3569 def ApplyStringOpsMom(A,SGData,SSGData,Mom): 3570 3570 '''Applies string operations to modulated magnetic moment components used in drawing 3571 Drawing matches Bilbao MVISUALIZE 3571 3572 ''' 3572 3573 SGOps = SGData['SGOps'] … … 3583 3584 NA += len(SGOps) 3584 3585 M,T = SGOps[nA] 3585 newMom = np.inner(Mom,M).T*SGData[' MagMom'][NA+nC]3586 newMom = np.inner(Mom,M).T*SGData['SpnFlp'][NA+nC] 3586 3587 if SSGData is not None: 3587 3588 if SSGData['SSGCen'][iAx//100][3]: #flip spin for BNS centered atoms
Note: See TracChangeset
for help on using the changeset viewer.