* Disable hyperpriors default setting so that the estimated priors

  on the means just use the maximum likelihood estimates.
* Added four more iterations on the Gauss-Siedel iterations of the
  template updates.
* Changed the tolerance for assessing convergence of the overall
  algorithm.

spm_mb_fit.m

@@ -167,7 +167,7 @@
         end
         fprintf('\n');
 
-        if it0==nit_max || (oE-E/nvox(dat) < sett.tol && it0>=nit_zm0)
+        if it0==nit_max || (oE-E/nvox(dat) < sett.tol*4 && it0>=nit_zm0)
             break;
         else
             % Compute deformations from velocities (unless this is to be done

spm_mb_shape.m

@@ -349,7 +349,6 @@
 % This could be made more efficient.
 
 % Parse function settings
-accel = sett.accel;
 accel = 0.99;
 B     = sett.B;
 d     = sett.ms.d;
@@ -474,7 +473,7 @@
 update_settings = [mu_settings(1:3) mu_settings(4:end)*(1-1/(K+1)) 2 2];
 if accel>0, s   = softmax(mu); end
 dmu = zeros(size(mu),'like',mu);
-for it=1:16
+for it=1:20
     for k=1:K
 
         % Diagonal elements of Hessian
@@ -536,7 +535,6 @@
 function dat = update_simple_affines(dat,mu,sett)
 
 % Parse function settings
-accel = sett.accel;
 accel = 0.99;
 B     = sett.B;
 if ~isempty(B)
@@ -634,7 +632,6 @@
 function dat = update_velocities(dat,mu,sett)
 
 % Parse function settings
-accel = sett.accel;
 accel = 0.99; % Speed up the registration
 nw    = get_num_workers(sett,4*sett.K+4*3);
 

tbx_cfg_mb.m

@@ -166,7 +166,8 @@
 pr_upd.tag     = 'hyperpriors';
 pr_upd.name    = 'Optimise';
 pr_upd.labels  = {'Yes','No'};
-pr_upd.values  = {{'b0_priors',{1000,10}}, []};
+%pr_upd.values = {{'b0_priors',{1000,10}}, []};
+pr_upd.values  = {{}, []};
 pr_upd.val     = {pr_upd.values{1}};
 pr_upd.help    = {['Specify whether the Gaussian-Wishart priors be updated at each iteration. '...
                    'Enabling this can slow down convergence if there are small numbers of subjects. '...