Changes in src/sas/sascalc/calculator/resolution_calculator.py [9a5097c:463e7ffc] in sasview

File:

• src/sas/sascalc/calculator/resolution_calculator.py (modified) (12 diffs)

src/sas/sascalc/calculator/resolution_calculator.py

@@ -12 +12 @@
 from math import sqrt
 import math
-import numpy as np
+import numpy
 import sys
 import logging
+
+logger = logging.getLogger(__name__)
 
 #Plank's constant in cgs unit
@@ -393 +395 @@
         dx_size = (self.qx_max - self.qx_min) / (1000 - 1)
         dy_size = (self.qy_max - self.qy_min) / (1000 - 1)
-        x_val = np.arange(self.qx_min, self.qx_max, dx_size)
-        y_val = np.arange(self.qy_max, self.qy_min, -dy_size)
-        q_1, q_2 = np.meshgrid(x_val, y_val)
+        x_val = numpy.arange(self.qx_min, self.qx_max, dx_size)
+        y_val = numpy.arange(self.qy_max, self.qy_min, -dy_size)
+        q_1, q_2 = numpy.meshgrid(x_val, y_val)
         #q_phi = numpy.arctan(q_1,q_2)
         # check whether polar or cartesian
@@ -887 +889 @@
         x_value = x_val - x0_val
         y_value = y_val - y0_val
-        phi_i = np.arctan2(y_val, x_val)
+        phi_i = numpy.arctan2(y_val, x_val)
 
         # phi correction due to the gravity shift (in phi)
@@ -893 +895 @@
         phi_i = phi_i - phi_0 + self.gravity_phi
 
-        sin_phi = np.sin(self.gravity_phi)
-        cos_phi = np.cos(self.gravity_phi)
+        sin_phi = numpy.sin(self.gravity_phi)
+        cos_phi = numpy.cos(self.gravity_phi)
 
         x_p = x_value * cos_phi + y_value * sin_phi
@@ -908 +910 @@
         nu_value = -0.5 * (new_x * new_x + new_y * new_y)
 
-        gaussian = np.exp(nu_value)
+        gaussian = numpy.exp(nu_value)
         # normalizing factor correction
         gaussian /= gaussian.sum()
@@ -954 +956 @@
             nu_value *= nu_value
             nu_value *= -0.5
-            gaussian *= np.exp(nu_value)
+            gaussian *= numpy.exp(nu_value)
             gaussian /= sigma
             # normalize
@@ -1006 +1008 @@
             detector_offset = self.sample2detector_distance[1]
         except:
-            logging.error(sys.exc_value)
+            logger.error(sys.exc_value)
 
         # detector size in [no of pix_x,no of pix_y]
@@ -1026 +1028 @@
                                                            offset_x, offset_y)
         # distance [cm] from the beam center on detector plane
-        detector_ind_x = np.arange(detector_pix_nums_x)
-        detector_ind_y = np.arange(detector_pix_nums_y)
+        detector_ind_x = numpy.arange(detector_pix_nums_x)
+        detector_ind_y = numpy.arange(detector_pix_nums_y)
 
         # shif 0.5 pixel so that pix position is at the center of the pixel
@@ -1041 +1043 @@
         detector_ind_y = detector_ind_y * pix_y_size
 
-        qx_value = np.zeros(len(detector_ind_x))
-        qy_value = np.zeros(len(detector_ind_y))
+        qx_value = numpy.zeros(len(detector_ind_x))
+        qy_value = numpy.zeros(len(detector_ind_y))
         i = 0
 
@@ -1061 +1063 @@
 
         # p min and max values among the center of pixels
-        self.qx_min = np.min(qx_value)
-        self.qx_max = np.max(qx_value)
-        self.qy_min = np.min(qy_value)
-        self.qy_max = np.max(qy_value)
+        self.qx_min = numpy.min(qx_value)
+        self.qx_max = numpy.max(qx_value)
+        self.qy_min = numpy.min(qy_value)
+        self.qy_max = numpy.max(qy_value)
 
         # Appr. min and max values of the detector display limits
@@ -1088 +1090 @@
             from sas.sascalc.dataloader.data_info import Data2D
             output = Data2D()
-            inten = np.zeros_like(qx_value)
+            inten = numpy.zeros_like(qx_value)
             output.data = inten
             output.qx_data = qx_value
             output.qy_data = qy_value
         except:
-            logging.error(sys.exc_value)
+            logger.error(sys.exc_value)
 
         return output
@@ -1107 +1109 @@
         plane_dist = dx_size
         # full scattering angle on the x-axis
-        theta = np.arctan(plane_dist / det_dist)
-        qx_value = (2.0 * pi / wavelength) * np.sin(theta)
+        theta = numpy.arctan(plane_dist / det_dist)
+        qx_value = (2.0 * pi / wavelength) * numpy.sin(theta)
         return qx_value