def compute_pulsation_at_node(self):
i, j, k = self.get_puls_node()
compute_sem_node_hist(self._c_init_data, self.c_stat_data,
self.c_trans_data, self._ts, self._num_ts_tot,
self.c_node_hist, i, j, k)
self._vel_puls = extract_pulsations_from_node_hist(self.c_node_hist)
free_node_hist(self.c_node_hist)
def compute_pulsation_at_node(self):
i, j, k = self.get_puls_node()
compute_smirnov_node_hist(init_data=self._c_init_data,
data=self.c_data,
ts=self._ts,
num_ts=self._num_ts_tot,
node_hist=self.c_node_hist,
i=i,
j=j,
k=k)
self._vel_puls = extract_pulsations_from_node_hist(self.c_node_hist)
free_node_hist(self.c_node_hist)
def test_compute_Smirnov_field_node(self):
data = smirnov.STG_SmirnovData()
node_hist = smirnov.STG_VelNodeHist()
smirnov.compute_smirnov_data(self.init_data, self.num_modes, data)
num_ts = len(self.time_arr) - 1
ts = self.time_arr[1] - self.time_arr[0]
smirnov.compute_smirnov_node_hist(self.init_data, data, ts, num_ts, node_hist, 0, 0, 0)
u, v, w = common.extract_pulsations_from_node_hist(node_hist)
self.assertNotEqual(u[0], u[1])
self.assertNotEqual(u[0], v[1])
self.assertNotEqual(w[0], v[1])
self.assertNotEqual(u[0], 0)
self.assertIsNotNone(u[0])
smirnov.free_smirnov_data(data)
smirnov.free_node_hist(node_hist)
def test_compute_davidson_node_hist(self):
ts = self.time_arr[1] - self.time_arr[0]
stat_data = davidson.STG_DavidsonData_Stationary()
davidson.compute_davidson_stat_data(self.init_data, self.num_modes, self.dissip_rate,
self.visc, ts, stat_data)
trans_data = davidson.STG_DavidsonData_Transient()
davidson.alloc_davidson_trans_data(self.init_data, self.num_modes, len(self.time_arr) - 1, trans_data)
davidson.compute_davidson_trans_data(stat_data, len(self.time_arr) - 1, trans_data)
node_hist = common.STG_VelNodeHist()
davidson.compute_davidson_node_hist(
self.init_data, stat_data, ts, len(self.time_arr) - 1, trans_data, node_hist, 0, 0, 0)
u, v, w = common.extract_pulsations_from_node_hist(node_hist)
self.assertNotEqual(u[0], u[1])
self.assertNotEqual(u[0], v[1])
self.assertNotEqual(w[0], v[1])
self.assertNotEqual(u[0], 0)
self.assertIsNotNone(u[0])
common.free_node_hist(node_hist)
davidson.free_davidson_stat_data(stat_data)
davidson.free_davidson_trans_data(trans_data)