QMOMKit: A C++ library for quadrature-based moment methods

QMOMKit is a C++ library implementing quadrature-based moment methods to represent and evolve number density functions in phase space. It reconstructs quadratures, sets of weighted nodes, from the raw moments of a number density and evolves them in phase space: it solves the phase-space part of a population balance equation at one location. Transport in physical space, the mesh, the time loop and the coupling to a host code belong to the host application. QMOMKit provides tools to: reconstruct a quadrature from raw moments, with the quadrature method of moments (QMOM) [McGraw, 1997] for one phase-space coordinate and the conditional quadrature method of moments (CQMOM) [Yuan & Fox, 2011] for several; let the quadrature reconstruction algorithm choose the number of nodes adaptively and bound the quadrature abscissas to an interval; choose the CQMOM conditioning order without changing the definition of the phase space; advance a state in phase space by operator splitting (Strang or Lie) of node advection and moment sources, including exchange across phase-space boundaries; check whether moments are realizable on a given support; rebuild states from moments that the host has transported in physical space; set up a case in C++ or from a .qmom configuration file, with moments given explicitly or generated from analytical distributions, mixtures and independent products; write moment and quadrature histories; compare with finite-sample Monte Carlo references from the optional companion package.