Changeset 4179

Timestamp:

Oct 16, 2019 9:19:14 AM (4 years ago)

Author:

vondreele

Message:

a bit more speedup for polarization correction

File:

• trunk/GSASIIpwd.py (modified) (1 diff)

trunk/GSASIIpwd.py

@@ -234 +234 @@
 
     """
-    pola = ((1.0-Pola)*npcosd(Azm)**2+Pola*npsind(Azm)**2)*npcosd(Tth)**2+   \
-        (1.0-Pola)*npsind(Azm)**2+Pola*npcosd(Azm)**2
-    dpdPola = -npsind(Tth)**2*(npsind(Azm)**2-npcosd(Azm)**2)
+    cazm = npcosd(Azm)
+    sazm = npsind(Azm)
+    pola = ((1.0-Pola)*cazm**2+Pola*npsind(Azm)**2)*npcosd(Tth)**2+   \
+        (1.0-Pola)*sazm**2+Pola*cazm**2
+    dpdPola = -npsind(Tth)**2*(sazm**2-cazm**2)
     return pola,dpdPola