  We present a new module of the parallel N-Body code P-GADGET3 for
cosmological simulations of light bosonic non-thermal dark matter, often
referred as Fuzzy Dark Matter (FDM). The dynamics of the FDM features a highly
non-linear Quantum Potential (QP) that suppresses the growth of structures at
small scales. Most of the previous attempts of FDM simulations either evolved
suppressed initial conditions, completely neglecting the dynamical effects of
QP throughout cosmic evolution, or resorted to numerically challenging
full-wave solvers. The code provides an interesting alternative, following the
FDM evolution without impairing the overall performance. This is done by
computing the QP acceleration through the Smoothed Particle Hydrodynamics (SPH)
routines, with improved schemes to ensure precise and stable derivatives. As an
extension of the P-GADGET3 code, it inherits all the additional physics modules
implemented up to date, opening a wide range of possibilities to constrain FDM
models and explore its degeneracies with other physical phenomena. Simulations
are compared with analytical predictions and results of other codes, validating
the QP as a crucial player in structure formation at small scales.
