  Non-commutative polynomial optimization is a powerful technique with numerous
applications in quantum nonlocality, quantum key distribution, causal
inference, many-body physics, amongst others. The standard approach is to
reduce such optimizations to a hierarchy of semi-definite programs, which can
be solved numerically using well-understood interior-point methods. A key, but
computationally costly, step is the formulation of moment matrices, whose size
(and hence cost) grows exponentially with the depth of the hierarchy. It is
therefore essential to have highly-optimized software to construct moment
matrices. Here, we introduce Moment: a toolkit that produces moment matrix
relaxations from the specification of a non-commutative optimization problem.
In order to obtain the absolute best performance, Moment is written in C++, and
for convenience of use provides an interface via MATLAB. We benchmark Moment's
performance, and see that it can be up to four orders of magnitude faster than
current software with similar functionality.
