We generalize the Thomas–Fermi approach to galaxy structure to include central supermassive black holes and find, self-consistently and non-linearly, the gravitational potential of the galaxy plus the central black hole (BH) system. This approach naturally incorporates the quantum pressure of the fermionic warm dark matter (WDM) particles and shows its full power and clearness in the presence of supermassive black holes. We find the main galaxy and central black hole magnitudes as the halo radius rh, halo mass Mh, black hole mass MBH, velocity dispersion σ, and phase space density, with their realistic astrophysical values, masses and sizes over a wide galaxy range. The supermassive black hole masses arise naturally in this framework. Our extensive numerical calculations and detailed analytic resolution of the Thomas–Fermi equations show that in the presence of the central BH, both DM regimes—classical (Boltzmann dilute) and quantum (compact)—do necessarily co-exist generically in any galaxy, from the smaller and compact galaxies to the largest ones. The ratio R(r) of the particle wavelength to the average interparticle distance shows consistently that the transition, R≃1, from the quantum to the classical region occurs precisely at the same point rA where the chemical potential vanishes. A novel halo structure with three regions shows up: in the vicinity of the BH, WDM is always quantum in a small compact core of radius rA and nearly constant density; in the region rA<r<ri until the BH influence radius ri, WDM is less compact and exhibits a clear classical Boltzmann-like behavior; for r>ri, the WDM gravity potential dominates, and the known halo galaxy shows up with its astrophysical size. DM is a dilute classical gas in this region. As an illustration, three representative families of galaxy plus central BH solutions are found and analyzed: small, medium and large galaxies with realistic supermassive BH masses of 105M⊙, 107M⊙ and 109M⊙, respectively. In the presence of the central BH, we find a minimum galaxy size and mass Mhmin≃107M⊙, larger (2.2233×103 times) than the one without BH, and reached at a minimal non-zero temperature Tmin. The supermassive BH heats up the DM and prevents it from becoming an exactly degenerate gas at zero temperature. Colder galaxies are smaller, and warmer galaxies are larger. Galaxies with a central black hole have large masses Mh>107M⊙>Mhmin; compact or ultracompact dwarf galaxies in the range 104M⊙<Mh<107M⊙ cannot harbor central BHs. We find novel scaling relations MBH=DMh38 and rh=CMBH43, and show that the DM galaxy scaling relations Mh=bΣ0rh2 and Mh=aσh4/Σ0 hold too in the presence of the central BH, Σ0 being the constant surface density scale over a wide galaxy range. The galaxy equation of state is derived: pressure P(r) takes huge values in the BH vicinity region and then sharply decreases entering the classical region, following consistently a self-gravitating perfect gas P(r)=σ2ρ(r) behavior.