Large-scale ferrofluid simulations on graphics processing units