The structure of magnetic domains, i.e. regions of uniform magnetization separated by domain walls, depends on the balance of competing interactions present in ferromagnetic (or ferrimagnetic) materials. When these interactions change then domain configurations also change as a result. Magnetite provides a good test bench to study these effects, as its magnetocrystalline anisotropy varies significantly with temperature. Using spin-polarized electron microscopy to map the micromagnetic domain structure in the (001) surface of a macroscopic magnetite crystal (~1 cm size) shows complex domain patterns with characteristic length-scales in the micrometer range and highly temperature dependent domain geometries. Although heating above the Curie temperature erases the domain patterns completely, cooling down reproduces domain patterns not only in terms of general characteristics: instead, complex microscopic domain geometries are reproduced in almost perfect fidelity between heating cycles. A possible explanation of the origin of the high-fidelity reproducibility is suggested to be a combination of the presence of hematite inclusions that lock bulk domains, together with the strong effect of the first order magnetocrystalline anisotropy which competes with the shape anisotropy to give rise to the observed complex patterns