  In this paper a new hp-adaptive strategy for elliptic problems based on
refinement history is proposed, which chooses h-, p- or hp-refinement on
individual elements according to a posteriori error estimate, as well as
smoothness estimate of the solution obtained by comparing the actual and
expected error reduction rate. Numerical experiments show that exponential
convergence can be achieved with this strategy.
