zhoupc / OASIS_matlab

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

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.