dp_s = np.sign(np.diff(p_s))