Ignore:
Timestamp:
Jul 24, 2015 2:55:09 PM (6 years ago)
Author:
vondreele
Message:

do export of pole figures & inverse pole figures as csv files (can be read by Excel).
fix some errors in texture calcs. - mostly numpy/nonnumpy issues

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIlattice.py

    r1945 r1947  
    11291129    import pytexture as ptx
    11301130    if SGLaue in ['m3','m3m']:
    1131         if 'float' in str(type(phi)):
     1131        if 'array' in str(type(phi)) and np.any(phi.shape):
    11321132            Kcl = np.zeros_like(phi)
    11331133        else:
     
    11351135        for j in range(0,L+1,4):
    11361136            im = j/4
    1137             if 'float' in str(type(phi)):
     1137            if 'array' in str(type(phi)) and np.any(phi.shape):
    11381138                pcrs = ptx.pyplmpsi(L,j,len(phi),phi)[0]
    11391139            else:
    1140                 pcrs,dum = ptx.pyplmpsi(L,j,1,phi)
    1141             Kcl += BOH['L='+str(L)][N-1][im]*pcrs*cosd(j*beta)       
     1140                pcrs = ptx.pyplmpsi(L,j,1,phi)[0]
     1141            Kcl += BOH['L=%d'%(L)][N-1][im]*pcrs*cosd(j*beta)       
    11421142    else:
    1143         if 'float' in str(type(phi)):
     1143        if 'array' in str(type(phi)) and np.any(phi.shape):
    11441144            pcrs = ptx.pyplmpsi(L,N,len(phi),phi)[0]
    11451145        else:
    1146             pcrs,dum = ptx.pyplmpsi(L,N,1,phi)
     1146            pcrs = ptx.pyplmpsi(L,N,1,phi)[0]
    11471147        pcrs *= RSQ2PI
    11481148        if N:
     
    11631163    'needs doc string'
    11641164    import pytexture as ptx
    1165     if 'float' in str(type(psi)):
     1165    if 'array' in str(type(psi)) and np.any(psi.shape):
     1166        psrs,dpdps = ptx.pyplmpsi(L,M,len(psi),psi)
     1167    else:
    11661168        psrs,dpdps = ptx.pyplmpsi(L,M,1,psi)
    1167     else:
    1168         psrs,dpdps = ptx.pyplmpsi(L,M,len(psi),psi)
    11691169    psrs *= RSQ2PI
    11701170    dpdps *= RSQ2PI
     
    11951195        Kcl = 0.0
    11961196        for j in range(0,L+1,4):
    1197             im = j/4+1
     1197            im = j/4
    11981198            pcrs,dum = ptx.pyplmpsi(L,j,1,phi)
    1199             Kcl += BOH['L='+str(L)][N-1][im-1]*pcrs*cosd(j*beta)       
     1199            Kcl += BOH['L=%d'%(L)][N-1][im]*pcrs*cosd(j*beta)       
    12001200    else:
    12011201        pcrs,dum = ptx.pyplmpsi(L,N,1,phi)
     
    12501250            Kcl = 0.0
    12511251            for j in range(0,l+1,4):
    1252                 im = j/4+1
     1252                im = j/4
    12531253                pcrs,dum = ptx.pyplmpsi(l,j,1,phi)
    1254                 Kcl += BOH['L='+str(l)][n-1][im-1]*pcrs*cosd(j*beta)       
     1254                Kcl += BOH['L='+str(l)][n-1][im]*pcrs*cosd(j*beta)       
    12551255        else:                #all but cubic
    12561256            pcrs,dum = ptx.pyplmpsi(l,n,1,phi)
     
    13071307                Kcl = 0.0
    13081308                for j in range(0,l+1,4):
    1309                     im = j/4+1
     1309                    im = j/4
    13101310                    pcrs,dum = ptx.pyplmpsi(l,j,len(beta),phi)
    1311                     Kcl += BOH['L='+str(l)][n-1][im-1]*pcrs*cosd(j*beta)       
     1311                    Kcl += BOH['L=%d'%(l)][n-1][im]*pcrs*cosd(j*beta)       
    13121312            else:                #all but cubic
    13131313                pcrs,dum = ptx.pyplmpsi(l,n,len(beta),phi)
Note: See TracChangeset for help on using the changeset viewer.