Parametrized nucleon density distributions are widely employed for the calculation of the properties of atomic nuclei and dense inhomogeneous matter in compact stars within the Thomas—Fermi method and its extensions. We show that the use of insufficiently smooth parametrizations may deteriorate the accuracy of this method. We discuss and clarify the smoothness condition using the example of the so-called "nuclear pasta" in the neutron star mantle.