joakim/MinPakke.jl
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions 1

Fingers crossed

Browse source
This commit is contained in:
Joakim Skogholt 2023-05-13 21:24:50 +02:00
parent 7f7660286e
commit cec7b0a9d7
1 changed files with 2 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
src/TR.jl
View file

@ -133,7 +133,7 @@ function TRPress(TR::TRSVD, y, lambdas)
denom = broadcast(.+, broadcast(./, lambdas, TR.s'), TR.s')'; denom = broadcast(.+, broadcast(./, lambdas, TR.s'), TR.s')';
resid = broadcast(.-, y, TR.U * broadcast(./, TR.s .* (TR.U'*y), denom)); resid = broadcast(.-, y, TR.U * broadcast(./, TR.s .* (TR.U'*y), denom));
H = broadcast(.+, TR.U.^2 * broadcast(./, TR.s, denom), 1/TR.n); H = broadcast(.+, TR.U.^2 * broadcast(./, TR.s, denom), 1/TR.n);
@inline press = TRPress(TR, y, lambdas, denom, H, resid); @inline press = TRPress(TR, y, lambdas, H, resid);
return press return press
end end
@ -158,7 +158,7 @@ function TRGCV(TR::TRSVD, y, lambdas)
resid = broadcast(.-, y, TR.U * broadcast(./, TR.s .* (TR.U'*y), denom)); resid = broadcast(.-, y, TR.U * broadcast(./, TR.s .* (TR.U'*y), denom));
H = broadcast(.+, TR.U.^2 * broadcast(./, TR.s, denom), 1/TR.n); H = broadcast(.+, TR.U.^2 * broadcast(./, TR.s, denom), 1/TR.n);
@inline GCV = TRGCV(TR, y, lambdas, denom, H, resid); @inline GCV = TRGCV(TR, y, lambdas, H, resid);
return GCV return GCV
end end