Magnetic domains have been the subject of much scientific investigation since their theoretical existence was first postulated by P.-E. Weiss over a century ago. Up to now, the three-dimensional (3D) domain structure of bulk magnets has never been observed owing to the lack of appropriate experimental methods. Domain analysis in bulk matter thus remains one of the most challenging tasks in research on magnetic materials. All current domain observation methods are limited to studying surface domains or thin magnetic films. As the properties of magnetic materials are strongly affected by their domain structure, the development of a technique capable of investigating the shape, size and distribution of individual domains in three dimensions is of great importance. Here, we show that the novel technique of Talbot-Lau neutron tomography with inverted geometry enables direct imaging of the 3D network of magnetic domains within the bulk of FeSi crystals.