We present a free open source package for high performance simulation and numerical investigation of boson samplers and, more generally, multi-photon interferometry. Our package is written in Julia, allowing C-like performance with easy notations and fast, high-level coding. Underlying building blocks can easily be modified without complicated low-level language modifications. We present a great variety of routines for tasks related to boson sampling, such as statistical tools, optimization methods and classical samplers. Special emphasis is put on validation of experiments, where we present novel algorithms. This package goes beyond the boson sampling paradigm, allowing for the investigation of new interferometric behaviours such as bosonic bunching.