# Resolution (knot-spacing) of warps in mm

--warpres=20,16,14,12,10,6,4,4,4

# Subsampling level (a value of 2 indicates that a 2x2x2 neighbourhood is collapsed to 1 voxel)

--subsamp=1,1,1,1,1,1,1,1,1

# FWHM of gaussian smoothing

--fwhm=8,6,4,3,3,2,1,0,0

# Maximum number of iterations

--miter=5,5,5,5,5,10,10,20,20

# Relative weight of regularisation

--lambda=0.0005,0.0001,0.00001,0.0000015,0.0000005,0.0000005,0.00000005,0.0000000005,0.00000000001

# If set to 1 lambda is multiplied by the current average squared difference

--ssqlambda=1

# Regularisation model

--regmod=bending_energy

# If set to 1 movements are estimated along with the field

--estmov=1,1,1,1,1,0,0,0,0

# 0=Levenberg-Marquardt, 1=Scaled Conjugate Gradient

--minmet=0,0,0,0,0,1,1,1,1

# Quadratic or cubic splines

--splineorder=3

# Precision for calculation and storage of Hessian

--numprec=double

# Linear or spline interpolation

--interp=spline

# If set to 1 the images are individually scaled to a common mean intensity 

--scale=1