Discrete conservation of nonnegativity for elliptic problems solved by the hp-FEM