def test_dec_gardner_knopoff(self):
# Testing the Gardner and Knopoff algorithm
config = {'time_distance_window': GardnerKnopoffWindow(),
'fs_time_prop': 1.0}
# Instantiate the declusterer and process the sample catalogue
dec = GardnerKnopoffType1()
vcl, flagvector = dec.decluster(self.cat, config)
print('vcl:', vcl)
print('flagvector:', flagvector, self.cat.data['flag'])
np.testing.assert_allclose(flagvector, self.cat.data['flag'])
def test_dec_gardner_knopoff_time_cutoff(self):
"""
Testing the Gardner and Knopoff algorithm using a cutoff
time of 100 days
"""
config = {'time_distance_window' : GardnerKnopoffWindow(),
'fs_time_prop' : 1.0,
'time_cutoff' : 100}
# Instantiate the declusterer and process the sample catalogue
dec = GardnerKnopoffType1()
vcl, flagvector = dec.decluster(self.cat, config)
print('vcl:', vcl)
catalog_flag = self.cat.data['flag']
catalog_flag[4] = 0 # event becomes mainshock when time_cutoff = 100
print('flagvector:', catalog_flag)
np.testing.assert_allclose(flagvector,self.cat.data['flag'])