We demonstrate a Monte Carlo algorithm for efficiently simulating ferroflui
ds. By identifying particle clusters and evolving them as single units, we
reduce correlation times by more than two orders of magnitude. This method
enables accurate calculations of ferrofluid thermodynamics in the limit of
strong magnetic coupling that would be impossible by conventional means. We
apply the method to study magnetic anisotropy in dilute thin films.