anisotropy is now computed at each RK4 step