We modify the open-source Sky3D code, which implements a time-dependent Hartree-Fock method with a Skyrme functional, so that it uses the Barcelona-Catania-Paris-Madrid (BCPM) energy density functional instead. We run static ground state simulations and dynamical giant quadrupole resonance simulations to demonstrate the program’s capabilities. Static simulations of nuclei in their ground state are benchmarked against a reference code. The binding energies of static ground state nuclei show differences below 0.3% with respect to the reference. We also simulate quadrupole resonances employing the time-dependent code and find results that are stable and compatible with previous Sky3D simulations.