Changeset d522352 in sasmodels

Timestamp:

Mar 28, 2019 6:29:49 AM (6 years ago)

Author:

awashington

Branches:

master

Children:

0c7b8d8

Parents:

298d2d4

Message:

Divide by 2π, not multiply.

Also, the average of a and b is ½(a+b), not a+b. Hence (2 π)² is
about 39.2, the constant fact that I was off by.

File:

• sasmodels/sesans.py (modified) (2 diffs)

sasmodels/sesans.py

@@ -73 +73 @@
         q = np.hstack([[0], q])
 
-        H0 = np.pi * (q[1:]**2 - q[:-1]**2)
+        H0 = (q[1:]**2 - q[:-1]**2) / (2 * np.pi) / 2
 
         # repq = np.tile(q, (SElength.size, 1)).T
@@ -80 +80 @@
         H *= q.reshape((-1, 1))
         H = H[1:] - H[:-1]
-        H *= 2 * np.pi / SElength
+        H /= 2 * np.pi * SElength
 
         lam = np.asarray(lam, dtype=np.float32)