We propose a method to approximate the solution to Schrodinger Bridge and multi marginal momentum Schrodinger Bridge problems with discrete or continuous Gaussian Mixture Model marginal distributions.