An automated parameter domain decomposition approach for gravitational wave surrogates using hp-greedy refinement