Skip to content
Snippets Groups Projects
Commit 462d033e authored by Andreas Weissenbrunner's avatar Andreas Weissenbrunner
Browse files

made the uncertainty calculation work again

parent 0917e3e5
No related branches found
No related tags found
No related merge requests found
...@@ -623,10 +623,9 @@ class mainPanel(wx.Panel): ...@@ -623,10 +623,9 @@ class mainPanel(wx.Panel):
def OnClicked_btn_uncertainty(self,e): def OnClicked_btn_uncertainty(self,e):
# calculate the uncertainty for the flow meter # calculate the uncertainty for the flow meter
if self.pathint.any(): if self.pathint.any():
self.fm_mean = np.zeros_like(self.pathint)
self.fm_std = np.zeros_like(self.pathint)
if flow.case == "SingleElbow": if flow.case == "SingleElbow":
self.fm_mean = np.zeros((len(flow.Rk),len(flow.dist),len(flow.phi)))
self.fm_std = np.zeros((len(flow.Rk),len(flow.dist),len(flow.phi)))
for i,rc in enumerate(flow.Rk): for i,rc in enumerate(flow.Rk):
for j,dist in enumerate(flow.dist): for j,dist in enumerate(flow.dist):
for k,phi in enumerate(flow.phi[:,0]): for k,phi in enumerate(flow.phi[:,0]):
...@@ -638,6 +637,8 @@ class mainPanel(wx.Panel): ...@@ -638,6 +637,8 @@ class mainPanel(wx.Panel):
self.regint_mean = interpol.RegularGridInterpolator((flow.Rk,flow.dist,flow.phi[:,0]),self.fm_mean,bounds_error = False,fill_value=None) self.regint_mean = interpol.RegularGridInterpolator((flow.Rk,flow.dist,flow.phi[:,0]),self.fm_mean,bounds_error = False,fill_value=None)
self.regint_std = interpol.RegularGridInterpolator((flow.Rk,flow.dist,flow.phi[:,0]),self.fm_std,bounds_error = False,fill_value=None) self.regint_std = interpol.RegularGridInterpolator((flow.Rk,flow.dist,flow.phi[:,0]),self.fm_std,bounds_error = False,fill_value=None)
else: else:
self.fm_mean = np.zeros((len(flow.Rk),len(flow.dl),len(flow.dist),len(flow.phi)))
self.fm_std = np.zeros((len(flow.Rk),len(flow.dl),len(flow.dist),len(flow.phi)))
for i,rc in enumerate(flow.Rk): for i,rc in enumerate(flow.Rk):
for j,dl in enumerate(flow.dl): for j,dl in enumerate(flow.dl):
for k,dist in enumerate(flow.dist): for k,dist in enumerate(flow.dist):
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment