Gaussian wave packet transform based numerical scheme for the semi-classical Schrödinger equation with random inputs