Using excitation energies to isobaric analog states (IAS) and charge invariance, we extract nuclear symmetry coefficients, from a mass formula, on a nucleus-by-nucleus basis. Consistently with charge invariance, the coefficients vary weakly across an isobaric chain. However, they change strongly with nuclear mass and range from a_a~10 MeV at mass A~10 to a_a~22 MeV at A~240. Following the considerations of a Hohenberg-Kohn functional for nuclear systems, we determine how to find in practice the symmetry coefficient using neutron and proton densities, even when those densities are simultaneously affected by significant symmetry-energy and Coulomb effects. These results facilitate extracting the symmetry coefficients from Skyrme-Hartree-Fock (SHF) calculations, that we carry out using a variety of Skyrme parametrizations in the literature. For the parametrizations, we catalog novel short-wavelength instabilities. In comparing the SHF and IAS results for the symmetry coefficients, we arrive at narrow (+-2.4 MeV) constraints on the symmetry energy values S(rho) at 0.04<rho<0.13 fm^-3. Towards normal density the constraints significantly widen, but the normal value of energy a_a^V and the slope parameter L are found to be strongly correlated. To narrow the constraints, we reach for the measurements of asymmetry skins and arrive at a_a^V=(30.2-33.7) MeV and L=(35-70) MeV, with those values being again strongly positively correlated along the diagonal of their combined region. Inclusion of the skin constraints allows to narrow the constraints on S(rho), at 0.04<rho<0.13 fm^-3, down to +-1.1 MeV. Several microscopic calculations, including variational, Bruckner-Hartree-Fock and Dirac-Bruckner-Hartree-Fock, are consistent with our constraint region on S(rho).