Changeset 3340 for trunk/GSASIIstrIO.py


Ignore:
Timestamp:
Apr 11, 2018 1:32:49 PM (4 years ago)
Author:
vondreele
Message:

make testDeriv py3 compatible
fix problem with Reload draw atoms - mag moments
add |Mag| to lst output of magnetic moments
correct errors in mag moment structure factor & derivatives

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIstrIO.py

    r3311 r3340  
    16561656           
    16571657    def PrintMomentsAndSig(General,Atoms,atomsSig):
     1658        cell = General['Cell'][1:7]
     1659        G = G2lat.fillgmat(cell)
     1660        ast = np.sqrt(np.diag(G))
     1661        GS = G/np.outer(ast,ast)
    16581662        pFile.write('\n Magnetic Moments:\n')    #add magnitude & angle, etc.? TBD
    1659         line = '   name      Mx        My        Mz'
     1663        line = '   name       Mx        My        Mz       |Mag|'
    16601664        cx,ct,cs,cia = General['AtomPtrs']
    16611665        cmx = cx+4
     
    16631667        pFile.write(line+'\n')
    16641668        pFile.write(135*'-'+'\n')
    1665         fmt = {0:'%7s',ct:'%7s',cmx:'%10.5f',cmx+1:'%10.5f',cmx+2:'%10.5f'}
     1669        fmt = {0:'%7s',ct:'%7s',cmx:'%10.3f',cmx+1:'%10.3f',cmx+2:'%10.3f'}
    16661670        noFXsig = {cmx:[10*' ','%10s'],cmx+1:[10*' ','%10s'],cmx+2:[10*' ','%10s']}
    16671671        for i,at in enumerate(Atoms):
     
    16771681                    else:
    16781682                        sigstr += noFXsig[ind][1]%(noFXsig[ind][0])
     1683                mag = np.array(at[cmx:cmx+3])
     1684                Mag = np.sqrt(np.inner(mag,np.inner(mag,GS)))
     1685                valstr += '%10.3f'%Mag
     1686                sigstr += 10*' '
    16791687                pFile.write(name+'\n')
    16801688                pFile.write(valstr+'\n')
Note: See TracChangeset for help on using the changeset viewer.