Computing a discrete Morse gradient from a watershed decomposition