Abstract We study the time-sliced thawed Gaussian propagation method, which was recently proposed for solving time-dependent Schrödinger equation. introduce a triplet of quadrature-based analysis, synthesis and re-initialization operators to give rigorous mathematical formulation method. Further, we derive combined error bounds discretization wave packet transform time-propagation basis functio...