DARTEL processing options used for each iteration. The penalizing energy term (linear elastic energy), the number of inner iterations (3), the Levenberg-Marquardt regularization (0.01), the number of cycles for the full multi-grid matrix solver (3) and the number of relaxations in each multi-grid cycle (3) as well as the third regularization parameter id (0.000001) were kept constant over iterations.