Efficient Multigrid based solvers for isogeometric analysis