# HG changeset patch # User jpicarbajal # Date 1332754110 0 # Node ID 4d562bed89951ae07c7b4953226163be0de4ed4f # Parent 0a825a74b5b605bf1879f991e4013d6ca8c4d504 linear-algebra: continue adapting nmf_bpas diff -r 0a825a74b5b6 -r 4d562bed8995 main/linear-algebra/devel/nmf_bpas.m --- a/main/linear-algebra/devel/nmf_bpas.m Sun Mar 25 23:59:38 2012 +0000 +++ b/main/linear-algebra/devel/nmf_bpas.m Mon Mar 26 09:28:30 2012 +0000 @@ -139,18 +139,19 @@ clear parser val_type ### PARSING TYPE +# TODO add callbacks here to use during main loop. See [1] % for regularized/sparse case salphaI = sqrt (par.alpha) * eye (k); zerokm = zeros (k,m); if strcmpi (par.type, 'regularized') - sbetaI = sqrt (par.beta) * eye (k); - zerokn = zeros (k,n); + sbetaI = sqrt (par.beta) * eye (k); + zerokn = zeros (k,n); elseif strcmpi (par.type, 'sparse') - sbetaE = sqrt (par.beta) * ones (1,k); - betaI = par.beta * ones (k,k); - zero1n = zeros (1,n); + sbetaE = sqrt (par.beta) * ones (1,k); + betaI = par.beta * ones (k,k); + zero1n = zeros (1,n); end ### @@ -193,7 +194,7 @@ SCconv = 0; SC_COUNT = 3; -#TODO: Replace with callbacks avoid switching each time +#TODO: [1] Replace with callbacks avoid switching each time for iter=1:par.max_iter switch par.type case 'plain'