Changeset 4571 for trunk/GSASIIpwd.py

Timestamp:

Sep 17, 2020 3:50:24 PM (3 years ago)

Author:

vondreele

Message:

Improvements to polarization stuff

File:

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

trunk/GSASIIpwd.py

@@ -237 +237 @@
 
     """
-    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)
+    cazm = npcosd(Azm)**2
+    sazm = npsind(Azm)**2
+    pola = ((1.0-Pola)*cazm+Pola*sazm)*npcosd(Tth)**2+(1.0-Pola)*sazm+Pola*cazm
+    dpdPola = -npsind(Tth)**2*(sazm-cazm)
     return pola,dpdPola