def test_rt_migration_true(self):
migrator = RtMigrator(self.wo)
nsta = migrator.nsta
ntr=len(self.obs_split[0])
#########################
# start loops
#########################
# loop over segments (simulate real-time data)
for itr in xrange(ntr):
# update all input streams
# loop over stations
data_list=[]
for ista in xrange(nsta):
tr = self.obs_split[ista][itr]
data_list.append(tr)
# update data
migrator.updateData(data_list)
# update stacks
migrator.updateStacks()
# update max
migrator.updateMax()
#########################
# end loops
#########################
# check we find the same absolute origin time
#migrator.max_out.plot()
st=migrator.max_out.stats.starttime+45
ed=migrator.max_out.stats.starttime+55
#max_out=migrator.max_out.slice(st,ed)
#x_out=migrator.x_out.slice(st,ed)
#y_out=migrator.y_out.slice(st,ed)
#z_out=migrator.z_out.slice(st,ed)
#plotMaxXYZ(max_out, x_out, y_out, z_out, 'test_out.png')
max_trace=migrator.max_out.data
tmax=np.argmax(max_trace)*self.dt
tdiff=(migrator.max_out.stats.starttime + tmax)-(self.starttime + self.ot)
self.assertEquals(tdiff,0)
def test_rt_migration_true(self):
migrator = RtMigrator(self.wo)
nsta = migrator.nsta
ntr = len(self.obs_split[0])
#########################
# start loops
#########################
# loop over segments (simulate real-time data)
for itr in xrange(ntr):
# update all input streams
# loop over stations
data_list = []
for ista in xrange(nsta):
tr = self.obs_split[ista][itr]
data_list.append(tr)
# update data
migrator.updateData(data_list)
# update stacks
migrator.updateStacks()
# update max
migrator.updateMax()
#########################
# end loops
#########################
# check we find the same absolute origin time
#migrator.max_out.plot()
st = migrator.max_out.stats.starttime + 45
ed = migrator.max_out.stats.starttime + 55
#max_out=migrator.max_out.slice(st,ed)
#x_out=migrator.x_out.slice(st,ed)
#y_out=migrator.y_out.slice(st,ed)
#z_out=migrator.z_out.slice(st,ed)
#plotMaxXYZ(max_out, x_out, y_out, z_out, 'test_out.png')
max_trace = migrator.max_out.data
tmax = np.argmax(max_trace) * self.dt
tdiff = (migrator.max_out.stats.starttime + tmax) - (self.starttime +
self.ot)
self.assertEquals(tdiff, 0)
for itr in xrange(ntr):
# update all input streams
# loop over stations
data_list = []
for ista in xrange(nsta):
tr = obs_split[ista][itr]
data_list.append(tr)
# update data
t0 = time.time()
migrator.updateData(data_list)
t_update_data.append(time.time() - t0)
# update stacks
t0 = time.time()
migrator.updateStacks()
t_update_stacks.append(time.time() - t0)
# update max
t0 = time.time()
migrator.updateMax()
t_update_max.append(time.time() - t0)
#########################
# end loops
#########################
tac = time.time()
print "Time taken for four-minute 100Hz real-data test for %d stations on %d test-points : %.2f s" % (
nsta, n_test, tac - tic)
for itr in xrange(ntr):
# update all input streams
# loop over stations
data_list=[]
for ista in xrange(nsta):
tr = obs_split[ista][itr]
data_list.append(tr)
# update data
t0=time.time()
migrator.updateData(data_list)
t_update_data.append(time.time()-t0)
# update stacks
t0=time.time()
migrator.updateStacks()
t_update_stacks.append(time.time()-t0)
# update max
t0=time.time()
migrator.updateMax()
t_update_max.append(time.time()-t0)
#########################
# end loops
#########################
tac=time.time()
print "Time taken for four-minute 100Hz real-data test for %d stations on %d test-points : %.2f s"%(nsta, n_test, tac-tic)
print "Times taken for updating data : %.2f s" % np.sum(t_update_data)
