Skip to content

Commit

Permalink
explicit type casting
Browse files Browse the repository at this point in the history
NPY_PROMOTION_STATE=weak_and_warn reports that several variables changed from int16 to int64. Casting to int64 addresses the issue.
  • Loading branch information
tompollard committed Oct 9, 2024
1 parent 56b7d5c commit d86a79b
Showing 1 changed file with 6 additions and 3 deletions.
9 changes: 6 additions & 3 deletions wfdb/io/convert/edf.py
Original file line number Diff line number Diff line change
Expand Up @@ -402,23 +402,26 @@ def read_edf(
temp_sig_data = np.fromfile(edf_file, dtype=np.int16)
temp_sig_data = temp_sig_data.reshape((-1, sum(samps_per_block)))
temp_all_sigs = np.hsplit(temp_sig_data, np.cumsum(samps_per_block)[:-1])

for i in range(n_sig):
# Check if `samps_per_frame` has all equal values
if samps_per_frame.count(samps_per_frame[0]) == len(samps_per_frame):
sig_data[:, i] = (
temp_all_sigs[i].flatten() - baseline[i]
(temp_all_sigs[i].flatten() - baseline[i]).astype(np.int64)
) / adc_gain_all[i]
else:
temp_sig_data = temp_all_sigs[i].flatten()

if samps_per_frame[i] == 1:
sig_data[:, i] = (temp_sig_data - baseline[i]) / adc_gain_all[i]
sig_data[:, i] = (temp_sig_data - baseline[i]).astype(
np.int64) / adc_gain_all[i]
else:
for j in range(sig_len):
start_ind = j * samps_per_frame[i]
stop_ind = start_ind + samps_per_frame[i]
sig_data[j, i] = np.mean(
(temp_sig_data[start_ind:stop_ind] - baseline[i])
/ adc_gain_all[i]
.astype(np.int64) / adc_gain_all[i]
)

# This is the closest I can get to the original implementation
Expand Down

0 comments on commit d86a79b

Please sign in to comment.