  We present and distribute a parallel finite-element toolbox written in the
free software FreeFem for computing the Bogoliubov-de Gennes (BdG) spectrum of
stationary solutions to one- and two-component Gross-Pitaevskii (GP) equations,
in two or three spatial dimensions. The parallelization of the toolbox relies
exclusively upon the recent interfacing of FreeFem with the PETSc library. The
latter contains itself a wide palette of state-of-the-art linear algebra
libraries, graph partitioners, mesh generation and domain decomposition tools,
as well as a suite of eigenvalue solvers that are embodied in the SLEPc
library. Within the present toolbox, stationary states of the GP equations are
computed by a Newton method. Branches of solutions are constructed using an
adaptive step-size continuation algorithm. The combination of mesh adaptivity
tools from FreeFem with the parallelization features from PETSc makes the
toolbox efficient and reliable for the computation of stationary states. Their
BdG spectrum is computed using the SLEPc eigenvalue solver. We perform
extensive tests and validate our programs by comparing the toolbox's results
with known theoretical and numerical findings that have been reported in the
literature.
