We introduce three numerical methods for characterizing the topological phases of three-dimensional multiband Hubbard models based on twisted boundary conditions, Wilson loops, as well as the local topological marker. We ...