Interplay between discretization and algebraic computation in adaptive numerical solutionof elliptic PDE problems