A GPU-based Multi-level Subspace Decomposition Scheme for Hierarchical Tensor Product Bases