Ferrofluid Simulations with the Barnes-Hut Algorithm on Graphics Processing Units