Hi John (and list),
I am trying to understand the code in spm_warp.c, and I think I may have
stumbled upon a bug:
in file: $spm2/toolbox/Deformations/spm_warp.c
At lines 798 - 801, it says:
its = mxGetPr(prhs[6])[0];
lambda = mxGetPr(prhs[6])[1];
epsilon = mxGetPr(prhs[6])[3];
meth = mxGetPr(prhs[6])[4];
but I believe it should read:
its = mxGetPr(prhs[6])[0];
lambda = mxGetPr(prhs[6])[1];
epsilon = mxGetPr(prhs[6])[2];
meth = mxGetPr(prhs[6])[3];
because in this case, prhs[6] has only 4 elements, and in C the indexing
should go 0,1,2,3, correct?
I believe the same bug is also in $spm5/toolbox/HDW/spm_warp.c
--
Matthew L. Senjem
Analyst Programmer - Imaging Systems
Mayo Clinic and Foundation
Mail Code: RO_CE_09_IMAGING
200 First St. SW
Rochester, MN 55901
[log in to unmask]
Fax: 507-284-9778
Phone: 507-538-0764
|