def test_determine_peaks_only_series_with_nonchanging_values():
values = np.array([0, 1, 1, -3, -5, 0]) # constant then reverse
cum_abs_delta_values = np.sum(np.abs(np.diff(values)))
expected_sum = cum_abs_delta_values / 2
peaks_only = eqsig.determine_pseudo_cyclic_peak_only_series(values)
cum_peaks = np.sum(peaks_only)
assert np.isclose(cum_peaks, expected_sum), cum_peaks
values = np.array([0, 1, 1, 3, -5, 0]) # constant the no reverse
cum_abs_delta_values = np.sum(np.abs(np.diff(values)))
expected_sum = cum_abs_delta_values / 2
peaks_only = eqsig.determine_pseudo_cyclic_peak_only_series(values)
cum_peaks = np.sum(peaks_only)
assert np.isclose(cum_peaks, expected_sum), cum_peaks
def test_determine_pseudo_cyclic_peak_only_series_with_triangle_series():
values = [0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0]
cum_abs_delta_values = np.sum(np.abs(np.diff(values)))
expected_sum = cum_abs_delta_values / 2
peaks_only = eqsig.determine_pseudo_cyclic_peak_only_series(values)
cum_peaks = np.sum(np.abs(peaks_only))
assert np.isclose(cum_peaks, expected_sum)
def test_determine_pseudo_cyclic_peak_only_series_with_non_zero_end():
end_value = 1.
values = np.array([0, 2, -1, 2, 0, end_value])
cum_abs_delta_values = np.sum(np.abs(np.diff(values)))
expected_sum = cum_abs_delta_values / 2 + end_value / 2
peaks_only = eqsig.determine_pseudo_cyclic_peak_only_series(values)
cum_peaks = np.sum(peaks_only)
assert np.isclose(cum_peaks, expected_sum)
def test_determine_pseudo_cyclic_peak_only_series_with_ground_motion():
record_path = TEST_DATA_DIR
record_filename = 'test_motion_dt0p01.txt'
rec = np.loadtxt(record_path + record_filename, skiprows=2)
cum_abs_delta_values = np.sum(np.abs(np.diff(rec)))
expected_sum = cum_abs_delta_values / 2
peaks_only = eqsig.determine_pseudo_cyclic_peak_only_series(rec)
cum_peaks = np.sum(peaks_only)
assert np.isclose(cum_peaks, expected_sum)
def test_determine_pseudo_cyclic_peak_only_series_with_sine_wave():
time = np.arange(99)
values = np.sin(time)
values[-1] = 0
cum_abs_delta_values = np.sum(np.abs(np.diff(values)))
expected_sum = cum_abs_delta_values / 2
peaks_only = eqsig.determine_pseudo_cyclic_peak_only_series(values)
cum_peaks = np.sum(np.abs(peaks_only))
assert np.isclose(cum_peaks, expected_sum)
def test_determine_pseudo_cyclic_peak_only_series_with_a_double_peak_and_offset(
):
values = np.array([0, 2, 1, 2, 0, 1, 0, -1, 0, 1, 0]) + 4
cum_abs_delta_values = np.sum(np.abs(np.diff(values)))
expected_sum = cum_abs_delta_values / 2
peaks_only = eqsig.determine_pseudo_cyclic_peak_only_series(values)
cum_peaks = np.sum(peaks_only)
expected_series = np.array([0, 2, -1, 2, 0, 1, 0, 1, 0, 1, 0])
assert np.sum(np.abs(peaks_only - expected_series)) == 0.0
assert np.isclose(cum_peaks, expected_sum)
