The algorithm simulates open quantum systems coupled to Gaussian environments valid for any configuration and coupling strength, and is applicable to problems with strongly coupled, or non-Markovian, ...