Skip to content

a potential type or bug in update_g #4

Description

@neurostar

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.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions