def filter_ts9(input):
import numpy as np
output = np.zeros(input.shape, dtype=np.float64)
inputfilter = DoubleInPointerFilter(input, False)
inputfilter.set_input_sampling_rate(sample_rate)
infilter = DoubleChamberlinFilter()
infilter.set_input_sampling_rate(sample_rate)
infilter.set_input_port(0, inputfilter, 0)
infilter.select(0)
infilter.set_attenuation(1)
infilter.set_cut_frequency(1000)
overfilter = DoubleOversampling6points5order_4Filter()
overfilter.set_input_sampling_rate(sample_rate)
overfilter.set_output_sampling_rate(sample_rate * 4)
overfilter.set_input_port(0, infilter, 0)
overdrivefilter = DoubleTS9OverdriveFilter()
overdrivefilter.set_input_sampling_rate(sample_rate * 4)
overdrivefilter.set_input_port(0, overfilter, 0)
overdrivefilter.set_drive(0.9)
lowpassfilter = DoubleButterworthLowPassFilter()
lowpassfilter.set_input_sampling_rate(sample_rate * 4)
lowpassfilter.set_cut_frequency(sample_rate)
lowpassfilter.set_order(5)
lowpassfilter.set_input_port(0, overdrivefilter, 0)
decimationfilter = DoubleDecimationFilter(1)
decimationfilter.set_input_sampling_rate(sample_rate * 4)
decimationfilter.set_output_sampling_rate(sample_rate)
decimationfilter.set_input_port(0, lowpassfilter, 0)
outfilter = DoubleOutPointerFilter(output, False)
outfilter.set_input_sampling_rate(sample_rate)
outfilter.set_input_port(0, decimationfilter, 0)
outfilter.process(input.shape[1])
return output
def filter_4(input): import numpy as np output = np.zeros(input.shape, dtype=np.float64) inputfilter = DoubleInPointerFilter(input, False) inputfilter.input_sampling_rate = sample_rate overfilter = DoubleOversampling6points5order_4Filter() overfilter.input_sampling_rate = sample_rate overfilter.output_sampling_rate = sample_rate * 4 overfilter.set_input_port(0, inputfilter, 0) overdrivefilter = DoubleTS9OverdriveFilter() overdrivefilter.input_sampling_rate = sample_rate * 4 overdrivefilter.set_input_port(0, overfilter, 0) overdrivefilter.drive = 0.9 lowpassfilter = DoubleButterworthLowPassFilter() lowpassfilter.input_sampling_rate = sample_rate * 4 lowpassfilter.cut_frequency = sample_rate lowpassfilter.order = 5 lowpassfilter.set_input_port(0, overdrivefilter, 0) decimationfilter = DoubleDecimationFilter(1) decimationfilter.input_sampling_rate = sample_rate * 4 decimationfilter.output_sampling_rate = sample_rate decimationfilter.set_input_port(0, lowpassfilter, 0) outfilter = DoubleOutPointerFilter(output, False) outfilter.input_sampling_rate = sample_rate outfilter.set_input_port(0, decimationfilter, 0) outfilter.process(input.shape[1]) return output