Changeset 9a5097c in sasview for src/sas/sasgui/perspectives/calculator/gen_scatter_panel.py
- Timestamp:
- Mar 26, 2017 11:33:16 PM (8 years ago)
- Branches:
- master, ESS_GUI, ESS_GUI_Docs, ESS_GUI_batch_fitting, ESS_GUI_bumps_abstraction, ESS_GUI_iss1116, ESS_GUI_iss879, ESS_GUI_iss959, ESS_GUI_opencl, ESS_GUI_ordering, ESS_GUI_sync_sascalc, costrafo411, magnetic_scatt, release-4.2.2, ticket-1009, ticket-1094-headless, ticket-1242-2d-resolution, ticket-1243, ticket-1249, ticket885, unittest-saveload
- Children:
- ed2276f
- Parents:
- 9146ed9
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/sas/sasgui/perspectives/calculator/gen_scatter_panel.py
r0f7c930 r9a5097c 7 7 import sys 8 8 import os 9 import numpy 9 import numpy as np 10 10 #import math 11 11 import wx.aui as aui … … 741 741 marker = 'o' 742 742 m_size = 3.5 743 sld_tot = (n umpy.fabs(sld_mx) + numpy.fabs(sld_my) + \744 n umpy.fabs(sld_mz) + numpy.fabs(output.sld_n))743 sld_tot = (np.fabs(sld_mx) + np.fabs(sld_my) + \ 744 np.fabs(sld_mz) + np.fabs(output.sld_n)) 745 745 is_nonzero = sld_tot > 0.0 746 746 is_zero = sld_tot == 0.0 … … 757 757 pix_symbol = output.pix_symbol[is_nonzero] 758 758 # II. Plot selective points in color 759 other_color = n umpy.ones(len(pix_symbol), dtype='bool')759 other_color = np.ones(len(pix_symbol), dtype='bool') 760 760 for key in color_dic.keys(): 761 761 chosen_color = pix_symbol == key 762 if n umpy.any(chosen_color):762 if np.any(chosen_color): 763 763 other_color = other_color & (chosen_color != True) 764 764 color = color_dic[key] … … 767 767 markeredgecolor=color, markersize=m_size, label=key) 768 768 # III. Plot All others 769 if n umpy.any(other_color):769 if np.any(other_color): 770 770 a_name = '' 771 771 if output.pix_type == 'atom': … … 795 795 draw magnetic vectors w/arrow 796 796 """ 797 max_mx = max(n umpy.fabs(sld_mx))798 max_my = max(n umpy.fabs(sld_my))799 max_mz = max(n umpy.fabs(sld_mz))797 max_mx = max(np.fabs(sld_mx)) 798 max_my = max(np.fabs(sld_my)) 799 max_mz = max(np.fabs(sld_mz)) 800 800 max_m = max(max_mx, max_my, max_mz) 801 801 try: … … 812 812 unit_z2 = sld_mz / max_m 813 813 # 0.8 is for avoiding the color becomes white=(1,1,1)) 814 color_x = n umpy.fabs(unit_x2 * 0.8)815 color_y = n umpy.fabs(unit_y2 * 0.8)816 color_z = n umpy.fabs(unit_z2 * 0.8)814 color_x = np.fabs(unit_x2 * 0.8) 815 color_y = np.fabs(unit_y2 * 0.8) 816 color_z = np.fabs(unit_z2 * 0.8) 817 817 x2 = pos_x + unit_x2 * max_step 818 818 y2 = pos_y + unit_y2 * max_step 819 819 z2 = pos_z + unit_z2 * max_step 820 x_arrow = n umpy.column_stack((pos_x, x2))821 y_arrow = n umpy.column_stack((pos_y, y2))822 z_arrow = n umpy.column_stack((pos_z, z2))823 colors = n umpy.column_stack((color_x, color_y, color_z))820 x_arrow = np.column_stack((pos_x, x2)) 821 y_arrow = np.column_stack((pos_y, y2)) 822 z_arrow = np.column_stack((pos_z, z2)) 823 colors = np.column_stack((color_x, color_y, color_z)) 824 824 arrows = Arrow3D(panel, x_arrow, z_arrow, y_arrow, 825 825 colors, mutation_scale=10, lw=1, … … 880 880 if self.is_avg or self.is_avg == None: 881 881 self._create_default_1d_data() 882 i_out = n umpy.zeros(len(self.data.y))882 i_out = np.zeros(len(self.data.y)) 883 883 inputs = [self.data.x, [], i_out] 884 884 else: 885 885 self._create_default_2d_data() 886 i_out = n umpy.zeros(len(self.data.data))886 i_out = np.zeros(len(self.data.data)) 887 887 inputs = [self.data.qx_data, self.data.qy_data, i_out] 888 888 … … 989 989 :Param input: input list [qx_data, qy_data, i_out] 990 990 """ 991 out = n umpy.empty(0)991 out = np.empty(0) 992 992 #s = time.time() 993 993 for ind in range(len(input[0])): … … 998 998 inputi = [input[0][ind:ind + 1], [], input[2][ind:ind + 1]] 999 999 outi = self.model.run(inputi) 1000 out = n umpy.append(out, outi)1000 out = np.append(out, outi) 1001 1001 else: 1002 1002 if ind % 50 == 0 and update != None: … … 1006 1006 input[2][ind:ind + 1]] 1007 1007 outi = self.model.runXY(inputi) 1008 out = n umpy.append(out, outi)1008 out = np.append(out, outi) 1009 1009 #print time.time() - s 1010 1010 if self.is_avg or self.is_avg == None: … … 1027 1027 self.npts_x = int(float(self.npt_ctl.GetValue())) 1028 1028 self.data = Data2D() 1029 qmax = self.qmax_x #/ n umpy.sqrt(2)1029 qmax = self.qmax_x #/ np.sqrt(2) 1030 1030 self.data.xaxis('\\rm{Q_{x}}', '\AA^{-1}') 1031 1031 self.data.yaxis('\\rm{Q_{y}}', '\AA^{-1}') … … 1048 1048 qstep = self.npts_x 1049 1049 1050 x = n umpy.linspace(start=xmin, stop=xmax, num=qstep, endpoint=True)1051 y = n umpy.linspace(start=ymin, stop=ymax, num=qstep, endpoint=True)1050 x = np.linspace(start=xmin, stop=xmax, num=qstep, endpoint=True) 1051 y = np.linspace(start=ymin, stop=ymax, num=qstep, endpoint=True) 1052 1052 ## use data info instead 1053 new_x = n umpy.tile(x, (len(y), 1))1054 new_y = n umpy.tile(y, (len(x), 1))1053 new_x = np.tile(x, (len(y), 1)) 1054 new_y = np.tile(y, (len(x), 1)) 1055 1055 new_y = new_y.swapaxes(0, 1) 1056 1056 # all data reuire now in 1d array 1057 1057 qx_data = new_x.flatten() 1058 1058 qy_data = new_y.flatten() 1059 q_data = n umpy.sqrt(qx_data * qx_data + qy_data * qy_data)1059 q_data = np.sqrt(qx_data * qx_data + qy_data * qy_data) 1060 1060 # set all True (standing for unmasked) as default 1061 mask = n umpy.ones(len(qx_data), dtype=bool)1061 mask = np.ones(len(qx_data), dtype=bool) 1062 1062 # store x and y bin centers in q space 1063 1063 x_bins = x 1064 1064 y_bins = y 1065 1065 self.data.source = Source() 1066 self.data.data = n umpy.ones(len(mask))1067 self.data.err_data = n umpy.ones(len(mask))1066 self.data.data = np.ones(len(mask)) 1067 self.data.err_data = np.ones(len(mask)) 1068 1068 self.data.qx_data = qx_data 1069 1069 self.data.qy_data = qy_data … … 1084 1084 :warning: This data is never plotted. 1085 1085 residuals.x = data_copy.x[index] 1086 residuals.dy = n umpy.ones(len(residuals.y))1086 residuals.dy = np.ones(len(residuals.y)) 1087 1087 residuals.dx = None 1088 1088 residuals.dxl = None … … 1091 1091 self.qmax_x = float(self.qmax_ctl.GetValue()) 1092 1092 self.npts_x = int(float(self.npt_ctl.GetValue())) 1093 qmax = self.qmax_x #/ n umpy.sqrt(2)1093 qmax = self.qmax_x #/ np.sqrt(2) 1094 1094 ## Default values 1095 1095 xmax = qmax 1096 1096 xmin = qmax * _Q1D_MIN 1097 1097 qstep = self.npts_x 1098 x = n umpy.linspace(start=xmin, stop=xmax, num=qstep, endpoint=True)1098 x = np.linspace(start=xmin, stop=xmax, num=qstep, endpoint=True) 1099 1099 # store x and y bin centers in q space 1100 1100 #self.data.source = Source() 1101 y = n umpy.ones(len(x))1102 dy = n umpy.zeros(len(x))1103 dx = n umpy.zeros(len(x))1101 y = np.ones(len(x)) 1102 dy = np.zeros(len(x)) 1103 dx = np.zeros(len(x)) 1104 1104 self.data = Data1D(x=x, y=y) 1105 1105 self.data.dx = dx … … 1171 1171 state = None 1172 1172 1173 n umpy.nan_to_num(image)1173 np.nan_to_num(image) 1174 1174 new_plot = Data2D(image=image, err_image=data.err_data) 1175 1175 new_plot.name = model.name + '2d' … … 1640 1640 for key in sld_list.keys(): 1641 1641 if ctr_list[0] == key: 1642 min_val = n umpy.min(sld_list[key])1643 max_val = n umpy.max(sld_list[key])1644 mean_val = n umpy.mean(sld_list[key])1642 min_val = np.min(sld_list[key]) 1643 max_val = np.max(sld_list[key]) 1644 mean_val = np.mean(sld_list[key]) 1645 1645 enable = (min_val == max_val) and \ 1646 1646 sld_data.pix_type == 'pixel' … … 1733 1733 npts = -1 1734 1734 break 1735 if n umpy.isfinite(n_val):1735 if np.isfinite(n_val): 1736 1736 npts *= int(n_val) 1737 1737 if npts > 0: … … 1770 1770 ctl.Refresh() 1771 1771 return 1772 if n umpy.isfinite(s_val):1772 if np.isfinite(s_val): 1773 1773 s_size *= s_val 1774 1774 self.sld_data.set_pixel_volumes(s_size) … … 1787 1787 try: 1788 1788 sld_data = self.parent.get_sld_from_omf() 1789 #nop = (nop * n umpy.pi) / 61789 #nop = (nop * np.pi) / 6 1790 1790 nop = len(sld_data.sld_n) 1791 1791 except:
Note: See TracChangeset
for help on using the changeset viewer.