The FLEUR code family is a program package for calculating ground-state as well as excited-state properties of solids. It is based on the full-potential linearized augmented-plane-wave (FLAPW) method [1-4]. The strength of the FLEUR code [5,6] lies in applications to bulk, semi-infinite, two- and one-dimensional solids [7], solids of nearly all chemical elements of the periodic table, solids with complex open structures, low symmetry, with complex non-collinear magnetism [8] in combination with spin-orbit interaction [9,10], external electric fields, and the treatment of spin-dependent transport properties [11,12]. It is an all-electron method and thus treats core and valence electrons and can deal with hyperfine properties. The inclusion of local orbitals allows a systematic extension of the LAPW basis that enables a precise treatment of semicore states [13], unoccupied states [14,15], and an elimination of the linearization error in general [16]. A large variety of local and semi-local (GGA) exchange and correlation functionals are implemented, including the LDA+U approach. In recent years the code has been developed further to make contact to electronically complex materials. Hybrid functionals [17,18] and the optimized-effective-potential (OEP) method [15,19] have been implemented. Wannier functions [20] can be constructed to make contact to realistic model Hamiltonians. Excitations can be treated on the basis of the GW approximation [21,22] and ladder diagrams are included to compute spin-wave excitations [23]. The Hubbard U can be calculated in the constrained random phase approximation (cRPA) [24].