  We derive an equilibrium statistical theory for the macroscopic description
of a ferromagnetic material at positive finite temperatures. Our formulation
describes the most-probable equilibrium macrostates that yield a coherent
deterministic large-scale picture varying at the size of the domain, as well as
it captures the effect of random spin fluctuations caused by the thermal noise.
We discuss connections of the proposed formulation to the Landau-Lifschitz
theory and to the studies of domain formation based on Monte Carlo lattice
simulations.
