A fully adaptive interpolated stochastic sampling method for random PDEs