k = (dp_s[:-1] != dp_s[1:]).sum()  # numpeaks