Changeset 0c7b8d8 in sasmodels

Timestamp:

Mar 28, 2019 6:33:05 AM (6 years ago)

Author:

awashington

Branches:

master

Children:

3f3df6c

Parents:

d522352

Message:

Use 64 bit floats in SesansTransform?

The new log spacing creates small matrices, so there's no longer a
reason not to use full double precision values.

File:

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

sasmodels/sesans.py

@@ -59 +59 @@
     def _set_hankel(self, SElength, lam, zaccept, Rmax):
         # type: (np.ndarray, float, float) -> None
-        # Force float32 arrays, otherwise run into memory problems on
-        # some machines
-        SElength = np.asarray(SElength, dtype='float32')
+        SElength = np.asarray(SElength, dtype='float64')
 
         # Rmax = #value in text box somewhere in FitPage?
@@ -70 +68 @@
         #                      dtype=np.float32))
         q = np.exp(np.linspace(np.log(q_min), np.log(q_max), 10*SElength.size,
-                               dtype=np.float32))
+                               dtype=np.float64))
         q = np.hstack([[0], q])
 
@@ -82 +80 @@
         H /= 2 * np.pi * SElength
 
-        lam = np.asarray(lam, dtype=np.float32)
+        lam = np.asarray(lam, dtype=np.float64)
         reptheta = np.outer(q[1:], lam)
-        reptheta /= np.float32(2*np.pi)
+        reptheta /= np.float64(2*np.pi)
         np.arcsin(reptheta, out=reptheta)
         # reptheta = np.arcsin(repq*replam/2*np.pi)