Computing multiparameter persistent homology through a discrete Morse-based approach