LinApart: Optimizing the univariate partial fraction decomposition

We present LinApart, a routine designed for efficiently performing the univariate partial fraction decomposition of large symbolic expressions. Our method is based on an explicit closed formula for the decomposition of rational functions with fully factorized denominators. We provide an implementation in the Wolfram Mathematica language, which can lead to very significant performance gains over the built-in Apart command. Furthermore, a C language library implementing the core functionality and suitable for interfacing with other software is also provided. Both codes are made available at https://github.com/fekeshazy/LinApart.