a potential type or bug in update_g
neurostar opened this issue · comments
I got a warning from MATLAB.
After review the code, It does look worth checking.
Line 35 of update_g.m
tmp_h = exp(log(g)*(0:maxl)'); % response kernel
tmp_hh = cumsum(h.*h); % hh(k) = h(1:k)'*h(1:k)
"h" looks to be tmp_h instead.
actual h is coming from rss_g function.
It may just work, since they are same.
Please check.