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)