small correction in the documentation
In the section "Computing local quantities: densities and currents", the function field_direction(pos,r0,delta) is defined with the condition:
if r==0: m_i=[0,0,1] This should be m_i =[0,0,-1]
This explains why we get some strange pattern in the present figures on the online documentation. With this minor correction, the result will be better and the topological hall effect clearly visible.