3-D parallel frequency-domain visco-acoustic wave modelling based on a hybrid direct/iterative solver