Reequilibration processes are often encountered within solids and they have long been described mathematically. However, computation of reequilibration data over entire processes is often difficult. An algorithm has been developed specifically for this purpose. It allows a fast and efficient computation of reequilibration parameters. Anisotropic diffusion can, moreover, be taken into account.