A multilevel variational method for Au= Bu on composite Grids