def test_resp_from_paz_loading_vs_evalresp(self):
zeros = [0., 0.]
poles = [-8.443 + 1.443j, -8.443 - 1.443j]
filename = os.path.join(self.data_dir,
'RESP.XX.NS306..SHZ.GS13.1.2180')
resp_er = read_inventory(filename)[0][0][0].response
loaded_resp = resp_er.get_evalresp_response(.1, 2**6, output='VEL')
# The optional kwargs are the same as those being set in the RESP file.
resp = Response.from_paz(zeros=zeros,
poles=poles,
stage_gain=22.0,
stage_gain_frequency=5.0,
normalization_frequency=5.0,
normalization_factor=1.070401)
paz_resp = resp.get_evalresp_response(.1, 2**6, output='VEL')
np.testing.assert_allclose(paz_resp, loaded_resp)
def test_resp_from_paz_setting_zeros_poles_and_sensitivity(self):
poles = [1 + 2j, 1 - 2j, 2 + 1j, 2 - 1j]
zeros = [0, 0, 5]
sensitivity = 1201. * (2**26 / 40.)
# a0 normalization factor at 1Hz for these values
normalization = np.prod(2 * pi * 1j - np.array(zeros))
normalization /= np.prod(2 * pi * 1j - np.array(poles))
normalization = np.abs(normalization)
resp = Response.from_paz(zeros, poles, sensitivity)
r_zeros = resp.response_stages[0].zeros
r_poles = resp.response_stages[0].poles
r_stage_gain = resp.response_stages[0].stage_gain
r_sens = resp.instrument_sensitivity.value
np.testing.assert_array_equal(r_zeros, zeros)
np.testing.assert_array_equal(r_poles, poles)
np.testing.assert_equal(r_stage_gain, sensitivity)
np.testing.assert_equal(r_sens / normalization, sensitivity)