Abstract We propose a greedy variational method for decomposing non-negative multivariate signal as weighted sum of Gaussians, which, borrowing the terminology from statistics, we refer to Gaussian mixture model. Notably, our has following features: (1) It accepts signals, i.e., sampled functions, histograms, time series, images, etc., input. (2) The can handle general (i.e., ellipsoidal) Gauss...