  Nuclear implementation of the density functional theory (DFT) is at present
the only microscopic framework applicable to the whole nuclear landscape. The
extension of DFT to superfluid systems in the spirit of the Kohn-Sham approach,
the superfluid local density approximation (SLDA) and its extension to
time-dependent situations, time-dependent superfluid local density
approximation (TDSLDA), have been extensively used to describe various static
and dynamical problems in nuclear physics, neutron star crust, and cold atom
systems. In this paper, we present the codes that solve the static and
time-dependent SLDA equations in three-dimensional coordinate space without any
symmetry restriction. These codes are fully parallelized with the message
passing interface (MPI) library and take advantage of graphic processing units
(GPU) for accelerating execution. The dynamic codes have checkpoint/restart
capabilities and for initial conditions one can use any generalized Slater
determinant type of wave function. The code can describe a large number of
physical problems: nuclear fission, collisions of heavy ions, the interaction
of quantized vortices with nuclei in the nuclear star crust, excitation of
superfluid fermion systems by time dependent external fields, quantum shock
waves, domain wall generation and propagation, the dynamics of the
Anderson-Bogoliubov-Higgs mode, dynamics of fragmented condensates, vortex
rings dynamics, generation and dynamics of quantized vortices, their crossing
and recombinations and the incipient phases of quantum turbulence.
