h-p discretization