Changeset 1452

Timestamp:

Aug 1, 2014 2:14:23 PM (9 years ago)

Author:

toby

Message:

fix bug that shows up in numpy 1.8

File:

• trunk/GSASIIlattice.py (modified) (2 diffs)

trunk/GSASIIlattice.py

@@ -1203 +1203 @@
     
 def polfcal(ODFln,SamSym,psi,gam):
-    'needs doc string'
+    '''Perform a pole figure computation.
+    Note that the the number of gam values must either be 1 or must
+    match psi. Updated for numpy 1.8.0
+    '''
     import pytexture as ptx
     RSQPI = 0.5641895835478
     SQ2 = 1.414213562373
-    PolVal = np.ones_like(gam)
+    PolVal = np.ones_like(psi)
     for term in ODFln:
         if abs(ODFln[term][1]) > 1.e-3:
@@ -1222 +1225 @@
                 else:
                     Ksl = RSQPI*psrs/SQ2
-            PolVal += ODFln[term][1]*Ksl 
+            PolVal += ODFln[term][1]*Ksl
     return PolVal