Ignore:
Timestamp:
Jul 23, 2015 3:36:38 PM (6 years ago)
Author:
vondreele
Message:

fix another new histo in new phase error
some cleanup in texture calculations - unneeded list generation

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIlattice.py

    r1820 r1945  
    11291129    import pytexture as ptx
    11301130    if SGLaue in ['m3','m3m']:
    1131         if phi.shape:
     1131        if 'float' in str(type(phi)):
    11321132            Kcl = np.zeros_like(phi)
    11331133        else:
    11341134            Kcl = 0.
    11351135        for j in range(0,L+1,4):
    1136             im = j/4+1
    1137             if phi.shape:
    1138                 pcrs = np.array([ptx.pyplmpsi(L,j,1,p)[0] for p in phi]).flatten()
     1136            im = j/4
     1137            if 'float' in str(type(phi)):
     1138                pcrs = ptx.pyplmpsi(L,j,len(phi),phi)[0]
    11391139            else:
    11401140                pcrs,dum = ptx.pyplmpsi(L,j,1,phi)
    1141             Kcl += BOH['L='+str(L)][N-1][im-1]*pcrs*cosd(j*beta)       
     1141            Kcl += BOH['L='+str(L)][N-1][im]*pcrs*cosd(j*beta)       
    11421142    else:
    1143         if phi.shape:
    1144             pcrs = np.array([ptx.pyplmpsi(L,N,1,p)[0] for p in phi]).flatten()
     1143        if 'float' in str(type(phi)):
     1144            pcrs = ptx.pyplmpsi(L,N,len(phi),phi)[0]
    11451145        else:
    11461146            pcrs,dum = ptx.pyplmpsi(L,N,1,phi)
     
    11661166        psrs,dpdps = ptx.pyplmpsi(L,M,1,psi)
    11671167    else:
    1168         psrs = np.array([ptx.pyplmpsi(L,M,1,p) for p in psi])
    1169         psrs,dpdps = np.reshape(psrs.flatten(),(-1,2)).T
     1168        psrs,dpdps = ptx.pyplmpsi(L,M,len(psi),psi)
    11701169    psrs *= RSQ2PI
    11711170    dpdps *= RSQ2PI
Note: See TracChangeset for help on using the changeset viewer.