In any quantum dynamics method that approximates wave functions as a linearly combined basis set, non-orthogonality can be is a problem. It has been proven in previous studies that, by using the most standard form of Matching Pursuit in combination with a Gaussian wave packet ansatz, exact quantum-mechanical correspondence can be obtained for particle tunneling in one and two dimensions. This study is an attempt to prove that this approach can be generally applicable to systems of arbitrary dimension propagating with an an-harmonic potential, and that adaptive initial state sampling can be used to make the method even more computationally efficient.