  First results from a fully self-consistent, temperature-dependent equation of
state that spans the whole density range of neutron stars and supernova cores
are presented. The equation of state (EoS) is calculated using a mean-field
Hartree-Fock method in three dimensions (3D). The nuclear interaction is
represented by the phenomenological Skyrme model in this work, but the EoS can
be obtained in our framework for any suitable form of the nucleon-nucleon
effective interaction. The scheme we employ naturally allows effects such as
(i) neutron drip, which results in an external neutron gas, (ii) the variety of
exotic nuclear shapes expected for extremely neutron heavy nuclei, and (iii)
the subsequent dissolution of these nuclei into nuclear matter. In this way,
the equation of state is calculated across phase transitions without recourse
to interpolation techniques between density regimes described by different
physical models. EoS tables are calculated in the wide range of densities,
temperature and proton/neutron ratios on the ORNL NCCS XT3, using up to 2000
processors simultaneously.
