Classical simulation of non-Gaussian bosonic circuits | Synapse