Make Matrix<Scalar>::multiply_with_vector parallel.
Tune both multigrid and multiscale.